hdu 3483 A Very Simple Problem

两种构造的方式都是正确的；

1.

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 60
#define ll long long
using namespace std;
int x;
ll M,n;
struct matrix
{
    int len_x;
    int len_y;
    ll data[maxn][maxn];
    void ini()
    {
        len_x=;
        len_y=;
        memset(data,,sizeof data);
    }
};

matrix mat_big;
ll f[maxn][maxn];
matrix s;

matrix mat_mul(matrix mat,matrix t)
{
    matrix ans;
    ans.ini();
    ans.len_x=mat.len_x;
    ans.len_y=t.len_y;
    for(int i=; i<mat.len_x; i++)
    {
        for(int j=; j<t.len_y; j++)
        {
            for(int k=; k<mat.len_y; k++)
            {
                ans.data[i][j]+=(mat.data[i][k]*t.data[k][j])%M;
                ans.data[i][j]%=M;
            }
        }
    }
    return ans;
}

int main()
{
    f[][]=;
    for(int i=; i<maxn; i++)
        for(int j=; j<=i; j++)
        {
            f[i][j]=f[i-][j]+f[i-][j-];
        }
    while(scanf("%lld%d%lld",&n,&x,&M)&&n!=-)
    {
        mat_big.ini();
        mat_big.len_x=x+;
        mat_big.len_y=x+;
        for(int i=; i<=x; i++)
            for(int j=; j<=i; j++)
            {
                mat_big.data[i][j]=x*f[i][j];
                if(mat_big.data[i][j]>M)mat_big.data[i][j]%=M;
            }
//        for(int i=0; i<=x; i++)
//            mat_big.data[x+1][i]=mat_big.data[x][i];
        mat_big.data[x+][x+]=;
        mat_big.data[x+][x]=;
        s.ini();
        s.len_y=;
        s.len_x=x+;
        s.data[][]=;
        s.data[x+][]=;
//        for(int i=0;i<=x+1;i++)
//            s.data[i][0]=1;
        matrix ret;
        ret.ini();
        ret.len_x=x+;
        ret.len_y=x+;
        for(int i=; i<ret.len_x; i++)
            ret.data[i][i]=;
        n++;
        while(n)
        {
            if(n&)ret=mat_mul(ret,mat_big);
            n>>=;
            mat_big=mat_mul(mat_big,mat_big);
        }
        matrix ans=mat_mul(ret,s);
        printf("%lld\n",ans.data[x+][]);
    }
    return ;
}

2.

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 60
#define ll long long
using namespace std;
int x;
ll M,n;
struct matrix
{
    int len_x;
    int len_y;
    ll data[maxn][maxn];
    void ini()
    {
        len_x=;
        len_y=;
        memset(data,,sizeof data);
    }
};

matrix mat_big;
ll f[maxn][maxn];
matrix s;

matrix mat_mul(matrix mat,matrix t)
{
    matrix ans;
    ans.ini();
    ans.len_x=mat.len_x;
    ans.len_y=t.len_y;
    for(int i=; i<mat.len_x; i++)
    {
        for(int j=; j<t.len_y; j++)
        {
            for(int k=; k<mat.len_y; k++)
            {
                ans.data[i][j]+=(mat.data[i][k]*t.data[k][j])%M;
                ans.data[i][j]%=M;
            }
        }
    }
    return ans;
}

int main()
{
    f[][]=;
    for(int i=; i<maxn; i++)
        for(int j=; j<=i; j++)
        {
            f[i][j]=f[i-][j]+f[i-][j-];
        }
    while(scanf("%lld%d%lld",&n,&x,&M)&&n!=-)
    {
        mat_big.ini();
        mat_big.len_x=x+;
        mat_big.len_y=x+;
        for(int i=; i<=x; i++)
            for(int j=; j<=i; j++)
            {
                mat_big.data[i][j]=x*f[i][j];
                if(mat_big.data[i][j]>M)mat_big.data[i][j]%=M;
            }
        for(int i=; i<=x; i++)
            mat_big.data[x+][i]=mat_big.data[x][i];
        mat_big.data[x+][x+]=;
//        mat_big.data[x+1][x]=1;
        s.ini();
        s.len_y=;
        s.len_x=x+;
        s.data[][]=;
        s.data[x+][]=;
//        for(int i=0;i<=x+1;i++)
//            s.data[i][0]=1;
        matrix ret;
        ret.ini();
        ret.len_x=x+;
        ret.len_y=x+;
        for(int i=; i<ret.len_x; i++)
            ret.data[i][i]=;
        while(n)
        {
            if(n&)ret=mat_mul(ret,mat_big);
            n>>=;
            mat_big=mat_mul(mat_big,mat_big);
        }
        matrix ans=mat_mul(ret,s);
        printf("%lld\n",ans.data[x+][]);
    }
    return ;
}