結果

問題 No.2971 無理積分
ユーザー 👑 p-adicp-adic
提出日時 2023-08-25 13:44:29
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 38 ms / 1,000 ms
コード長 541 bytes
コンパイル時間 581 ms
コンパイル使用メモリ 82,324 KB
実行使用メモリ 54,300 KB
最終ジャッジ日時 2024-08-16 01:00:21
合計ジャッジ時間 3,180 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 36 ms
52,776 KB
testcase_01 AC 35 ms
53,252 KB
testcase_02 AC 36 ms
54,300 KB
testcase_03 AC 36 ms
54,128 KB
testcase_04 AC 35 ms
52,704 KB
testcase_05 AC 36 ms
52,916 KB
testcase_06 AC 36 ms
53,016 KB
testcase_07 AC 34 ms
53,260 KB
testcase_08 AC 35 ms
53,392 KB
testcase_09 AC 35 ms
53,148 KB
testcase_10 AC 33 ms
53,292 KB
testcase_11 AC 35 ms
53,084 KB
testcase_12 AC 35 ms
53,808 KB
testcase_13 AC 35 ms
53,188 KB
testcase_14 AC 36 ms
54,080 KB
testcase_15 AC 35 ms
53,360 KB
testcase_16 AC 36 ms
53,404 KB
testcase_17 AC 35 ms
53,048 KB
testcase_18 AC 38 ms
53,628 KB
testcase_19 AC 38 ms
54,020 KB
testcase_20 AC 36 ms
53,820 KB
testcase_21 AC 36 ms
53,376 KB
testcase_22 AC 36 ms
53,000 KB
testcase_23 AC 35 ms
53,564 KB
testcase_24 AC 35 ms
52,704 KB
testcase_25 AC 35 ms
52,984 KB
testcase_26 AC 34 ms
53,248 KB
testcase_27 AC 35 ms
53,084 KB
testcase_28 AC 34 ms
53,052 KB
testcase_29 AC 36 ms
54,228 KB
testcase_30 AC 34 ms
53,004 KB
testcase_31 AC 35 ms
53,264 KB
testcase_32 AC 35 ms
53,196 KB
testcase_33 AC 35 ms
53,180 KB
testcase_34 AC 36 ms
53,268 KB
testcase_35 AC 38 ms
53,472 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

R,O=range,print
T=10**12
D=[1]*14
A=D[:]
B=D[:]
for i in R(1,14):D[i],A[i],B[i]=D[i-1]*(3*i+1)*i*2,A[i-1]*(3-2*i),B[i-1]*i*2
N=int(input())
if N<1:O(0.4),exit()
if N<2:O(10/9),exit()
n=0
while pow(N,3*n-5)<=T:n+=1
n-=1
l,r=0,N
while l+1<r:
	m=(l+r)//2
	if m*m>N:r=m
	else:l=m
s,M=N*N*sum((D[n-1]//(3*i+1)//B[i])*A[i]*N**(3*(n-1-i))for i in R(n))*T,N**(3*n-2)
l,r=l*s//D[n-1]//M,r*s//max(1,D[n-1]-1)//max(1,M-1)
s=N*s**2
while l+1<r:
	m=(l+r)//2
	if (m*D[n-1]*M)**2>s:r=m
	else:l=m
a=str(l)
a="0"*(13-len(a))+a
a=a[:-12]+"."+a[-12:]
O(a[:-6])
0