結果

問題 No.2971 無理積分
ユーザー 👑 p-adicp-adic
提出日時 2023-08-25 13:39:03
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 31 ms / 1,000 ms
コード長 541 bytes
コンパイル時間 375 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 11,008 KB
最終ジャッジ日時 2024-08-16 01:00:18
合計ジャッジ時間 2,910 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 28 ms
10,752 KB
testcase_01 AC 31 ms
10,752 KB
testcase_02 AC 27 ms
10,880 KB
testcase_03 AC 30 ms
10,880 KB
testcase_04 AC 30 ms
10,880 KB
testcase_05 AC 28 ms
10,752 KB
testcase_06 AC 28 ms
10,880 KB
testcase_07 AC 30 ms
10,752 KB
testcase_08 AC 28 ms
10,752 KB
testcase_09 AC 28 ms
10,880 KB
testcase_10 AC 28 ms
10,752 KB
testcase_11 AC 28 ms
11,008 KB
testcase_12 AC 30 ms
10,752 KB
testcase_13 AC 29 ms
10,880 KB
testcase_14 AC 29 ms
10,880 KB
testcase_15 AC 27 ms
10,880 KB
testcase_16 AC 28 ms
10,880 KB
testcase_17 AC 27 ms
10,752 KB
testcase_18 AC 27 ms
10,752 KB
testcase_19 AC 27 ms
11,008 KB
testcase_20 AC 27 ms
10,880 KB
testcase_21 AC 28 ms
10,880 KB
testcase_22 AC 29 ms
10,752 KB
testcase_23 AC 29 ms
10,880 KB
testcase_24 AC 28 ms
10,752 KB
testcase_25 AC 27 ms
10,752 KB
testcase_26 AC 29 ms
10,752 KB
testcase_27 AC 28 ms
10,752 KB
testcase_28 AC 27 ms
10,880 KB
testcase_29 AC 27 ms
10,880 KB
testcase_30 AC 29 ms
10,752 KB
testcase_31 AC 29 ms
10,752 KB
testcase_32 AC 29 ms
11,008 KB
testcase_33 AC 28 ms
10,752 KB
testcase_34 AC 28 ms
11,008 KB
testcase_35 AC 27 ms
10,752 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