結果
問題 | No.144 エラトステネスのざる |
ユーザー | Navier_Boltzmann |
提出日時 | 2023-11-15 18:40:23 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 800 ms / 2,000 ms |
コード長 | 1,433 bytes |
コンパイル時間 | 219 ms |
コンパイル使用メモリ | 81,700 KB |
実行使用メモリ | 101,004 KB |
最終ジャッジ日時 | 2023-11-15 18:40:32 |
合計ジャッジ時間 | 8,667 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 115 ms
97,164 KB |
testcase_01 | AC | 122 ms
97,164 KB |
testcase_02 | AC | 122 ms
97,164 KB |
testcase_03 | AC | 119 ms
97,164 KB |
testcase_04 | AC | 120 ms
97,164 KB |
testcase_05 | AC | 120 ms
97,164 KB |
testcase_06 | AC | 158 ms
100,364 KB |
testcase_07 | AC | 158 ms
100,364 KB |
testcase_08 | AC | 158 ms
100,364 KB |
testcase_09 | AC | 159 ms
100,364 KB |
testcase_10 | AC | 159 ms
100,364 KB |
testcase_11 | AC | 158 ms
100,364 KB |
testcase_12 | AC | 158 ms
100,364 KB |
testcase_13 | AC | 775 ms
101,004 KB |
testcase_14 | AC | 800 ms
101,004 KB |
testcase_15 | AC | 795 ms
101,004 KB |
testcase_16 | AC | 773 ms
101,004 KB |
testcase_17 | AC | 784 ms
101,004 KB |
testcase_18 | AC | 756 ms
101,004 KB |
testcase_19 | AC | 781 ms
101,004 KB |
ソースコード
from collections import * from itertools import * from functools import * from heapq import * import sys,math N,p = input().split() N = int(N) p = float(p) class prime_factorize(): def __init__(self,M=10**6): self.sieve = [-1]*(M+1) self.sieve[1] = 1 self.p = [False]*(M+1) self.mu = [1]*(M+1) for i in range(2,M+1): if self.sieve[i] == -1: self.p[i] = True i2 = i**2 for j in range(i2,M+1,i2): self.mu[j] = 0 for j in range(i,M+1,i): self.sieve[j] = i self.mu[j] *= -1 def factors(self,x): tmp = [] while self.sieve[x] != x: tmp.append(self.sieve[x]) x //= self.sieve[x] tmp.append(self.sieve[x]) return tmp def is_prime(self,x): return self.p[x] def mobius(self,x): return self.mu[x] pf = prime_factorize() ans = 0 mem = defaultdict(lambda:-1) mem[0] = 1 def f(x): if mem[x]!=-1: return mem[x] mem[x] = f(x-1)*(1-p) return mem[x] for i in range(2,N+1): C = Counter(pf.factors(i)) t = 1 for v in C.values(): t *= (1+v) # print(i,t) ans += f(t-2) print(ans)