結果
| 問題 | No.781 円周上の格子点の数え上げ | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 2021-05-25 23:35:33 | 
| 言語 | PyPy3 (7.3.15) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 356 ms / 2,000 ms | 
| コード長 | 892 bytes | 
| コンパイル時間 | 168 ms | 
| コンパイル使用メモリ | 82,432 KB | 
| 実行使用メモリ | 216,448 KB | 
| 最終ジャッジ日時 | 2024-10-14 14:19:46 | 
| 合計ジャッジ時間 | 3,557 ms | 
| ジャッジサーバーID (参考情報) | judge4 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| other | AC * 21 | 
ソースコード
import sys
input = lambda : sys.stdin.readline().rstrip()
sys.setrecursionlimit(2*10**5+10)
write = lambda x: sys.stdout.write(x+"\n")
debug = lambda x: sys.stderr.write(x+"\n")
writef = lambda x: print("{:.12f}".format(x))
x,y = list(map(int, input().split()))
eps = 10**-7
def isqrt(n):
    """x*x<=nなる最大のx
    """
    x,y = n, (n+1)//2
    while y<x:
        x,y = y, (y + n//y) // 2
    return x
def sub(x):
    ans = 0
    for v in range(x):
        vv = x - v*v
        if vv<=0:
            break
        if (int((vv**0.5)+eps)**2 - vv) == 0:
            ans += 4
#             print(v,vv)
    return ans
vals = [0]*(y+1)
for v in range(4000):
    if v*v>=y:
        break
    v2 = v*v
    for u in range(4000):
        u2 = u*u
        if v2 + u2>y:
            break
        vals[v2+u2] += 4
#         if v2+u2==25:
#             print(v,u)
ans = max(vals[x:])
print(ans)
            
            
            
        