线性筛

线性筛法

#define M 100000 
LL prime[M/3];
bool flag[M];
inline void get_prime()
{
	LL i,j;
      	memset(flag,false,sizeof(flag));
    	for(i=2;i<M;i++){
         	if(!flag[i])                            
         	prime[cnt++]=i;
         	for(j=0;j<cnt&&i*prime[j]<M;j++){
             		flag[i*prime[j]]=true;            
             		if(i%prime[j]==0)             
                 		break;
        	}
     	}
}