結果
問題 | No.2264 Gear Coloring |
ユーザー | vwxyz |
提出日時 | 2023-05-25 02:42:16 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 493 ms / 2,000 ms |
コード長 | 1,168 bytes |
コンパイル時間 | 293 ms |
コンパイル使用メモリ | 82,248 KB |
実行使用メモリ | 66,560 KB |
最終ジャッジ日時 | 2024-06-06 10:23:31 |
合計ジャッジ時間 | 3,074 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 42 ms
53,504 KB |
testcase_01 | AC | 40 ms
53,632 KB |
testcase_02 | AC | 41 ms
53,760 KB |
testcase_03 | AC | 45 ms
60,032 KB |
testcase_04 | AC | 42 ms
53,888 KB |
testcase_05 | AC | 66 ms
66,560 KB |
testcase_06 | AC | 47 ms
59,776 KB |
testcase_07 | AC | 65 ms
66,304 KB |
testcase_08 | AC | 46 ms
60,032 KB |
testcase_09 | AC | 44 ms
59,264 KB |
testcase_10 | AC | 46 ms
59,392 KB |
testcase_11 | AC | 47 ms
60,416 KB |
testcase_12 | AC | 47 ms
59,776 KB |
testcase_13 | AC | 47 ms
61,184 KB |
testcase_14 | AC | 493 ms
66,048 KB |
testcase_15 | AC | 80 ms
62,336 KB |
testcase_16 | AC | 411 ms
66,560 KB |
testcase_17 | AC | 54 ms
62,208 KB |
testcase_18 | AC | 55 ms
63,744 KB |
testcase_19 | AC | 66 ms
62,336 KB |
testcase_20 | AC | 59 ms
63,872 KB |
testcase_21 | AC | 53 ms
62,336 KB |
ソースコード
import sys readline=sys.stdin.readline import math from collections import defaultdict def Divisors(N): divisors=[] for i in range(1,N+1): if i**2>=N: break elif N%i==0: divisors.append(i) if i**2==N: divisors+=[i]+[N//i for i in divisors[::-1]] else: divisors+=[N//i for i in divisors[::-1]] return divisors def Factorize(N): assert N>=1 factors=defaultdict(int) for p in range(2,N): if p**2>N: break while N%p==0: factors[p]+=1 N//=p if N!=1: factors[N]+=1 return factors def LCM(n,m): if n or m: return abs(n)*abs(m)//math.gcd(n,m) return 0 N,M=map(int,readline().split()) A=list(map(int,readline().split())) mod=998244353 L=1 for a in A: L=LCM(L,a) D=Divisors(L) cnt=defaultdict(int) for d in D: cnt[d]=L//d for p in Factorize(L): for d in D: if d%p==0: cnt[d//p]-=cnt[d] ans=0 for d,c in cnt.items(): p=1 for a in A: g=math.gcd(a,d) p*=pow(M,g,mod) p%=mod ans+=p*c ans%=mod ans*=pow(L,mod-2,mod) ans%=mod print(ans)