http://acm.hdu.edu.cn/showproblem.php?pid=2256&&构造矩阵求值

#include<iostream> #include<string.h> #include<algorithm> #define N 1024 using namespace std; typedef long long L; typedef struct str { L s[2][2];//注意数据范围 }Node; Node a,b; Node ceil(Node p,Node q) { Node c; memset(c.s,0,sizeof(c.s)); for(int i=0;i<2;++i) for(int j=0;j<2;++j) for(int t=0;t<2;++t) c.s[i][j]=(c.s[i][j]+p.s[i][t]*q.s[t][j])%N; return c; } Node doit(int k) { Node p=a,q=b; while(k) { if(1&k) p=ceil(p,q); q=ceil(q,q); k=k>>1; } return p; } int main() { int T; scanf("%d",&T); while(T--) { a.s[0][0]=a.s[1][1]=1; a.s[0][1]=a.s[1][0]=0; b.s[0][0]=b.s[1][1]=5; b.s[0][1]=12;b.s[1][0]=2; int n; scanf("%d",&n); if(n==1) { printf("9\n"); continue; } n--; Node c=doit(n); int sum=(5*c.s[0][0]+2*c.s[0][1])%N; printf("%d\n",(2*sum-1)%N);//注意结果是对1024取余 }return 0; } /* 思路：构造矩阵，然后进行矩阵二分幂运算。 构造过程： (sqrt(2)+sqrt(3))^2=5+2*sqrt(6); (5+2*sqrt(6))^n=An+Bn*sqrt(6); Bn*sqrt(6)=An-(5-2*sqrt(6))^n；//因为此项等于(0.101...)^n 约等于0 (5+2*sqrt(6))^n=2*An-Bn*sqrt(6); An+Bn*sqrt(6)=(A(n-1)+B(n-1)*sqrt(6))*(5+2*sqrt(6))=(5*A(n-1)+12B(n-1))+(5*B(n-1)+2*A(n-1))*sqrt(6) 所以有： An=5*A(n-1)+12B(n-1) Bn=2*A(n-1)+5*B(n-1) (5+2*sqrt(6))^n=An+Bn*sqrt(6)+(5-2*sqrt(6))^n=An-Bn*sqrt(6)=2An-(0.101...)^n 所以(5+2*sqrt(6))^n=2An-1 最后的答案ans=2An-1; 而矩阵为 | 5 12| | 2 5| */