51nod 1126 矩阵快速幂 水

有一个序列是这样定义的：f(1) = 1, f(2) = 1, f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7.

给出A，B和N，求f(n)的值。

Input

输入3个数：A,B,N。数字之间用空格分割。(-10000 <= A, B <= 10000, 1 <= N <= 10^9)

Output

输出f(n)的值。

Input示例

3 -1 5

Output示例

6
题意：f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7
题解：矩阵快速幂处理 
 1  1   a  1
      *
 0  0   b  0

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<vector>
 #define ll __int64
 using namespace std;
 ll a,b,n;
 struct matrix
 {
     ll m[][];
 }ans,exm;
 struct matrix matrix_mlit(struct matrix aa,struct matrix bb)
 {
     struct matrix there;
     for(int i=;i<;i++)
     {
         for(int j=;j<;j++)
         {
             there.m[i][j]=;
             for(int u=;u<;u++)
             there.m[i][j]=(there.m[i][j]+aa.m[i][u]*bb.m[u][j]%)%;
         }
     }
     return there;
 }
 ll  matrix_quick(ll aa,ll bb,ll gg)
 {
     exm.m[][]=aa;exm.m[][]=;
     exm.m[][]=bb;exm.m[][]=;
     ans.m[][]=;ans.m[][]=;
     ans.m[][]=;ans.m[][]=;
     gg-=;
     while(gg)
     {
         if(gg&)
         ans=matrix_mlit(ans,exm);
         exm=matrix_mlit(exm,exm);
          gg>>=;
     }
     return ans.m[][];
 }
 int main()
 {
     scanf("%I64d %I64d %I64d",&a,&b,&n);
     if(n==){
         printf("1\n");
         return ;
         }
     if(n==){
         printf("1\n");
         return ;
     }
     printf("%I64d\n",(matrix_quick(a,b,n)+)%);
     return ;
 }