結果
問題 | No.1274 楽しい格子点 |
ユーザー | mkawa2 |
提出日時 | 2021-09-01 17:26:51 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 34 ms / 2,000 ms |
コード長 | 914 bytes |
コンパイル時間 | 145 ms |
コンパイル使用メモリ | 12,288 KB |
実行使用メモリ | 10,880 KB |
最終ジャッジ日時 | 2024-05-05 21:40:59 |
合計ジャッジ時間 | 3,309 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 33 ms
10,752 KB |
testcase_01 | AC | 31 ms
10,752 KB |
testcase_02 | AC | 31 ms
10,752 KB |
testcase_03 | AC | 30 ms
10,880 KB |
testcase_04 | AC | 31 ms
10,752 KB |
testcase_05 | AC | 32 ms
10,624 KB |
testcase_06 | AC | 31 ms
10,752 KB |
testcase_07 | AC | 30 ms
10,624 KB |
testcase_08 | AC | 30 ms
10,752 KB |
testcase_09 | AC | 30 ms
10,624 KB |
testcase_10 | AC | 31 ms
10,752 KB |
testcase_11 | AC | 30 ms
10,624 KB |
testcase_12 | AC | 31 ms
10,752 KB |
testcase_13 | AC | 29 ms
10,624 KB |
testcase_14 | AC | 31 ms
10,752 KB |
testcase_15 | AC | 30 ms
10,752 KB |
testcase_16 | AC | 31 ms
10,624 KB |
testcase_17 | AC | 30 ms
10,624 KB |
testcase_18 | AC | 30 ms
10,752 KB |
testcase_19 | AC | 31 ms
10,752 KB |
testcase_20 | AC | 30 ms
10,752 KB |
testcase_21 | AC | 31 ms
10,624 KB |
testcase_22 | AC | 32 ms
10,752 KB |
testcase_23 | AC | 32 ms
10,752 KB |
testcase_24 | AC | 31 ms
10,752 KB |
testcase_25 | AC | 32 ms
10,624 KB |
testcase_26 | AC | 32 ms
10,752 KB |
testcase_27 | AC | 31 ms
10,752 KB |
testcase_28 | AC | 31 ms
10,624 KB |
testcase_29 | AC | 31 ms
10,624 KB |
testcase_30 | AC | 31 ms
10,624 KB |
testcase_31 | AC | 30 ms
10,624 KB |
testcase_32 | AC | 31 ms
10,752 KB |
testcase_33 | AC | 31 ms
10,752 KB |
testcase_34 | AC | 30 ms
10,752 KB |
testcase_35 | AC | 31 ms
10,752 KB |
testcase_36 | AC | 31 ms
10,752 KB |
testcase_37 | AC | 31 ms
10,624 KB |
testcase_38 | AC | 31 ms
10,752 KB |
testcase_39 | AC | 29 ms
10,880 KB |
testcase_40 | AC | 30 ms
10,752 KB |
testcase_41 | AC | 30 ms
10,624 KB |
testcase_42 | AC | 31 ms
10,624 KB |
testcase_43 | AC | 30 ms
10,880 KB |
testcase_44 | AC | 30 ms
10,624 KB |
testcase_45 | AC | 30 ms
10,624 KB |
testcase_46 | AC | 30 ms
10,752 KB |
testcase_47 | AC | 30 ms
10,880 KB |
testcase_48 | AC | 31 ms
10,752 KB |
testcase_49 | AC | 34 ms
10,752 KB |
testcase_50 | AC | 31 ms
10,624 KB |
testcase_51 | AC | 30 ms
10,624 KB |
testcase_52 | AC | 31 ms
10,624 KB |
testcase_53 | AC | 31 ms
10,752 KB |
testcase_54 | AC | 33 ms
10,624 KB |
testcase_55 | AC | 30 ms
10,624 KB |
testcase_56 | AC | 33 ms
10,624 KB |
testcase_57 | AC | 31 ms
10,752 KB |
ソースコード
# g=gcd(A,B)とすると、xの移動距離dx=g*i、yの移動距離dy=g*jとなる # a=A/g,b=B/g(このときaとbは互いに素)とすると # a,bともに奇数の場合・・・i,jの偶奇は等しい必要がある # a,bの一方が奇数、もう一方が偶数の場合・・・任意のi,jがとれる from math import gcd A, B = map(int, input().split()) if A==0 and B==0: print(0.25) exit() g = gcd(A, B) a, b = A//g, B//g ans = 0 if a & 1 and b & 1: for i in range(20): x = 1+g*i if x > 20: break for j in range(20): if (i & 1) ^ (j & 1): continue y = 1+g*j if y > 20: break ans += pow(x+y, -x-y) else: for i in range(20): x = 1+g*i if x > 20: break for j in range(20): y = 1+g*j if y > 20: break ans += pow(x+y, -x-y) print("{:.7f}".format(ans))