結果
問題 | No.1529 Constant Lcm |
ユーザー | kohei2019 |
提出日時 | 2022-06-03 01:07:45 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,425 bytes |
コンパイル時間 | 305 ms |
コンパイル使用メモリ | 82,576 KB |
実行使用メモリ | 272,012 KB |
最終ジャッジ日時 | 2024-09-21 02:15:11 |
合計ジャッジ時間 | 34,068 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 43 ms
55,312 KB |
testcase_01 | AC | 78 ms
74,628 KB |
testcase_02 | AC | 42 ms
55,304 KB |
testcase_03 | AC | 41 ms
54,192 KB |
testcase_04 | AC | 42 ms
54,508 KB |
testcase_05 | AC | 42 ms
54,772 KB |
testcase_06 | AC | 43 ms
54,640 KB |
testcase_07 | AC | 46 ms
55,112 KB |
testcase_08 | AC | 42 ms
55,440 KB |
testcase_09 | AC | 42 ms
54,116 KB |
testcase_10 | AC | 2,711 ms
258,024 KB |
testcase_11 | AC | 1,019 ms
148,132 KB |
testcase_12 | AC | 111 ms
80,968 KB |
testcase_13 | AC | 2,315 ms
231,448 KB |
testcase_14 | AC | 1,380 ms
181,404 KB |
testcase_15 | AC | 1,467 ms
189,368 KB |
testcase_16 | AC | 1,124 ms
158,048 KB |
testcase_17 | AC | 611 ms
121,648 KB |
testcase_18 | AC | 659 ms
125,320 KB |
testcase_19 | AC | 1,521 ms
184,764 KB |
testcase_20 | TLE | - |
testcase_21 | TLE | - |
testcase_22 | TLE | - |
testcase_23 | TLE | - |
testcase_24 | TLE | - |
testcase_25 | TLE | - |
ソースコード
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)