hdu 2855

这道题要用到整体的思想，原因在于有组合数，组合数C(N,K)和C(n,k+1)之间没有直接的关系。受到昨天的启发，利用二项式公式，就把fib(n)拆成直接的表达式，也就是带根号5的内个表达式，之后整理，发现g(n) = 3 * g(n-1) -g(n-2), 附代码：

#include <iostream>
#include <cstdio>

#define N 40
#define ss(a) scanf("%d",&a)
#define cl(a) memset(a,0,sizeof(a))

using namespace std;

int a[N][3][3],r[3],mod;

void init()
{
    cl(a);cl(r);
    a[0][1][1]=3;a[0][1][2]=mod-1;
    a[0][2][1]=1;a[0][2][2]=0;
    r[1]=1;r[2]=0;
}

void rmul(int x)
{
    int i,j,k;
    for (i=1;i<=2;i++)
        for (j=1;j<=2;j++)
        {
            int w=0;
            for (k=1;k<=2;k++) w=(w+a[x-1][i][k]*a[x-1][k][j])%mod;
            a[x][i][j]=w;
        }
}

void zmul(int x)
{
    int y1,y2;
    y1=(a[x][1][1]*r[1]+a[x][1][2]*r[2])%mod;
    y2=(a[x][2][1]*r[1]+a[x][2][2]*r[2])%mod;
    r[1]=y1;
    r[2]=y2;
}

void quick(int n)
{
    int i=0;
    while (n>0)
    {
        if (n&1) zmul(i);
        i++;
        rmul(i);
        n>>=1;
    }
}

int main()
{
    int n,T;
    ss(T);
    while (T--)
    {
        ss(n);ss(mod);
        if (n<=1) cout<<n<<endl;
        else
        {
            init();
            quick(n-1);
            cout<<r[1]<<endl;
        }
    }
    return 0;
}