UVa10539

Posted 2020-08-09 SilverNebula

http://vjudge.net/problem/UVA-10539

 

先打出来sqrt(n)以内的素数表，然后对于每个素数x，他对答案的贡献就是最大的p使x^p<=n，即log(x,n)。注意精度误差。 
用1..r的减去1..l-1的就是答案。

 

 1 /*by SilverN*/
 2 #include<iostream>
 3 #include<algorithm>
 4 #include<cstring>
 5 #include<cstdio>
 6 #include<cmath>
 7 #define LL long long
 8 using namespace std;
 9 const int mxn=1200000;
10 int pri[mxn],cnt=0;
11 bool vis[mxn];
12 int num[mxn];
13 void Pri(){
14     cnt=0;int i,j;
15     for(i=2;i<mxn;i++){
16         if(!vis[i])pri[++cnt]=i;
17         for(j=1;j<=cnt && i*pri[j]<mxn;j++){
18             vis[i*pri[j]]=1;
19             if(i%pri[j]==0)break;
20         }
21     }
22     return;
23 }
24 LL query(LL n){
25     LL ans=0;
26     int i,j;
27     for(i=1;i<=cnt && (LL)pri[i]*pri[i]<=n;i++){
28         ans+=log(n+0.1)/log(pri[i])-1;
29     }
30     return ans;
31 }
32 int n;
33 LL l,r;
34 int main(){
35     Pri();
36     int T;
37     scanf("%d",&T);
38     while(T--){
39         scanf("%lld%lld",&l,&r);
40         printf("%lld\n",query(r)-query(l-1));
41     }
42     return 0;
43 }