結果
問題 | No.2669 Generalized Hitting Set |
ユーザー |
|
提出日時 | 2024-01-16 21:06:51 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 3,152 ms / 4,000 ms |
コード長 | 1,144 bytes |
コンパイル時間 | 326 ms |
コンパイル使用メモリ | 81,588 KB |
実行使用メモリ | 215,036 KB |
最終ジャッジ日時 | 2024-09-28 02:58:54 |
合計ジャッジ時間 | 34,025 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 41 |
ソースコード
import sysreadline = sys.stdin.readlinepopcnt = [0] * (2 ** 12)for t in range(1, 2 ** 12):popcnt[t] = 1 + popcnt[t & (t - 1)]n, m, k = map(int, readline().split())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]f = [0] * (2 ** n)for i in range(m):# reversed, but it's oks = int(readline(), 2)f[s] += 1a = [0] * (n + 1)for l in range(k, n + 1):a[l] = (-1) ** (l - k) * binom[l - 1][l - k]# supset zetab = 1while b < 2 ** n:for l in range(0, 2 ** n, 2 * b):for p in range(l, l + b):f[p] += f[p + b]b *= 2for t in range(2 ** n):f[t] *= a[popcnt[t >> 12] + popcnt[t & 0xfff]]# subset zetab = 1while b < 2 ** n:for l in range(0, 2 ** n, 2 * b):for p in range(l, l + b):f[p + b] += f[p]b *= 2g = [n + 1] * (m + 1)for t in range(2 ** n):g[f[t]] = min(g[f[t]], popcnt[t >> 12] + popcnt[t & 0xfff])for p in reversed(range(m)):g[p] = min(g[p], g[p + 1])print('\n'.join(map(str, g[1:])))