結果

問題 No.826 連絡網
ユーザー takakintakakin
提出日時 2020-07-14 22:07:58
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
WA  
実行時間 -
コード長 706 bytes
コンパイル時間 165 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 24,068 KB
最終ジャッジ日時 2024-11-18 22:22:09
合計ジャッジ時間 11,253 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

diff # 

import sys
input=lambda: sys.stdin.readline().rstrip()
n,p=map(int,input().split())
pf={}
for i in range(2,int(p**0.5)+1):
  while p%i==0:
    pf[i]=pf.get(i,0)+1
    p//=i
if p>1:
  pf[p]=1

def primes(n):
  is_prime=[True]*(n+1)
  is_prime[0]=False
  is_prime[1]=False
  for i in range(2,int(n**0.5)+1):
    if not is_prime[i]:
      continue
    for j in range(i*2,n+1,i):
      is_prime[j]=False
  return [i for i in range(n+1) if is_prime[i]]

P=primes(n)
PF=set(pf.keys())
for pp in P:
	if p*pp<=n:
		PF.add(pp)
	else:
		break
Chk=[0]*n
for i in range(n):
	if Chk[i]==1:
		continue
	else:
		if i+1 in PF:
			Chk[i]=1
			mul=2
			while mul*(i+1)<=n:
				Chk[mul*(i+1)-1]=1
				mul+=1
print(sum(Chk))
0