結果

問題 No.1653 Squarefree
ユーザー titia
提出日時 2021-08-21 04:11:40
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 901 ms / 2,000 ms
コード長 803 bytes
コンパイル時間 410 ms
コンパイル使用メモリ 82,036 KB
実行使用メモリ 246,272 KB
最終ジャッジ日時 2024-10-14 12:50:59
合計ジャッジ時間 32,384 ms
ジャッジサーバーID
(参考情報)
judge3 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 38
権限があれば一括ダウンロードができます

ソースコード

diff #

L0,R0=map(int,input().split())

if L0==1:
    ANS=1
    L0+=1
else:
    ANS=0
A=[]

for i in range(L0,R0+1):
    A.append(i)

import math

x=10**7+100

L=math.floor(math.sqrt(x)) # 平方根を求める

Primelist=[i for i in range(x+1)]
Primelist[1]=0 # 1は素数でないので0にする.
 
for i in Primelist:
    if i>L:
        break
    if i==0:
        continue
    for j in range(2*i,x+1,i):
        Primelist[j]=0

Primes=[Primelist[j] for j in range(x+1) if Primelist[j]!=0]

for p in Primes:
    for i in range((L0+p-1)//p*p,(R0+p)//p*p,p):
        x=i-L0

        if A[x]%(p*p)==0:
            A[x]=0

        if A[x]%p==0:
            A[x]//=p

for a in A:
    if a==0:
        continue
    sq=int(math.sqrt(a))
    if a==1 or sq*sq!=a:
        ANS+=1

print(ANS)
            
        



0