求一千以内的孪生素数

这个程序是我自己写的，但可能用的方法回合其他人相符合，但绝对原创，这个是有这自行说的

#include <stdio.h>
#include <math.h>
main()
{
    int i,j;
    printf ("........................The prime number is :.......................\n");
    for (j = 0,i =3;i <= 999;i += 2 )
    if (prime (i))
        if (prime (i+2))
            printf (j%4?"%4d:%4d+%4d":"%4d:%4d+%4d\n",++j,i,i+2);

}
   上面这是主程序

#include <stdio.h>
#include <math.h>
prime (i)
        int i;
    {
        int j;
        if (i <= 1)         return (0);
        if (i == 2)         return (1);
        if (!(i%2))         return (0);
        j = sqrt((double)i) + 1;


        for (j = 3;j <= (int)(sqrt((double)i)+1);j += 2)
            if (!(i%j))     return (0);
        return (1);


    }

这是子程序，也是用来判别是否为素数的程序