Coprimes（从 i 位置往后第一个与其互质的数的位置）

 bitset维护每个质因子出现的位置

从后往前跑，每次找到一个位置使该数的所有质因子都没有出现过。

把该数的所有质因子的bitset或起来，然后找从当前位置往后的第一个0位置。

取反后用bitset的_Find_first函数可以解决。

 

#include <bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i <= (b); i++)
#define per(i, a, b) for(int i = (a); i >= (b); i--)
using namespace std;
const int N = 51000;
int n,m;

int prime[10*N]; //素数表下标从1开始
bool mark[10*N]; //mark[i]==1表示i不是素数
int tot;
void get_list(int n){
    for(int i = 2;i <= n;i ++){
        if(!mark[i]) prime[++tot] = i;
        for(int j = 1;j<=tot&&prime[j]*i<=n;j ++){
            mark[i*prime[j]] = 1;
            if(!(i%prime[j])) break;
        }
    }
}
int s[N],dp[N];
bitset<51234>bt[42000],re;
int main() {
    //freopen("a.txt","r",stdin);
    ios::sync_with_stdio(0);
    get_list(500000);
    cin>>n>>m;
    rep(i, 1, n) cin>>s[i];
    dp[n+1] = n+1;
    rep(i, 0, tot) bt[i][n+1] = 1;
    re[n+1] = 1;
    per(i, n, 1) {
        int x = s[i];
        bitset<51234> sta;
        for(int j = 1; prime[j]*prime[j]<=x; j++){
            int p = prime[j];
            if(x%p==0) sta |= bt[j],bt[j][i] = 1;
            while(x%p==0) x/=p;
        }
        if(x>1) rep(j, 1, tot) if(prime[j]==x) {sta|=bt[j];bt[j][i]=1;break;}
        sta^=re;
        re[i] = 1;
        dp[i] = sta._Find_first();
        dp[i] = min(dp[i],dp[i+1]);
    }
    rep(i, 1, m){
        int l,r;
        cin>>l>>r;
        if(dp[l]<=r) cout<<"S"<<endl;
        else cout<<"N"<<endl;
    }
    return 0;
}