素数判定
对于表达式n^2+n+41,当n在(x,y)范围内取整数值时(包括x,y)(-39<=x<y<=50),判定该表达式的值是否都为素数。
Input
输入数据有多组,每组占一行,由两个整数x,y组成,当x=0,y=0时,表示输入结束,该行不做处理。
Output
对于每个给定范围内的取值,如果表达式的值都为素数,则输出"OK",否则请输出“Sorry”,每组输出占一行。
Sample Input
0 1 0 0
Sample Output
OK
这个题应该算作素数里面最简单的了,但是同样关于素数还是有很多知识。。。(以下内容摘自学长PPT)
如何判断素数呢?
•试除法:给定n,如果它不能整除[2, sqrt(n)]内所有的数,它就是素数。
埃式筛法(Eratosthenes sieve)
n初始队列{2、3,4,5,6,7,8,9,10,11,12,13,...,n},操作步骤:
(1)输出最小的素数2,然后划掉2的倍数,
剩下{3,5,7,9,11,13,...}
(2)输出最小的素数3,然后划掉3的倍数,
剩下{5,7,11,13,...}
(3)输出最小的素数5,然后划掉5的倍数,
剩下{7,11、13,...}
继续以上步骤,直到队列为空。
埃式筛法虽然还不错,但其实做了一些无用功,某个数会被筛到好几次,比如12,被2和3筛了2次。
还有一种欧拉筛选法,自己还没有看..
这题自己用的试除法。也看到过网上有人用欧拉筛法,感觉自己还没达到那个高度..hhh
#include<stdio.h>
#include<cmath>
int main()
{
int x,y,n,N;
while(scanf("%d%d",&x,&y)!=EOF)
{
if(x == 0 && y == 0)
break;
N=1;
for(int i = x;i<=y;i++)
{
n=i*i+i+41;
for(int j = 2; j*j<=n;j++)
if(n%j == 0)
{
N = 0;
break;
}
if(N == 0)break;
}
if(!N)
printf("Sorry\n");
else
printf("OK\n");
}
return 0;
}
代码中for()中,我用的是j*j<n,也是在PPT中看到的提示,相比于开根号,精确度应该更高一点。