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

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