結果
問題 | No.1727 [Cherry 3rd Tune] Stray |
ユーザー | vwxyz |
提出日時 | 2023-05-25 02:01:20 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 31 ms / 6,000 ms |
コード長 | 993 bytes |
コンパイル時間 | 413 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 10,752 KB |
最終ジャッジ日時 | 2024-12-24 00:47:57 |
合計ジャッジ時間 | 1,208 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 27 ms
10,752 KB |
testcase_01 | AC | 28 ms
10,752 KB |
testcase_02 | AC | 27 ms
10,752 KB |
testcase_03 | AC | 31 ms
10,752 KB |
testcase_04 | AC | 30 ms
10,624 KB |
testcase_05 | AC | 27 ms
10,752 KB |
testcase_06 | AC | 28 ms
10,752 KB |
testcase_07 | AC | 30 ms
10,624 KB |
ソースコード
import sys readline=sys.stdin.readline 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 T=int(readline()) for t in range(T): N,C=map(int,readline().split()) D=Divisors(N) fact=Factorize(N) cnt=defaultdict(int) for d in D: cnt[d]+=N//d for p in fact: for d in D: if d%p==0: cnt[d//p]-=cnt[d] ans=pow(C,N)*N for d,c in cnt.items(): ans+=pow(C,d*2)*c ans//=2*N print(ans)