結果

問題 No.2829 GCD Divination
ユーザー 👑 p-adicp-adic
提出日時 2024-06-20 15:34:38
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 124 ms / 2,000 ms
コード長 530 bytes
コンパイル時間 307 ms
コンパイル使用メモリ 82,764 KB
実行使用メモリ 77,216 KB
最終ジャッジ日時 2024-06-20 23:01:55
合計ジャッジ時間 3,414 ms
ジャッジサーバーID
(参考情報)
judge1 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 44 ms
60,604 KB
testcase_01 AC 45 ms
61,792 KB
testcase_02 AC 45 ms
60,840 KB
testcase_03 AC 50 ms
62,168 KB
testcase_04 AC 60 ms
68,748 KB
testcase_05 AC 124 ms
77,112 KB
testcase_06 AC 121 ms
77,216 KB
testcase_07 AC 121 ms
77,140 KB
testcase_08 AC 122 ms
76,932 KB
testcase_09 AC 119 ms
77,124 KB
testcase_10 AC 49 ms
63,064 KB
testcase_11 AC 47 ms
62,988 KB
testcase_12 AC 50 ms
64,116 KB
testcase_13 AC 48 ms
63,032 KB
testcase_14 AC 51 ms
63,208 KB
testcase_15 AC 51 ms
64,112 KB
testcase_16 AC 49 ms
62,416 KB
testcase_17 AC 49 ms
63,132 KB
testcase_18 AC 51 ms
63,048 KB
testcase_19 AC 49 ms
63,516 KB
testcase_20 AC 47 ms
62,288 KB
testcase_21 AC 51 ms
63,772 KB
testcase_22 AC 44 ms
62,108 KB
testcase_23 AC 52 ms
65,048 KB
testcase_24 AC 48 ms
62,868 KB
testcase_25 AC 47 ms
62,704 KB
testcase_26 AC 48 ms
62,744 KB
testcase_27 AC 50 ms
63,728 KB
testcase_28 AC 48 ms
63,404 KB
testcase_29 AC 48 ms
62,700 KB
testcase_30 AC 49 ms
63,996 KB
testcase_31 AC 48 ms
62,072 KB
testcase_32 AC 53 ms
62,136 KB
testcase_33 AC 51 ms
63,276 KB
testcase_34 AC 49 ms
62,668 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