L1-028. 判断素数（神奇的判定素数方法）

本题的目标很简单，就是判断一个给定的正整数是否素数。

输入格式：

输入在第一行给出一个正整数N（<=10），随后N行，每行给出一个小于2³¹的需要判断的正整数。

输出格式：

对每个需要判断的正整数，如果它是素数，则在一行中输出“Yes”，否则输出“No”。

输入样例：

2
11
111

输出样例：

Yes
No

解题思路：

很神奇的判断素数的方法：http://blog.csdn.net/huang_miao_xin/article/details/51331710

也就是方法3

下面讲一下这种更快速的判断方法；
首先看一个关于质数分布的规律：大于等于5的质数一定和6的倍数相邻。例如5和7，11和13,17和19等等；

证明：令x≥1，将大于等于5的自然数表示如下：

······ 6x-1，6x，6x+1，6x+2，6x+3，6x+4，6x+5，6(x+1），6(x+1)+1 ······

可以看到，不在6的倍数两侧，即6x两侧的数为6x+2，6x+3，6x+4，由于2(3x+1)，3(2x+1)，2(3x+2)，所以它们一定不是素数，再除去6x本身，显然，素数要出现只可能出现在6x的相邻两侧。这里有个题外话，关于孪生素数，有兴趣的道友可以再另行了解一下，由于与我们主题无关，暂且跳过。这里要注意的一点是，在6的倍数相邻两侧并不是一定就是质数。

根据以上规律，判断质数可以6个为单元快进，即将方法（2）循环中i++步长加大为6，加快判断速度，代码如下：

int isp(long long num)  
{  
    if(num<=1) return 0;
    if(num==2||num==3) return 1;//两个较小数另外处理 
    if(num%6!=1&&num%6!=5) return 0;//不在6的倍数两侧的一定不是质数   
    //在6的倍数两侧的也可能不是质数  
    for(long long i=5;i*i<=num;i+=6)
        if(num%i==0||num%(i+2)==0) return 0;
    return 1;//排除所有，剩余的是质数 
} 

代码稍作修改，补充上判断0和1的情况。

AC代码：

#include<stdio.h>

int isp(long long num)  
{  
	if(num<=1) return 0;
    if(num==2||num==3) return 1;//两个较小数另外处理 
    if(num%6!=1&&num%6!=5) return 0;//不在6的倍数两侧的一定不是质数   
    //在6的倍数两侧的也可能不是质数  
    for(long long i=5;i*i<=num;i+=6)
		if(num%i==0||num%(i+2)==0) return 0;
    return 1;//排除所有，剩余的是质数 
}  

int main()
{
	int t;
	long long n;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%lld",&n);
		printf("%s\n",isp(n)==1?"Yes":"No");
	}
	return 0;
}

方法2的AC代码：

（在数据量较少的情况下体现不出方法3的效果，但是方法3把步长从1提高到6，肯定是有优势的）

#include<stdio.h>
int isp(long long n)
{
	if(n<=1) return 0;
	else
	{
		for(long long i=2;i*i<=n;i++)
		{
			if(n%i==0) return 0;
		}
	}
	return 1;
}

int main()
{
	int t;
	long long n;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%lld",&n);
		printf("%s\n",isp(n)==1?"Yes":"No");
	}
	return 0;
}