結果
問題 |
No.2318 Phys Bone Maker
|
ユーザー |
|
提出日時 | 2023-05-27 01:07:36 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,288 bytes |
コンパイル時間 | 1,033 ms |
コンパイル使用メモリ | 82,344 KB |
実行使用メモリ | 79,996 KB |
最終ジャッジ日時 | 2024-12-25 12:38:57 |
合計ジャッジ時間 | 29,666 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 TLE * 1 |
other | AC * 44 TLE * 1 |
ソースコード
def make_divisors(n): lower_divisors , upper_divisors = [], [] i = 1 while i*i <= n: if n % i == 0: lower_divisors.append(i) if i != n // i: upper_divisors.append(n//i) i += 1 return lower_divisors + upper_divisors[::-1] def factorization(n): arr = dict() temp = n for i in range(2, int(-(-n**0.5//1))+1): if temp%i==0: cnt=0 while temp%i==0: cnt+=1 temp //= i arr[i]=cnt if temp!=1: arr[temp]=1 if len(arr)==0: arr[n]=1 return arr MOD = 998244353 N = int(input()) Y = make_divisors(N) F = [] for y in Y: F.append(factorization(y)) dp = [0 for _ in range(len(Y))] dp[0] = 1 for i in range(len((Y))): for j in range(i+1, len(Y)): #Y[i] -> Y[j] temp = dp[i] for p in F[i]: if p==1: continue if p not in F[j]: temp=0 break else: ej = F[j][p] ei = F[i][p] if ei>ej: temp=0 break elif ei==ej: temp=temp*(ej+1)%MOD dp[j] = (dp[j]+temp)%MOD print(dp[-1])