結果
問題 | No.1653 Squarefree |
ユーザー | asayia |
提出日時 | 2021-08-21 11:56:06 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 202 ms / 2,000 ms |
コード長 | 1,595 bytes |
コンパイル時間 | 1,473 ms |
コンパイル使用メモリ | 168,744 KB |
実行使用メモリ | 20,992 KB |
最終ジャッジ日時 | 2024-10-15 03:16:17 |
合計ジャッジ時間 | 7,590 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 189 ms
20,856 KB |
testcase_01 | AC | 2 ms
5,248 KB |
testcase_02 | AC | 2 ms
5,248 KB |
testcase_03 | AC | 2 ms
5,248 KB |
testcase_04 | AC | 2 ms
5,248 KB |
testcase_05 | AC | 2 ms
5,248 KB |
testcase_06 | AC | 2 ms
5,248 KB |
testcase_07 | AC | 2 ms
5,248 KB |
testcase_08 | AC | 2 ms
5,248 KB |
testcase_09 | AC | 3 ms
5,248 KB |
testcase_10 | AC | 2 ms
5,248 KB |
testcase_11 | AC | 87 ms
19,200 KB |
testcase_12 | AC | 88 ms
19,280 KB |
testcase_13 | AC | 183 ms
20,736 KB |
testcase_14 | AC | 185 ms
20,736 KB |
testcase_15 | AC | 202 ms
20,864 KB |
testcase_16 | AC | 202 ms
20,864 KB |
testcase_17 | AC | 188 ms
20,736 KB |
testcase_18 | AC | 202 ms
20,992 KB |
testcase_19 | AC | 193 ms
20,864 KB |
testcase_20 | AC | 183 ms
20,864 KB |
testcase_21 | AC | 187 ms
20,864 KB |
testcase_22 | AC | 193 ms
20,864 KB |
testcase_23 | AC | 190 ms
20,736 KB |
testcase_24 | AC | 189 ms
20,728 KB |
testcase_25 | AC | 181 ms
20,864 KB |
testcase_26 | AC | 185 ms
20,736 KB |
testcase_27 | AC | 162 ms
20,736 KB |
testcase_28 | AC | 183 ms
20,736 KB |
testcase_29 | AC | 181 ms
20,740 KB |
testcase_30 | AC | 180 ms
20,608 KB |
testcase_31 | AC | 175 ms
20,736 KB |
testcase_32 | AC | 190 ms
20,736 KB |
testcase_33 | AC | 196 ms
20,864 KB |
testcase_34 | AC | 186 ms
20,864 KB |
testcase_35 | AC | 174 ms
20,736 KB |
testcase_36 | AC | 20 ms
5,248 KB |
testcase_37 | AC | 21 ms
5,248 KB |
testcase_38 | AC | 20 ms
5,248 KB |
testcase_39 | AC | 21 ms
5,376 KB |
testcase_40 | AC | 20 ms
5,376 KB |
ソースコード
#include <bits/stdc++.h> #define maxx(a, b, c) max(max(a, b), c) #define minn(a, b, c) min(min(a, b), c) #define between(x, l, r) (x >= l && x <= r) #define point(a, b) ((a - 1) * m + b) #define y second #define x first using namespace std; const int N = 1000010; typedef long long LL; LL primes[N], cnt; bool st[N]; LL p[N], check[N]; void get_primes(int n) // 线性筛质数 { for (int i = 2; i <= n; i ++ ) { if (!st[i]) primes[cnt ++ ] = i; for (int j = 0; primes[j] <= n / i; j ++ ) { st[primes[j] * i] = true; if (i % primes[j] == 0) break; } } } signed main() { LL l, r; cin >> l >> r; //auto now = clock(); LL d = pow(r, 1.0 / 3) + 2; // 筛根号三次内的数 get_primes(d); // O(n) for (int i = 0; i <= r - l; i ++ ) p[i] = l + i; // O(n) int res = r - l + 1; for (int i = 0; i < cnt; i ++ ) // O(nloglogn * logn) { LL start = (primes[i] - l % primes[i]) % primes[i]; //if(i == 1) cout << start << '\n'; for (int j = start; j <= r - l; j += primes[i]) { if(check[j]) continue; int tmp = 0; while (p[j] % primes[i] == 0) { p[j] /= primes[i]; tmp ++ ; } if(tmp >= 2 && !check[j]) check[j] = 1, res -- ; } } for (int i = 0; i <= r - l; i ++ ) { if(check[i] || p[i] == 1) continue; LL k = sqrt(p[i]); if(k * k == p[i]) res -- ; } cout << res << '\n'; return 0; } /* 607148 */