CF955C Sad powers

题目大意

给你 \(q\) 个询问，每次询问 \([l,r]\) 这个区间内满足 \(x=a^p(a>0,p>1)\) 的 \(x\) 的数量。

\(1⩽q⩽10^5\)，\(1\leqslant l\leqslant r\leqslant 10^{18}\)。

解题思路

显然，\(\sqrt[2]{10^{18}}=10^9\)，\(\sqrt[3]{10^{18}}=10^6\)，\(\sqrt[18]{10^{18}}=10\)，且区间 \([l,r]\) 里的指数为偶数的数为 \(\sqrt[2]{r}-\sqrt[2]{l-1}\)。

故可以对于 \(x \in [2,10^6],y \in [2,10^6],2 \nmid y\) 预处理出所有满足 \(x^y\) 不是平方数的数，然后排序去重即可。

每次询问二分再加上 \(\sqrt[2]{r}-\sqrt[2]{l-1}\) 即可。

CODE

#include <bits/stdc++.h>

using namespace std;

#define int long long

const int N = 3e6 + 7;

const int MX = 1e18;

int q, cnt;

int l, r, num[N];

inline int solve(int x)
{
	int idx = lower_bound(num + 1, num + cnt + 1, x) - num;
	if((idx <= cnt && num[idx] > x) || idx > cnt) idx--;
	return idx + (int)sqrt(x);
}

signed main()
{
	for(int i = 2; i <= 1000000; ++i)
	{
		long long k = i * i;
		for(; k <= MX / i; )
		{
			k *= i;
			int sqrtk = sqrt(k);
			if(sqrtk * sqrtk != k) num[++cnt] = k;
		}
	}
	sort(num + 1, num + cnt + 1);
	cnt = unique(num + 1, num + cnt + 1) - num - 1;
	cin >> q;
	while(q--)
	{
		cin >> l >> r;
		cout << solve(r) - solve(l - 1) << "\n";
	}
	return 0;
}

标签：10,idx,int,18,sqrt,CF955C,num,Sad,powers
来源： https://www.cnblogs.com/orzz/p/15503237.html