ihe 1110310415 linyiming QQ：545931195

/*愧疚了，好久没写什么有意思的题了。 这道基于二分的矩阵连乘，方法很经典。 求幂二分，求和再二分。经典。 求幂的时候，若奇数，A^n=A^(n/2)*A^(n/2)*A. 偶数时 A^n=A^(n/2)*A^(n/2)。 求和的时候，若为奇数，A^1+A^2+...+A^(2*m+1)=(A^1+...+A^m)+(A^1+...+A^m)*A^m*A, 若为偶数 A^1+A^2+...+A^

#include #include #include using namespace std; int tmp[21][21],ans[21][21],b[21][21]; int N; void mul(int a[][21],int b[][21]) { memset(tmp,0,sizeof(tmp)); for(int i=1; i<=N; i++)

#include #include #include #define mod 1000000000 using namespace std; typedef long long LL; LL a[4][4],tmp[4][4],ans[4][4]; void mul(LL a[][4],LL b[][4]) { for(int i=1; i<=4; i++) for(

#include #include #include using namespace std; typedef long long LL; LL tmp[11][11],a[11][11],b[11][11]; int mod; void mul(LL a[][11],LL b[][11]) { for(int i=0; i<10; i++) for(int j=0;

#include #include #define mod 9999997 long long tmp[3][3],ans[3][3],a[3][3]; void mul(long long a[][3],long long b[][3]) { for(int i=1; i<=2; i++) for(int j=1; j<=2; j++) {

#include #include #include #define mod 10000000 using namespace std; typedef long long LL; LL a[4][4],tmp[4][4],ans[4][4]; void mul(LL a[][4],LL b[][4]) { for(int i=1; i<=3; i++) for(in

#include #include #include #include #define c (sqrt(5.0)+1.0)/2.0 #define mod 10000 using namespace std; typedef long long LL; LL tmp[3][3],a[3][3],ans[3][3]; void mul(LL a[][3],LL b[][3]) { f

/*题目所给的a，b和f的初始值很重要。 一次来确定矩阵的元素和幂。 [a b1 0]×[f(n−1)f(n−2)]=[f(n)f(n−1)] [a b1 0]n−2×[11]=[f(n)f(n−1)] 因为f[1]=1,f[2]=1, 所以递推到n-2次幂，间接矩阵为1,1. 如果初始条件改变的话，次幂也间接矩阵也会随之改变。 要灵活应用。*/ #include #include const

#include #include const int mod=10000; int A[3][3],ans[3][3]; void mul(int a[][3],int b[][3]) { int tmp[3][3]; for(int i=1; i<=2; i++) for(int j=1; j<=2; j++) { t

构造矩阵如下： Ai*bi AX*BX AX*BY AY*BX AY*BY 0 a(i-1)*b(i-1) Ai 0 AX 0 AY 0 a(i-1) Bi