HDU5104 Primes Problem && BestCoder Round #18 1001

首先将10000内的素数筛出来(约1000个），(p1,p2,p3)枚举三元组前两个p1,p2，可知若存在p3满足条件，必有p3=n−p1−p2，故令t=n−p1−p2。若t为不小于p2的素数，则t满足p3的条件。则答案加一。

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;

const int maxn=10010;
int pr[maxn]={0};
int prime[1500];

int main()
{
    int i,j,t;
    t=sqrt(maxn*1.0)+1;
    for(i=2;i<=t;i++)
        for(j=2;j<=maxn/i;j++)
            pr[j*i]=1;
    j=0;
    for(i=2;i<=maxn;i++)
        if(pr[i]==0)
            prime[j++]=i;
    //cout<<j<<endl;
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        int p,q,ans=0;
        double m=n/3;
        for(i=0;i<=j;i++)
        {
            if(prime[i]>m)
                break;
            for(p=i;p<=j;p++)
            {
                if(prime[p]>=n)
                    break;
                q=n-prime[i]-prime[p];
                if(q<prime[p])
                    continue;
                if(pr[q]==0)
                    ans++;//cout<<prime[i]<<endl;cout<<prime[p]<<endl;}
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}