正文

洛谷p3383模板线性筛素数(miller_rabin)

自为  自为  2022-10-07  153

关键词：

题目描述

如题，给定一个范围N，你需要处理M个某数字是否为质数的询问（每个数字均在范围1-N内）

输入输出格式

输入格式：

第一行包含两个正整数N、M，分别表示查询的范围和查询的个数。

接下来M行每行包含一个不小于1且不大于N的整数，即询问该数是否为质数。

输出格式：

输出包含M行，每行为Yes或No，即依次为每一个询问的结果。

输入输出样例

输入样例#1：复制

输出样例#1：复制

Yes
Yes
No
No
Yes

说明

时空限制：500ms 128M

数据规模：

对于30%的数据：N<=10000，M<=10000

对于100%的数据：N<=10000000，M<=100000

样例说明：

N=100，说明接下来的询问数均不大于100且不小于1。

所以2、3、97为质数，4、91非质数。

故依次输出Yes、Yes、No、No、Yes。

补一发Miller Rabin

注意开long long

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#define LL long long 
using namespace std;
const LL MAXN=2*1e7+10;
const LL INF=1e7+10;
inline char nc()
{
    static char buf[MAXN],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,MAXN,stdin),p1==p2)?EOF:*p1++;
}
inline LL read()
{
    char c=nc();LL x=0,f=1;
    while(c<‘0‘||c>‘9‘){if(c==‘-‘)f=-1;c=nc();}
    while(c>=‘0‘&&c<=‘9‘){x=x*10+c-‘0‘;c=nc();}
    return x*f;
}
LL fastpow(LL a,LL p,LL mod)
{
    LL base=1;
    while(p)
    {
        if(p&1) base=(base*a)%mod;
        a=(a*a)%mod;
        p>>=1;
    }
    return base;
}
LL num[]= {2,3,5,7,11,13,17,19};
bool Miller_Rabin(LL n)
{
    if (n==2) return 1;
    if((n&1)==0||n==1) return false;
    for (LL i=0; i<8; i++) if (n==num[i]) return 1;
    
    LL temp=n-1,t=0,nxt;
    while((temp&1)==0) temp>>=1,t++;
    
    for(LL i=0;i<8;i++)
    {
        LL a=num[i];
        LL now=fastpow(a,temp,n);
        nxt=now;
        for(LL j=1;j<=t;j++)
        {
            nxt=(now*now)%n;
            if(nxt==1&&now!=n-1&&now!=1) return false;
            now=nxt;
        }
        if(now!=1) return false;
    }
    return true;
}
int main()
{
    #ifdef WIN32
    freopen("a.in","r",stdin);
    #else
    #endif 
    LL N=read(),M=read();
    while(M--)
    {
        LL opt=read();
        if(Miller_Rabin(opt))    printf("Yes
");
        else printf("No
");
    }
    return 0;
}

线性筛洛谷p3383线性筛模板

思路：如果我们要筛出[1,n]内的所有素数，使用[1,√n]内的素数去筛就可以了设bool型数组a，a[i]表示i是否被某个素数筛过从2开始枚举每个数i：若a[i]=false，表示i没有更小的素因子，从而知道i是素数。枚举i的所有倍数j，令a[j]=1这... 查看详情

洛谷p3383模板线性筛素数

题目链接：https://www.luogu.org/problem/show?pid=3383题目描述如题，给定一个范围N，你需要处理M个某数字是否为质数的询问（每个数字均在范围1-N内）输入输出格式输入格式：第一行包含两个正整数N、M，分别表示查询的范围和查询的... 查看详情

p3383模板线性筛素数洛谷

https://www.luogu.org/problem/show?pid=3383#sub题目描述如题，给定一个范围N，你需要处理M个某数字是否为质数的询问（每个数字均在范围1-N内）输入输出格式输入格式：第一行包含两个正整数N、M，分别表示查询的范围和查询的个数... 查看详情

洛谷p3383模板线性筛素数

题目描述如题，给定一个范围N，你需要处理M个某数字是否为质数的询问（每个数字均在范围1-N内）输入输出格式输入格式：第一行包含两个正整数N、M，分别表示查询的范围和查询的个数。接下来M行每行包含一个不小于1且... 查看详情

洛谷p3383模板线性筛素数

洛谷p3383模板线性筛素数(miller_rabin)

[洛谷p3383]线性筛素数-欧拉筛法

Description如题，给定一个范围N，你需要处理M个某数字是否为质数的询问（每个数字均在范围1-N内）Input&OutputInput第一行包含两个正整数N、M，分别表示查询的范围和查询的个数。接下来M行每行包含一个不小于1且不大于N的整数... 查看详情

p3383模板线性筛素数(代码片段)

P3383【模板】线性筛素数欧拉筛O(n)#include<iostream>#include<cstdio>usingnamespacestd;intn,m,cnt,prime[10000002],v[10000002];//prime:素数表v:存某数的最小质因数intmain()scanf("%d%d",&n,&m);for(inti=2;i< 查看详情

p3383模板线性筛素数

题目描述如题，给定一个范围N，你需要处理M个某数字是否为质数的询问（每个数字均在范围1-N内）输入输出格式输入格式：第一行包含两个正整数N、M，分别表示查询的范围和查询的个数。接下来M行每行包含一个不小于1且不大... 查看详情

p3383模板线性筛素数

p3383模板线性筛素数(代码片段)

//P3383【模板】线性筛素数#include<bits/stdc++.h>usingnamespacestd;intis_prime[10000005];voidFind_prime(intn)memset(is_prime,1,sizeof(is_prime));is_prime[1]=0;for(inti=2;i*i<=n;i++)if(is_prime[i])for(intj=i*i;j<=n;j+=i)is_prime[j]=0;intmain()intn,m;scanf("%d%d",&n,&m);Fin... 查看详情

p3383模板线性筛素数(代码片段)

P3383【模板】线性筛素数埃氏筛->欧拉筛普通埃氏筛（O（nlognlogn））for(inti=2;i<=n;i++)//注意终止条件if(!notpri[i])for（intj=2;j*i<=n;j++）notpri[i*j]=true;优化for(inti=2;i<=n;i++)//到根号if(!notpri[i])for（intj=i;j*i<=n;j++）//j从i开始，因为... 查看详情

p3383模板线性筛素数(代码片段)

普及组模板——线性筛素数

题目：【模板】线性筛素数(洛谷_3383)#include<iostream>#include<cstdio>#include<algorithm>#include<cstring>#include<cmath>usingnamespacestd;inlineintread(){intt=1,num=0;charc=getchar();w 查看详情

（模板）线性筛素数

————————————————————————————&mda 查看详情

线性筛素数模板

传送门：线性筛素数 Prime：1#include<cstdio>23constintMAXN=10000100;4intPrime[MAXN],n,m,Size;5boolVis[MAXN]={1,1};67intmain()8{9scanf("%d%d",&n,&m);10for(inti=2;i<n;i++)11{12if(!Vis[i])13Pr 查看详情

模板线性筛素数

#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>usingnamespacestd;intn,m,vis[10000010],prime[5000010],tot;voidinit(){for(inti=2;i<=n;i++){if(!vis[i])pr 查看详情

模板线性筛求素数(代码片段)

线性筛求素数1#include<iostream>2#include<cstring>3#include<algorithm>4#include<cstdio>56constintMAX_N=10000000+10;78boolvis[MAX_N];9intPrime[MAX_N>>1];10inlinevoidGetPrime(intn 查看详情