题目描述
Emiya 是个擅长做菜的高中生,他共掌握 �n 种烹饪方法,且会使用 �m 种主要食材做菜。为了方便叙述,我们对烹饪方法从 1∼�1∼n 编号,对主要食材从 1∼�1∼m 编号。
Emiya 做的每道菜都将使用恰好一种烹饪方法与恰好一种主要食材。更具体地,Emiya 会做 ��,�ai,j 道不同的使用烹饪方法 �i 和主要食材 �j 的菜 (1≤�≤�,1≤�≤�)(1≤i≤n,1≤j≤m),这也意味着 Emiya 总共会做 ∑�=1�∑�=1���,�i=1∑nj=1∑mai,j 道不同的菜。
Emiya 今天要准备一桌饭招待 Yazid 和 Rin 这对好朋友,然而三个人对菜的搭配有不同的要求,更具体地,对于一种包含 �k 道菜的搭配方案而言:
-
Emiya 不会让大家饿肚子,所以将做至少一道菜,即 �≥1k≥1;
-
Rin 希望品尝不同烹饪方法做出的菜,因此她要求每道菜的烹饪方法互不相同;
-
Yazid 不希望品尝太多同一食材做出的菜,因此他要求每种主要食材至多在一半的菜(即 ⌊�2⌋⌊2k⌋ 道菜)中被使用(这里的 ⌊�⌋⌊x⌋ 为下取整函数,表示不超过 �x 的最大整数)。
这些要求难不倒 Emiya,但他想知道共有多少种不同的符合要求的搭配方案。两种方案不同,当且仅当存在至少一道菜在一种方案中出现,而不在另一种方案中出现。
Emiya 找到了你,请你帮他计算,你只需要告诉他符合所有要求的搭配方案数对质数 998244353998244353 取模的结果。
输入格式
从文件 meal.in
中读入数据。
第 11 行两个用单个空格隔开的整数 �,�n,m。
第 22 行至第 �+1n+1 行,每行 �m 个用单个空格隔开的整数,其中第 �+1i+1 行的 �m 个数依次为 ��,1,��,2,…,��,�ai,1,ai,2,…,ai,m。
输出格式
输出到文件 meal.out
中。
仅一行一个整数,表示所求方案数对 998244353998244353 取模的结果。
输入数据 1
2 3
1 0 1
0 1 1
输出数据 1
3
样例说明 1
由于在这个样例中,对于每组 �,�i,j,Emiya 都最多只会做一道菜,因此我们直接通过给出烹饪方法、主要食材的编号来描述一道菜。
符合要求的方案包括:
- 做一道用烹饪方法 11、主要食材 11 的菜和一道用烹饪方法 22、主要食材 22 的菜
- 做一道用烹饪方法 11、主要食材 11 的菜和一道用烹饪方法 22、主要食材 33 的菜
- 做一道用烹饪方法 11、主要食材 33 的菜和一道用烹饪方法 22、主要食材 22 的菜
因此输出结果为 3 mod 998244353=33mod998244353=3。 需要注意的是,所有只包含一道菜的方案都是不符合要求的,因为唯一的主要食材在超过一半的菜中出现,这不满足 Yazid 的要求。
输入数据 2
3 3
1 2 3
4 5 0
6 0 0
输出数据 2
190
样例说明 2
Emiya 必须至少做 22 道菜。
做 22 道菜的符合要求的方案数为 100100。
做 33 道菜的符合要求的方案数为 9090。
因此符合要求的方案数为 100+90=190100+90=190。
输入数据 3
5 5
1 0 0 1 1
0 1 0 1 0
1 1 1 1 0
1 0 1 0 1
0 1 1 0 1
输出数据 3
742
样例 4
样例 5
数据范围与提示
测试点编号 | �=n= | �=m= | ��,�<ai,j< |
---|---|---|---|
11 | 22 | 22 | 22 |
22 | 33 | ||
33 | 55 | 22 | |
44 | 33 | ||
55 | 1010 | 22 | |
66 | 33 | ||
77 | 1010 | 22 | 103103 |
88 | 33 | ||
9∼129∼12 | 4040 | 22 | |
13∼1613∼16 | 33 | ||
17∼2117∼21 | 500500 | ||
22∼2522∼25 | 100100 | 2×1032×103 | 998244353998244353 |
对于所有测试点,保证 1≤�≤1001≤n≤100,1≤�≤2×1031≤m≤2×103,0≤��,�<9982443530≤ai,j<998244353。
代码示例
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod = 998244353;
int n, m;
ll a[110][2010];
ll sum[110];
ll f[110][2010];
ll ans;
ll read() {
ll xx = 0; char ch = getchar();
while (ch < '0'|| ch > '9') ch = getchar();
while (ch >= '0' && ch <= '9') xx = (xx << 1) + (xx << 3) + ch - '0', ch = getchar();
return xx;
}
inline int h(int x) { return x + n + 5; }
int main() {
n = read(); m = read();
for (int i = 1; i <= n; i++)
for (int j = 1; j <= m; j++) {
a[i][j] = read();
sum[i] = (sum[i] + a[i][j]) % mod;
}
ans = 1;
for (int i = 1; i <= n; i++)
ans = (ans * (sum[i] + 1)) % mod;
ans--;
for (int j = 1; j <= m; j++) {
memset(f, 0, sizeof f);
f[1][h(0)] = 1;
f[1][h(1)] = a[1][j];
f[1][h(-1)] = (sum[1] + mod - a[1][j]) % mod;
for (int i = 2; i <= n; i++)
for (int k = -n; k <= n; k++)
f[i][h(k)] = ( (f[i - 1][h(k)] + a[i][j] * f[i - 1][h(k - 1)])
+ (f[i - 1][h(k + 1)] * (sum[i] + mod - a[i][j])) ) % mod;
for (int k = 1; k <= n; k++)
ans = (ans + mod - f[n][h(k)]) % mod;
}
printf("%lld\n", ans);
return 0;
}