結果

問題 No.2829 GCD Divination
ユーザー 👑 p-adicp-adic
提出日時 2024-06-20 15:34:49
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 177 ms / 2,000 ms
コード長 530 bytes
コンパイル時間 89 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 11,136 KB
最終ジャッジ日時 2024-06-20 23:01:56
合計ジャッジ時間 3,026 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 37 ms
10,880 KB
testcase_01 AC 37 ms
10,880 KB
testcase_02 AC 33 ms
10,880 KB
testcase_03 AC 34 ms
10,880 KB
testcase_04 AC 52 ms
10,880 KB
testcase_05 AC 177 ms
10,880 KB
testcase_06 AC 174 ms
11,008 KB
testcase_07 AC 158 ms
11,136 KB
testcase_08 AC 153 ms
11,136 KB
testcase_09 AC 146 ms
11,136 KB
testcase_10 AC 37 ms
10,880 KB
testcase_11 AC 32 ms
10,880 KB
testcase_12 AC 39 ms
10,752 KB
testcase_13 AC 34 ms
11,008 KB
testcase_14 AC 34 ms
11,008 KB
testcase_15 AC 39 ms
10,752 KB
testcase_16 AC 34 ms
11,008 KB
testcase_17 AC 38 ms
10,880 KB
testcase_18 AC 41 ms
10,880 KB
testcase_19 AC 36 ms
10,880 KB
testcase_20 AC 34 ms
10,880 KB
testcase_21 AC 43 ms
11,008 KB
testcase_22 AC 33 ms
10,880 KB
testcase_23 AC 51 ms
10,880 KB
testcase_24 AC 38 ms
10,880 KB
testcase_25 AC 35 ms
11,008 KB
testcase_26 AC 38 ms
11,008 KB
testcase_27 AC 38 ms
10,880 KB
testcase_28 AC 36 ms
10,752 KB
testcase_29 AC 35 ms
10,880 KB
testcase_30 AC 36 ms
10,880 KB
testcase_31 AC 33 ms
10,880 KB
testcase_32 AC 34 ms
11,008 KB
testcase_33 AC 42 ms
10,880 KB
testcase_34 AC 38 ms
11,008 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

R=range
L=9**4
P=[]
C=[0]*L
for i in R(2,L):
	if C[i]<1:
		P+=[i];j=i*i
		while j<L:C[j],j=1,j+i
E={}
def e(n):
	if n in E:return E[n]
	a=m=n
	for p in P:
		c=0
		while n%p<1:n//=p;c+=1
		if c:a=a//p*(p-1)
	if n>1:a=a//n*(n-1)
	E[m]=a
	return a
F={1:0}
def f(n):
	if n in F:return F[n]
	D=[1];m=n
	for p in P:
		c=0
		while n%p<1:n//=p;c+=1
		r=len(D)
		for i in R(c):D+=[d*p for d in D[-r:]]
	r=len(D)
	if n>1:D+=[d*n for d in D[-r:]];r*=2
	F[m]=(m+sum(f(m/D[d])*e(D[d])for d in R(1,r)))/(m-1)
	return F[m]
print(f(int(input())))
0