Almost Prime Numbers UVA - 10539

问题

分析

比较笨的方法，对素数打表后，再对每一个素数的指数进行打表，然后二分查找

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <vector>
using namespace std;
typedef long long LL;
const LL maxn=1000005,Inf=0x3f3f3f3f;
int isprime[maxn],prime[maxn],pn=0;
vector<LL> table[maxn];
void getPrime(){
    fill(isprime,isprime+maxn,1);
    for(LL i=2;i<maxn;++i){
        if(isprime[i]) prime[pn++]=i;
        for(int j=0;j<pn && i*prime[j]<maxn;++j){
            isprime[i*prime[j]]=0;
            if(i%prime[j]==0) break;
        }
    }
}

void getTable(){
    for(long long i=0;i<pn;++i){
        long long x=prime[i],t=x;
        x*=t;
        while(x<1000000000000LL){
            table[i].push_back(x);
            x*=t;
        }
    }
}

int main(void){
    getPrime();
    getTable();
    LL n,a,b;
    scanf("%lld",&n);
    while(n--){
        scanf("%lld%lld",&a,&b);
        LL lim=sqrt(b+0.5);
        LL ans=0;
        for(int i=0;i<pn && prime[i]<=lim;++i){
            int t1=upper_bound(table[i].begin(),table[i].end(),b)-table[i].begin();
            int t2=lower_bound(table[i].begin(),table[i].end(),a)-table[i].begin();
            ans+=t1-t2;
        }
        printf("%lld\n",ans);
    }
    return 0;
}

使用log进行计算
参照：https://www.cnblogs.com/SilverNebula/p/5891134.html
在使用log的时候还是不够熟练，我一开始去想如何取整消除误差去了，应该直接往上加一个0.1然后取整就能很快地消除误差了，后面减1是为了消除质数的一次幂

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <vector>
using namespace std;
typedef long long LL;
const LL maxn=1000005,Inf=0x3f3f3f3f;
const double eps=1e-5;
int isprime[maxn],prime[maxn],pn=0;
vector<LL> table[maxn];
void getPrime(){
    fill(isprime,isprime+maxn,1);
    for(LL i=2;i<maxn;++i){
        if(isprime[i]) prime[pn++]=i;
        for(int j=0;j<pn && i*prime[j]<maxn;++j){
            isprime[i*prime[j]]=0;
            if(i%prime[j]==0) break;
        }
    }
}

int main(void){
    getPrime();
    LL n,a,b;
    scanf("%lld",&n);
    while(n--){
        scanf("%lld%lld",&a,&b);
        LL lim=sqrt(b+0.5);
        LL ans=0;
        for(int i=0;i<pn && prime[i]<=lim;++i){
            int t1=log(b+0.1)/log(prime[i])-1; //减1是为了去掉质数的一次方
            ans+=t1;
        }
        lim=sqrt(a+0.5);
        for(int i=0;i<pn && prime[i]<=lim;++i){
            int t1=log(a-1+0.1)/log(prime[i])-1; //减1是为了去掉质数
            ans-=t1;
        }
        printf("%lld\n",ans);
    }
    return 0;
}