首页 > 其他 > 详细

P1835 素数密度_NOI导刊2011提高(04)

时间:2018-06-26 22:24:18      阅读:219      评论:0      收藏:0      [点我收藏+]

题目描述

给定区间[L,R](L≤R≤2147483647,R-L≤1000000),请计算区间中素数的个数。

输入输出格式

输入格式:

 

两个数L和R。

 

输出格式:

 

一行,区间中素数的个数。

 

输入输出样例

输入样例#1: 
2 11
输出样例#1: 
5

 

Solution:

  本题比较模板()。

  因为$int_{max}$范围内的合数只需用$\sqrt {int_{max}}\leq 50000$内的质数就能筛出。

  于是,预处理出$50000$内的质数,然后埃式筛掉$[l,r]$范围内的合数就好了,随便瞎搞一下就好了。

  最后统计一下质数个数并输出。

代码:

 1 #include<bits/stdc++.h>
 2 #define il inline
 3 #define ll long long
 4 #define For(i,a,b) for(ll (i)=(a);(i)<=(b);(i)++)
 5 #define Bor(i,a,b) for(int (i)=(b);(i)>=(a);(i)--)
 6 using namespace std;
 7 const int N=50005,M=1000005;
 8 int l,r,prime[N+5],cnt;
 9 bool isp[N+5],isprime[M+5];
10 
11 il void init(){
12     For(i,2,N){
13         if(!isp[i]) prime[++cnt]=i;
14         for(int j=1;j<=cnt&&prime[j]*i<=N;j++){
15             isp[prime[j]*i]=1;
16             if(i%prime[j]==0)break;
17         }
18     }
19 }
20 
21 il void solve(){
22     ll p,q;
23     For(i,1,cnt){
24         if(prime[i]>r)break;
25         p=ceil(l*1.0/prime[i]),q=ceil(r*1.0/prime[i]);
26         if(p==1)p=2;
27         for(int j=p;j<=q&&prime[i]*j<=r;j++)
28             isprime[prime[i]*j-l+1]=1;
29     }
30     int tot=0;
31     For(i,1,r-l+1) if(!isprime[i]) tot++;
32     cout<<tot;
33 }
34 
35 int main(){
36     init();
37     cin>>l>>r;
38     solve();
39     return 0;
40 }

 

P1835 素数密度_NOI导刊2011提高(04)

原文:https://www.cnblogs.com/five20/p/9231418.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!