区间数字因子的个数为质数的和

Lucky Numbers

#include <iostream>
#include <algorithm>
#include <map>
#include <stack>
#include <queue>
#include <set>
#include <vector>
#include <cstdlib>
#include <cmath>
#include <string>
#include <cstdio>
#include <cstring>
#include <functional>
#define  all(x) x.begin(),x.end()
#define  lowbit(x) (x&(-x))
#define  ll long long
#define  ull unsigned long long
#define  max(a,b) (a>b?a:b)
#define  min(a,b) (a>b?b:a)
#define  PI acos(-1.0)
#define  fir first
#define  sec second
#define  inf 0x3f3f3f3f
#define  E 2.7182818284590452353602874713527
using namespace std;

/**********************************************Head-----Template****************************************/
bool Finish_read;
template<class T>inline void read(T &x){Finish_read=0;x=0;int f=1;char ch=getchar();while(!isdigit(ch)){if(ch=='-')f=-1;if(ch==EOF)return;ch=getchar();}while(isdigit(ch))x=x*10+ch-'0',ch=getchar();x*=f;Finish_read=1;}
template<class T>inline void print(T x){if(x/10!=0)print(x/10);putchar(x%10+'0');}
template<class T>inline void writeln(T x){if(x<0)putchar('-');x=abs(x);print(x);putchar('\n');}
template<class T>inline void write(T x){if(x<0)putchar('-');x=abs(x);print(x);}
ll gcd(ll a,ll b){return b==0?a:gcd(b,a%b);}
ll lcm(ll a,ll b){return a/gcd(a,b)*b;}
/********************************Head----Template**********************************************/

const int maxn = 1e5+7;
int vis[maxn];
vector<int>prime;

const ll limit=1e10+7;

vector<pair<ll,int> >vec;

void init(){
    memset(vis,0,sizeof(vis));
    for(int i=2;i<maxn;i++){
        if(!vis[i]){
            prime.push_back(i);
            for(int j=i+i;j<maxn;j+=i) vis[j]=1;
        }
    }
    for(int x:prime){
        ll tmp=1ll*x;
        for(int i=2;;i++){
            tmp*=x;
            if(tmp>=1&&tmp<=limit) {
                if(!vis[i+1]) vec.push_back({tmp,i+1});
            }
            else break;
        }
    }
    sort(all(vec));
    int sz=vec.size();
    for(int i=1;i<sz;i++) vec[i].sec+=vec[i-1].sec;
}

ll solve(ll l,ll r){
    int retr=lower_bound(all(vec),make_pair(r+1,-1))-vec.begin();
    int retl=lower_bound(all(vec),make_pair(l,-1))-vec.begin();
    ll ans=0;
    if(retr) ans+=vec[retr-1].sec;//本身是对r+1进行二分，因子找到的是比r+1大的位置，这样就去除了相等后位置+1的情况;
    if(retl) ans-=vec[retl-1].sec;
    return ans;
}

int main (){
    init();
    int q;read(q);
    while(q--){
        ll l,r;read(l),read(r);
        writeln(solve(l,r));
    }
    return 0;
}
/**
链接：https://www.codechef.com/problems/UEMALGO5
题意：区间满足条件的数字的因子个数和  Q:1e5  L,R 1e10
条件: 1. 该数本身是非质数
      2. 该数因子个数是质数
分析:{4,9} 满足条件的数只能表示成n = p^(k) 且p为质数,k+1必须为质数
因此我们可以进行打表 
对于2的贡献 只能是{2^2,3},{2^4,5},{2^6,7}....{2^k,k+1};//最多30几个
对于3的贡献 只能是{3^2,3},{3^4,5},{3^6,7}....{2^k,k+1};
.......
因此可以得到全部的满足条件的二元组，存vector进行排序二分查找，注意细节;

  */