結果
| 問題 |
No.2082 AND OR XOR
|
| コンテスト | |
| ユーザー |
taiga0629kyopro
|
| 提出日時 | 2022-09-12 17:53:30 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 226 ms / 2,000 ms |
| コード長 | 3,148 bytes |
| コンパイル時間 | 249 ms |
| コンパイル使用メモリ | 82,428 KB |
| 実行使用メモリ | 82,108 KB |
| 最終ジャッジ日時 | 2024-11-29 20:11:16 |
| 合計ジャッジ時間 | 5,347 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 29 |
ソースコード
from _collections import defaultdict
mod=998244353
def sol1(n,a,b,c):
def op(x,y):
return (a*(x&y)+b*(x|y)+c*(x^y))%4
ans=0
for bit in range(4**n):
A=[]
for i in range(n):
A.append(bit%4)
bit//=4
res1=A[0]
res2=A[n-1]
for i in range(n-1):
res1=op(res1,A[i+1])
res2=op(res2,A[n-2-i])
if res1==res2:
ans+=1
else:
print(A,res1,res2)
return ans%mod
def sol2(n,a,b,c):
def op(x,y):
return (a*(x&y)+b*(x|y)+c*(x^y))%4
def mp(f,g):
ff=[]
gg=[]
for i in range(4):
ff.append(f%4)
gg.append(g%4)
f//=4
g//=4
res=[0]*4
for i in range(4):
res[i]=ff[gg[i]]
return res[0]+4*res[1]+16*res[2]+64*res[3]
def tof(f):
ff=[]
for i in range(4):
ff.append(f%4)
f//=4
return ff
func=[0]*4
for i in range(4):
fs=[0]*4
for x in range(4):fs[x]=op(x,i)
func[i]=fs[0]+4*fs[1]+16*fs[2]+64*fs[3]
dp=[defaultdict(int) for i in range(n+1)]
for i in range(4):
dp[1][func[i],func[i]]+=1
for i in range(1,n):
for f,g in dp[i]:
dp[i][f,g]%=mod
for j in range(4):
dp[i+1][mp(func[j],f),mp(g,func[j])]+=dp[i][f,g]
ans=0
for i in range(4):
for j in range(4):
for f,g in dp[n-2]:
res1=tof(f)[i]
res1=op(res1,j)
res2=tof(g)[j]
res2=op(res2,i)
if res1==res2:
ans+=dp[n-2][f,g]
return ans%mod
def sol3(n,a,b,c):
def op(x,y):
return (a*(x&y)+b*(x|y)+c*(x^y))%4
def mp(f,g):
ff=[]
gg=[]
for i in range(4):
ff.append(f%4)
gg.append(g%4)
f//=4
g//=4
res=[0]*4
for i in range(4):
res[i]=ff[gg[i]]
return res[0]+4*res[1]+16*res[2]+64*res[3]
def tof(f):
ff=[]
for i in range(4):
ff.append(f%4)
f//=4
return ff
func=[0]*4
for i in range(4):
fs=[0]*4
for x in range(4):fs[x]=op(x,i)
func[i]=fs[0]+4*fs[1]+16*fs[2]+64*fs[3]
dp=[defaultdict(int) for i in range(n+1)]
for i in range(4):
dp[1][i,func[i]]+=1
for i in range(1,n):
for x,F in dp[i]:
dp[i][x,F]%=mod
for j in range(4):
dp[i+1][op(x,j),mp(F,func[j])]+=dp[i][x,F]
ans=0
for an in range(4):
for x,F in dp[n-1]:
res1=op(x,an)
res2=tof(F)[an]
if res1==res2:ans+=dp[n-1][x,F]
return ans%mod
n,a,b,c=map(int,input().split())
print(sol3(n,a,b,c))
exit()
cnt=0
for n in range(N,N+1):
for a in range(4):
for b in range(4):
for c in range(4):
ans=sol2(n,a,b,c)
if ans!= 4**n and ans!=(4**n)//2:
cnt+=1
print(a,b,c,cnt)
print(ans,4**n)
taiga0629kyopro