HDU 5015 233 Matrix(矩阵快速幂）

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<string>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#include<bitset>
using namespace std;
#define REPF( i , a , b ) for ( int i = a ; i <= b ; ++ i )
#define REP( i , n ) for ( int i = 0 ; i < n ; ++ i )
#define CLEAR( a , x ) memset ( a , x , sizeof a )
typedef long long LL;
typedef pair<int,int>pil;
const int INF = 0x3f3f3f3f;
const int MOD=10000007;
LL n,m;
LL num[20];
struct Matrix{
    LL mat[15][15];
    void Clear()
    {
        CLEAR(mat,0);
    }
};
Matrix mult(Matrix m1,Matrix m2)
{
    Matrix ans;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
        {
            ans.mat[i][j]=0;
            for(int k=0;k<n;k++)
                ans.mat[i][j]=(ans.mat[i][j]+m1.mat[i][k]*m2.mat[k][j])%MOD;
        }
    return ans;
}
Matrix Pow(Matrix m1,LL b)
{
    Matrix ans;ans.Clear();
    for(int i=0;i<n;i++)
        ans.mat[i][i]=1;
    while(b)
    {
        if(b&1)
            ans=mult(ans,m1);
        b>>=1;
        m1=mult(m1,m1);
    }
    return ans;
}

int main()
{
    while(~scanf("%lld%lld",&n,&m))
    {
        Matrix A,B;A.Clear();B.Clear();
        B.mat[0][0]=23;B.mat[n+1][0]=3;
        REPF(i,1,n)  scanf("%lld",&B.mat[i][0]);
        n+=2;//n==5
        for(int i=0;i<n-1;i++)  A.mat[i][0]=10;
        for(int i=0;i<n;i++)  A.mat[i][n-1]=1;
        for(int i=1;i<n-1;i++)
            for(int j=1;j<=i;j++)  A.mat[i][j]=1;
        A=Pow(A,m);
        B=mult(A,B);
        printf("%lld\n",B.mat[n-2][0]);
    }
    return 0;
}