BZOJ 4026: dC Loves Number Theory（套路主席树）

题目
发现就是区间乘积/${\sum }_{pr出现在[l,r]的质因数中}(pr-1)/pr$
关于出现过就加入这一条件。
有一个套路：
把每一个数$y$与他前面第一个和他相等的数$x$形成一个二元组$(loc[x],loc[y])$，第1个往0连。那么如果$loc[x]<landl<=loc[y]<=r$那么$y$所代表的数在$[l,r]$出现过。且只有一组$loc[x]<landl<=loc[y]<=r$，其他的$y$所代表的数都不会造成贡献。
那么就是一个二维数点问题，直接上主席树。
AC Code:

#include<bits/stdc++.h>
#define maxn 500005
#define maxpt maxn * 30
#define mod 1000777
using namespace std;

int n,q,a[maxn],sum[maxn],pre[mod],inv[mod]={1,1};
int rt[maxn],lc[maxpt],rc[maxpt],val[maxpt],tot;
int cnt = 0;
vector<pair<int,int> >G[maxn];

void Insert(int &now,int l,int r,int pos,int v){
	val[++tot]=max(1,val[now])*1ll*v%mod,lc[tot]=lc[now],rc[tot]=rc[now];
	now = tot;
	if(l == r) return;
	int mid = (l+r) >> 1;
	if(pos <= mid) Insert(lc[now],l,mid,pos,v);
	else Insert(rc[now],mid+1,r,pos,v);
}

int query(int now,int l,int r,int ql,int qr){
	if(l>qr || r<ql || !now) return 1;
	if(ql<=l && r<=qr) return val[now];
	int mid = (l+r) >> 1;
	return 1ll * query(lc[now],l,mid,ql,qr) * query(rc[now],mid+1,r,ql,qr) % mod;
}

int main(){
	scanf("%d%d",&n,&q);sum[0]=1;
	for(int i=1;i<=n;i++) scanf("%d",&a[i]),sum[i]=1ll*sum[i-1]*a[i]%mod;
	for(int i=2;i<mod;i++) inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;
	for(int i=1;i<=n;i++){
		for(int j=2;j*j<=a[i];j++)
			if(a[i] % j == 0){
				G[pre[j]].push_back(make_pair(i,1ll*(j-1)*inv[j]%mod));
				pre[j] = i;
				while(a[i] % j == 0) a[i] /= j;
			}
		if(a[i] > 1) 
			G[pre[a[i]]].push_back(make_pair(i,1ll*(a[i]-1)*inv[a[i]]%mod)) , pre[a[i]]=i;
	}
	for(int i=0;i<n;i++){
		if(i) rt[i] = rt[i-1];
		for(int j=0,siz=G[i].size();j<siz;j++)
			Insert(rt[i],1,n,G[i][j].first,G[i][j].second);
	}
	for(int l,r,la=0;q--;){
		scanf("%d%d",&l,&r);
		l^=la,r^=la;
		la=sum[r] * 1ll * inv[sum[l-1]] % mod * query(rt[l-1],1,n,l,r) % mod;
		printf("%d\n",la=(la+mod)%mod);
	}
}