結果
問題 | No.2669 Generalized Hitting Set |
ユーザー |
|
提出日時 | 2024-01-14 02:49:54 |
言語 | PyPy3 (7.3.15) |
結果 |
RE
|
実行時間 | - |
コード長 | 1,688 bytes |
コンパイル時間 | 433 ms |
コンパイル使用メモリ | 82,460 KB |
実行使用メモリ | 203,500 KB |
最終ジャッジ日時 | 2024-09-28 02:41:22 |
合計ジャッジ時間 | 10,102 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | RE * 41 |
ソースコード
from typing import Listimport sysdef supset_zeta(f: List[int]):n = len(f)b = 1while b < n:l = 0while l < n:for p in range(l, l + b):f[p] += f[p + b]l += 2 * bb *= 2def subset_zeta(f: List[int]):n = len(f)b = 1while b < n:l = 0while l < n:for p in range(l, l + b):f[p + b] += f[p]l += 2 * bb *= 2N_MAX = 24popcnt = [0] * (1 << N_MAX)for t in range(1, 1 << N_MAX):popcnt[t] = 1 + popcnt[t & (t - 1)]def calc_f(n: int, m: int, k: int, s: List[int]):binom = [[0] * (n + 1) for _ in range(n + 1)]for i in range(n + 1):binom[i][0] = 1for j in range(1, i + 1):binom[i][j] = binom[i - 1][j - 1] + binom[i - 1][j]a = [0] * (n + 1)for l in range(k, n + 1):a[l] = 1for i in range(k, l):a[l] -= a[i] * binom[l][i]f = [0] * (1 << n)for u in s:f[u] += 1supset_zeta(f)for t in range(1 << n):f[t] *= a[popcnt[t]]subset_zeta(f)return fdef calc_g(n: int, m: int, f: List[int]):g = [n + 1] * (m + 1)for t in range(1 << n):g[f[t]] = min(g[f[t]], popcnt[t])for p in reversed(range(m)):g[p] = min(g[p], g[p + 1])return gif __name__ == '__main__':readline = sys.stdin.readlinen, m, k = map(int, readline.split())s = [0] * mfor i in range(m):si = readline()for j in range(n):s[i] |= (si[j] == '1') << jf = calc_f(n, m, k, s)g = calc_g(n, m, f)print('\n'.join(map(str, g[1:])))