[CSP-S2019] Emiya 家今天的饭

---

Day 6+7 DP（那一年的day2题1，完全被支配了┭┮﹏┭┮）

---

代码如下（示例）：

#include<bits/stdc++.h>
#define mod 998244353
using namespace std;
typedef long long ll;
const int maxn = 105, maxm = 2005;
int a[maxn][maxm],sum[maxn][maxm];
ll f[maxn][maxn*2],g[maxn][maxn];
int main()
{   int n,m;
	cin >>n>>m;
	for(int i=1;i<=n;i++)
	    for(int j=1;j<=m;j++)
	    {
	        cin>>a[i][j];
	        sum[i][0] = (sum[i][0]+a[i][j])%mod;
		}
    for(int i = 1; i<=n; i++)
        for(int j = 1; j<=m; j++)
            sum[i][j] = (sum[i][0]-a[i][j]+mod)%mod;//除去需要非法列的a值后剩余总和 
    ll ans = 0;
    for(int col = 1; col<=m; col++)
    {
        memset(f,0,sizeof(f));
        f[0][n] = 1;//求非法列 
        for(int i=1;i<=n;i++)
            for(int j = n-i; j<=n+i; j++) 
                f[i][j] = (f[i-1][j]+f[i-1][j-1]*a[i][col]%mod+f[i-1][j+1]*sum[i][col]%mod)%mod;
        for(int j = 1; j<=n; j++)
            ans = (ans+f[n][n+j])%mod;
	}
	g[0][0] = 1;
	for(int i = 1; i<=n; i++)
	    for(int j = 0; j<=n; j++) 
		    g[i][j] = (g[i-1][j]+(j>0?g[i-1][j-1]*sum[i][0]%mod:0))%mod;//总方案数 
    for(int j = 1; j<=n; j++)
	    ans = (ans-g[n][j]+mod)%mod;  
	cout << ans*(mod-1)%mod << endl;
	return 0;
}