最好的输出一个范围内素数的方法，你有更好的吗?

#include "stdafx.h"
#include "math.h"
#include <iostream>
using namespace std;
#define M 100000
int is(int n)
{
if (2==n)
{return 0;
}
double N=sqrt((double)(n+1));
for (int i=3;i<N;i+=2)
if (n%i==0)
{
return 0;
}
return 1;

}

int _tmain(int argc, _TCHAR* argv[])
{
for (int j=1;j<=M;j+=2)
{
if (is(j))
{
cout<<j<<',';
}
}
return 0;
}