#CSPS2019D. Emiya 家今天的饭

题目描述

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∑n​j=1∑m​ai,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

详见附加文件 meal4.in、meal4.ans。

样例 5

详见附加文件 meal5.in、meal4.ans。

数据范围与提示

测试点编号	�=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;
}