結果

問題 No.2552 Not Coprime, Not Divisor
ユーザー konchakoncha
提出日時 2023-11-25 15:26:48
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 722 bytes
コンパイル時間 323 ms
コンパイル使用メモリ 81,700 KB
実行使用メモリ 280,420 KB
最終ジャッジ日時 2023-11-25 15:26:55
合計ジャッジ時間 6,910 ms
ジャッジサーバーID
(参考情報)
judge11 / judge12
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 TLE -
testcase_01 -- -
testcase_02 -- -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
権限があれば一括ダウンロードができます

ソースコード

diff #

from collections import deque

def primes(n):
    p = []
    flgs = [False]*(n+1)
    
    for i in range(2, n+1):
        if not flgs[i]:
            p.append(i)
            for j in range(i*2, n+1, i):
                flgs[j] = True
    return p


N = int(input())
answer = 0
ps = primes(N)

queue = deque([(p, i, 0) for i, p in enumerate(ps)])

while queue:
    now, end, count = queue.popleft()
    
    cnt = N//now
    
    if count%2 == 0:
        answer += cnt*(cnt-1) // 2
    else:
        answer -= cnt*(cnt-1) // 2
        
    for j, p in enumerate(ps[end+1:], start=end+1):
        if now * p > N: break
        queue.append((now*p, j, count+1))

for i in range(2, N+1):
    answer -= (N//i)-1

print(answer)
0