結果
問題 | No.2578 Jewelry Store |
ユーザー | 👑 p-adic |
提出日時 | 2023-12-07 15:26:14 |
言語 | PyPy3 (7.3.15) |
結果 |
RE
|
実行時間 | - |
コード長 | 828 bytes |
コンパイル時間 | 474 ms |
コンパイル使用メモリ | 82,916 KB |
実行使用メモリ | 83,016 KB |
最終ジャッジ日時 | 2024-09-27 02:11:31 |
合計ジャッジ時間 | 9,548 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | RE | - |
testcase_01 | RE | - |
testcase_02 | RE | - |
testcase_03 | RE | - |
testcase_04 | RE | - |
testcase_05 | RE | - |
testcase_06 | RE | - |
testcase_07 | RE | - |
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 | RE | - |
testcase_17 | RE | - |
testcase_18 | RE | - |
testcase_19 | RE | - |
testcase_20 | RE | - |
testcase_21 | RE | - |
testcase_22 | RE | - |
testcase_23 | RE | - |
testcase_24 | RE | - |
testcase_25 | RE | - |
testcase_26 | RE | - |
testcase_27 | RE | - |
testcase_28 | RE | - |
testcase_29 | RE | - |
testcase_30 | RE | - |
testcase_31 | RE | - |
testcase_32 | RE | - |
testcase_33 | RE | - |
testcase_34 | RE | - |
testcase_35 | RE | - |
testcase_36 | RE | - |
testcase_37 | TLE | - |
testcase_38 | -- | - |
testcase_39 | -- | - |
testcase_40 | -- | - |
testcase_41 | -- | - |
testcase_42 | -- | - |
testcase_43 | -- | - |
testcase_44 | -- | - |
testcase_45 | -- | - |
testcase_46 | -- | - |
testcase_47 | -- | - |
testcase_48 | -- | - |
testcase_49 | -- | - |
testcase_50 | -- | - |
testcase_51 | -- | - |
testcase_52 | -- | - |
testcase_53 | -- | - |
ソースコード
R=range I=lambda:map(int,input().split()) T,m=I() M=[m] P=set() 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() if n<99: i=2 while i*i<=n: if n%i<1 and not i in P:r(i) while n%i<1:n//=i i+=1 if n>1:r(n) else: for t in R(99): x,d,c=t,1,0 y=(x*x+1)%n while d<2 and c*c<n:x,y,d,c=(x*x+1)%n,(y**4+2*y*y+2)%n,G(abs(x-y),n),c+1 if max(d,c*c)<n: M+=[n//d,d] break else:r(n) assert(0) 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)