CF 955C Sad Power 枚举(Perfect Powers)

题意:Q次询问,每次询问[L,R]中有多少个数x 满足 x=a^p (a>0,p>1)
Q<=1e5. 1<=L<=R<=1e18.


上界为1e18 那么p最大不会超过61.
现在计算出f[1,n]有多少个数x 满足x=a^p .那么答案就是f[r]-f[l-1]
枚举p,二分此时a^p<=n的最大a为多少 然后发现重复的处理不来....wa...

可以先暴力预处理出1e18内,p>=3的所有Perfect Number. O(1e6).

然后对于给定的n,p==2的有sqrt(n)个.然后p>=3时,在vector中二分中最后一个<=n的位置即可.

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=2e5+5;
vector<ll> v;
void init()
{
	ll n=1e18;
	for(ll i=2;i*i*i<=n;i++) //when p>=3  enumerate a
	{
		ll s=i*i;
		while(s<=n/i)
		{
			s*=i;
			ll t=sqrt(s);
			if(t*t!=s) // 
				v.push_back(s);
		}
	}
	sort(v.begin(),v.end());
	v.erase(unique(v.begin(),v.end()),v.end());
}
ll calc(ll n)
{
	ll idx=lower_bound(v.begin(),v.end(),n)-v.begin();	
	if(idx<v.size()&&v[idx]>n)
		idx--;
	if(idx==v.size())
		idx--;
	//cout<<n<<' '<<sqrt(n)<<' '<<idx<<'\n';
	return ll(sqrt(n))+idx;
} 
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	init(); 
	ll Q,l,r;
	cin>>Q;
	while(Q--)
	{
		cin>>l>>r;
		cout<<calc(r)-calc(l-1)<<'\n';
	}
	return 0;
}