素数打表

-素数打表

最初学c语言的时候，我们学了两种求素数的方法：
例如：从键盘上输入一个大于1的正整数m，判断它是不是素数,是输出yes，不是输出no
方法1：

for(i=2;i<=m-1;i++)
    if(m%i==0) 
           break;
if(i==m)
    printf("yes\n")；
else
    printf("no\n");

方法2：

k=sqrt(m);
for(i=2;i<=k;i++)
    if(m%i==0)
            break;
if(i==k+1)
    printf("yes\n")；
else
    printf("no\n");

由上述我们可以知道如果m非常大的话，方法2明显要比方法1少循环了很多次，但是仍然不是很理想。于是就出现了更高效的素数打表方法。
假设输入连续输入多个数x（每个数据一行）（1<a<=1000000）判断其是否为素数部分代码如下：

for(j=2;j<=1000000;j++)
    if(a[j]==0)
            for(i=j+j;i<=1000000;i+=j)
                        a[i]=1;

定义一个数组a[1000001]（清零），由上述可知从2开始循环2，3，4，5，6，7，8，9，10…，那么a[2]=0,进入到下一个for循环里，执行后，a[4]=1,a[6]=1,a[8]=1…也就是说，我们通过这一次循环，标记了1000000以内所有2的倍数，之后的3，5，7也是同理，所以判断一个数是不是素数，只需要判断a[x]的值是0还是1。
进一步分析上面这个算法，你有没有发现还是有很多重复的情况，比如说10这个数，它是2的倍数，也是5的倍数，当j等于2的时候，把它标记了一次，当j等于5的时候，又把它标记了一次，所以，上面这个算法，我们可以进一步优化：

for(j=2;j<=1000000;j++)
    if(a[j]==0)
            for(i=j*j;i<=1000000;i+=j)
                        a[i]=1;

把i=j+j改成i=j*j，又会少了很多循环的次数，还是刚才10那个数，现在这个代码，当j等于5的时候，就不会再把它标记一次了，下一重循环直接从25开始。
如有错误，欢迎指正。