北邮12月份月赛 D 矩阵二分幂

     是一道数学题吧，用矩阵二分幂可以轻松解决。以前在NYOJ上做过类似的题，所以比赛时直接贴的以前的代码，，水过。。。。。数据太水啊。。。。。题目：

D  零零漆的跳槽

Accept:42    Submit:117Time Limit:1000MS    Memory Limit:65536KB

Description

世上没有铁饭碗,　也没有永远的老板, 跳槽风也刮到了组织, 尽管组织极力想挽留零零漆, 但他还是依然要跳到另一个有关部门. 凭他的经验以及高超的杀猪功力, 他顺利的通过了有关部门前面的测试, 来到了算法测试关. 给他的问题很简单----给两个整数n,m, 求斐波纳契数fib[n] % m...

算卖肉钱久了, 零零漆还真想不起怎么去计算斐波纳契数了, 但是他在考场竟然能通过安装在皮鞋里的电话和你通信, 他只能寄希望于你了...

当然你还是知道fib[n]的定义的:

/ fib[0]=0

| fib[1]=1

\ fib[n]=fib[n-1]+fib[n-2] (n>1)

Input

输入第一行是一个整数 c, 0 < c <= 5000, 表示要计算多少个fib[n]

接下来的c行, 每行有两个整数n,m, 0 <= n <= 2147483647, 0 < m < 32767, 意义如前所述

Output

对于每对n,m, 对应输出单独的一行, 包含一个整数 r = fib[n] % m

Sample Input

8

42 8468

6335 6501

19170 5725

11479 9359

26963 4465

5706 8146

23282 6828

9962 492

Sample Output

3712

3547

1210

5683

1502

5894

5113

1

ac代码：

#include <stdio.h>//int f(int x,int y,int a,int b)//{//return a*x+b*y;//}int main(){int n,m,i,j,k,t;scanf("%d",&t);while(t--){scanf("%d%d",&n,&m);int a[2][2]={1,0,1,0},b[2][2]={1,1,1,0},t[2][2]; while(n){if(n&1){t[0][0]=(a[0][0]*b[0][0]+a[0][1]*b[1][0])%m;t[0][1]=(a[0][0]*b[0][1]+a[0][1]*b[1][1])%m;t[1][0]=(a[1][0]*b[0][0]+a[1][1]*b[1][0])%m;t[1][1]=(a[1][0]*b[0][1]+a[1][1]*b[1][1])%m;a[0][0]=t[0][0];a[0][1]=t[0][1];a[1][0]=t[1][0];a[1][1]=t[1][1];}t[0][0]=(b[0][0]*b[0][0]+b[0][1]*b[1][0])%m;t[0][1]=(b[0][0]*b[0][1]+b[0][1]*b[1][1])%m;t[1][0]=(b[1][0]*b[0][0]+b[1][1]*b[1][0])%m;t[1][1]=(b[1][0]*b[0][1]+b[1][1]*b[1][1])%m;b[0][0]=t[0][0];b[0][1]=t[0][1];b[1][0]=t[1][0];b[1][1]=t[1][1];n=n>>1;}printf("%d\n",a[0][1]);}}