250-like Number

题意

给定一个整数 \(n\)，求有多少小于等于 \(n\) 的满足以下条件的整数 \(k\)：

\(k\) 可以被表示为 \(k = p \times q^3\)，其中 \(p \lt q\)，并且 \(p, q\) 均为质数。

数据范围

\(1 \leqslant n \leqslant 10^{18}\)，\(n\) 是整数。

思路

首先，我们发现这个式子中有一个立方，这个肯定就是突破口，因为 \(p < q\)，所以 \(p < q \leqslant 10^6\)。

因为要求 \(p,q\) 均为质数，所以先要筛出所有 \(10^6\) 以内的质数。

根据突破口，我们可以发现枚举 \(q\) 是不会TLE的，那么考虑第一步枚举 \(q\)，并且随着 \(q\) 的增大，\(p\) 的取值范围是不会变小的！

那么做法就出来了，先预处理出所有 \(10^6\) 以内的质数，然后用双指针求出答案即可。

复杂度

时间：\(O(V\log\log V + n)\)，说实话有点险，不建议在比赛时这么做
空间：\(O(n)\)

Code

点击查看代码

#include <iostream>
#include <vector>
#include <algorithm>

using namespace std;
using ll = long long;

const int N = 1e6 + 10;

int f[N];
ll n, ans;
vector<int> v;

void ES () { // 埃氏筛筛出所有质数
  f[1] = 1;
  for (int i = 2; i <= 1e6; i++) {
    if (f[i]) {
      continue;
    }
    v.push_back(i);
    for (int j = 2; i * j <= 1e6; j++) {
      f[i * j] = 1;
    }
  }
}

int main () {
  ios::sync_with_stdio(0), cin.tie(0);
  ES();
  cin >> n;
  for (int i = 1, k = v.size() - 1; i < v.size(); i++) { // 枚举 q, 双指针求答案
    ll j = 1ll * v[i] * v[i] * v[i]; // q 的立方
    if (j > n) { // 超出范围
      break;
    }
    for (; k >= 0 && 1ll * v[k] * j > n; k--) { // 移动 p 的范围
    }
    ans += min(i, k + 1); // 求答案，且 p 不能超过 q
  }
  cout << ans;
  return 0;
}