1. 程式人生 > >【UOJ#38】【清華集訓2014】奇數國

【UOJ#38】【清華集訓2014】奇數國

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); } 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); } } }

【UOJ#38】【清華集訓2014】奇數國