[NWERC 2017]Factor-Free Tree

最新推荐文章于 2022-07-03 21:03:51 发布

StaroForgin

最新推荐文章于 2022-07-03 21:03:51 发布

阅读量241

点赞数 1

分类专栏： # 树 ------构造------

本文链接：https://blog.csdn.net/Tan_tan_tann/article/details/113828277

版权

树同时被 2 个专栏收录

57 篇文章 0 订阅

订阅专栏

------构造------

33 篇文章 0 订阅

订阅专栏

这是一篇关于算法的博客，主要讨论如何找到区间内的根节点，使得区间内的每个数与根节点互质。通过质因数分解预处理，可以有效地判断数之间的互质关系，并采用双指针方法确保线性时间复杂度。文章提供了详细的代码实现。

摘要由CSDN通过智能技术生成

Factor-Free Tree

题解

挺有趣的一道题。

我们首先应该很容易发现，对于一段区间，这段区间中任意一个与区间中其他数互质的数都可以作为这段区间的根。
因为如果这个数做根时会使其子区间无解，即使其他数做了根这个子区间无解的地方依旧会被传递下去，最后也一定会产生无解的状况。
所以，对于一段区间，我们如果找到了一个与该区间内其他数互质的数，我们可以直接将其提出来做根。

对于如何判断一个数是否与一段区间内其他数互质，我们可以先预处理出这个数所有的质因数，很明显，对于一个小于 $10^7$ 的数，不同的质因数不会超过10个。
然后我们往前往后找第一个与其有相同质因数，即第一个不与它互质的数。
很明显，如果当前区间在找出来的这个区间中，那就可以。
至于这个区间，我们可以将其预处理出来。

注意扫的时候我们不能从一边扫过去，要从两边扫向中间，这样才能保证时间复杂度严格 $O(nlog\, n)$ 。
如果从一边扫的话是可以被卡到 $O\left(n^2\right)$ 的。

源码

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
#define MAXN 1000005
#define reg register
typedef long long LL;
typedef unsigned long long uLL;
typedef pair<int,int> pii;
const int INF=0x7f7f7f7f;
template<typename _T>
_T Fabs(_T x){return x<0?-x:x;}
template<typename _T>
void read(_T &x){
	_T f=1;x=0;char s=getchar();
	while(s>'9'||s<'0'){if(s=='-')f=-1;s=getchar();}
	while('0'<=s&&s<='9'){x=(x<<3)+(x<<1)+(s^48);s=getchar();}
	x*=f;
}
int n,prime[MAXN],cntp,pre[MAXN*10],fa[MAXN],idx,ls[MAXN],rs[MAXN],vis[MAXN];
bool oula[MAXN*10];
struct Number{
	int x,sta[10],stak;
	void work(){while(x>1){if(sta[stak]!=pre[x])sta[++stak]=pre[x];x/=prime[pre[x]];}}
}a[MAXN];
void init(){
	for(int i=2;i<=1e7;i++){
		if(!oula[i])prime[++cntp]=i,pre[i]=cntp;
		for(int j=1;j<=cntp&&1ll*i*prime[j]<=1e7;j++){
			oula[i*prime[j]]=1;pre[i*prime[j]]=j;
			if(i%prime[j]==0)continue;
		}
	}
}
void sakura(int l,int r,int father){
	int L=l,R=r;if(l==r){fa[l]=father;return ;}if(l>r)return ;
	while(L<=R){
		if(ls[L]<l&&rs[L]>r){sakura(l,L-1,L);sakura(L+1,r,L);fa[L]=father;return ;}L++;
		if(ls[R]<l&&rs[R]>r){sakura(R+1,r,R);sakura(l,R-1,R);fa[R]=father;return ;}R--;
	}
	puts("impossible");exit(0);
}
signed main(){
	read(n);init();for(int i=1;i<=n;i++)read(a[i].x),a[i].work(),rs[i]=n+1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=a[i].stak;j++)
			ls[i]=max(ls[i],vis[a[i].sta[j]]),vis[a[i].sta[j]]=i;
	for(int i=1;i<=cntp;i++)vis[i]=n+1;
	for(int i=n;i>0;i--)
		for(int j=1;j<=a[i].stak;j++)
			rs[i]=min(rs[i],vis[a[i].sta[j]]),vis[a[i].sta[j]]=i;
	sakura(1,n,0);for(int i=1;i<=n;i++)printf("%d ",fa[i]);
	return 0;
}