結果
| 問題 |
No.2669 Generalized Hitting Set
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-04-02 20:26:38 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 2,536 ms / 4,000 ms |
| コード長 | 1,304 bytes |
| コンパイル時間 | 217 ms |
| コンパイル使用メモリ | 82,352 KB |
| 実行使用メモリ | 217,204 KB |
| 最終ジャッジ日時 | 2024-09-30 23:18:44 |
| 合計ジャッジ時間 | 28,818 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 41 |
ソースコード
import sys
input = sys.stdin.readline
def popcount32(x):
x = x - ((x >> 1) & 0x55555555)
x = (x & 0x33333333) + ((x >> 2) & 0x33333333)
x = (x + (x >> 4)) & 0x0F0F0F0F
x += x >> 8
x += x >> 16
return x & 0x0000003F
n, m, k = map(int, input().split())
dp = [0] * (1 << n)
for _ in range(m):
S = int(input(), 2)
dp[S] += 1
nCk = [[0] * (n + 1) for _ in range(n + 1)]
nCk[0][0] = 1
for i in range(1, n + 1):
nCk[i][0] = 1
for j in range(1, i + 1):
nCk[i][j] = nCk[i - 1][j - 1] + nCk[i - 1][j]
C = [0] * (n + 1)
C[k] = 1
for i in range(k + 1, n + 1):
C[i] = 1
for j in range(i):
C[i] -= C[j] * nCk[i][j]
for i in range(n):
dbit = 1 << i
for s in range(0, 1 << n, 2 * dbit):
for j in range(s, s + dbit):
dp[j] += dp[dbit + j]
for bit in range(1 << n):
p = popcount32(bit)
dp[bit] *= C[p]
for i in range(n):
dbit2 = 1 << i
for s in range(0, 1 << n, 2 * dbit2):
for j in range(s, s + dbit2):
dp[dbit2 + j] += dp[j]
ans = [n] * m
for bit in range(1 << n):
c = dp[bit] - 1
if c >= 0:
p = popcount32(bit)
ans[c] = min(ans[c], p)
else:
pass
for i in range(m - 2, -1, -1):
ans[i] = min(ans[i], ans[i + 1])
print(*ans, sep="\n")