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;
}