【數論】8.30題解-prime素數密度

prime

洛谷p835


題目描述

給定區間[L, R](L <= R <= 2147483647, R-L <= 1000000),請計算區間中
素數的個數。

輸入輸出

輸入
兩個數 L 和 R。
輸出
一行,區間中素數的個數。

樣例

樣例輸入

2 11

樣例輸出

5

說明

時空限制
時間限制 1s/testcase
空間限制 32MB

思路

R-L<=1000000
L <= R <= 2147483647
時間上用質數篩能過
求出2~45000的所有質數(sqrt(2147483647)大約是4300+)
然后將所有是這些質數的倍數的數刪掉,剩下的就是質數
空間上只有32mb,數組不能開太大,考慮到R-L<=1000000
完全可以只開到100005,之后對于大于100000的L和R,在存質數時
減去L即可。

代碼

/*
b[i]==1  i不是質數
b[i]==0  i是質數
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int i,j,k,l,t,n,m,ans;
int zhi[100005];
bool b[100005],z[1000005];
inline int read() {
    int x=0,w=1;
    char ch=0;
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') w=-1,ch=getchar();
    while(ch>='0'&&ch<='9')x=x*10+ch-48,ch=getchar();
    return x*w;
}
int main() {
    freopen("prime.in","r",stdin);
    freopen("prime.out","w",stdout);
    for(i=2; i<=sqrt(100000); i++) {
        if(!b[i]) {
            for(j=2; j<=100000/i; j++) {
                b[i*j]=1;
            }
        }
    }
    for(i=2; i<=100000; i++)if(!b[i])zhi[++zhi[0]]=i;
    n=read();
    m=read();
    for(i=1; i<=zhi[0]; i++) {
        for(j=max(n/zhi[i],1); j<=m/zhi[i]; j++) {
            if(j==1)continue;
            if(j*zhi[i]<n)continue;
            z[zhi[i]*j-n]=1;                   //平移數組,控制空間大小
        }
    }
    for(i=n-n; i<=m-n; i++) {
        if(!z[i])++ans;
    }
    printf("%d\n",ans);
    return 0;
}
最后編輯于
?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發布,文章內容僅代表作者本人觀點,簡書系信息發布平臺,僅提供信息存儲服務。

推薦閱讀更多精彩內容

  • 個人學習批處理的初衷來源于實際工作;在某個迭代版本有個BS(安卓手游模擬器)大需求,從而在測試過程中就重復涉及到...
    Luckykailiu閱讀 4,782評論 0 11
  • 題目描述:讀入N名學生的成績,將獲得某一給定分數的學生人數輸出。 輸入:測試輸入包含若干測試用例,每個測試用例的格...
    我是阿喵醬閱讀 785評論 0 5
  • 生活大爆炸版石頭剪刀布 題目描述 石頭剪刀布是常見的猜拳游戲:石頭勝剪刀,剪刀勝布,布勝石頭。如果兩個人出拳一樣,...
    bbqub閱讀 487評論 0 0
  • 北上廣不相信安全套 我承認,這又是一個該死的標題黨。我有罪…...
    鈞法崇是閱讀 291評論 0 0
  • 因為不知道見了這次什么時候還可以再見,媽媽給我拍了很多照片,有我的,她的,還有她抱著我的。媽媽在鎮上住了兩天就走了...
    一塵不染靜女子閱讀 234評論 0 0