Простые числа. Решето Эратосфена: различия между версиями
Перейти к навигации
Перейти к поиску
Ctrlalt (обсуждение | вклад) |
Ctrlalt (обсуждение | вклад) |
||
| Строка 28: | Строка 28: | ||
===Факторизация числа (получение списка простых делителей), O(sqrt(N))=== | ===Факторизация числа (получение списка простых делителей), O(sqrt(N))=== | ||
{| width="100%" | |||
| width=50% | | |||
vector<int> factorize(int n) { | vector<int> factorize(int n) { | ||
vector<int> factorization; | vector<int> factorization; | ||
| Строка 42: | Строка 44: | ||
return factorization; | return factorization; | ||
} | } | ||
| width="10px" | | |||
| width=50% | | |||
map<int, int> factorize(int n) { | |||
map<int, int> factorization; | |||
for (long long d = 2; d * d <= n; d++) { | |||
while (n % d == 0) { | |||
factorization[d]++; | |||
n /= d; | |||
} | |||
} | |||
if (n != 1) | |||
factorization[n]++; | |||
return factorization; | |||
} | |||
|} | |||
===Решето Эратосфена, O(NloglogN)=== | ===Решето Эратосфена, O(NloglogN)=== | ||
Версия от 03:54, 15 ноября 2022
Проверка числа на простоту, O(sqrt(N))
bool isPrime(int n) {
if (n < 2)
return 0;
for (long long d = 2; d * d <= n; d++)
if (n % d == 0)
return 0;
return 1;
}
Получение списка всех делителей числа, O(sqrt(N))
vector<int> getDivisors(int n) {
vector<int> divisors;
for (long long d = 1; d * d <= n; d++) {
if (n % d == 0) {
divisors.push_back(d);
if (d * d != n)
divisors.push_back(n / d);
}
}
sort(divisors.begin(), divisors.end());
return divisors;
}
Факторизация числа (получение списка простых делителей), O(sqrt(N))
vector<int> factorize(int n) {
vector<int> factorization;
for (long long d = 2; d * d <= n; d++) {
while (n % d == 0) {
factorization.push_back(d);
n /= d;
}
}
if (n != 1)
factorization.push_back(n);
return factorization;
}
|
map<int, int> factorize(int n) {
map<int, int> factorization;
for (long long d = 2; d * d <= n; d++) {
while (n % d == 0) {
factorization[d]++;
n /= d;
}
}
if (n != 1)
factorization[n]++;
return factorization;
}
|
Решето Эратосфена, O(NloglogN)
vector<int> getPrimes() {
const int N = 1e6;
vector<int> isPrime(N, 1), primes;
for (int i = 2; i < N; i++) {
if (isPrime[i]) {
primes.push_back(i);
for (long long j = 1LL * i * i; j < N; j += i)
isPrime[j] = 0;
}
}
return primes;
}
Решето Эратосфена, O(N)
vector<int> getPrimes() {
const int N = 1e6;
vector<int> minDivisor(N), primes;
for (int i = 2; i < N; i++) {
if (!minDivisor[i]) {
minDivisor[i] = i;
primes.push_back(i);
}
for (int j = 0; j < primes.size() && primes[j] <= minDivisor[i] && 1LL * i * primes[j] < N; j++)
minDivisor[i * primes[j]] = primes[j];
}
return primes;
}
Ссылки
Теория:
Код:
Задачи: