結果
問題 | No.1653 Squarefree |
ユーザー |
|
提出日時 | 2021-08-20 23:07:30 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,330 ms / 2,000 ms |
コード長 | 871 bytes |
コンパイル時間 | 224 ms |
コンパイル使用メモリ | 82,788 KB |
実行使用メモリ | 84,072 KB |
最終ジャッジ日時 | 2024-10-14 08:37:23 |
合計ジャッジ時間 | 39,239 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 38 |
ソースコード
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)) 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 l,r = list(map(int, input().split())) ok = [1]*(r-l+1) # ok[i] : l+i for i in range(2, 10**6+100): ii = i*i ll = ((l+ii-1)//ii)*ii for k in range(ll, r+1, ii): ok[k-l] = 0 for i in range(1,10**6+100): ir = r//i k = isqrt(ir) kk = k*k if kk<=10**6: continue if k>1 and l<=i*kk<=r: while i*kk-l>=0: ok[i*kk-l] = 0 i -= 1 for v in range(max(2,l), r+1): vv = isqrt(v) if vv*vv==v: ok[v-l] = 0 ans = sum(ok) print(ans)