結果

問題 No.1653 Squarefree
ユーザー 👑 KazunKazun
提出日時 2021-08-21 00:24:49
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 255 ms / 2,000 ms
コード長 1,967 bytes
コンパイル時間 359 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 97,236 KB
最終ジャッジ日時 2024-04-22 09:47:07
合計ジャッジ時間 9,110 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 242 ms
93,184 KB
testcase_01 AC 61 ms
74,368 KB
testcase_02 AC 67 ms
74,368 KB
testcase_03 AC 75 ms
79,232 KB
testcase_04 AC 69 ms
74,752 KB
testcase_05 AC 68 ms
74,368 KB
testcase_06 AC 72 ms
74,624 KB
testcase_07 AC 76 ms
74,624 KB
testcase_08 AC 71 ms
74,240 KB
testcase_09 AC 70 ms
74,624 KB
testcase_10 AC 69 ms
74,240 KB
testcase_11 AC 180 ms
89,088 KB
testcase_12 AC 174 ms
89,344 KB
testcase_13 AC 234 ms
93,056 KB
testcase_14 AC 240 ms
97,236 KB
testcase_15 AC 230 ms
93,752 KB
testcase_16 AC 243 ms
93,440 KB
testcase_17 AC 231 ms
93,656 KB
testcase_18 AC 252 ms
97,152 KB
testcase_19 AC 227 ms
94,976 KB
testcase_20 AC 253 ms
97,036 KB
testcase_21 AC 231 ms
92,736 KB
testcase_22 AC 255 ms
96,640 KB
testcase_23 AC 237 ms
93,152 KB
testcase_24 AC 231 ms
92,544 KB
testcase_25 AC 229 ms
92,856 KB
testcase_26 AC 251 ms
96,768 KB
testcase_27 AC 226 ms
93,312 KB
testcase_28 AC 235 ms
94,080 KB
testcase_29 AC 231 ms
92,864 KB
testcase_30 AC 247 ms
96,896 KB
testcase_31 AC 226 ms
93,056 KB
testcase_32 AC 243 ms
95,232 KB
testcase_33 AC 243 ms
94,488 KB
testcase_34 AC 231 ms
93,608 KB
testcase_35 AC 245 ms
93,184 KB
testcase_36 AC 62 ms
74,624 KB
testcase_37 AC 69 ms
74,624 KB
testcase_38 AC 70 ms
74,496 KB
testcase_39 AC 70 ms
74,368 KB
testcase_40 AC 69 ms
74,880 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#エラトステネスの篩
def Sieve_of_Eratosthenes(N,mode=False):
    """Nまでのエラトステネスの篩を実行

    N:自然数
    mode:False->素数のリスト,True->素数かどうかのリスト
    (False->[2,3,5,...],True->[0,0,1,1,0,1,...])
    """

    if N==0:
        return [0]

    T=[1]*(N+1)
    T[0]=T[1]=0

    for x in range(4,N+1,2): T[x]=0
    for x in range(9,N+1,3): T[x]=0

    a=5
    Flag=0
    while a*a<=N:
        if T[a]:
            b=a*a
            c=2*a
            while b<=N:
                T[b]=0
                b+=c
        a+=2+2*Flag
        Flag^=1

    if mode:
        return T
    else:
        return [k for k in range(N+1) if T[k]]

#floor(a^(1/k)) を求める.
def Floor_Root(a,k):
    """floor(a^(1/k)) を求める.

    a:非負整数
    k:正の整数
    """
    assert 0<=a and 0<k
    if a==0: return 0
    if k==1: return a

    #大体の値を求める.
    x=int(pow(a,1/k))

    #増やす
    while pow(x+1,k)<=a:
        x+=1

    #減らす
    while pow(x,k)>a:
        x-=1
    return x

def kth_Power(a,k):
    """ 整数 a が k 乗数かどうかを求め, そうならば, b^k=a を満たす k を返す.

    [Input]
    a:int
    k:int (k>0)

    [Output]
    存在しない  : None
    存在する    : b^k=a を満たす b
    """

    a_abs=abs(a)
    if a: sgn=a//a_abs
    else: sgn=0

    b=Floor_Root(a_abs,k)
    if pow(sgn*b,k)==a:
        return sgn*b
    else:
        return None
#==================================================
L,R=map(int,input().split())
P=Sieve_of_Eratosthenes(10**6)

X=list(range(L,R+1))
for p in P:
    p2=p*p
    q=p*((L+(p-1))//p)
    for x in range(q,R+1,p):
        if X[x-L]>0:
            if X[x-L]%p2==0:
                X[x-L]=0
            elif X[x-L]%p==0:
                X[x-L]//=p

Y=0
for x in range(L,R+1):
    if X[x-L]==1:
        Y+=1
        continue

    if X[x-L]>0 and kth_Power(X[x-L],2)==None:
        Y+=1

print(Y)
0