結果

問題 No.781 円周上の格子点の数え上げ
ユーザー convexineqconvexineq
提出日時 2021-03-16 09:46:11
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 709 bytes
コンパイル時間 188 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 409,372 KB
最終ジャッジ日時 2024-04-25 10:07:15
合計ジャッジ時間 6,786 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 TLE -
testcase_01 TLE -
testcase_02 TLE -
testcase_03 TLE -
testcase_04 TLE -
testcase_05 TLE -
testcase_06 TLE -
testcase_07 TLE -
testcase_08 TLE -
testcase_09 TLE -
testcase_10 TLE -
testcase_11 TLE -
testcase_12 TLE -
testcase_13 TLE -
testcase_14 TLE -
testcase_15 TLE -
testcase_16 TLE -
testcase_17 TLE -
testcase_18 TLE -
testcase_19 TLE -
testcase_20 TLE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

def Eratosthenes(N): #N以下の素数のリストを返す
    #iが素数のときis_prime_list[i]=1，それ以外は0
    N+=1
    is_prime_list = [True]*N
    m = int(N**0.5)+1
    for i in range(3,m,2):
        if is_prime_list[i]:
            is_prime_list[i*i::2*i]=[False]*((N-i*i-1)//(2*i)+1)
    return [2] + [i for i in range(3,N,2) if is_prime_list[i]]

n = 10**7
r = [1]*(n+1)
r[0] = 0
for p in Eratosthenes(n):
    if p==2: continue
    for i in range(p,n+1,p):
        j = i//p
        c = 1
        while j%p==0:
            j //= p
            c += 1
        if p%4==1:
            r[i] *= c+1
        else:
            if c%2: r[i] = 0

x,y = map(int,input().split())
print(max(r[x:y+1])*4)
0