Isprime

筛法的思想是去除要求范围内所有的合数，剩下的就是素数了，而任何合数都可以表示为素数的乘积，因此如果已知一个数为素数，则它的倍数都为合数。

1. #include"cstdio"  
2. #include"cstring"  
3. using namespace std;  
4. #define MAX 100000//求MAX范围内的素数  
5. long long su[MAX],cnt;  
6. bool isprime[MAX];  
7. void prime()  
8. {  
9.     cnt=1;  
10.     memset(isprime,1,sizeof(isprime));//初始化认为所有数都为素数  
11.     isprime[0]=isprime[1]=0;//0和1不是素数  
12.     for(long long i=2;i<=MAX;i++)  
13.     {  
14.         if(isprime[i])  
15.             su[cnt++]=i;//保存素数i  
16.         for(long long j=1;j<cnt&&su[j]*i<MAX;j++)  
17.         {  
18.             isprime[su[j]*i]=0;//筛掉小于等于i的素数和i的积构成的合数  
19.         }  
20.     }  
21. }  
22. int main()  
23. {  
24.     prime();  
25.     for(long long i=1;i<cnt;i++)  
26.         printf("%d  ",su[i]);  
27.     return 0;  
28. }