hdu 5015 推导+矩阵快速幂

题目链接：https://vjudge.net/problem/HDU-5015
转自：https://blog.csdn.net/dragon60066/article/details/60339236
题意：有一个矩阵横排第一行初始是0,233,2333,23333…，然后给出计算公式a[i][j] = a[i - 1][j] + a[i][j - 1]，然后给出了n和m和第一列的各值要求计算a[n][m]的值。
思路：推导+矩阵快速幂。

#include <bits/stdc++.h>
using namespace std;
const int MOD=10000007;
struct mat
{
    long long a[15][15];
    mat()
    {
        memset(a,0,sizeof(a));
    }
};
int n,m;
mat mat_mul(mat x,mat y)
{
    mat res;
    memset(res.a,0,sizeof(res.a));
    for(int i=1; i<=n+2; i++)
    {
        for(int j=1; j<=n+2; j++)
        {
            for(int k=1; k<=n+2; k++)
            {
                res.a[i][j]=(res.a[i][j]+x.a[i][k]*y.a[k][j])%MOD;
            }

        }
    }
    return res;
}
mat mat_pow(mat a,int c)
{
    mat ans;
    for(int i=1; i<=n+2; i++)
    {
        ans.a[i][i]=1;
    }
    while(c)
    {
        if(c&1)
        {
            ans=mat_mul(a,ans);
        }

        a=mat_mul(a,a);
        c=c>>1;
    }
    return ans;
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        mat A,B;
        A.a[1][1]=23;
        for(int i=1; i<=n; i++)
        {
            scanf("%lld",&A.a[i+1][1]);
        }
        A.a[n+2][1]=3;
        for(int i=1; i<=n+1; i++)
        {
            B.a[i][1]=10;
        }
        for(int i=1; i<=n+2; i++)
        {
            B.a[i][n+2]=1;
        }
        for(int i=1; i<n+2; i++)
        {
            for(int j=2; j<=i; j++)
            {
                B.a[i][j]=1;
            }
        }
        B=mat_pow(B,m);
        A=mat_mul(B,A);
        printf("%lld\n",A.a[n+1][1]);
    }
    return 0;
}