結果

問題 No.1946 ロッカーの問題
ユーザー 👑 SPD_9X2SPD_9X2
提出日時 2022-05-20 22:47:23
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 577 ms / 3,000 ms
コード長 1,987 bytes
コンパイル時間 1,664 ms
コンパイル使用メモリ 81,320 KB
実行使用メモリ 151,204 KB
最終ジャッジ日時 2023-10-20 13:37:40
合計ジャッジ時間 6,714 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 37 ms
53,332 KB
testcase_01 AC 39 ms
53,332 KB
testcase_02 AC 383 ms
151,204 KB
testcase_03 AC 577 ms
151,008 KB
testcase_04 AC 390 ms
151,020 KB
testcase_05 AC 297 ms
127,548 KB
testcase_06 AC 313 ms
133,112 KB
testcase_07 AC 277 ms
122,308 KB
testcase_08 AC 211 ms
103,412 KB
testcase_09 AC 479 ms
143,648 KB
testcase_10 AC 400 ms
132,680 KB
testcase_11 AC 96 ms
77,364 KB
testcase_12 AC 216 ms
109,768 KB
testcase_13 AC 38 ms
53,344 KB
testcase_14 AC 37 ms
53,344 KB
testcase_15 AC 72 ms
77,112 KB
testcase_16 AC 41 ms
53,344 KB
testcase_17 AC 290 ms
113,428 KB
testcase_18 AC 563 ms
151,020 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from sys import stdin

def Sieve(n): #n以下の素数全列挙(O(nloglogn)) retは素数が入ってる。divlisはその数字の素因数が一つ入ってる

    ret = []
    divlis = [-1] * (n+1) #何で割ったかのリスト(初期値は-1)
    
    flag = [True] * (n+1)
    flag[0] = False
    flag[1] = False

    ind = 2
    while ind <= n:

        if flag[ind]:
            ret.append(ind)

            ind2 = ind ** 2

            while ind2 <= n:
                flag[ind2] = False
                divlis[ind2] = ind
                ind2 += ind

        ind += 1

    return ret,divlis


# return Prime factorization dict of l,(l+1),...,(r-1)
# plis must include all primes under M**0.5 ( M=max(r**0.5 , r-l) )
# dict style: ret[18] = [2,3,3]
def primeFactorizationFactorial(l,r,plis):

    a = {}
    ret = {}
    for i in range(l,r):
        a[i] = i
        ret[i] = []
    
    for p in plis:
        if p**2 >= r:
            continue
        divfi = (l+p-1) // p * p
        for i in range(divfi,r,p):
            while a[i] % p == 0:
                a[i] //= p
                ret[i].append(p)
    for i in a:
        if a[i] != 1:
            ret[i].append(a[i])
    return ret



N,M = map(int,stdin.readline().split())

plis,_ = Sieve(N+10)
dic = primeFactorizationFactorial(1,N+10,plis)

#print (dic)

A = list(map(int,stdin.readline().split()))

lis = [0] * (N+1)
ANS = 0

for i in A:
    lis[i] = 1

for i in range(N,0,-1):

    if lis[i] == 1:

        d2 = {}
        for v in dic[i]:
            if v not in d2:
                d2[v] = 1
            else:
                d2[v] += 1

        lisA = [1]
        for p in d2:
            newlis = []
            wow = p
            for j in range(d2[p]):
                for k in lisA:
                    newlis.append(wow * k)
                wow *= p

            for j in newlis:
                lisA.append(j)

        for j in lisA:
            lis[j] ^= 1


    else:
        ANS += 1

print (ANS)
0