hdu 3483

来源：互联网发布：三维码生成软件编辑：程序博客网时间：2024/04/27 11:03

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;}