void init()
{
    int m=sqrt(maxn+0.5);
    memset(vis,0,sizeof(vis));
    for(int i=2;i<=m;i++)
    {
        if(!vis[i])
        for(int j=i*i;j<=maxn;j+=i)
        {
            vis[j]=1;
        }
    }
}