結果
問題 | No.1973 Divisor Sequence |
ユーザー |
👑 ![]() |
提出日時 | 2022-06-10 21:45:20 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 203 ms / 2,000 ms |
コード長 | 1,029 bytes |
コンパイル時間 | 230 ms |
コンパイル使用メモリ | 82,212 KB |
実行使用メモリ | 76,112 KB |
最終ジャッジ日時 | 2024-09-21 07:28:11 |
合計ジャッジ時間 | 2,785 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
#素因数分解def Prime_Factorization(N):if N==0:return [[0,1]]if N<0:R=[[-1,1]]else:R=[]N=abs(N)if N&1==0:C=0while N&1==0:N>>=1C+=1R.append([2,C])if N%3==0:C=0while N%3==0:N//=3C+=1R.append([3,C])k=5Flag=0while k*k<=N:if N%k==0:C=0while N%k==0:C+=1N//=kR.append([k,C])k+=2+2*FlagFlag^=1if N!=1:R.append([N,1])return R#==================================================N,M=map(int,input().split())F=Prime_Factorization(M)X=1Mod=10**9+7D={}for p,e in F:if e not in D:DP=[0]*(e+1); DP[0]=1for n in range(1,N+1):for i in range(1,e+1):DP[i]+=DP[i-1]DP[i]%=ModDP.reverse()DP[e]=sum(DP)%ModX*=DP[e]X%=Modprint(X)