結果

問題 No.1304 あなたは基本が何か知っていますか？私は知っています．
ユーザー vwxyz
提出日時 2023-03-20 02:10:39
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
TLE  
実行時間 -
コード長 1,889 bytes
コンパイル時間 505 ms
コンパイル使用メモリ 12,160 KB
実行使用メモリ 34,976 KB
最終ジャッジ日時 2025-06-22 03:08:15
合計ジャッジ時間 6,272 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 2 TLE * 1 -- * 71
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
readline=sys.stdin.readline

def Hadamard(polynomial,n,mod=0,inverse=False):
    polynomial_=[x for x in polynomial]+[0]*((1<<n)-len(polynomial))
    for bit in range(n):
        for i in range(1<<n):
            ii=i^(1<<bit)
            if i>ii:
                continue
            polynomial_[i],polynomial_[ii]=polynomial_[i]+polynomial_[ii],polynomial_[i]-polynomial_[ii]
            if mod:
                polynomial_[i]%=mod
                polynomial_[ii]%=mod
    if inverse:
        if mod:
            inve_2=pow((mod+1)//2,n)
            for i in range(1<<n):
                polynomial_[i]*=inve_2
                polynomial_[i]%=mod
        else:
            pow_2=pow(2,n)
            for i in range(1<<n):
                polynomial_[i]/=pow_2
    return polynomial_

def XOR_Convolution(polynomial0,polynomial1,mod=0):
    n=(max(len(polynomial0),len(polynomial1))-1).bit_length()
    Hadamard_polynomial0=Hadamard(polynomial0,n,mod=mod)
    Hadamard_polynomial1=Hadamard(polynomial1,n,mod=mod)
    if mod:
        convolution=[x*y%mod for x,y in zip(Hadamard_polynomial0,Hadamard_polynomial1)]
    else:
        convolution=[x*y for x,y in zip(Hadamard_polynomial0,Hadamard_polynomial1)]
    convolution=Hadamard(convolution,n,mod=mod,inverse=True)
    return convolution

N,K,X,Y=map(int,readline().split())
Y+=1
A=list(map(int,readline().split()))
mod=998244353
dp=[[0]*(1<<10) for i in range(1<<10)]
for a in A:
    dp[a][a]=1
for n in range(1,N):
    prev=dp
    dp=[[0]*(1<<10) for i in range(1<<10)]
    cnt=[0]*(1<<10)
    for i in range(1<<10):
        for j in range(1<<10):
            cnt[j]+=prev[i][j]
            cnt[j]%=mod
    for a in A:
        for j in range(1<<10):
            dp[a][a^j]+=cnt[j]-prev[a][j]
            dp[a][a^j]%mod
ans=0
for i in range(1<<10):
    for j in range(X,Y):
        ans+=dp[i][j]
        ans%=mod
print(ans)
0