結果
問題 |
No.1529 Constant Lcm
|
ユーザー |
![]() |
提出日時 | 2021-06-06 13:29:48 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 2,857 ms / 3,000 ms |
コード長 | 1,955 bytes |
コンパイル時間 | 278 ms |
コンパイル使用メモリ | 82,432 KB |
実行使用メモリ | 306,924 KB |
最終ジャッジ日時 | 2024-11-22 20:02:22 |
合計ジャッジ時間 | 29,660 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 24 |
ソースコード
import math from collections import defaultdict def get_sieve_of_eratosthenes(n): if not isinstance(n, int): raise TypeError('n is int type.') if n < 2: raise ValueError('n is more than 2') #nが1の場合はエラーを返す。 prime = [2] if n == 2: return prime limit = int(n**0.5) #ルートn以上で割れることはない data = [i + 1 for i in range(2, n, 2)] #最初から3スタートの奇数にしておく while True: p = data[0] if limit < p: #ここでイコールが入るとlimitが素数の時にlimit^2が残ってダメ。 return prime + data prime.append(p) #ここで一番最初に取り出した素数を加える data = [e for e in data if e % p != 0] #残りでその素数で割れない数字を取得。 def main(): n = int(input()); MOD = 998244353 if n == 3: print(2);exit() if n == 5: print(12);exit() Prime = get_sieve_of_eratosthenes(n) m = n//2 mx = m*2 #数列の長さ*2 L = [defaultdict(int) for _ in range(m)] ans_dic = defaultdict(int) #print(Prime) for x in Prime: val = x #cnt = 1 while val < n: idx = val #print(idx) while idx < mx: if idx <= m: loc = idx - 1 else: loc = n - 1 - idx #print("x",x,loc,idx) L[loc][x] += 1 if n%2 == 0 and loc == m - 1: L[loc][x] += 1 ans_dic[x] = max(ans_dic[x], L[loc][x]) idx += val #print(idx) #print(L) val *= x #cnt += 1 #print(ans_dic) ans = 1 for x in ans_dic.keys(): val = ans_dic[x] ans *= pow(x,val,MOD) ans %= MOD print(ans) if __name__ == '__main__': main()