結果
問題 | No.1653 Squarefree |
ユーザー | 👑 emthrm |
提出日時 | 2021-08-20 22:57:30 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 163 ms / 2,000 ms |
コード長 | 2,073 bytes |
コンパイル時間 | 2,565 ms |
コンパイル使用メモリ | 201,612 KB |
実行使用メモリ | 19,584 KB |
最終ジャッジ日時 | 2024-04-22 09:18:22 |
合計ジャッジ時間 | 7,872 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 155 ms
19,564 KB |
testcase_01 | AC | 11 ms
7,716 KB |
testcase_02 | AC | 10 ms
7,720 KB |
testcase_03 | AC | 10 ms
7,724 KB |
testcase_04 | AC | 11 ms
7,716 KB |
testcase_05 | AC | 11 ms
7,844 KB |
testcase_06 | AC | 11 ms
7,712 KB |
testcase_07 | AC | 10 ms
7,844 KB |
testcase_08 | AC | 10 ms
7,716 KB |
testcase_09 | AC | 11 ms
7,844 KB |
testcase_10 | AC | 11 ms
7,840 KB |
testcase_11 | AC | 142 ms
19,444 KB |
testcase_12 | AC | 134 ms
19,440 KB |
testcase_13 | AC | 151 ms
19,500 KB |
testcase_14 | AC | 153 ms
19,464 KB |
testcase_15 | AC | 148 ms
19,436 KB |
testcase_16 | AC | 146 ms
19,380 KB |
testcase_17 | AC | 151 ms
19,516 KB |
testcase_18 | AC | 146 ms
19,436 KB |
testcase_19 | AC | 151 ms
19,364 KB |
testcase_20 | AC | 162 ms
19,572 KB |
testcase_21 | AC | 145 ms
19,456 KB |
testcase_22 | AC | 153 ms
19,568 KB |
testcase_23 | AC | 163 ms
19,320 KB |
testcase_24 | AC | 151 ms
19,356 KB |
testcase_25 | AC | 151 ms
19,420 KB |
testcase_26 | AC | 157 ms
19,332 KB |
testcase_27 | AC | 152 ms
19,412 KB |
testcase_28 | AC | 145 ms
19,476 KB |
testcase_29 | AC | 145 ms
19,416 KB |
testcase_30 | AC | 152 ms
19,584 KB |
testcase_31 | AC | 150 ms
19,404 KB |
testcase_32 | AC | 155 ms
19,548 KB |
testcase_33 | AC | 145 ms
19,500 KB |
testcase_34 | AC | 141 ms
19,424 KB |
testcase_35 | AC | 144 ms
19,408 KB |
testcase_36 | AC | 10 ms
7,716 KB |
testcase_37 | AC | 10 ms
7,844 KB |
testcase_38 | AC | 10 ms
7,908 KB |
testcase_39 | AC | 10 ms
7,716 KB |
testcase_40 | AC | 10 ms
7,716 KB |
ソースコード
#define _USE_MATH_DEFINES #include <bits/stdc++.h> using namespace std; #define FOR(i,m,n) for(int i=(m);i<(n);++i) #define REP(i,n) FOR(i,0,n) #define ALL(v) (v).begin(),(v).end() using ll = long long; constexpr int INF = 0x3f3f3f3f; constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL; constexpr double EPS = 1e-8; constexpr int MOD = 1000000007; // constexpr int MOD = 998244353; constexpr int dy[] = {1, 0, -1, 0}, dx[] = {0, -1, 0, 1}; constexpr int dy8[] = {1, 1, 0, -1, -1, -1, 0, 1}, dx8[] = {0, -1, -1, -1, 0, 1, 1, 1}; template <typename T, typename U> inline bool chmax(T &a, U b) { return a < b ? (a = b, true) : false; } template <typename T, typename U> inline bool chmin(T &a, U b) { return a > b ? (a = b, true) : false; } struct IOSetup { IOSetup() { std::cin.tie(nullptr); std::ios_base::sync_with_stdio(false); std::cout << fixed << setprecision(20); } } iosetup; std::vector<int> prime_sieve(int n, bool get_only_prime) { std::vector<int> prime, smallest_prime_factor(n + 1); std::iota(smallest_prime_factor.begin(), smallest_prime_factor.end(), 0); for (int i = 2; i <= n; ++i) { if (smallest_prime_factor[i] == i) prime.emplace_back(i); for (int p : prime) { if (i * p > n || p > smallest_prime_factor[i]) break; smallest_prime_factor[i * p] = p; } } return get_only_prime ? prime : smallest_prime_factor; } int main() { constexpr int M = 1000000; ll l, r; cin >> l >> r; vector<int> is_squarefree(r - l + 1, true); vector<ll> tmp(r - l + 1); iota(ALL(tmp), l); for (int p : prime_sieve(M, true)) { assert(p >= 1); for (ll i = (l + p - 1) / p * p; i <= r; i += p) { const int idx = i - l; if (!is_squarefree[idx] || tmp[idx] % p > 0) continue; tmp[idx] /= p; if (tmp[idx] % p == 0) is_squarefree[idx] = false; } } int ans = 0; REP(i, r - l + 1) { if (!is_squarefree[i]) continue; if (tmp[i] == 1) { ++ans; } else { ll sq = llround(sqrt(tmp[i])); ans += sq * sq != tmp[i]; } } cout << ans << '\n'; return 0; }