結果
問題 | No.1653 Squarefree |
ユーザー |
|
提出日時 | 2021-08-20 22:14:21 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 440 ms / 2,000 ms |
コード長 | 1,741 bytes |
コンパイル時間 | 2,013 ms |
コンパイル使用メモリ | 196,520 KB |
最終ジャッジ日時 | 2025-01-23 23:37:01 |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 38 |
ソースコード
#include <bits/stdc++.h>using namespace std;#define rep(i, n) for(int i = 0; i < n; i++)#define rep2(i, x, n) for(int i = x; i <= n; i++)#define rep3(i, x, n) for(int i = x; i >= n; i--)#define each(e, v) for(auto &e: v)#define pb push_back#define eb emplace_back#define all(x) x.begin(), x.end()#define rall(x) x.rbegin(), x.rend()#define sz(x) (int)x.size()using ll = long long;using pii = pair<int, int>;using pil = pair<int, ll>;using pli = pair<ll, int>;using pll = pair<ll, ll>;const int MOD = 1000000007;//const int MOD = 998244353;const int inf = (1<<30)-1;const ll INF = (1LL<<60)-1;template<typename T> bool chmax(T &x, const T &y) {return (x < y)? (x = y, true) : false;};template<typename T> bool chmin(T &x, const T &y) {return (x > y)? (x = y, true) : false;};struct io_setup{io_setup(){ios_base::sync_with_stdio(false);cin.tie(NULL);cout << fixed << setprecision(15);}} io_setup;int main(){ll L, R; cin >> L >> R;int K = R-L;vector<bool> ok(K+1, true);ll M = 1000000;rep2(i, 1, M){ll l = 1, r = inf; //(l,r]while(r-l > 1){ll m = (l+r)/2;ll t = m*m;if(t > (L-1)/i) r = m;else l = m;}for(ll t = r; t*t*i <= R; t++){ok[t*t*i-L] = false;}}vector<bool> ps(M+1, true);for(ll i = 2; i <= M; i++){if(!ps[i]) continue;for(ll t = ((L+i-1)/i)*i; t <= R; t += i){if(t%(i*i) == 0) ok[t-L] = false;}for(int j = i+i; j <= M; j += i){ps[j] = false;}}int ans = 0;rep(i, K+1){if(ok[i]) ans++;}cout << ans << '\n';}