欧拉筛

来源:互联网 发布:网络连接异常 编辑:程序博客网 时间:2024/06/11 18:22
#define mx 10000005#define mod 1000000007#define ll long longusing namespace std;bool isprime[mx];int prime[mx],cnt;int ans[mx];void get_prime (){    mem(isprime);    mem(prime);    cnt=0;    for(int i=2;i<mx;++i)     {        if(!isprime[i]) prime[cnt++] = i;        for(int j=0;j<cnt&&i*prime[j]<mx;++j)        {           isprime[i*prime[j]] = true;           if(!i%prime[j]) break;        }     }}
原创粉丝点击