結果
| 問題 | 
                            No.1426 Got a Covered OR
                             | 
                    
| コンテスト | |
| ユーザー | 
                             nok0
                         | 
                    
| 提出日時 | 2021-01-17 11:11:24 | 
| 言語 | Python3  (3.13.1 + numpy 2.2.1 + scipy 1.14.1)  | 
                    
| 結果 | 
                             
                                WA
                                 
                             
                            
                         | 
                    
| 実行時間 | - | 
| コード長 | 1,473 bytes | 
| コンパイル時間 | 235 ms | 
| コンパイル使用メモリ | 12,800 KB | 
| 実行使用メモリ | 32,528 KB | 
| 最終ジャッジ日時 | 2024-11-29 10:36:03 | 
| 合計ジャッジ時間 | 8,059 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge1 / judge2 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | WA * 3 | 
| other | WA * 24 | 
ソースコード
import math
import sys
read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines
def popcnt(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
mod = 1_000_000_007
MAX = 10 ** 5 + 10
fact = [1, 1]
factinv = [1, 1]
inv = [0, 1]
for i in range(2, MAX + 1):
    fact.append((fact[-1] * i) % mod)
    inv.append((-inv[mod % i] * (mod // i)) % mod)
    factinv.append((factinv[-1] * inv[-1]) % mod)
def COM(n, r):
    if (r < 0) or (n < r):
        return 0
    return fact[n] * factinv[r] * factinv[n - r] % mod
n = int(readline())
b = list(map(int, readline().split()))
base, cont, res = 0, 0, 1
print(b)
for v in b:
    if v == -1:
        cont += 1
    else:
        if (v | (base ^ v)) != v:
            print(0)
            exit()
        a = popcnt(base)
        b = popcnt(v)
        x = b - a
        tmp = 0
        for i in range(0, x + 1):
            tmp += COM(x, i) * pow((pow(2, b - i, mod) - 1), cont + 1, mod) * (-1 if i % 2 else 1)
            tmp %= mod
        res *= tmp
        res %= mod
        cont, base = 0, v
print(res)
            
            
            
        
            
nok0