HDU5950 (2016沈阳现场赛)

给一个方程p[i]=p[i-1]+p[i-2]*2+i^4。

i的范围是2^31

很显然这里需要矩阵快速幂,可是这个i^4需要考虑如何用矩阵表示（用比i小的形式，或者说需要之前的数表达）。

这里可以考虑将其二项式展开，变成(x-1)^4+4*(x-1)^3+6*(x-1)^2+4*(x-1)^1+1

如果化成这样，就可以用矩阵的方法去表达x^4,x^3,x^2,x^1了

|1   2   1   4   6   4   1|     |   f(n+1)   |           |    f(n+2)    |


|1   0   0   0   0   0   0|     |     f(n)     |           |    f(n+1)    |


|0   0   1   4   6   4   1|     | (n+1)^4  |           |  (n+2)^4   |


|0   0   0   1   3   3   1|  * | (n+1)^3  |     =    |  (n+2)^3   |


|0   0   0   0   1   2   1|     | (n+1)^2  |           |  (n+2)^2  |


|0   0   0   0   0   1   1|     |    n+1     |           |     n+2      |


|0   0   0   0   0   0   1|     |      1       |           |       1        |

这个是化成的矩阵方式，那么就可以用矩阵快速幂进行求解了。

这里有个坑点wa了好几发，就是这个mod数，其是大于int类型了， 需要long long 去存储

#include <iostream>
#include <algorithm>
#include <vector>
#include<stdio.h>
#include <string.h>
#include <stack>
#include <queue>
using namespace std;
typedef long long ll;
const ll mod=2147493647;
struct ttt{
    ll map1[8][8];
};
int n;
ttt mul(ttt &a,ttt &b){
    ttt c;
    int i,j,k;
    for(i=1;i<=n;i++){
        for(j=1;j<=n;j++){
            c.map1[i][j]=0;
            for(k=1;k<=n;k++){
                c.map1[i][j]+=a.map1[i][k]*b.map1[k][j];
                c.map1[i][j]%=mod;
            }
        }
    }
    return c;
}
ttt mul_self(ttt &a){
    ttt b;
    int i,j,k;
    for(i=1;i<=n;i++)
    for(j=1;j<=n;j++){
        b.map1[i][j]=0;
        for(k=1;k<=n;k++){
            b.map1[i][j]+=a.map1[i][k]*a.map1[k][j];
            b.map1[i][j]%=mod;
        }
    }
    return b;
}
ttt pow1(ttt &a,ll k){
    ttt b;
    int i,j;
    memset(b.map1,0,sizeof(b.map1));
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
        if(i==j)
        b.map1[i][j]=1;
    while(k){
        if(k&1)
        b=mul(b,a);
        k=k>>1;
        a=mul_self(a);
    }
    return b;
}
ttt A;
ll aa[7][7]={{1,2,1,4,6,4,1},{1,0,0,0,0,0,0},{0,0,1,4,6,4,1},{0,0,0,1,3,3,1},{0,0,0,0,1,2,1},{0,0,0,0,0,1,1},{0,0,0,0,0,0,1}};
int main(){
    int i,j,k,f1,f2,f3,f4,t1,t2,t3,t4;
    //freopen("in.txt","r",stdin);
    //freopen("out2.txt","w",stdout);
    int T;
    while(scanf("%d",&T)==1){
    ll g1,g2,g3,g4;
    n=7;
    while(T--){
        scanf("%lld %lld %lld",&g1,&g2,&g3);
        g2%=mod;g3%=mod;
        if(g1==1){
        printf("%lld\n",g2);continue;
        }else if(g1==2){
        printf("%lld\n",g3);continue;
        }
        ttt B;
        for(i=1;i<=7;i++)
            for(j=1;j<=7;j++){
                A.map1[i][j]=aa[i-1][j-1];
            }
        B=pow1(A,g1-2);
    printf("%lld\n",(B.map1[1][1]*g3%mod+B.map1[1][2]*g2+B.map1[1][3]*16ll%mod+B.map1[1][4]*8ll%mod+B.map1[1][5]*4ll+
    B.map1[1][6]*2ll+B.map1[1][7])%mod);
    }
    }
    return 0;
}