結果
問題 | No.1529 Constant Lcm |
ユーザー | kohei2019 |
提出日時 | 2022-06-03 01:07:45 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,425 bytes |
コンパイル時間 | 193 ms |
コンパイル使用メモリ | 81,768 KB |
実行使用メモリ | 271,204 KB |
最終ジャッジ日時 | 2023-10-21 01:25:10 |
合計ジャッジ時間 | 32,728 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge9 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 40 ms
55,788 KB |
testcase_01 | AC | 74 ms
73,844 KB |
testcase_02 | AC | 44 ms
55,632 KB |
testcase_03 | AC | 40 ms
55,788 KB |
testcase_04 | AC | 40 ms
55,788 KB |
testcase_05 | AC | 40 ms
55,788 KB |
testcase_06 | AC | 41 ms
55,788 KB |
testcase_07 | AC | 41 ms
55,788 KB |
testcase_08 | AC | 41 ms
55,788 KB |
testcase_09 | AC | 40 ms
55,788 KB |
testcase_10 | AC | 2,454 ms
257,588 KB |
testcase_11 | AC | 981 ms
147,592 KB |
testcase_12 | AC | 105 ms
80,356 KB |
testcase_13 | AC | 2,175 ms
231,076 KB |
testcase_14 | AC | 1,280 ms
180,968 KB |
testcase_15 | AC | 1,383 ms
188,928 KB |
testcase_16 | AC | 1,064 ms
157,684 KB |
testcase_17 | AC | 592 ms
121,364 KB |
testcase_18 | AC | 627 ms
124,728 KB |
testcase_19 | AC | 1,418 ms
184,152 KB |
testcase_20 | AC | 2,975 ms
271,204 KB |
testcase_21 | AC | 2,960 ms
271,200 KB |
testcase_22 | AC | 2,781 ms
270,672 KB |
testcase_23 | AC | 2,828 ms
271,200 KB |
testcase_24 | TLE | - |
testcase_25 | AC | 2,995 ms
270,680 KB |
ソースコード
def primeset(N): #N以下の素数をsetで求める.エラトステネスの篩O(√Nlog(N)) lsx = [1]*(N+1) for i in range(2,int(-(-N**0.5//1))+1): if lsx[i] == 1: for j in range(i,N//i+1): lsx[j*i] = 0 setprime = set() for i in range(2,N+1): if lsx[i] == 1: setprime.add(i) return setprime def factorization_all_n(n):#n以下の自然数すべてをを素因数分解 lspn = [[] for i in range(n+1)] lsnum = [i for i in range(n+1)] lsp = list(primeset(n)) lsp.sort() for p in lsp: for j in range(1,n//p+1): cnt = 0 while lsnum[p*j]%p==0: lsnum[p*j] //= p cnt += 1 lspn[j*p].append((p,cnt)) return lspn def modPow(a,n,mod):#繰り返し二乗法 a**n % mod if n==0: return 1 if n==1: return a%mod if n & 1: return (a*modPow(a,n-1,mod)) % mod t = modPow(a,n>>1,mod) return (t*t)%mod import collections N = int(input()) mod = 998244353 keym = collections.Counter() lspn = factorization_all_n(N-1) for i in range(1,N): v = N-i a = collections.Counter() for p,cnt in lspn[i]: a[p] += cnt for p,cnt in lspn[v]: a[p] += cnt for p in a.keys(): keym[p] = max(keym[p],a[p]) ans = 1 for p,cnt in keym.items(): ans *= modPow(p, cnt, mod) ans %= mod print(ans)