n最大为2000000000(不知为什么OJ上是1000),若为判断2000000000是素数,则必有一个素数在sqrt(n)内,求出这个范围 的所有素数,其比最大数据小的n'的sqrt(n')也在这个范围 内。
#include#include #include #include #include using namespace std;const int Max=50000;bool prime[Max+1];int isprime[Max];int main(){ memset(prime,true,sizeof(prime)); prime[0]=prime[1]=false; int e=(int)sqrt(1.0*Max); int cp=0; for(int i=2;i<=e;i++){ if(prime[i]){ for(int j=i+i;j<=Max;j+=i) prime[j]=false; } } for(int i=0;i<=Max;i++) if(prime[i]) isprime[cp++]=i; int n; while(scanf("%d",&n)!=EOF){ if(n==1) { printf("NO\n"); continue; } int c=0; int lim; lim=sqrt(n*1.0); bool flag=true; for(int i=0;i