hdu 4602 Partition（快速幂）

推公式+快速幂

公式有很多形式，可以写矩阵

1、前n-1项和的两倍+2的(n-2)次方，这个写不出啥

2、递推式：f(n)=2*f(n-1)+2的(n-3)次方

3、公式：2的（n-k-2）次方*（n-k+1）+2的（n-k-1）

代码什么的看他的吧http://blog.csdn.net/liuledidai/article/details/9449301

第一次写矩阵就不献丑了

#include<stdio.h>

const int mod=1e9+7;

#define LL __int64

LL p[2][2];
LL q[1][2];

LL bb(LL a,LL b,LL c,LL d)
{
    return (a*c%mod+b*d%mod)%mod;
}

LL aa(int n)
{
    LL a,b,c,d;
    LL x,y;
    p[0][0]=2;
    p[0][1]=0;
    p[1][0]=1;
    p[1][1]=2;

    q[0][0]=2;
    q[0][1]=1;
    while(n)
    {
        a=p[0][0];b=p[0][1];c=p[1][0];d=p[1][1];
        x=q[0][0];y=q[0][1];
        if(n&1){
            q[0][0]=bb(x,y,a,c);
            q[0][1]=bb(x,y,b,d);
        }
        p[0][0]=bb(a,b,a,c);
        p[0][1]=bb(a,b,b,d);
        p[1][0]=bb(c,d,a,c);
        p[1][1]=bb(c,d,b,d);

        n>>=1;
    }
    return (q[0][0]%mod);
}

int main()
{
    int T,n,m;
    scanf("%d",&T);

    while(T--)
    {
        scanf("%d%d",&n,&m);
        if(n<m)
            printf("0
");
        else if(n-m==0)
            printf("1
");
        else if(n-m==1)
            printf("2
");
        else 
            printf("%I64d
",aa(n-m-1));
    }
    return 0;
}