symons

摘要： 本体有2种方法，一种是矩阵快速幂，一种是找规律。我用的是矩阵快速幂。f(1)=1; f(2)=1; f(n)=( A*(f(n-1)) + B*(f(n-2)) )mod7;[ f(n) ] = [ A B ] * [ f(n-1) ] => [ f(n) ] = [ A B ] (n-2） * [ f(1) ][ f(n-1) ] [ 1 0 ] [ f(n-2) ] => [ f(n-1) ] [ 1 0 ] [ f(2) ] 1 #include 2 #include 3 #include 4 #include 5 using namespace std; 6 struct 阅读全文

摘要： 1 //*代表的是矩阵乘的运算,res是结果矩阵2 while(n){3 if(N&1)4 res=res*A;5 n>>=1;6 A=A*A;7 } 阅读全文