結果
問題 | No.1731 Product of Subsequence |
ユーザー |
👑 ![]() |
提出日時 | 2021-11-05 22:05:43 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,168 bytes |
コンパイル時間 | 292 ms |
コンパイル使用メモリ | 82,396 KB |
実行使用メモリ | 87,812 KB |
最終ジャッジ日時 | 2024-11-06 12:56:40 |
合計ジャッジ時間 | 22,055 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 28 TLE * 3 |
ソースコード
#素因数分解 def Prime_Factorization(N): if N<0: R=[[-1,1]] else: R=[] N=abs(N) if N&1==0: C=0 while N&1==0: N>>=1 C+=1 R.append([2,C]) if N%3==0: C=0 while N%3==0: N//=3 C+=1 R.append([3,C]) k=5 Flag=0 while k*k<=N: if N%k==0: C=0 while N%k==0: C+=1 N//=k R.append([k,C]) k+=2+2*Flag Flag^=1 if N!=1: R.append([N,1]) return R #================================================== N,K=map(int,input().split()) A=list(map(int,input().split())) Mod=10**9+7 if K==1: exit(print((pow(2,N,Mod)-1)%Mod)) P=[(p,e) for p,e in Prime_Factorization(K)] m=len(P) DP={(0,)*m:1} v=[0]*m for a in A: E=DP.copy() f=[0]*m for i in range(m): p,_=P[i] while a%p==0: f[i]+=1 a//=p for u in E: for i in range(m): v[i]=min(u[i]+f[i],P[i][1]) DP[tuple(v)]=(DP.get(tuple(v),0)+E[u])%Mod end=tuple([e for p,e in P]) print(DP.get(end,0))