結果
| 問題 |
No.826 連絡網
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2020-06-03 17:35:11 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 12 ms / 2,000 ms |
| コード長 | 2,035 bytes |
| コンパイル時間 | 1,579 ms |
| コンパイル使用メモリ | 172,292 KB |
| 実行使用メモリ | 7,636 KB |
| 最終ジャッジ日時 | 2024-11-25 05:14:23 |
| 合計ジャッジ時間 | 2,665 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 30 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
#define rep(i, n) for(int i = 0; i < (int)n; ++i)
#define FOR(i, a, b) for(int i = a; i < (int)b; ++i)
#define rrep(i, n) for(int i = ((int)n - 1); i >= 0; --i)
typedef long long ll;
typedef long double ld;
const int Inf = 1e9;
const double EPS = 1e-9;
const int MOD = 1e9 + 7;
class LinearSieve {
public:
vector<int> lp, pr;
LinearSieve() {}
LinearSieve(int n) {
lp.resize(n + 1);
lp[0] = -1, lp[1] = -1;
init(n);
}
void init(int x) {
FOR (i, 2, x + 1) {
if (lp[i] == 0) {
lp[i] = i;
pr.push_back(i);
}
for (int j = 0; j < (int)pr.size() && pr[j] <= lp[i] && (ll)i * pr[j] <= x; ++j) {
lp[i * pr[j]] = pr[j];
}
}
}
bool isPrime(int x) {
return lp[x] == x;
}
vector<int> factorList(int x) {
vector<int> res;
while (x != 1) {
res.push_back(lp[x]);
x /= lp[x];
}
return res;
}
vector<pair<int, int> > factorize(int x) {
vector<int> p = factorList(x);
if (p.size() == 0) return {};
vector<pair<int, int> > res(1, make_pair(p[0], 0));
rep (i, p.size()) {
if (res.back().first == p[i]) {
res.back().second++;
} else {
res.emplace_back(p[i], 1);
}
}
return res;
}
};
int main() {
cin.tie(nullptr);
ios::sync_with_stdio(0);
int n, p;
cin >> n >> p;
LinearSieve ls = LinearSieve(n + 10);
if (p == 1 || (p > n / 2 && ls.isPrime(p))) {
cout << 1 << endl;
return 0;
}
int cnt = 1;
FOR (i, n / 2 + 1, n + 1) {
if (ls.isPrime(i)) cnt++;
}
cout << n - cnt << endl;
return 0;
}