結果
| 問題 |
No.1289 RNG and OR
|
| コンテスト | |
| ユーザー |
ayaoni
|
| 提出日時 | 2020-11-14 07:31:43 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 578 ms / 2,000 ms |
| コード長 | 1,067 bytes |
| コンパイル時間 | 257 ms |
| コンパイル使用メモリ | 82,168 KB |
| 実行使用メモリ | 138,104 KB |
| 最終ジャッジ日時 | 2024-07-22 22:42:22 |
| 合計ジャッジ時間 | 3,092 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 21 |
ソースコード
import sys
def I(): return int(sys.stdin.readline().rstrip())
def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))
def popcount(n):
c = (n & 0x5555555555555555) + ((n >> 1) & 0x5555555555555555)
c = (c & 0x3333333333333333) + ((c >> 2) & 0x3333333333333333)
c = (c & 0x0f0f0f0f0f0f0f0f) + ((c >> 4) & 0x0f0f0f0f0f0f0f0f)
c = (c & 0x00ff00ff00ff00ff) + ((c >> 8) & 0x00ff00ff00ff00ff)
c = (c & 0x0000ffff0000ffff) + ((c >> 16) & 0x0000ffff0000ffff)
c = (c & 0x00000000ffffffff) + ((c >> 32) & 0x00000000ffffffff)
return c
N = I()
A = LI()
s = sum(A)
mod = 998244353
dp = [[0]*(N+1) for _ in range(1<<N)]
for i in range(1<<N):
for j in range(N+1):
if j == 0:
dp[i][j] = A[i]
continue
dp[i][j] += dp[i][j-1]
if i&(1<<(j-1)):
dp[i][j] += dp[i&~(1<<(j-1))][j-1]
ans = 0
for i in range(1,1<<N):
x = popcount(i)
y = ((1<<N)-1)^i
a = s*pow(s-dp[y][-1],mod-2,mod)
if x % 2:
ans += a
else:
ans -= a
ans %= mod
print(ans)
ayaoni