結果

問題 No.2973 シュニレルマン積分入門
ユーザー ecotteaecottea
提出日時 2024-12-02 21:59:13
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 967 bytes
コンパイル時間 229 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 89,984 KB
最終ジャッジ日時 2024-12-02 21:59:32
合計ジャッジ時間 18,078 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 463 ms
89,984 KB
testcase_01 AC 480 ms
89,724 KB
testcase_02 AC 479 ms
89,600 KB
testcase_03 AC 511 ms
89,320 KB
testcase_04 AC 484 ms
89,620 KB
testcase_05 AC 511 ms
89,216 KB
testcase_06 AC 483 ms
89,216 KB
testcase_07 AC 481 ms
89,572 KB
testcase_08 AC 489 ms
89,216 KB
testcase_09 AC 492 ms
89,400 KB
testcase_10 AC 513 ms
89,472 KB
testcase_11 AC 486 ms
89,552 KB
testcase_12 AC 514 ms
89,812 KB
testcase_13 AC 485 ms
89,472 KB
testcase_14 AC 483 ms
89,404 KB
testcase_15 AC 511 ms
89,600 KB
testcase_16 AC 487 ms
89,564 KB
testcase_17 AC 511 ms
89,472 KB
testcase_18 AC 485 ms
89,216 KB
testcase_19 AC 485 ms
89,620 KB
testcase_20 AC 512 ms
89,824 KB
testcase_21 WA -
testcase_22 AC 498 ms
89,708 KB
testcase_23 WA -
testcase_24 AC 487 ms
89,344 KB
testcase_25 WA -
testcase_26 AC 510 ms
89,472 KB
testcase_27 WA -
testcase_28 WA -
testcase_29 AC 493 ms
89,600 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import cmath
import math
from fractions import Fraction

N = int(input())

M = 1 * 3 * 4 * 5 * 6 * 7 * 8 * 9 * 10

val = 0
for k in range(M):
    num = cmath.exp(2 * k * N * math.pi * 1j / M)
    dnm = cmath.exp(4 * k * math.pi * 1j / M)
    dnm += cmath.exp(2 * k * math.pi * 1j / M)
    dnm += 10
    val += num / dnm
val /= M
# print(val)

# これもだめ
dnm = 10**max(-N+1,0)
num = round(val.real * dnm)
print(num, dnm, sep='/')

# Fraction で有理数近似しようと思ったけど精度が足りない
# x = Fraction(val.real)
# y = x.limit_denominator(10**(-N+1))
# # print(y)

# y_num = y.numerator
# y_dnm = y.denominator

# tmp = y_dnm

# c2 = 0
# c5 = 0

# while tmp % 2 == 0:
#     c2 += 1
#     tmp //= 2
# while tmp % 5 == 0:
#     c5 += 1
#     tmp //= 5
# # print(c2, c5, tmp)

# if c2 > c5:
#     y_num *= 5**(c2-c5)
#     y_dnm *= 5**(c2-c5)
# elif c2 < c5:
#     y_num *= 2**(c5-c2)
#     y_dnm *= 2**(c5-c2)

# print(y_num, y_dnm, sep='/')
0