結果

問題 No.2971 無理積分
ユーザー 👑 p-adic
提出日時 2023-08-25 13:39:03
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 36
権限があれば一括ダウンロードができます

ソースコード

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