Coprimes Gym - 101492C (bitset)

题目

https://cn.vjudge.net/problem/Gym-101492C

题意

给你n个数 每次询问一个区间内是否存在互质的数

思路

记录i位置后面第一个与其互质的数的位置

bitset对应每个素数出现的位置

异或后即这些素数出现的位置 即都不互质

再异或 即可得到互质的位置 bitset_Find_frist 查找第一个可行的位置

代码

#include <bits/stdc++.h>
using namespace std;
const int maxn = 5e4+100;
typedef long long ll;

int dp[maxn*10],a[maxn];
int prime[maxn*10],p[maxn*10],isprime[maxn*10];
bitset<maxn>dap[maxn],ok,ret;
int cnt;
int n,m;
void getpri()
{
    cnt = 0;
    for(int i = 2;i < maxn*10;i++)
    {
        if(isprime[i] == 0)
        {
            p[i] = cnt;
            prime[cnt++] = i;
        }
        for(int j = 0;j < cnt&&i*prime[j]<maxn*10;j++)
        {
            isprime[i*prime[j]] = 1;
            if(i%prime[j]==0)
                break;
        }
    }
}


int main()
{
    ios::sync_with_stdio(false);
    getpri();
    cin>>n>>m;
    dp[n+1] = n+1;
    for(int i = 1;i <= n;i++)
        cin>>a[i];
    for(int i = n;i >= 1;i--)
    {
        ok[i] = 1;
        ret.reset(); //
        int x = a[i];
        for(int j = 0;j < cnt&&prime[j]*prime[j] <= x;j++)
        {
            if(x % prime[j] == 0)
            {
                dap[p[prime[j]]][i] = 1;
                ret|=dap[p[prime[j]]];
                while(x%prime[j]==0)
                    x /= prime[j];
            }
        }
        if(x != 1)
            dap[p[x]][i] = 1,ret|=dap[p[x]];
        dp[i] = (ret^ok)._Find_first();
        dp[i] = min(dp[i],dp[i+1]);
    }
    while(m--)
    {
        int l,r;
        cin>>l>>r;
        if(dp[l] <= r)
            cout<<"S"<<endl;
        else cout<<"N"<<endl;
    }
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值