Простые числа. Решето Эратосфена: различия между версиями
Перейти к навигации
Перейти к поиску
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; }
Ссылки
Теория:
Код:
Задачи: