[CC-ANUGCD]Maximum number, GCD condition

[CC-ANUGCD]Maximum number, GCD condition

题目大意：

一个长度为\(n(n\le10^5)\)的数列\(A(A_i\le10^5)\)，\(m(m\le10^5)\)次询问，每次询问\(l\sim r\)中不与\(g\)互质的数中的最大数以及最大数的个数。

思路：

对于每个质数维护一棵线段树，记录区间内包含这个质数的数的和。询问时将\(g\)分解质因数，在线段树上寻找最大值。

统计个数时将所有数以数值为第一关键字、下标为第二关键字排序后二分即可。

源代码：

#include<cstdio>
#include<cctype>
#include<algorithm>
inline int getint() {
	register char ch;
	while(!isdigit(ch=getchar()));
	register int x=ch^'0';
	while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
	return x;
}
const int N=1e5+1,P=9593;
int p[P];
bool vis[N];
std::pair<int,int> v[N];
inline void sieve() {
	for(register int i=2;i<N;i++) {
		if(!vis[i]) p[++p[0]]=i;
		for(register int j=1;j<=p[0]&&i*p[j]<N;j++) {
			vis[i*p[j]]=true;
			if(i%p[j]==0) break;
		}
	}
}
class SegmentTree {
	#define mid ((b+e)>>1)
	private:
		struct Node {
			int val;
			Node *left,*right;
			Node() {
				val=-1;
				left=right=NULL;
			}
		};
	public:
		Node *root;
		void modify(Node* &p,const int &b,const int &e,const int &x,const int &y) {
			if(!p) p=new Node();
			p->val=std::max(p->val,y);
			if(b==e) return;
			if(x<=mid) modify(p->left,b,mid,x,y);
			if(x>mid) modify(p->right,mid+1,e,x,y);
		}
		int query(const Node* const &p,const int &b,const int &e,const int &l,const int &r) const {
			if(!p) return -1;
			if(b==l&&e==r) return p->val;
			int ret=-1;
			if(l<=mid) ret=std::max(ret,query(p->left,b,mid,l,std::min(mid,r)));
			if(r>mid) ret=std::max(ret,query(p->right,mid+1,e,std::max(mid+1,l),r));
			return ret;
		}
	#undef mid
};
SegmentTree t[P];
int main() {
	sieve();
	const int n=getint(),m=getint();
	for(register int i=1,b;i<=n;i++) {
		v[i].first=b=getint();
		v[i].second=i;
		for(register int j=1;p[j]*p[j]<=b;j++) {
			if(b%p[j]!=0) continue;
			while(b%p[j]==0) b/=p[j];
			t[j].modify(t[j].root,1,n,i,v[i].first);
		}
		if(b!=1) {
			const int j=std::lower_bound(&p[1],&p[p[0]]+1,b)-p;
			t[j].modify(t[j].root,1,n,i,v[i].first);
		}
	}
	std::sort(&v[1],&v[n]+1);
	for(register int i=0;i<m;i++) {
		int g=getint();
		const int l=getint(),r=getint();
		int max=-1;
		for(register int i=1;p[i]*p[i]<=g;i++) {
			if(g%p[i]!=0) continue;
			while(g%p[i]==0) g/=p[i];
			max=std::max(max,t[i].query(t[i].root,1,n,l,r));
		}
		if(g!=1) {
			const int &i=std::lower_bound(&p[1],&p[p[0]]+1,g)-p;
			max=std::max(max,t[i].query(t[i].root,1,n,l,r));
		}
		const int cnt=std::upper_bound(&v[1],&v[n]+1,std::make_pair(max,r))-std::lower_bound(&v[1],&v[n]+1,std::make_pair(max,l));
		printf("%d %d\n",max,cnt?:-1);
	}
	return 0;
}