#include<bits/stdc++.h>
using namespace std;

using ll=long long;
using pii=pair<int,int>;
#define all(a) a.begin(),a.end()
#define sz(a) ((int)a.size())
#define pb push_back

const int N=105;
int mod;

int add(int x, int y){
    x+=y; if(x>=mod) x-=mod;
    return x;
}
int sub(int x, int y){
    x-=y; if(x<0) x+=mod;
    return x;
}
int mul(int x, int y){
    return ((ll)x)*y%mod;
}
int Pow(int x, ll y){
    int res=1;
    for(; y; x=mul(x,x),y>>=1) if(y&1) res=mul(res,x);
    return res;
}

int n,m,C[N][N],dp[N][N*N*2],sum[N][N*N*2],dp2[N][N][N],sum2[N][N][N];

void build_dp(int mx){
    memset(dp,0,sizeof dp);
    memset(sum,0,sizeof sum);
    dp[0][0]=1;
    for(int i=0; i<N*N*2; ++i) sum[0][i]=1;
    for(int i=1; i<=n; ++i){
        for(int j=0; j<=mx*i; ++j){
            dp[i][j]=sum[i-1][j];
            if(j-mx-1>=0) dp[i][j]=sub(dp[i][j],sum[i-1][j-mx-1]);
        }
        sum[i][0]=dp[i][0];
        for(int j=1; j<=mx*(i+1); ++j){
            sum[i][j]=add(sum[i][j-1],dp[i][j]);
        }
    }
}

signed main(){
    ios_base::sync_with_stdio(0),cin.tie(0);
    cin >> n >> m >> mod;
    C[0][0]=1;
    for(int i=1; i<N; ++i){
        for(int j=0; j<=i; ++j){
            if(j==0||j==i) C[i][j]=1;
            else C[i][j]=add(C[i-1][j-1],C[i-1][j]);
        }
    }
    int res=Pow(2*m+1,n);
    build_dp(m*2);
    for(int mx=0; mx<=n; ++mx){
        dp2[mx][0][0]=1;
        for(int i=0; i<N; ++i) sum2[mx][0][i]=1;
        for(int i=1; i<=n; ++i){
            for(int j=0; j<=n; ++j){
                dp2[mx][i][j]=sum2[mx][i-1][j];
                if(j-mx-1>=0) dp2[mx][i][j]=sub(dp2[mx][i][j],sum2[mx][i-1][j-mx-1]);
            }
            sum2[mx][i][0]=dp2[mx][i][0];
            for(int j=1; j<=n; ++j){
                sum2[mx][i][j]=add(sum2[mx][i][j-1],dp2[mx][i][j]);
            }
        }
    }
    for(int mx=n; mx>=1; --mx){
        for(int i=0; i<N; ++i) for(int j=0; j<N; ++j) dp2[mx][i][j]=sub(dp2[mx][i][j],dp2[mx-1][i][j]);
    }
    res=sub(res,dp[n][n*m]);
    //cout << res << "\n";
    for(int x=1; x<=m; ++x){
        int y=x;
        while(y+1<=m&&m/(y+1)==m/x) y++;
        int tot=0;
        for(int mx=0; mx<=m/x; ++mx){
            for(int i=0; i<n; ++i){
                for(int j=0; j<=n; ++j) if(j<n-i){
                    if(j+mx+2>n-i){
                        tot=add(tot,mul(C[n][i],dp2[mx][i][j]));
                        //cout << mx << ' ' << i << ' ' << j << ' ' << dp2[mx][i][j] << "\n";
                    }
                }
            }
        }
        res=sub(res,mul(tot,2*(y-x+1)));
        x=y;
    }
    cout << res << "\n";
}