結果
| 問題 |
No.1274 楽しい格子点
|
| コンテスト | |
| ユーザー |
mkawa2
|
| 提出日時 | 2021-09-01 17:26:51 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
AC
|
| 実行時間 | 33 ms / 2,000 ms |
| コード長 | 914 bytes |
| コンパイル時間 | 614 ms |
| コンパイル使用メモリ | 12,288 KB |
| 実行使用メモリ | 10,880 KB |
| 最終ジャッジ日時 | 2024-11-27 23:03:29 |
| 合計ジャッジ時間 | 3,879 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 57 |
ソースコード
# 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))
mkawa2