SPOJ_1_PRIME1

USACO快刷完了,從今天開始寫SPOJ的題解了。

這一到題不難。

高斯篩法+埃拉托色尼篩法,先求出1-33000的所有素數,再用這些素數去求得區間內所有素數

/*
 *NAME: prime1
 *LANG: C++
 */
#include <cstdio>
#include <cstdlib>
#include <cstring>
using namespace std;

int prime[7000];
bool check[103000];
int n,m;

int main(){
    //先求出1-33000的所有素數
    memset(check,true,sizeof(check));
    int tot=0;int N=32000;
    for (int i=2;i<=N;++i){
	if (check[i]) prime[tot++]=i;
	for (int j=0;j<tot;++j){
	    if (i*prime[j]>N) break;
	    check[i*prime[j]]=false;
	    if (i%prime[j] == 0)break;
	}
    }
    
    //solve
    int k;scanf("%d",&k);
    for (int t=1;t<=k;++t){
	scanf("%d %d",&n,&m);if (n==1) n=2;
	memset(check,true,sizeof(check));
	for (int i=0;i<tot;++i)
	    for (int j=n/prime[i];j<=m/prime[i];++j)
		if (prime[i]*j>=n && j!=1) check[prime[i]*j-n]=false;
	for (int i=0;i<=m-n;++i) 
	    if (check[i]) printf("%d\n",i+n);
	printf("\n");
    }
    return 0;
}


發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章