Gym - 101492C Coprimes [bitset+质因素分解]

最新推荐文章于 2019-08-22 08:37:13 发布

ACTerminate

最新推荐文章于 2019-08-22 08:37:13 发布

阅读量493

点赞数

分类专栏： Codeforces bitset 文章标签： bitset 分解素因子

本文链接：https://blog.csdn.net/ACTerminate/article/details/79264341

版权

Codeforces 同时被 2 个专栏收录

76 篇文章 0 订阅

订阅专栏

bitset

3 篇文章 0 订阅

订阅专栏

题意：给你n个数，m个询问，每次询问[L,R]，若L，R中有互质的数输出S,否则输出N。

题解：对每个数质因数分解，然后bitset维护每个素因子包含的位置，通过位操作与_Fnd_first,获取每个位置的数右边最近的互质数的位置，然后O（1）判断询问。

AC代码：

#include<stdio.h>
#include<bitset>
#include<string.h>
using namespace std;
bitset<50005>bt[42005],ret,an,OK;
int a[50005],prim[42005],pos[500005],Right[50005],dp[500005],tot;
bool f[500005];
int main()
{
	memset(f,true,sizeof(f));
	f[0]=f[1]=false;
	for(int i=2;i<=500000;i++)
		if(f[i])
		{
			prim[tot++]=i;
			for(int j=i*2;j<=500000;j+=i)
				f[j]=false;
		}
	for(int i=0;i<tot;i++)
		pos[prim[i]]=i;
	int n,m;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		scanf("%d",&a[i]);
	dp[n+1]=50005;
	for(int i=n;i>=1;i--)
	{
		int now=a[i];
		OK[i]=1;
		ret=OK;an.reset();
		for(int j=0;prim[j]*prim[j]<=now;j++)
			if(now%prim[j]==0)
			{
				bt[pos[prim[j]]][i]=1;
				an|=bt[pos[prim[j]]];
				while(now%prim[j]==0)now/=prim[j];
			}
		if(now!=1)bt[pos[now]][i]=1,an|=bt[pos[now]];
		Right[i]=(ret^an)._Find_first();
		dp[i]=Right[i];
		dp[i]=min(dp[i],dp[i+1]);
	}
	while(m--)
	{
		int l,r;
		scanf("%d%d",&l,&r);
		if(dp[l]<=r)printf("S\n");
		else printf("N\n");
	}
}