ACM之旅——HDU2012

素数判定

 

对于表达式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中看到的提示，相比于开根号，精确度应该更高一点。