結果
| 問題 |
No.1653 Squarefree
|
| コンテスト | |
| ユーザー |
titia
|
| 提出日時 | 2021-08-21 04:12:01 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 803 bytes |
| コンパイル時間 | 99 ms |
| コンパイル使用メモリ | 12,672 KB |
| 実行使用メモリ | 100,352 KB |
| 最終ジャッジ日時 | 2024-10-14 12:51:27 |
| 合計ジャッジ時間 | 6,508 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | -- * 3 |
| other | TLE * 1 -- * 37 |
ソースコード
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**6+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)
titia