题目描述:
1<=a,b<=10^18,p,q都是素数 2<=p,q<=10^9;
求在[a,b]内可以表示为 x*p^k*q^m k > m 的数的个数
分析:
由于要小于b,因此m一定小于 log10(b)/log10(p*q);
因此我们可以枚举m,中间计数的时候需要用到容斥。
具体看代码:
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
LL mypow(LL a,int b)
{
LL ans = 1;
while(b){
if(b&1){
ans=ans*a;
b--;
}
b>>=1;
a=a*a;
}
return ans;
}
int main()
{
LL a,b,p,q;
while(~scanf("%lld%lld%lld%lld",&a,&b,&p,&q)){
int mmax = log10(b*1.0)/log10(p*q*1.0)+1;
LL ans = 0;
for(int i=0;i<=mmax;i++){
if(mypow(p,i+1)>b*1.0/mypow(q,i))//防止爆long long
break;
for(int j=i+1;j<64;j++){
if(mypow(p,j)>b*1.0/mypow(q,i)) break;//防止爆long long
LL tmp=mypow(p,j)*mypow(q,i);
LL cnt1 = b/tmp,cnt2=(a-1)/tmp;//由于是闭区间 因此要用a-1;
ans += cnt1;
ans -= cnt2;
ans -= cnt1/p;
ans -= cnt1/q;
ans += cnt1/p/q;
ans += cnt2/p;
ans += cnt2/q;
ans -=cnt2/p/q;
}
}
printf("%lld\n",ans);
}
return 0;
}