HUST 1214 Cubic-free numbers II(区间n=x^3*k的n的个数、枚举x容斥)

题目链接；
HUST 1214 Cubic-free numbers II
题意：
求区间$n\in [L,R]$，且$n$不可以表示成$n=x^{3}*k（x是大于1的正数）$的形式的$n$的个数。
数据范围：大概$10^7$
分析：
其实这题和HDU 2204 Eddy’s 爱好
差不多。
HDU那道题我们是枚举指数，然后容斥，因为考虑到数据范围。
同样的道理这里我们枚举底数，但是这样子依然会重复，例如$2^{3}*64$和$8^3*1$就被计算了两次，同时我们可以发现多计算的一定是底数的倍数，而且我们需要倒着处理。
这样子可以得到区间$[1,r]$中的解，那么区间$[L,R]$中的解减一下就好了。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <climits>
#include <cmath>
#include <ctime>
#include <cassert>
#define IOS ios_base::sync_with_stdio(0); cin.tie(0);
using namespace std;
typedef long long ll;
const double eps = 1e-8;
const int MAX_N = 20000100;

ll num[MAX_N];

ll solve(ll n)
{
    if(n <= 0) return 0;
    ll top = (ll)(pow(n * 1.0, 1.0 / 3.0) + eps);
    ll ans = 0;
    for(ll i = top; i >= 2; --i) {
        num[i] = n / (i * i * i);
        for(ll j = 2; j * i <= top; ++j) {
            num[i] -= num[j * i];
        }
        ans += num[i];
    }
    return n - ans;
}

int main()
{
    int T;
    ll L, R;
    scanf("%d", &T);
    while(T--) {
        scanf("%lld%lld", &L, &R);
        printf("%lld\n", solve(R - 1) - solve(L - 1));
    }
    return 0;
}