数数题

我的思路：
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;
}