题解: 由题意我们可以得到这个式子 ans=f(ai*...*aj)=f(pi^ki.....pj*kj)=(pi-1)*pi^(ki-1).....(pj-1)*pj^(kj-1) 所以ans=pi^ki....*pj*kj*((pi-1)/pi......*(pj-1)/pj) 然后前面部分我们可以通过前缀积得到 后面的可以通过主席树+永久化标记维护
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <stack>
#include <queue>
#include <cmath>
#include <set>
#include <map>
#define mp make_pair
#define pb push_back
#define pii pair<int,int>
#define link(x) for(edge *j=h[x];j;j=j->next)
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,r,l) for(int i=r;i>=l;i--)
const int MAXN=5e4+10;
const int NM=1e5+10;
const double eps=1e-8;
const int mod=1e6+777;
#define ll long long
using namespace std;
struct edge{int t,v;edge*next;}e[MAXN<<1],*h[MAXN],*o=e;
void add(int x,int y,int vul){o->t=y;o->v=vul;o->next=h[x];h[x]=o++;}
ll read(){
ll x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return x*f;
}
int n,q,a[MAXN],rt[MAXN],cnt,last[1000005];
ll sum[MAXN];
bool flag;
ll ksm(ll x,ll y){
ll res=1;
while(y){
if(y&1)res=res*x%mod;
x=x*x%mod;y=y>>1;
}
return res;
}
typedef struct node{
int l,r;ll tag;
}node;
node d[NM*101];
void update1(int &x,int y,int l,int r,int ql,int qr,ll vul){
//cout<<ql<<" "<<qr<<" "<<vul<<endl;
x=++cnt;d[x]=d[y];
if(!d[x].tag)d[x].tag=1;
if(ql<=l&&r<=qr){d[x].tag=(1LL*d[x].tag*vul)%mod;return ;}
int mid=(l+r)>>1;
if(ql<=mid)update1(d[x].l,d[y].l,l,mid,ql,qr,vul);
if(qr>mid)update1(d[x].r,d[y].r,mid+1,r,ql,qr,vul);
}
ll ans;
void querty(int x,int l,int r,int t){
//cout<<l<<" "<<r<<" "<<d[x].tag<<endl;
if(d[x].tag)ans=(1ll*ans*d[x].tag)%mod;
if(l==r)return ;
int mid=(l+r)>>1;
if(t<=mid)querty(d[x].l,l,mid,t);
else querty(d[x].r,mid+1,r,t);
}
int main(){
//cout<<4*ksm(5,mod-2)%mod<<endl;
n=read();q=read();int t;
sum[0]=1;
inc(i,1,n){
rt[i]=++cnt;d[rt[i]]=d[rt[i-1]];
t=read();sum[i]=sum[i-1];
for(int j=2;j<=sqrt(t);j++){
if(t%j==0){
int num=0;
while(t%j==0)num++,t/=j;
sum[i]=(sum[i]*ksm(j,num))%mod;
update1(rt[i],rt[i],1,n,last[j]+1,i,1LL*(j-1)*ksm(j,mod-2)%mod);
last[j]=i;
}
}
if(t!=1){
sum[i]=(sum[i]*t)%mod;
update1(rt[i],rt[i],1,n,last[t]+1,i,1LL*(t-1)*ksm(t,mod-2)%mod);
last[t]=i;
}
}
ll res=0;int l,r;
while(q--){
l=read()^res;r=read()^res;
ans=1;querty(rt[r],1,n,l);
printf("%lld
",res=1LL*sum[r]*ksm(sum[l-1],mod-2)%mod*ans%mod);
}
return 0;
}
4026: dC Loves Number Theory
Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 803 Solved: 240
[Submit][Status][Discuss]
Description
dC 在秒了BZOJ 上所有的数论题后，感觉萌萌哒，想出了这么一道水题，来拯救日益枯
竭的水题资源。
给定一个长度为 n的正整数序列A，有q次询问，每次询问一段区间内所有元素乘积的
φ（φ(n)代表1~n 中与n互质的数的个数） 。由于答案可能很大，所以请对答案 mod 10^6 +
777。 （本题强制在线，所有询问操作的l,r都需要 xor上一次询问的答案 lastans，初始时，
lastans = 0）
Input
第一行，两个正整数，N，Q，表示序列的长度和询问的个数。
第二行有N 个正整数，第i个表示Ai.
下面Q行，每行两个正整数，l r，表示询问[l ^ lastans,r ^ lastans]内所有元素乘积的φ
Output
Q行，对于每个询问输出一个整数。
Sample Input
5 10
3 7 10 10 5
3 4
42 44
241 242
14 9
1201 1201
0 6
245 245
7 7
6 1
1203 1203
3 7 10 10 5
3 4
42 44
241 242
14 9
1201 1201
0 6
245 245
7 7
6 1
1203 1203
Sample Output
40
240
12
1200
2
240
4
4
1200
4
240
12
1200
2
240
4
4
1200
4
HINT
1 <= N <= 50000
1 <= Q <= 100000
1 <= Ai <= 10^6