【UOJ#38】【清华集训2014】奇数国

考虑欧拉函数的性质，60很小，压位存下线段树每个节点出现质数。

#include<bits/stdc++.h>
const int N=;
const int yql=;
typedef long long ll;
using namespace std;
int prime[],n,m,vis[];
ll val,p,rev[];
struct Segment_Tree{
#define lson (o<<1)
#define rson (o<<1|1)
    ll sumv[N<<],pri[N<<];
    inline void pushup(int o){
        sumv[o]=(sumv[lson]*sumv[rson])%yql;
        pri[o]=pri[lson]|pri[rson];
    }
    inline void build(int o,int l,int r){
        if(l==r){sumv[o]=;pri[o]=;return;}
        int mid=(l+r)>>;
        build(lson,l,mid);build(rson,mid+,r);
        pushup(o);
    }
    inline void change(int o,int l,int r,int q,int v){
        if(l==r){
            sumv[o]=v;pri[o]=;
            for(int i=;i<;i++)if(v%prime[i]==)pri[o]|=1LL<<i;
            return;
        }
        int mid=(l+r)>>;
        if(q<=mid)change(lson,l,mid,q,v);
        else change(rson,mid+,r,q,v);
        pushup(o);
    }
    inline void query(int o,int l,int r,int ql,int qr){
        if(ql<=l&&r<=qr){val=val*sumv[o]%yql;p|=pri[o];return;}
        int mid=(l+r)>>;
        if(ql<=mid)query(lson,l,mid,ql,qr);
        if(qr>mid)query(rson,mid+,r,ql,qr);
    }
}T;
inline ll fpow(ll x,ll p){
    ll ans=;
    for(;p;p>>=,x=x*x%yql)if(p&)ans=ans*x%yql;
    return ans;
}
inline ll read(){
    ll f=,x=;char ch;
    do{ch=getchar();if(ch=='-')f=-;}while(ch<''||ch>'');
    do{x=x*+ch-'';ch=getchar();}while(ch>=''&&ch<='');
    return f*x;
}
int main(){
    int cnt=;
    for(int i=;cnt<;i++)if(!vis[i]){prime[cnt++]=i;for(int j=i;j<=;j+=i)vis[j]=;}
    for(int i=;i<;i++)rev[i]=fpow(prime[i],yql-)*1LL*(prime[i]-)%yql;
    T.build(,,N);int m=read();n=N;
    while(m--){
        int opt=read(),l=read(),r=read();
        if(opt==)T.change(,,n,l,r);
        else{
            val=;p=;T.query(,,n,l,r);
            for(int j=;j<;j++,p>>=)if(p&)val=1LL*val*rev[j]%yql;
            printf("%lld\n",val);
        }
    }
}