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

}

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

}