結果
問題 | No.106 素数が嫌い!2 |
ユーザー |
![]() |
提出日時 | 2020-04-05 12:19:44 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 45 ms / 5,000 ms |
コード長 | 2,590 bytes |
コンパイル時間 | 975 ms |
コンパイル使用メモリ | 115,440 KB |
実行使用メモリ | 20,288 KB |
最終ジャッジ日時 | 2024-07-03 08:06:39 |
合計ジャッジ時間 | 2,490 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 13 |
ソースコード
#include <algorithm>#include <bitset>#include <cassert>#include <cmath>#include <ctime>#include <cstring>#include <functional>#include <iostream>#include <iomanip>#include <limits>#include <list>#include <map>#include <queue>#include <set>#include <stack>#include <string>#include <regex>#include <vector>#define fix(n) cout<<fixed<<setprecision(n)#define rep(i,n) for (int i = 0; i < (n); ++i)#define all(a) (a).begin(), (a).end()#define sort(a) sort(all(a))#define uniq(a) sort(a);(a).erase(unique(all(a)), (a).end())#define reverse(a) reverse(all(a))#define ctos(c) string(1, (c))#define out(d) cout << (d)#define outl(d) std::cout<<(d)<<"\n"#define YES() cout << "YES" << endl#define NO() cout << "NO" << endl#define Yes() cout << "Yes" << endl#define No() cout << "No" << endl#define ceil(x, y) ((x + y - 1) / (y))#define debug(x) cerr << #x << ": " << (x) << '\n'#define debug2(x, y) cerr << #x << ": " << (x) << ", " << #y << ": " << (y) << '\n'#define debug3(x, y, z) cerr << #x << ": " << (x) << ", " << #y << ": " << (y) << ", " << #z << ": " << (z) << '\n'#define dbg(v) for (size_t _ = 0; _ < v.size(); ++_){ cerr << #v << "[" << _ << "] : " << v[_] << '\n'; }#define pb push_back#define fst first#define int long long#define INF __LONG_LONG_MAX__using namespace std;using ll = long long;using ld = long double;using P = pair<ll,ll>;const ll MOD = 1000000007; // 10^9 + 7const int dx[8] = {1, 0, -1, 0, 1, -1, -1, 1};const int dy[8] = {0, 1, 0, -1, 1, 1, -1, -1};// 正の整数n未満の素数を全て求めるvector<ll> sieve(ll n) {vector<bool> is_prime(n);rep(i, n) is_prime[i] = true;is_prime[1] = false;for (ll i = 4; i < n; i += 2) is_prime[i] = false;ll lim = ll(sqrt(n));for (ll i = 3; i <= lim; i += 2) {for (ll j = 3; i * j < n; j += 2) {is_prime[i * j] = false;}}vector<ll> ret;for (ll i = 1; i < n; i++) {if (is_prime[i]) ret.push_back(i);}return ret;}void solve() {const ll MAX = 2 * 1000000;vector<ll> p = sieve(MAX);vector<ll> nums(MAX+1, 0);for (ll pi : p)for (int i = pi; i <= MAX; i += pi) nums[i]++;ll N, K; cin >> N >> K;ll ans = 0;for (int i = 1; i <= N; i++) {if (nums[i] >= K) ans++;}outl(ans);}signed main() {cin.tie(0);ios::sync_with_stdio(false);srand((unsigned)time(NULL));fix(12);solve();}