P5664 Emiya 家今天的飯

Miku

這個題很顯然的可以從部分分推到正解

64
上去就是一個四維dp，dp[i][j][k][z]表示在第1行的時候第一行選了j個，第2行選了k個，第3行選了z個的

情況下的方案數，轉移用手就能推。

只是有個小細節

#include<iostream>
#include<cstdio> 
#include<algorithm>
#include<cstring>
#define int long long
using namespace std;
int n,m;
int ma[101][2001];
int mod=998244353;
int dp[101][41][41][41];
int ans; 
signed main(){
	scanf("%lld%lld",&n,&m);
	for(int i=1;i<=n;++i){
		for(int j=1;j<=m;++j){
			scanf("%lld",&ma[i][j]);
		}
	}
	dp[0][0][0][0]=1;
	//我是 小細節 
	for(int i=1;i<=n;++i){
		for(int j=0;j<=i;++j){
			for(int k=0;k<=i;++k){
				for(int z=0;z<=i;++z){
					dp[i][j][k][z]=(dp[i][j][k][z]+dp[i-1][j][k][z])%mod;
					if(j) dp[i][j][k][z]=(dp[i][j][k][z]+(dp[i-1][j-1][k][z]*ma[i][1])%mod)%mod;
					if(k) dp[i][j][k][z]=(dp[i][j][k][z]+(dp[i-1][j][k-1][z]*ma[i][2])%mod)%mod;
					if(z) dp[i][j][k][z]=(dp[i][j][k][z]+(dp[i-1][j][k][z-1])*ma[i][3]%mod)%mod; 
				}
			}
		}
	}
	for(int i=0;i<=n/2;++i){
		for(int j=0;j<=n/2;++j){
			for(int z=0;z<=n/2;++z){
				if((i+j)>=z&&(z+i)>=j&&(z+j)>=i){
					ans+=dp[n][i][j][z];
					ans%=mod;
				}
			}
		}
	}
	cout<<ans-1;
	//我也是 
	return 0;
}

64

要考慮更多的列，那何不壓縮一下，改成三維，前兩維一樣，但是第三維改成其餘行選了幾個

這樣的話，就要分開考慮每一列了。然後再運用一點點數學知識

顯然考慮合法方案很噁心，可是scz說說過“正難則反”%%%

求出所有解減去不合法的。

#include<iostream>
#include<cstdio> 
#include<algorithm>
#include<cstring>
#define int long long
using namespace std;
int n,m;
int ma[101][2001];
int mod=998244353;
int dp[101][51][80];
int sum[101];
int ans=1; 
int simex; 
signed main(){
	scanf("%lld%lld",&n,&m);
	for(int i=1;i<=n;++i){
		for(int j=1;j<=m;++j){
			scanf("%lld",&ma[i][j]);
			sum[i]+=ma[i][j];
			sum[i]%=mod;
		}
		//sum[i]%=mod;
		ans=(ans*(sum[i]+1)%mod)%mod;
	}
	//dp[0][0][0]=1;
	for(int q=1;q<=m;++q){
		memset(dp,0,sizeof(dp));
		dp[0][0][0]=1;
		for(int i=1;i<=n;++i){
			for(int j=0;j<=i;++j){
				for(int z=0;z<=i-j;++z){
					dp[i][j][z]=(dp[i][j][z]+dp[i-1][j][z])%mod;
					if(j){
						dp[i][j][z]=(dp[i][j][z]+(dp[i-1][j-1][z]*ma[i][q])%mod)%mod;
					} 
					if(z){
						dp[i][j][z]=(dp[i][j][z]+(dp[i-1][j][z-1]*(sum[i]-ma[i][q])%mod)%mod)%mod;
					}
				}
			}
		}
		for(int ii=1;ii<=n;++ii){
			for(int j=0;j<=n-ii;++j){
				if(ii>j){
					simex+=(dp[n][ii][j]);
					simex%=mod; 
				}
			}
			simex=(simex+mod)%mod;
		} 
	}
	cout<<(ans-simex-1+mod)%mod;
	return 0;
}

100

要想拿到這個，還要再壓一維。我們關心選了幾個嗎？不，正相反，這個差纔是關心的。

又少了一位，但是差可能是個負數，最小是-n，那麼就全部加上一個n。