杭电 1262 寻找素数对 【素数】

#include <stdio.h>
#include <math.h>
#define N 10000
 int p[N];   //判断p[i]中的i是否为素数
int pn[N];   //存储素数

void prim_num()
{
    int i,j,n;
    for(i=1; i<=N; i++)
        p[i]=1;
    n=(int)sqrt(N);
    for(i=2; i<=n; i++)
    {
        for(j=i+i; j<=N; j+=i)   //素数的合肯定不是素数，这就是判断哪些数不是素数
        {
            p[j]=0;
        }
    }
    j=1;
    for(i=1; i<=N; i++)  //把素数存储入pn[N],下表从1开始
    {
        if(p[i])
        {
            pn[j++]=i;
        }
    }
}

int main()
{
    int i,j,k;
    int n;
    prim_num();
    while(scanf("%d",&n)!=EOF)
    {
        int flag=1;
        for(i=2;i<N;i++)
        {
            if(n/2==pn[i])
            {
                k=i;
                break;
            }
            if(n/2>pn[i]&&n/2<pn[i+1])
            {
                k=i;
                break;
            }

        }
        for(i=k;i>=0&&flag;i--)
        {
            for(j=i;j<N;j++)
            {
                if(pn[i]+pn[j]==n)
                {
                    printf("%d %d\n",pn[i],pn[j]);
                    flag=0;

                }
            }
        }
    }
    return 0;
}