hdu 1005 Number Sequence（矩阵连乘+二分快速求幂）

题目：
http://acm.hdu.edu.cn/showproblem.php?pid=1005

代码：

 #include<iostream>
 #include<stdio.h>
 #include<string.h>
 #include<algorithm>
 using namespace std;
 const int mod=;
 struct matrix
 {
   int f[][];
 };

 matrix mul(matrix a,matrix b)
 {
     int i,j,k;
     matrix c;
     memset(c.f,,sizeof(c.f));
     for(i=;i<;i++)
     {
         for(j=;j<;j++)
         {
             for(k=;k<;k++)
             {
                 c.f[i][j]+=a.f[i][k]*b.f[k][j];
                 c.f[i][j]%=mod;
             }
         }
     }
     return c;
 }

 matrix pow_mod(matrix a,int b)
 {
     matrix s;
     int t;
     memset(s.f,,sizeof(s.f));
     for(int i=;i<;i++)
     {
         s.f[i][i]=;
     }
     while(b)
     {
         /*if(b&1)
         {
             s=mul(s,a);
         }
         a=mul(a,a);
         b=b>>1;*/
         t=b%;
         b/=;
         if(t!=)
             s=mul(s,a);
          a=mul(a,a);
     }
     return s;
 }
 int main()
 {
     int a,b,n;
     while(scanf("%d%d%d",&a,&b,&n))
     {
         if(a==&&b==&&n==)
             break;
          matrix e;
     e.f[][]=a;
     e.f[][]=;
     e.f[][]=b;
     e.f[][]=;
     e=pow_mod(e,n-);
     printf("%d\n",(e.f[][]+e.f[][])%mod);
     }

     return ;
 }