#216. REQ

讲道理刚看完题目就听到题解了233

直接离线。从左往右枚举r，用a[i]记录l=i时的答案，考虑每次r++后的影响：

对于r的每个质因数p，找到左边第一个被p整除的数i，a[i+1]到a[r]乘以φ(p)，a[1]到a[i]乘以p。

然而值得了六十分。哪里WA了呢？我也不知道。可能是因为我叫wamach。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdlib>
using namespace std;
#define rep(i,j,k) for(i=j;i<=k;++i)
#define per(i,j,k) for(i=j;i>=k;--i)
#define ll long long
#define db double
#define mkp(x,y) make_pair(x,y)
#define pii pair<int,int>
#define X first
#define Y second
const int N=200005,inf=1000001,mod=1000000007;
pii qj[N];bool vis[N];
int n,a[N],Q,sa[N],pri[inf],phi[inf],cnt,ans[N],b[inf];ll lz[N<<2];
bool cmp(int x,int y){
	return qj[x]<qj[y];
}
void getpri(){
	int i,j,y;
	for(i=2;i<inf;++i){
		if(!vis[i])phi[pri[++cnt]=i]=i-1;
		for(j=1;(y=pri[j]*i)<inf;++j){
			vis[y]=1;
			if(i%pri[j]==0){
				phi[y]=phi[i]*pri[j];break;
			}
			phi[y]=(pri[j]-1)*phi[i];
		}
	}
}
void build(int l,int r,int num){
	lz[num]=1;
	if(l==r)return;
	int mid=l+r>>1;
	build(l,mid,num<<1);
	build(mid+1,r,num<<1|1);
}
void pd(int num){
	lz[num<<1]=lz[num<<1]*lz[num]%mod;
	lz[num<<1|1]=lz[num<<1|1]*lz[num]%mod;
	lz[num]=1;
}
void add(int x,int L,int R,int l,int r,int num){
	if(L<=l&&r<=R){
		lz[num]=lz[num]*x%mod;return;
	}
	pd(num);
	int mid=l+r>>1;
	if(mid>=L)add(x,L,R,l,mid,num<<1);
	if(mid<R)add(x,L,R,mid+1,r,num<<1|1);
}
void modify(int i,int x,int y){
	if(x>2)add(phi[x],b[y]+1,i,1,n,1);
	if(b[y])add(x,1,b[y],1,n,1);
	b[y]=i;
}
ll query(int x,int l,int r,int num){
	if(l==r)return lz[num];
	pd(num);
	int mid=l+r>>1;
	return x>mid?query(x,mid+1,r,num<<1|1):query(x,l,mid,num<<1);
}
int main(){
	getpri();
	int i,j,x,y;
	scanf("%d",&n);
	rep(i,1,n)scanf("%d",&a[i]);
	scanf("%d",&Q);
	rep(i,1,Q){
		scanf("%d%d",&x,&y);
		qj[sa[i]=i]=mkp(y,x);
	}
	sort(sa+1,sa+1+Q,cmp);
	build(1,n,1);
	y=1;
	rep(i,1,n){
		for(j=1;pri[j]*pri[j]<=a[i];++j){
			if(a[i]%pri[j]==0){
				for(x=1;a[i]%pri[j]==0;x*=pri[j])a[i]/=pri[j];
				modify(i,x,pri[j]);
			}
		}
		if(a[i]>1)modify(i,a[i],a[i]);
		for(;y<=Q&&qj[sa[y]].X==i;++y)
			ans[sa[y]]=query(qj[sa[y]].Y,1,n,1);
	}
	rep(i,1,Q)printf("%d\n",ans[i]);
	return 0;
}