http://acm.hdu.edu.cn/showproblem.php?pid=3483
被别人黑了,做这题.....思路上卡了,在网上搜了题解后,恍然大悟……把我的代码该了下就行了,然后就AC了
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
ll a[52][52],b[52][52],c[52][52],mod,zu[52][52];
void cal(ll a1[][52],ll b1[][52],ll c1[][52],int size)
{
for(int i=0;i<size;i++)
for(int j=0,r;j<size;j++)
for(r=0,c1[i][j]=0;r<size;r++)
c1[i][j]=(c1[i][j]+a1[i][r]*b1[r][j])%mod;
}
ll quick_cal(int n,int size)
{
while(n)
{
if(n&1)
{
cal(a,c,b,size);
memcpy(c,b,sizeof(b));
}
n>>=1;
cal(a,a,b,size);
memcpy(a,b,sizeof(b));
}
return c[0][1]%mod;
ll ans=size-2;
for(int i=0;i<size;i++)
ans=(ans+c[0][i])%mod;
return ans;
}
int main()
{
int n,x;
zu[0][0]=1;
zu[1][0]=zu[1][1]=1;
for(int i=2;i<=50;i++)
{
zu[i][0]=zu[i][i]=1;
for(int j=1;j<i;j++)
zu[i][j]=zu[i-1][j-1]+zu[i-1][j];
}
while(scanf("%d%d%I64d",&n,&x,&mod)==3)
{
if(x<0&&x<0&&mod<0) break;
memset(a,0,sizeof(a));
a[0][0]=1;
for(int i=1;i<x+2;i++)
a[0][i]=x*zu[x][i-1]%mod;
for(int i=1;i<x+2;i++)
{
for(int j=1;j<=i;j++)
a[i][j]=x*zu[i-1][j-1]%mod;
}
memset(c,0,sizeof(c));
for(int i=0;i<x+2;i++)
c[i][i]=1;
printf("%I64d\n",quick_cal(n,x+2));
}
return 0;
}