結果
| 問題 | No.1529 Constant Lcm |
| コンテスト | |
| ユーザー |
ygd.
|
| 提出日時 | 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()
ygd.