BZOJ 4026: dC Loves Number Theory 可持久化线段树 + 欧拉函数 + 数学

Code: 

#include <bits/stdc++.h>
#define ll long long  
#define maxn 50207 
#define setIO(s) freopen(s".in","r",stdin) 
#define mod 1000777 
using namespace std;
struct Tree
{
    int tot; 
    int lson[maxn*270],rson[maxn*270]; 
    ll mul[maxn*270];  
    int ins(int x,int l,int r,int p,ll v) 
    {
        int o=++tot; 
        lson[o]=lson[x],rson[o]=rson[x]; 
        if(x) mul[o]=mul[x]*v%mod; 
        else mul[o]=v; 
        if(l==r) return o; 
        int mid=(l+r)>>1;  
        if(p<=mid) lson[o]=ins(lson[x],l,mid,p,v); 
        else rson[o]=ins(rson[x],mid+1,r,p,v); 
        return o;
    } 
    ll query(int x,int l,int r,int L,int R) 
    {
        if(l>=L&&r<=R) return mul[x]; 
        ll re=1; 
        int mid=(l+r)>>1; 
        if(L<=mid) re=re*query(lson[x],l,mid,L,R)%mod; 
        if(R>mid) re=re*query(rson[x],mid+1,r,L,R)%mod; 
        return re;  
    }
}tr; 
int cnt,n,m; 
int prime[1000004],vis[1000004],pre[1000006],rt[maxn],cur[maxn];  
ll arr[maxn],mul[maxn]; 
ll qpow(ll base,ll k) 
{
    ll re=1; 
    while(k) 
    {
        if(k&1) re=re*base%mod; 
        base=base*base%mod; 
        k>>=1; 
    }
    return re; 
}
void Initialize(int N) 
{
    for(int i=2;i<=N;++i) 
    {
        if(!vis[i]) prime[++cnt]=i; 
        for(int j=1;j<=cnt&&1ll*prime[j]*i<=N;++j) 
        {
            vis[i*prime[j]]=1; 
            if(i%prime[j]==0) break; 
        }
    } 
    for(int i=1;i<=n;++i) 
    {
        int x=arr[i],cc=0;   
        rt[i]=rt[i-1];  
        if(x==1) 
        {
            rt[i]=tr.ins(rt[i-1],1,n,i,1); 
            continue; 
        }
        cur[++cc]=rt[i];   
        for(int j=1;j<=cnt&&1ll*prime[j]*prime[j]<=x;++j) 
        {
            if(x%prime[j]==0) 
            { 
                if(pre[prime[j]])
                {
                    cur[cc+1]=tr.ins(cur[cc],1,n,pre[prime[j]],1ll*prime[j]*qpow((prime[j]-1), mod-2)%mod); 
                    ++cc; 
                }
                cur[cc+1]=tr.ins(cur[cc],1,n,pre[prime[j]]=i,1ll*(prime[j]-1)*qpow(prime[j],mod-2)%mod); 
                ++cc;     
                while(x%prime[j]==0) x/=prime[j];       
            }
        }
        if(x>1) 
        {
            if(pre[x]) 
            { 
                cur[cc+1]=tr.ins(cur[cc],1,n,pre[x],1ll*x*qpow((x-1), mod-2)%mod); 
                ++cc; 
            }
            cur[cc+1]=tr.ins(cur[cc],1,n,pre[x]=i,1ll*(x-1)*qpow(x,mod-2)%mod); 
            ++cc;     
        }
        rt[i]=cur[cc];  
        for(int i=1;i<=cc;++i) cur[i]=0; 
    }
} 
int main() 
{
    // setIO("input"); 
    scanf("%d%d",&n,&m); 
    int mx=0; 
    mul[0]=1;  
    for(int i=1;i<=n;++i) 
    {
        scanf("%lld",&arr[i]); 
        mul[i]=mul[i-1]*arr[i]%mod; 
        mx=max(mx, (int)arr[i]); 
    } 
    Initialize(mx);   
    ll lastans=0; 
    for(int i=1;i<=m;++i) 
    {
        // lastans=0;
        int l,r; 
        scanf("%d%d",&l,&r); 
        l^=lastans,r^=lastans;             
        ll rev=qpow(mul[l-1],mod-2); 
        lastans=mul[r]*rev%mod*tr.query(rt[r],1,n,l,r)%mod;    
        printf("%lld\n",lastans); 
    }
    return 0;
}