和Leo一起做爱线段树的好孩子之 HDU 4630 No Pain No Game

神仙线段树系列

陈老师的确出的题个个牛皮

首先我们知道如果一个因数在一段区间出现了两次那么这个因数是一个合法的GCD

注意n不是1e5说明可以暴力查找每个因数

然后维护一个当前因数的上一次出现位置

如果出现了那对上一次出现位置来说是合法的

当然这只是左端点

于是对右端点排序

离线处理

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int N=3e5+100;
#define lc (p<<1)
#define rc (p<<1|1)
inline void read(int &x){
	x=0;
	int f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){
		if(ch=='-')f=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9'){
		x=x*10+ch-'0';
		ch=getchar();
	}
	x*=f;
}
struct Node{
	int l,r,Id;
}A[N];
bool cmp (Node A,Node B){
	return A.r<B.r;
}
int sum[N]={};
int cnt=0;
void PushUp(int p){
	sum[p]=max(sum[lc],sum[rc]);
}
void Update(int p,int l,int r,int pos,int val){
	if(l==r){
		sum[p]=max(sum[p],val);	
		return;
	}
	int mid=(l+r)/2;
	if(pos<=mid)Update(lc,l,mid,pos,val);
	else Update(rc,mid+1,r,pos,val);
	PushUp(p);	
}
int Query(int p,int l,int r,int dl,int dr){
	if(dl<=l&&r<=dr){
		return sum[p];
	}
	int mid=(l+r)/2;
	int ret=0;
	if(dl<=mid)ret=max(ret,Query(lc,l,mid,dl,dr));
	if(mid< dr)ret=max(ret,Query(rc,mid+1,r,dl,dr));
	return ret;
}
int Pre[N]={};
int val[N]={};
int ans[N]={};
int n;
void Clr(){
	memset(Pre,0,sizeof(Pre));
	memset(sum,0,sizeof(sum));
}
int main(){
//	freopen("test.in","r",stdin);
	int Cas;
	read(Cas);
	while(Cas--){
		Clr();
		read(n);
		for(int i=1;i<=n;++i){
			read(val[i]);
		}
		int Q;
		read(Q);
		for(int i=1;i<=Q;++i){
			read(A[i].l);
			read(A[i].r);
			A[i].Id=i;
		}
		sort(A+1,A+1+Q,cmp);
		int now=1;
		for(int i=1;i<=n;++i){
			for(int j=1;j*j<=val[i];++j){
				if(val[i]%j==0){
					if(Pre[j]){
						Update(1,1,n,Pre[j],j);
					}
					if(j*j!=val[i]&&Pre[val[i]/j]){
						Update(1,1,n,Pre[val[i]/j],val[i]/j);
					}
					Pre[j]=Pre[val[i]/j]=i;					
				}
			}	
			while(A[now].r==i){
				ans[A[now].Id]=Query(1,1,n,A[now].l,i);
				now++;
			}		
		}
		for(int i=1;i<=Q;++i)cout<<ans[i]<<'\n';
	}
} 

 

转载于:https://www.cnblogs.com/Leo-JAM/p/10079133.html