結果
問題 | No.1653 Squarefree |
ユーザー | kwm_t |
提出日時 | 2021-08-21 16:33:15 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 136 ms / 2,000 ms |
コード長 | 2,072 bytes |
コンパイル時間 | 1,878 ms |
コンパイル使用メモリ | 169,932 KB |
実行使用メモリ | 12,024 KB |
最終ジャッジ日時 | 2024-04-23 07:26:36 |
合計ジャッジ時間 | 6,009 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 136 ms
11,860 KB |
testcase_01 | AC | 11 ms
6,940 KB |
testcase_02 | AC | 11 ms
6,944 KB |
testcase_03 | AC | 11 ms
6,944 KB |
testcase_04 | AC | 10 ms
6,944 KB |
testcase_05 | AC | 11 ms
6,940 KB |
testcase_06 | AC | 10 ms
6,940 KB |
testcase_07 | AC | 10 ms
6,940 KB |
testcase_08 | AC | 10 ms
6,940 KB |
testcase_09 | AC | 9 ms
6,944 KB |
testcase_10 | AC | 10 ms
6,940 KB |
testcase_11 | AC | 70 ms
11,844 KB |
testcase_12 | AC | 67 ms
11,896 KB |
testcase_13 | AC | 128 ms
11,840 KB |
testcase_14 | AC | 128 ms
11,712 KB |
testcase_15 | AC | 129 ms
11,844 KB |
testcase_16 | AC | 132 ms
11,876 KB |
testcase_17 | AC | 130 ms
11,964 KB |
testcase_18 | AC | 132 ms
11,924 KB |
testcase_19 | AC | 131 ms
12,024 KB |
testcase_20 | AC | 129 ms
11,964 KB |
testcase_21 | AC | 132 ms
11,952 KB |
testcase_22 | AC | 132 ms
11,916 KB |
testcase_23 | AC | 131 ms
11,840 KB |
testcase_24 | AC | 128 ms
11,960 KB |
testcase_25 | AC | 135 ms
11,708 KB |
testcase_26 | AC | 133 ms
11,772 KB |
testcase_27 | AC | 130 ms
11,904 KB |
testcase_28 | AC | 128 ms
11,660 KB |
testcase_29 | AC | 130 ms
11,904 KB |
testcase_30 | AC | 132 ms
11,840 KB |
testcase_31 | AC | 128 ms
11,756 KB |
testcase_32 | AC | 136 ms
11,900 KB |
testcase_33 | AC | 134 ms
11,876 KB |
testcase_34 | AC | 130 ms
11,952 KB |
testcase_35 | AC | 128 ms
11,836 KB |
testcase_36 | AC | 10 ms
6,940 KB |
testcase_37 | AC | 11 ms
6,940 KB |
testcase_38 | AC | 10 ms
6,940 KB |
testcase_39 | AC | 10 ms
6,940 KB |
testcase_40 | AC | 10 ms
6,940 KB |
ソースコード
#include "bits/stdc++.h" //#include "atcoder/all" using namespace std; //using namespace atcoder; //using mint = modint1000000007; //const int mod = 1000000007; //using mint = modint998244353; //const int mod = 998244353; const int INF = 1e9; //const long long LINF = 1e18; //const bool debug = false; #define rep(i, n) for (int i = 0; i < (n); ++i) #define rep2(i,l,r)for(int i=(l);i<(r);++i) #define rrep(i, n) for (int i = (n-1); i >= 0; --i) #define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i) #define all(x) (x).begin(),(x).end() #define allR(x) (x).rbegin(),(x).rend() #define endl "\n" #define P pair<int,int> template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; } template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; } //エラトステネスの篩(素数の数を返している) //11→2,3,5,7,11 vector<long long> sieve(long long n = 1000000) { vector<bool>p(n + 1, true); vector<long long>ret; p[0] = false; p[1] = false; for (int i = 2; i <= n; ++i) { if (p[i]) { ret.push_back(i); for (int j = 2 * i; j <= n; j += i) p[j] = false; } } return ret; } bool squarefree(long long n) { long long ng = 0; long long ok = INF; while (1 < abs(ok - ng)) { long long c = (ng + ok) / 2; if (n <= (c * c)) ok = c; else ng = c; } bool ret = false; if (n == ok * ok)ret = true; return ret; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); long long l, r; cin >> l >> r; auto p = sieve(); vector<long long>x(r - l + 1); rep(i, x.size()) x[i] = l + i; rep(i, p.size()) { long long pri = p[i]; long long fir = ((l + (pri - 1)) / pri) * pri; if (r < fir)continue; for (; fir <= r; fir += pri) { long long idx = fir - l; if (0 == (x[idx] % (pri * pri)))x[idx] = 0; else if (0 == (x[idx] % pri))x[idx] /= pri; } } int ans = 0; rep(i, x.size()) { if (0 == x[i])continue; else if (1 == x[i])ans++; else ans += 1 - squarefree(x[i]); } cout << ans << endl; return 0; }