题目信息
解题思路:dp
-
对于一种做菜方案,至多只有一种主要食材在超过[k/2]道菜里出现过。所以我们可以用排除法,用总方案数 - 不合法方案数即可。
-
我们需要枚举主要食材,dp[i][j][k]代表前i种方式中包含j种食材,有k种烹饪方式未使用的方案数。
-
对于任意x,若x > [k / 2],则2x + n - k > n。
-
因此,可以将使用过某种食材的菜看作使用两次此食材,并为每种烹饪方式加一种每种食材都用过一次的菜。
-
用过某食材 > n次 => 不合法方案。
代码实现
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod = 998244353;
const int N = 105, M = 2005;
ll dp[N][2 * N];
ll a[N][M], b[N];
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];
b[i] = (b[i] + a[i][j]) % mod;
}
}
ll res = 1;
for (int i = 1; i <= n; ++i) {
res = res * (b[i] + 1) % mod;
}
res = (res - 1 + mod) % mod;
for (int i = 1; i <= m; ++i) {
memset(dp, 0, sizeof dp);
dp[0][0] = 1;
for (int j = 1; j <= n; ++j) {
for (int k = 0; k <= (j - 1) << 1; ++k) {
dp[j][k] = (dp[j][k] + dp[j - 1][k] * (b[j] - a[j][i] + mod)) % mod;
dp[j][k + 1] = (dp[j][k + 1] + dp[j - 1][k]) % mod;
dp[j][k + 2] = (dp[j][k + 2] + dp[j - 1][k] * a[j][i]) % mod;
}
}
for (int j = n + 1; j <= n << 1; ++j) {
res = (res - dp[n][j] + mod) % mod;
}
}
cout << res << '\n';
return 0;
}