結果
問題 | No.1289 RNG and OR |
ユーザー |
![]() |
提出日時 | 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)