結果
問題 |
No.2318 Phys Bone Maker
|
ユーザー |
![]() |
提出日時 | 2023-05-26 22:37:34 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,367 bytes |
コンパイル時間 | 204 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 201,616 KB |
最終ジャッジ日時 | 2024-12-25 08:49:50 |
合計ジャッジ時間 | 45,570 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 TLE * 1 |
other | AC * 37 TLE * 8 |
ソースコード
class PrimeNumbers: def __init__(self,nmax): rootnmax = int(nmax**0.5) self.prime_judgement = [True]*(rootnmax+1) self.prime_judgement[0] = self.prime_judgement[1] = False for i in range(2,rootnmax+1): if self.prime_judgement[i]: for j in range(2,rootnmax//i+1): self.prime_judgement[i*j] = False self.prime_list = [] for i,flag in enumerate(self.prime_judgement): if flag: self.prime_list.append(i) def prime_factorization(self,n): return_list = [] for i in self.prime_list: if n==1 or i*i>n: break if n%i==0: return_list.append([i,0]) while n%i==0: return_list[-1][1] += 1; n //= i if n!=1: return_list.append([n,1]) return return_list def divisors_enumeration(self,n): divisors = [1] for p,v in self.prime_factorization(n): for i in range(len(divisors)): for j in range(v): divisors.append(divisors[i]*p**(j+1)) return set(divisors) from math import gcd n = int(input()) mod = 998244353 pn = PrimeNumbers(n) dn = pn.divisors_enumeration(n) dp = {v:1 for v in dn} for u in sorted(dn-{1}): for v in dn-pn.divisors_enumeration(u): l = u*v//gcd(u,v) dp[l] = (dp[l]+dp[u])%mod print(dp[n])