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

1007. 素数对猜想 (20)

时间限制
400 ms
内存限制
65536 kB
代码长度限制
8000 B
判题程序
Standard
作者
CHEN, Yue

让我们定义 dn 为:dn = pn+1 - pn,其中 pi 是第i个素数。显然有 d1=1 且对于n>1有 dn 是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。

现给定任意正整数N (< 105),请计算不超过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以内,至于这两个程序到底哪个跟准确,还有待考证,希望大神们能多给我一些指点^_^
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值