エラトステネスのふるい

#include <stdio.h>
#include <math.h>

int main(){
    int prm[1000],i,j,k,n;
    for(i=0;i<1000;i++) prm[i]=1;
    i=1;
    n=sqrt(1000.0);
    do{
        if(prm[i]!=0){
            /* エラトステネスのふるい */
            for (j=i+(i+1);j<1000;j=j+(i+1)){prm[j] = 0;}
        }
        i++;
    }while(n >= i);
    k=0;
    for(i=2;i<1000;i++){
        if(prm[i-1]!=0){printf("%3d,",i); k++;}
        if(k==10){printf("\n");k=0;}
    }
    printf("\n");
}