题意简述
多组数据。每次给定lll,rrr(l,rl,rl,r都<=1e18<=1e18<=1e18),求lll~rrr中有多少xxx满足x=apx=a^px=ap,其中aaa和ppp是两个整数满足a>0,p>1a>0,p>1a>0,p>1。
数据
输入:
6
1 4
9 9
5 7
12 29
137 591
1 1000000
输出:
2
1
0
3
17
1111
思路
会发现,当p>=3p>=3p>=3时,使得0<=ap<=1e180<=a^p<=1e180<=ap<=1e18的aaa并不多,p=3p=3p=3大概有1e61e61e6,然后后面更少。经计算,总共的个数(去重后)是100333210033321003332,就比1e61e61e6多一点。考虑暴力存下,排序去重后存着。p>=3p>=3p>=3的情况只要upperbound−lowerboundupperbound-lowerboundupperbound−lowerbound即可(lower_bound和upper_bound返回的都是指针,如果相减就珂以得到中间长度了。不用+1+1+1的原因是upper_bound找到的那个不珂以取)。记为a1a_1a1。
那如何处理p=2p=2p=2的情况?
相当于求有多少aaa满足l<=a2<=rl<=a^2<=rl<=a2<=r。
我们会发现,如果设S(n)S(n)S(n)为满足1<=a2<=n1<=a^2<=n1<=a2<=n的aaa个数,则我们要求的就是S(r)−S(l−1)S(r)-S(l-1)S(r)−S(l−1)。
由数学常识得,S(n)=nS(n)=\sqrt{n}S(n)=n。
故我们要求的就是r−l−1\sqrt{r}-\sqrt{l-1}r−l−1。记为a2a_2a2
然后只要把a1a_1a1和a2a_2a2加一起即可。
代码:
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int Max=1000000ll*1000000ll*1000000ll;//1e18
int Root(int x)
{
return (int)sqrt(x*1.0+0.5);
}
vector<int>tmp;
void Build()
{
for(int i=2;i<=1000000;i++)//枚举底数
{
int s=i*i*i;//至少三次方
while(s<=Max)
{
int rs=Root(s);
if (rs*rs<s) tmp.push_back(s);
//这个就是判s是否是平方数,为了避免和后面考虑p=2的情况重合
if (s>Max/i) break;//如果爆了赶快退出
else s*=i;//没爆继续乘
}
}
sort(tmp.begin(),tmp.end());
tmp.erase(unique(tmp.begin(),tmp.end()),tmp.end());//排序去重,方便lower_bound和upper_bound
//printf("%d\n",tmp.size());
//=1003332,如果要自己写(即不完全Copy我的),珂以先对一下这个答案
}
int l,r;
void Solve()
{
int a1=upper_bound(tmp.begin(),tmp.end(),r)-lower_bound(tmp.begin(),tmp.end(),l);
int a2=Root(r)-Root(l-1);
printf("%I64d\n",a1+a2);//上面说的a1和a2
}
main()
{
Build();
int T;scanf("%I64d",&T);
while(T--)
{
scanf("%I64d%I64d",&l,&r);
Solve();
}
return 0;
}
本文详细解析了Codeforces 955C题目,通过题意简述和数据介绍,阐述了解题思路。主要讨论了当p>=3和p=2时的不同处理方法,利用数学性质优化了暴力求解的过程,最终给出代码实现。
365

被折叠的 条评论
为什么被折叠?



