【UOJ#38】【清華集訓2014】奇數國
阿新 • • 發佈:2017-09-08
lin sum chang esp color 奇數 efi bsp urn
考慮歐拉函數的性質,60很小,壓位存下線段樹每個節點出現質數。
#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); } inlinevoid 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); } } }
【UOJ#38】【清華集訓2014】奇數國