hdu-1005 Number Sequence 解题报告

题面链接：hdu-1005 Number Sequence

题面：f(1) = 1, f(2) = 1, f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7.
给你这个公式 给你A B 让你输出第n个数的大小；
思路：找规律是不可能的 一辈子不可能的 有的规律你找不到的 找到了也T了 别找了
所以用矩阵快速幂来做 ；
如果你问矩阵的性质在这有什么用 那我就这样解释了
矩阵快速幂传送门（转载）
看懂了 估计就会了；

#include<iostream>
using namespace std;
struct Matrix{
	int m[2][2];
};
Matrix mul(Matrix a,Matrix b){
    Matrix ans;
    for(int i = 0; i < 2; i++){
        for(int j = 0; j < 2; j++){     
            ans.m[i][j] = 0;
            for(int k = 0; k < 2; k++){
                ans.m[i][j] = (ans.m[i][j] + a.m[i][k] * b.m[k][j] % 7) % 7;
            }

        }
    }
    return ans;
}
Matrix exp_mod(Matrix a,int n)
{
    Matrix t;
    t.m[0][0]=t.m[1][1]=1;
    if(n==0) return t;
    if(n==1) return a;
    t=exp_mod(a,n/2);
    t=mul(t,t);
    if((n&1)==1) t=mul(t,a);
    return t;
}
int main(){
	int A,B,n;
	while(~scanf("%d%d%d",&A,&B,&n))
	{
		if(A==0&&B==0&&n==0)break;
		if(n<=2){
			printf("1\n");
			continue;
		}
		Matrix a,ans;
		a.m[0][0]=A;
		a.m[0][1]=B;
		a.m[1][0]=1;
		a.m[1][1]=0;
		ans.m[0][0]=1;
		ans.m[1][0]=1;
		ans.m[0][1]=0;
		ans.m[1][1]=0;
		n-=2;
		ans=mul(exp_mod(a,n),ans);
		cout<<ans.m[0][0]<<endl;
	}
}