BZOJ4026: dC Loves Number Theory

题解: 由题意我们可以得到这个式子  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\n",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 MB
Submit: 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

Sample Output

40
240
12
1200
2
240
4
4
1200
4

HINT

 

 1 <= N <= 50000 

1 <= Q <= 100000 

1 <= Ai <= 10^6