agc039F Min Product Sum

agc039F Min Product Sum

Solution

• 考虑行和列分别按照最小值排序，那么$A_{i,j}$应该是这样的：

  111111111

  111111111

  111122222

  111122333

• 为若干个$L$型。每个位置乘在一起就是贡献。

• 考虑方案数为这样的：

  111122333

  111122333

  222222333

  333333333

  表示的是那个位置填的数要>=上面的矩阵中那个位置的数。

• 考虑按照$1-K$从左上角往下加入一行一列，但是并不能满足恰好的条件，所以需要容斥，也就是新加入的$t$有$a$行$b$列不能以$t$为最小值，容斥系数$(-1{)}^{a+b}$。

• 先加入行再加入列需要多记录一些东西，不是很好。

• 考虑先加入不容斥的行，不容斥的列，再先后加入容斥的行、容斥的列，转移起来就没有问题了。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define maxn 105
#define ll long long 
using namespace std;

int n,m,c,mo,i,j,k;
ll f[2][maxn][maxn];
ll mul[maxn][maxn*maxn],C[maxn][maxn];

ll ksm(ll x,ll y){
	ll s=1;
	for(;y;y/=2,x=x*x%mo) if (y&1)
		s=s*x%mo;
	return s;
}

int main(){
	freopen("ceshi.in","r",stdin);
	scanf("%d%d%d%d",&n,&m,&c,&mo);
	for(mul[0][0]=1,i=1;i<=c;i++){
		mul[i][0]=1;
		for(j=1;j<=n*m;j++) mul[i][j]=mul[i][j-1]*i%mo;
	}
	for(C[0][0]=1,i=1;i<=max(n,m);i++) for(C[i][0]=1,j=1;j<=i;j++)
		C[i][j]=(C[i-1][j]+C[i-1][j-1])%mo;
	int p=0,q=1;
	for(int x=1;x<=c;x++){
		p^=1,q^=1; memset(f[p],0,sizeof(f[p]));
		for(i=0;i<=max(n,m);i++) f[q][0][i]=f[q][i][0]=0; f[q][0][0]=1;
		for(i=n;i>=0;i--) for(j=m;j>=0;j--) if (f[q][i][j]) for(k=0;k<=n-i;k++)
			(f[p][i+k][j]+=f[q][i][j]*mul[x][(m-j)*k]%mo*mul[c-x+1][j*k]%mo*C[i+k][k])%=mo;
		p^=1,q^=1; memset(f[p],0,sizeof(f[p]));
		for(i=n;i>=0;i--) for(j=m;j>=0;j--) if (f[q][i][j]) for(k=0;k<=m-j;k++)
			(f[p][i][j+k]+=f[q][i][j]*mul[x][(n-i)*k]%mo*mul[c-x+1][i*k]%mo*C[j+k][k])%=mo;
		p^=1,q^=1; memset(f[p],0,sizeof(f[p]));
		for(i=n;i>=0;i--) for(j=m;j>=0;j--) if (f[q][i][j]) for(k=0;k<=n-i;k++)
			(f[p][i+k][j]+=f[q][i][j]*mul[x][(m-j)*k]%mo*mul[c-x][j*k]%mo*C[i+k][k]*((k&1)?-1:1))%=mo;
		p^=1,q^=1; memset(f[p],0,sizeof(f[p]));
		for(i=n;i>=0;i--) for(j=m;j>=0;j--) if (f[q][i][j]) for(k=0;k<=m-j;k++)
			(f[p][i][j+k]+=f[q][i][j]*mul[x][(n-i)*k]%mo*mul[c-x][i*k]%mo*C[j+k][k]*((k&1)?-1:1))%=mo;
	}
	printf("%lld\n",(f[p][n][m]%mo+mo)%mo);
}