bzoj3813 奇数国（欧拉函数+树状数组）

其实就是求区间积的欧拉函数。每个数最多分解为前60个质数的乘积，单点修改。开60个树状数组来搞就好了。或者线段树+bitset什么的都好吧qaq

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100010
#define mod 19961993
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n=100000,m,prime[100],tot=0,c[61][N],cnt[61][N];
bool notprime[300];
inline void init(){
    notprime[1]=1;
    for(int i=2;i<=281;++i){
        if(!notprime[i]) prime[++tot]=i;
        for(int j=1;prime[j]*i<=281;++j){
            notprime[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
}
inline void add(int id,int x,int val){
    for(;x<=n;x+=x&-x) c[id][x]+=val;
}
inline int ask(int id,int x){
    int res=0;for(;x;x-=x&-x) res+=c[id][x];return res;
}
inline int ksm(ll x,int k){
    ll res=1;for(;k;k>>=1,x=x*x%mod) if(k&1) res=res*x%mod;return res;
}
int main(){
//  freopen("a.in","r",stdin);
    m=read();init();
    for(int i=1;i<=n;++i) cnt[2][i]=1,add(2,i,1);
    while(m--){
        int op=read(),x=read(),y=read();
        if(op==0){
            ll res=1;
            for(int i=1;i<=60;++i){
                int d=ask(i,y)-ask(i,x-1);
                if(d) res=res*ksm(prime[i],d-1)%mod*(prime[i]-1)%mod;
            }printf("%lld\n",res);
        }else{
            for(int i=1;i<=60;++i){
                int d=0;while(y%prime[i]==0) y/=prime[i],++d;
                int dd=d-cnt[i][x];if(dd) add(i,x,dd);cnt[i][x]=d;
            }
        }
    }return 0;
}