PAT 1007. 素数对猜想经典方法，与少见方法之我见

1007. 素数对猜想 (20)

时间限制

400 ms

内存限制

65536 kB

代码长度限制

8000 B

判题程序

Standard

作者

CHEN, Yue

让我们定义 d_n 为：d_n = p_n+1 - p_n，其中 p_i 是第i个素数。显然有 d₁=1 且对于n>1有 d_n 是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。

现给定任意正整数N (< 10⁵)，请计算不超过N的满足猜想的素数对的个数。

输入格式：每个测试输入包含1个测试用例，给出正整数N。

输出格式：每个测试用例的输出占一行，不超过N的满足猜想的素数对的个数。

输入样例：

20

输出样例：

4

方法一：我原来的方法

#include<stdio.h>
#include<math.h>
int pansu(int n);
int main()
{
	int k=0;
	int a[100000];
	int i,n,N=0;
	scanf("%d",&n);
	for(i=1;i<=n;i++)
	{
		if(pansu(i)==1)
		{
			a[k]=i;
			k++;
		}
			
	}
	for(i=0;i<k-1;i++)
	{
		if(a[i+1]-a[i]==2)
		N++;
	}
	printf("%d",N);
	
 } 
 int pansu(int n)
 {
 	int i;
 	for(i=2;i<=sqrt(n);i++)
 	{
 		if(n%i==0)
 		{
 			return 0;
 			break;
		 }
	 }
	 if(i>=sqrt(n))
	 {
	 	return 1;
	 }
 }

这个方法一开始在最后一个测试点会过不去，显示段错误，是内存超限的原因，于是我将数组大小从1000扩大到了100000于是就可以满分了

方法二，网上大神的方法，大概意思是向素数数组中存储100个素数，然后判断方法是这个数能否被数组下标小于等于这个数字开根号的数组中的数值整除，如果能就不是素数。下面码上砖

#include <stdio.h>
int main()
{
    int N;
    scanf("%d", &N);
    /* record three successive numbers if they are prime numbers, start from 2, 3, 4 */
    int iPrimeMinus2 = 1, iPrimeMinus1 = 1, iPrime;
    /* record the prime numbers before sqrt(10^5) */
    int primes[100] = {2, 3};
    int paircount = 0;
    int primecount = 2;
    
    for (int i = 4; i <= N; i++)
    {
        iPrime = 1;
        /* test if i is a prime number */
        for(int j = 0; primes[j] * primes[j] <= i; j++) if(i % primes[j] == 0)
        {
            iPrime = 0;
            break;
        }
        /* i is a prime number */
        if(iPrime == 1)
        {
            if(primecount < 100)    primes[primecount++] = i;
            if(iPrimeMinus2 == 1)   paircount++;    /* a prime pair found */
        }
        /* change the flags */
        iPrimeMinus2 = iPrimeMinus1;
        iPrimeMinus1 = iPrime;
    }
    printf("%d", paircount);
    
    return 0;
}

这两个方法都可以得到满分，在10000的测试下得到的结果相同都是1224，但是还是有些不同就是在100,0000数据时，方法一得到的结果是8169，而方法二得到的结果是6195，方法二的结果偏下，但是测试点应该在10,0000以内，至于这两个程序到底哪个跟准确，还有待考证，希望大神们能多给我一些指点^_^