我的思路:
dp[i][j][k]表示⾛到i,j位置,和是k,枚举下⼀个⾛啥,然后把
另⼀边的⼀⾏或者⼀列⼲掉。因为这个时间复杂度太高,所以只有40分
#include<bits/stdc++.h>
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define IL inline
#define ls (t<<1)
#define rs (t<<1|1)
using namespace std;
typedef long long LL;
typedef pair<int,int> PII;
IL LL read()
{
LL x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
return x*f;
}
#define io read()
const int N=108,S=108;
int n,m,s,mod;
int mi[N];
int f[N][N][S];
IL void add_ans(int &x,int y)
{
x+=y;
if(x>=mod) x-=mod;
}
int main()
{
n=io;m=io;s=io;mod=io;
mi[0]=1;
for(int i=1;i<N;++i){
mi[i]=1ll*mi[i-1]*(s+1)%mod;
}
f[1][1][0]=1;
for(int i=1;i<=n;++i){
for(int j=1;j<=m;++j){
for(int k=0;k<=s;++k){
for(int t=0;k+t<=s;++t){
if(i==n){
add_ans(f[i][j+1][k+t],f[i][j][k]);
}
else{
add_ans(f[i][j+1][k+t],1llu*f[i][j][k]*(t+1)%mod*mi[n-i-1]%mod);
}
if(j==m){
add_ans(f[i+1][j][k+t],f[i][j][k]);
}
else{
add_ans(f[i+1][j][k+t],1llu*f[i][j][k]*t%mod*mi[m-j-1]%mod);
}
}
}
}
}
/*for(int i=1;i<=n;++i){
for(int j=1;j<=m;++j){
for(int k=0;k<=s;++k){
printf("[%d][%d][%d]:%d\n",i,j,k,f[i][j][k]);
}
}
}*/
/*for(int i=1;i<n;++i){
for(int k=0;k<=s;++k){
for(int t=0;k+t<=s;++t){
add_ans(f[i+1][m][k+t],f[i][m][k]);
}
}
}
for(int j=1;j<m;++j){
for(int k=0;k<=s;++k){
for(int t=0;k+t<=s;++t){
add_ans(f[n][j+1][k+t],f[n][j][k]);
}
}
}*/
printf("%d\n",f[n][m][s]);
//cout<<clock()<<'\n';
return 0;
}
考虑这条路径什么时候碰到右下边界。
可以算出往右⾛了多少步,往下⾛了多少步。
往右⾛⼀步要求下⾯的格⼦<=当前格⼦,往下要求右边的
格⼦<当前格⼦,剩下的格⼦可以随便选。
往右和往下分别dp,然后枚举往下的和和往右的和合并。
时间复杂度O((n+m)s^2(还是满可靠的)
#include<bits/stdc++.h>
using namespace std;
#define N 105
#define M 2505
int n,m,s,p,res,f[N][N][N],pn[M*M],f1[M][N],f2[M][N],g[N],g1[N],fac[M*M],ifac[M*M];
int pw(int a,int b){int r=1;for(;b;b>>=1,a=1ll*a*a%p)if(b&1)r=1ll*r*a%p;return r;}
int C(int a,int b){return a<0||b<0||a<b?0:1ll*fac[a]*ifac[b]%p*ifac[a-b]%p;}
int main()
{
scanf("%d%d%d%d",&n,&m,&s,&p);
for(int i=pn[0]=1;i<=n*m;i++)
pn[i]=1ll*pn[i-1]*(s+1)%p;
if(n==1||m==1||n<=40)
{
f[1][1][0]=1;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
for(int k=0;k<=s;k++)
if(f[i][j][k])
for(int l=0;k+l<=s;l++)
{
if(i!=n&&j!=m)
{
int x=i+1,y=j;
f[x][y][k+l]=(f[x][y][k+l]+1ll*f[i][j][k]*l%p*pn[m-j-1])%p;
x=i,y=j+1;
f[x][y][k+l]=(f[x][y][k+l]+1ll*f[i][j][k]*(l+1)%p*pn[n-i-1])%p;
}
else if(i!=n)
{
int x=i+1,y=j;
f[x][y][k+l]=(f[x][y][k+l]+f[i][j][k])%p;
}
else if(j!=m)
{
int x=i,y=j+1;
f[x][y][k+l]=(f[x][y][k+l]+f[i][j][k])%p;
}
}
}
}
printf("%d\n",f[n][m][s]);
return 0;
}
for(int i=f1[0][0]=1;i<=n;i++)
for(int j=0;j<=s;j++)
for(int k=0;j+k<=s;k++)
f1[i][j+k]=(f1[i][j+k]+1ll*f1[i-1][j]*k%p)%p;
for(int i=f2[0][0]=1;i<=m;i++)
for(int j=0;j<=s;j++)
for(int k=0;j+k<=s;k++)
f2[i][j+k]=(f2[i][j+k]+1ll*f2[i-1][j]*(k+1)%p)%p;
int lim=M*M;
for(int i=fac[0]=1;i<lim;i++)
fac[i]=1ll*fac[i-1]*i%p;
ifac[lim-1]=pw(fac[lim-1],p-2);
for(int i=lim-1;i;i--)
ifac[i-1]=1ll*ifac[i]*i%p;
for(int i=1;i<n;i++)
{
for(int j=0;j<=s;j++)
g[j]=g1[j]=0;
int x=i-1,y=m-2,c=C(x+y,x);
for(int j=0;j<=s;j++)
for(int k=0;j+k<=s;k++)
g[j+k]=(g[j+k]+1ll*f1[x][j]*f2[y][k])%p;
for(int j=0;j<=s;j++)
for(int k=0;j+k<=s;k++)
g1[j+k]=(g1[j+k]+1ll*g[j]*(k+1))%p;
for(int j=0;j<=s;j++)
{
int a=n*m-i-m-x-y-(n-i);
int r=1ll*g1[j]*c%p*pn[a]%p*C(s-j+n-i-1,n-i-1)%p;
res=(res+r)%p;
}
}
for(int i=1;i<m;i++)
{
for(int j=0;j<=s;j++)
g[j]=g1[j]=0;
int x=n-2,y=i-1,c=C(x+y,x);
for(int j=0;j<=s;j++)
for(int k=0;j+k<=s;k++)
g[j+k]=(g[j+k]+1ll*f1[x][j]*f2[y][k])%p;
for(int j=0;j<=s;j++)
for(int k=0;j+k<=s;k++)
g1[j+k]=(g1[j+k]+1ll*g[j]*k)%p;
for(int j=0;j<=s;j++)
{
int a=n*m-n-i-x-y-(m-i);
int r=1ll*g1[j]*c%p*pn[a]%p*C(s-j+m-i-1,m-i-1)%p;
res=(res+r)%p;
}
}
printf("%d\n",res);
return 0;
}