結果

問題 No.1751 Fortune Nim
ユーザー chineristACchineristAC
提出日時 2021-11-19 22:44:30
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 1,963 bytes
コンパイル時間 243 ms
コンパイル使用メモリ 82,256 KB
実行使用メモリ 109,744 KB
最終ジャッジ日時 2024-06-10 10:00:52
合計ジャッジ時間 4,908 ms
ジャッジサーバーID
(参考情報)
judge5 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 AC 60 ms
67,584 KB
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 AC 61 ms
67,456 KB
testcase_06 AC 59 ms
67,328 KB
testcase_07 AC 60 ms
67,584 KB
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 AC 105 ms
101,756 KB
testcase_17 AC 101 ms
96,896 KB
testcase_18 AC 108 ms
105,104 KB
testcase_19 AC 99 ms
99,584 KB
testcase_20 AC 114 ms
108,980 KB
testcase_21 AC 104 ms
98,560 KB
testcase_22 AC 114 ms
108,916 KB
testcase_23 AC 109 ms
104,268 KB
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 AC 114 ms
108,948 KB
testcase_29 AC 115 ms
108,816 KB
testcase_30 AC 112 ms
108,944 KB
testcase_31 AC 113 ms
108,792 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

N = 2*10**5
mod = 998244353
g1 = [1]*(N+1) # 元テーブル
g2 = [1]*(N+1) #逆元テーブル
inverse = [1]*(N+1) #逆元テーブル計算用テーブル

for i in range( 2, N + 1 ):
    g1[i]=( ( g1[i-1] * i ) % mod )
    inverse[i]=( ( -inverse[mod % i] * (mod//i) ) % mod )
    g2[i]=( (g2[i-1] * inverse[i]) % mod )
inverse[0]=0

def cmb(n, r, mod):
    if ( r<0 or r>n ):
        return 0
    r = min(r, n-r)
    return g1[n] * g2[r] * g2[n-r] % mod

def fwt(n,A):
    assert len(A) == 2**n
    for i in range(n):
        t = 2**i
        for j in range(2**n):
            if j&t==0:
                A[j] += A[j|t]
    return A

def ifwt(n,A):
    assert len(A) == 2**n
    for i in range(n):
        t = 2**i
        for j in range(2**n):
            if j&t==0:
                A[j] -= A[j|t]
    return A

inv = pow(1024,mod-2,mod)
 
def _fourier(f, inverse = False):
    f = f[:]
    n = (len(f) - 1).bit_length()
    for d in range(n):
        for U in range(1 << n):
            if not U >> d & 1:
                s, t = f[U], f[U | 1 << d]
                f[U], f[U | 1 << d] = (s + t)%mod, (s - t)%mod
    if inverse:
        f = [v *inv % mod for v in f]
    return f
 
def convolution(f, g):
    return _fourier([a * b  % mod for a, b in zip(_fourier(f), _fourier(g))], inverse = 1)

import sys,random,bisect
from collections import deque,defaultdict
from heapq import heapify,heappop,heappush
from itertools import permutations
from math import log,gcd

input = lambda :sys.stdin.readline()
mi = lambda :map(int,input().split())
li = lambda :list(mi())

N = int(input())
A = li()

S = sum(A)

G = 0
for a in A:
    G ^= a

assert G==0
if G==0:
    n = S
    res = (1-pow(inverse[3],n,mod)) * pow(2,mod-2,mod) % mod
    print(res)

else:
    n = 10**100
    for i in range(N):
        if A[i]^G <= A[i]:
            tmp = S - (A[i]-A[i]^G) + 1
            n = min(n,tmp)
    res = (pow(inverse[3],n,mod)+1) * pow(2,mod-2,mod) % mod
    print(res)
0