結果
| 問題 |
No.1426 Got a Covered OR
|
| コンテスト | |
| ユーザー |
👑 SPD_9X2
|
| 提出日時 | 2021-03-12 22:57:39 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 206 ms / 2,000 ms |
| コード長 | 1,961 bytes |
| コンパイル時間 | 153 ms |
| コンパイル使用メモリ | 82,304 KB |
| 実行使用メモリ | 108,672 KB |
| 最終ジャッジ日時 | 2024-10-14 13:31:56 |
| 合計ジャッジ時間 | 3,953 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 24 |
ソースコード
"""
累積ORがB
Bが最後に出たやつに包含されていない場合は0
X -> Y の区間について考える
長さをLとする
Xで既にあるbitは、自由にしてok
2^L
新たなbitは、1回登場すればおk
(2^L-1)
Aが正なのが嫌だ
0がある場合を省かなくてはいけない
X箇所0がある場合…
みたいに包除すればよい
"""
import sys
from sys import stdin
def modfac(n, MOD):
f = 1
factorials = [1]
for m in range(1, n + 1):
f *= m
f %= MOD
factorials.append(f)
inv = pow(f, MOD - 2, MOD)
invs = [1] * (n + 1)
invs[n] = inv
for m in range(n, 1, -1):
inv *= m
inv %= MOD
invs[m - 1] = inv
return factorials, invs
def modnCr(n,r,mod,fac,inv): #上で求めたfacとinvsを引数に入れるべし(上の関数で与えたnが計算できる最大のnになる)
if r == 0:
return 1
return fac[n] * inv[n-r] * inv[r] % mod
mod = 10**9+7
fac,inv = modfac(200000,mod)
def popcnt(X):
ret = 0
while X > 0:
if X % 2 == 1:
ret += 1
X //= 2
return ret
N = int(stdin.readline())
B = list(map(int,stdin.readline().split()))
ans = 1
last = 0
for i in B:
if i != -1:
if last | i != i:
print ("0")
sys.exit()
last = i
lastind = -1
last = 0
for i in range(N):
if B[i] != -1:
ready = popcnt(last)
new = popcnt(B[i] ^ last)
L = i-lastind
tans = 0
for znum in range(L+1):
if znum == L and new > 0:
continue
now = pow(pow(2,L-znum,mod),ready,mod) * pow(pow(2,L-znum,mod)-1,new,mod)
now *= modnCr(L,znum,mod,fac,inv)
if znum % 2 == 0:
tans += now
else:
tans -= now
ans *= tans
ans %= mod
lastind = i
last = B[i]
#print (i,ans)
print (ans)
SPD_9X2