結果
問題 | No.2578 Jewelry Store |
ユーザー | 👑 p-adic |
提出日時 | 2023-12-07 16:08:27 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 795 bytes |
コンパイル時間 | 474 ms |
コンパイル使用メモリ | 82,196 KB |
実行使用メモリ | 197,560 KB |
最終ジャッジ日時 | 2024-09-27 02:12:14 |
合計ジャッジ時間 | 26,148 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | WA | - |
testcase_01 | WA | - |
testcase_02 | WA | - |
testcase_03 | WA | - |
testcase_04 | WA | - |
testcase_05 | WA | - |
testcase_06 | WA | - |
testcase_07 | WA | - |
testcase_08 | WA | - |
testcase_09 | WA | - |
testcase_10 | WA | - |
testcase_11 | WA | - |
testcase_12 | WA | - |
testcase_13 | WA | - |
testcase_14 | WA | - |
testcase_15 | WA | - |
testcase_16 | WA | - |
testcase_17 | WA | - |
testcase_18 | WA | - |
testcase_19 | WA | - |
testcase_20 | WA | - |
testcase_21 | WA | - |
testcase_22 | WA | - |
testcase_23 | WA | - |
testcase_24 | WA | - |
testcase_25 | WA | - |
testcase_26 | WA | - |
testcase_27 | WA | - |
testcase_28 | WA | - |
testcase_29 | WA | - |
testcase_30 | WA | - |
testcase_31 | WA | - |
testcase_32 | WA | - |
testcase_33 | WA | - |
testcase_34 | WA | - |
testcase_35 | WA | - |
testcase_36 | WA | - |
testcase_37 | AC | 204 ms
77,092 KB |
testcase_38 | WA | - |
testcase_39 | WA | - |
testcase_40 | WA | - |
testcase_41 | WA | - |
testcase_42 | AC | 217 ms
77,492 KB |
testcase_43 | WA | - |
testcase_44 | AC | 228 ms
76,984 KB |
testcase_45 | WA | - |
testcase_46 | WA | - |
testcase_47 | TLE | - |
testcase_48 | -- | - |
testcase_49 | -- | - |
testcase_50 | -- | - |
testcase_51 | -- | - |
testcase_52 | -- | - |
testcase_53 | -- | - |
ソースコード
R=range I=lambda:map(int,input().split()) T,m=I() P=set() for i in R(2,min(m+1,31623)): if m%i<1:P.add(i) while m%i<1:m//=i M=[m] def G(a,b): while a:a,b=b%a,a return b def r(n): for j in R(len(M)): while M[j]%n<1:M[j]//=n P.add(n) while M: n=M.pop() c=int(n**0.25) if n<10**10:r(n) else: while c: x,d=c,1 y=(x*x+1)%n while d<2 and c:x,y,d,c=(x*x+1)%n,(y**4+2*y*y+2)%n,G(abs(x-y),n),c-1 if d<n and c: M+=[n//d,d] break else:r(n) P=list(P) L=len(P) Q=998244353 U=1<<L for t in R(T): N,B,C,D=I() c=[1]*U for a in I(): b=0 for i in R(L):b|=(m//a%P[i]<1)<<i c[b],B=c[b]*(1+B*(m%a<1))%Q,(C*B+D)%Q for i in R(L): for d in R(U): if d>>i&1<1:c[d]=c[d]*c[d|1<<i]%Q v=0 for d in R(U):v+=c[d]*[1,-1][len([1for i in R(L)if d>>i&1])%2] print((v-(m<2))%Q)