题目链接: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;
}