Fibonacci数列（codevs 1250）

题目描述 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<=10⁹, 1<=q<=30000)

输出描述 Output Description

文件包含T行，每行对应一个答案。

样例输入 Sample Input

6 2

7 3

7 11

样例输出 Sample Output

数据范围及提示 Data Size & Hint

1<=T<=10000

n<=10⁹, 1<=q<=30000

/*

  矩阵乘法求斐波那契数列

*/

#include<cstdio>

#include<iostream>

#define ll long long

using namespace std;

ll a[][],b[][],c[][],ans[][],n,mod;

void work()

{

    b[][]=ans[][]=;

    b[][]=ans[][]=;

    b[][]=ans[][]=;

    b[][]=ans[][]=;

    while(n)

    {

        if(n&)

        {

            for(int k=;k<=;k++)

              for(int i=;i<=;i++)

                for(int j=;j<=;j++)

                  c[i][j]=(c[i][j]+(ans[i][k]*b[k][j]%mod)%mod);

            for(int i=;i<=;i++)

              for(int j=;j<=;j++)

                ans[i][j]=c[i][j],c[i][j]=;

        }

        for(int k=;k<=;k++)

          for(int i=;i<=;i++)

            for(int j=;j<=;j++)

              c[i][j]=(c[i][j]+(b[i][k]*b[k][j]%mod)%mod);

        for(int i=;i<=;i++)

          for(int j=;j<=;j++)

            b[i][j]=c[i][j],c[i][j]=;

        n/=;

    }

    printf("%lld\n",(ans[][]+ans[][])%mod);

}

int main()

{

    int T;

    scanf("%d",&T);

    while(T--)

    {

        cin>>n>>mod;

        if(n<=)printf("%lld\n",n);

        else n--,work();

    }

    return ;

}

巴特西

Fibonacci数列（codevs 1250）

最新文章

热门文章