pku3070 （矩阵1）

http://162.105.81.212/JudgeOnline/problem?id=3070

题意很简单，给出一个数n，就是要求第n个斐波纳契数模10000（就是求最后4位数）。

分析：要是用递推式F[n] = F[n-1] + F[n-2]来计算，则复杂度是O(n)的，由于n的范围达到 0 ≤ n ≤ 1,000,000,000肯定会超时的。

按题目要求要用矩阵来做的，详情看原题。

那么对于一个矩阵怎么来取模呢？实际上这与整数的取模幂运算a^b mod n 是一样的，代码如下：

int modexp(int a, int b, int n) { int t = a, ret = 1; while(b != 0) { if(b % 2 == 1) ret *= t % n; t = t * t % n; b /= 2; } return ret; } 

根据上面的过程实现矩阵的模取幂运算。

#include<iostream> using namespace std; const int mod = 10000; struct matrix //定义矩阵2x2 { long a[2][2]; }; matrix mul(matrix x, matrix y, long mod) //定义返回的矩阵 { matrix ret; //按照矩阵相乘求ret每个元素的值； ret.a[0][0] = (x.a[0][0] * y.a[0][0] + x.a[0][1] * y.a[1][0]) % mod; ret.a[0][1] = (x.a[0][0] * y.a[0][1] + x.a[0][1] * y.a[1][1]) % mod; ret.a[1][0] = (x.a[1][0] * y.a[0][0] + x.a[1][1] * y.a[1][0]) % mod; ret.a[1][1] = (x.a[1][0] * y.a[0][1] + x.a[1][1] * y.a[1][1]) % mod; return ret; } matrix mypow(matrix x, long e, long mod) //求矩阵x的幂取模，e为指数； { matrix ret, temp; //返回矩阵和临时矩阵 if(e == 0) //指数为0时，返回单位矩阵 { ret.a[0][0] = 1; ret.a[0][1] = 0; ret.a[1][0] = 0; ret.a[1][1] = 1; return ret; } if(e == 1) return x; //指数为1时，返回原来的矩阵 temp = mypow(x, e>>1, mod); //x的e/2次方； ret = mul(temp, temp, mod); //ret = temp * temp; if(e & 1) //如果e为奇数，ret乘以x ret = mul(ret, x, mod); return ret; } int main() { int n; matrix ans; while(scanf("%d",&n) && n!=-1) { ans.a[0][0] = 1; ans.a[0][1] = 1; //初始化 ans.a[1][0] = 1; ans.a[1][1] = 0; if(n) { ans = mypow(ans, n-1, mod); //注意是n-1; printf("%ld/n",ans.a[0][0]); } else printf("0/n"); } return 0; }