C - Coprimes Gym - 101492C

题意：给你一个长度为n的序列，给你m次询问区间 [l,r] ,回答区间内是否有一对互质的数。

题解： 首先素数筛，开素数筛个bitset，bitset的每一位对应序列的每一位。从后往前循环，维护两个bitset，a清空，b每次循环都把后一位从0变为1，然后求当前数的素数因子，然后a把素数因子的bitset或起来，然后与b异或，使用_Find_first(),加1就是当前位的后面最近互质的数的位置，从后向前维护最小值minn数组，如果minn[l]<=r 输出S否则输出N。

代码：

#include <bits/stdc++.h>
using namespace std;
const int maxn=5e5+5;
bitset<maxn/10>p[50000],b,v;
int prime[maxn],vis[maxn];
int a[maxn],pos[maxn],minn[maxn];
void f_prim(){
    prime[0]=0;
    for(int i=2;i<maxn;i++){
        if(!prime[i]){
            prime[++prime[0]]=i;
            for(int j=2;j*i<maxn;j++){
                prime[j*i]=1;
            }
        }
    }
    for(int i=1;i<=prime[0];i++){
        p[i].reset();
        vis[prime[i]]=i;
    }
}

int main()
{
    f_prim();
    int n,m,x;
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
    }
     b.reset();
    for(int i=n;i>=1;i--){
        x=a[i]; v.reset();
        for(int j=1;prime[j]*prime[j]<=x&&j<=prime[0];j++){
            if(x%prime[j]==0){
                v=v|p[j];
                p[j][i-1]=1;
                while(x%prime[j]==0){
                    x/=prime[j];
                }
            }
        }
        if(x>1){
            v=v|p[vis[x]];
            p[vis[x]][i-1]=1;
        }
        b[i]=1;
        v=v^b;
        pos[i]=v._Find_first()+1;
    }
    minn[n+1]=n+1;
    for(int i=n;i>=1;i--){
        minn[i]=min(minn[i+1],pos[i]);
    }
    while(m--){
        int l,r;
        scanf("%d %d",&l,&r);
        if(r>=minn[l]){
            printf("S\n");
        }else{
            printf("N\n");
        }
    }
    return 0;
}