codevs 1250 Fibonacci数列

时间限制: 1 s
空间限制: 128000 KB
题目等级 : 钻石 Diamond
题解
题目描述 Description
定义：f0=f1=1, fn=fn-1+fn-2(n>=2)。{fi}称为Fibonacci数列。

输入n，求fn mod q。其中1<=q<=30000。

输入描述 Input Description
第一行一个数T(1<=T<=10000)。

以下T行，每行两个数，n，q(n<=109, 1<=q<=30000)

输出描述 Output Description
文件包含T行，每行对应一个答案。

样例输入 Sample Input
3

6 2

7 3

7 11

样例输出 Sample Output
1

0

10

数据范围及提示 Data Size & Hint
1<=T<=10000

n<=109, 1<=q<=30000

思路：斐波那契数列的相关题，求f(n)%mod的值，由于n比较大，用普通的求斐波那契数列的方法肯定会超时，那么我们就要用到一个神奇的东西——矩阵乘法了。
关于矩阵乘法怎么算，在这里我就不再赘述了，如果有不知道的童鞋请自行百度或查阅相关资料。
关于矩阵乘法，我们发现如下规律：

这么一来，我们就可以用矩阵的快速幂来快速处理f(n)%mod的值了。矩阵快速幂的原理与普通快速幂原理大致相同，核心都是二进制加速，不同的是普通快速幂的初始数值是1（即单位数字），而矩阵快速幂的初始矩阵是个2*2的单位矩阵，即

之后的方法同普通的快速幂。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
int T,n,mod;
struct matrix{
    int ma[3][3];
}a,b,c;
matrix multiply(matrix A,matrix B)
{
    matrix C;
    for(int i=1;i<=2;i++)
        for(int j=1;j<=2;j++)
            C.ma[i][j]=0;
    for(int i=1;i<=2;i++)
        for(int j=1;j<=2;j++)
            for(int k=1;k<=2;k++)
            {
                C.ma[i][j]=(C.ma[i][j]%mod+(A.ma[i][k]%mod*B.ma[k][j]%mod)%mod)%mod;
            }
    return C;
}
matrix ksm(matrix A,int y)
{
    matrix ret,base=A;
    ret.ma[1][1]=1,ret.ma[1][2]=0;
    ret.ma[2][1]=0,ret.ma[2][2]=1;
    while(y)
    {
        if(y&1) ret=multiply(ret,base);
        base=multiply(base,base);
        y>>=1;
    }
    return ret;
}
int main()
{
    a.ma[1][1]=1,a.ma[1][2]=1;
    a.ma[2][1]=1,a.ma[2][2]=0;
    b.ma[1][1]=1,b.ma[1][2]=0;
    b.ma[2][1]=1,b.ma[2][2]=0;
    scanf("%d",&T);
    int ans;
    while(T--)
    {
        scanf("%d%d",&n,&mod);
        c=multiply(ksm(a,n-1),b);
        ans=c.ma[1][1]%mod;
        printf("%d\n",ans);
    }
    return 0;
}