結果
| 問題 | No.1653 Squarefree |
| ユーザー |
 kusaf_
|
| 提出日時 | 2024-10-03 13:19:15 |
| 言語 | C++23 (gcc 12.3.0 + boost 1.83.0) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 1,790 bytes |
| コンパイル時間 | 3,950 ms |
| コンパイル使用メモリ | 259,240 KB |
| 実行使用メモリ | 347,888 KB |
| 最終ジャッジ日時 | 2024-10-03 13:19:26 |
| 合計ジャッジ時間 | 10,061 ms |
|
ジャッジサーバーID （参考情報） |
judge4 / judge3 |
（要ログイン）
テストケース
テストケース表示| 入力 | 結果 | 実行時間 実行使用メモリ |
|---|---|---|
| testcase_00 | TLE | - |
| testcase_01 | -- | - |
| testcase_02 | -- | - |
| testcase_03 | -- | - |
| testcase_04 | -- | - |
| testcase_05 | -- | - |
| testcase_06 | -- | - |
| testcase_07 | -- | - |
| testcase_08 | -- | - |
| testcase_09 | -- | - |
| testcase_10 | -- | - |
| testcase_11 | -- | - |
| testcase_12 | -- | - |
| testcase_13 | -- | - |
| testcase_14 | -- | - |
| testcase_15 | -- | - |
| testcase_16 | -- | - |
| testcase_17 | -- | - |
| testcase_18 | -- | - |
| testcase_19 | -- | - |
| testcase_20 | -- | - |
| testcase_21 | -- | - |
| testcase_22 | -- | - |
| testcase_23 | -- | - |
| testcase_24 | -- | - |
| testcase_25 | -- | - |
| testcase_26 | -- | - |
| testcase_27 | -- | - |
| testcase_28 | -- | - |
| testcase_29 | -- | - |
| testcase_30 | -- | - |
| testcase_31 | -- | - |
| testcase_32 | -- | - |
| testcase_33 | -- | - |
| testcase_34 | -- | - |
| testcase_35 | -- | - |
| testcase_36 | -- | - |
| testcase_37 | -- | - |
| testcase_38 | -- | - |
| testcase_39 | -- | - |
| testcase_40 | -- | - |
ソースコード
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
vector<ll> PrimeEnum(ll n) {
vector<bool> sieve(n / 3 + 1, true);
for(ll p = 5, d = 4, i = 1, sq = sqrt(n); p <= sq; p += d = 6 - d, i++) {
if(!sieve[i]) { continue; }
for(ll q = p * p / 3, r = d * p / 3 + (d * p % 3 == 2), s = 2 * p, qe = ssize(sieve); q < qe; q += r = s - r) { sieve[q] = false; }
}
vector<ll> ret = {2, 3};
for(ll p = 5, d = 4, i = 1; p <= n; p += d = 6 - d, i++) {
if(sieve[i]) { ret.emplace_back(p); }
}
while(!ret.empty() && ret.back() > n) { ret.pop_back(); }
return ret;
}
uint64_t Root(uint64_t a, int k = 2) {
if(a <= 1 || k == 1) { return a; }
if(k >= 64) { return 1; }
auto check = [&](__uint128_t n) {
__uint128_t x = 1, m = n;
for(ll p = k; p; p >>= 1, m *= m) {
if(p & 1) { x *= m; }
}
return x <= a;
};
uint64_t n = powl(a, (long double)(1.0) / k);
while(!check(n)) { --n; }
while(check(n + 1)) { ++n; }
return n;
}
vector<vector<pair<ll, ll>>> PrimeRange(ll l, ll r) {
int n = r - l + 1;
vector<ll> v(n);
iota(v.begin(), v.end(), l);
vector<vector<pair<ll, ll>>> p(n);
for(auto &i : PrimeEnum(Root(r))) {
for(ll j = (l + i - 1) / i * i; j <= r; j += i) {
int c = 0;
while(v[j - l] % i == 0) {
v[j - l] /= i;
c++;
}
if(c) { p[j - l].emplace_back(i, c); }
}
}
for(int i = 0; i < n; i++) {
if(v[i] != 1) { p[i].emplace_back(v[i], 1); }
}
return p;
}
int main() {
ios::sync_with_stdio(false);
cin.tie(nullptr);
ll L, R;
cin >> L >> R;
ll ans = 0;
for(auto &v : PrimeRange(L, R)) {
bool t = true;
for(auto &[p, c] : v) {
if(c >= 2) {
t = false;
break;
}
}
ans += t;
}
cout << ans << "\n";
}