题目：BZOJ3813、UOJ#38、洛谷P4140。
题目大意：某国家有100000个银行，初始每个银行只有3块钱。有两个操作：①修改某个银行存的钱；②求某段区间内所有钱的“总和”的欧拉函数值。
这些钱数的质因子只会包含前60个质数，这个国家的加法相当于我们的乘法（什么奇怪的国家(￣_,￣ )）。
解题思路：原题对第②种操作的描述是，总和为product，让你求编号1~product内，多少数number满足$number·x+product·y=1$。
我们将其转化为$number·xequiv 1(mod product)$，要使这个同余方程有解，当且仅当number和product互质。
那么这个问题就转化为[1,product]里与product互质的数有多少个，即$varphi(product)$。
然后利用公式$varphi(n)=n*prod_{p_i|n}(1-frac{1}{p_i})=n*prod_{p_i|n}frac{p_i-1}{p_i}$求解即可。
区间维护，显然用线段树，那么还要保存每段出现的质因子怎么办？由于最多出现60个质数，我们用一个long long的每一位表示一个质数，然后用或运算(or,|)即可实现相乘操作。
然后它又要模一个19961993，所以求欧拉函数时要用到这60个质数的逆元，而19961993是个质数，我们线性推一波逆元即可。
概括算法为：线段树+位运算+数论。
所以总时间复杂度为$O(tlog (100000+60))$。
这道那么玄妙的题目居然被我1A了，有点不可思议(⊙o⊙)
C++ Code：
#include<cstdio>
using namespace std;
#define mo 19961993
const int prime[61]={0,2,3,5,7,11,13,17,19,23,29,
31,37,41,43,47,53,59,61,67,71,73,79,83,89,97,101,103,107,109,113,
127,131,137,139,149,151,157,163,167,173,
179,181,191,193,197,199,211,223,227,229,
233,239,241,251,257,263,269,271,277,281};//60个素数
int inv[300];
struct node{
long long s,p;//s表示区间乘积，p储存该乘积包含哪些素数，第i个二进制位为1表示包含第i个素数（从1开始）
}d[400004],ans;
void bt(int l,int r,int o){//建树
if(l==r){
d[o].s=3;
d[o].p=2;
return;
}
int mid=l+r>>1;
bt(l,mid,o<<1);
bt(mid+1,r,o<<1|1);
d[o].s=d[o<<1].s*d[o<<1|1].s%mo;
d[o].p=2;
}
void change(int l,int r,int o,const int t,const int k){//单点修改
if(l==r){
d[o].s=k;
long long p=0;
for(int i=1;i<=60;++i){
if(!(k%prime[i]))p|=1LL<<(i-1);
}
d[o].p=p;
return;
}
int mid=l+r>>1;
if(t<=mid)change(l,mid,o<<1,t,k);else
change(mid+1,r,o<<1|1,t,k);
d[o].s=d[o<<1].s*d[o<<1|1].s%mo;
d[o].p=d[o<<1].p|d[o<<1|1].p;
}
void query(int l,int r,int o,int L,int R){//查询
if(L<=l&&r<=R){
ans.s=ans.s*d[o].s%mo;
ans.p|=d[o].p;
return;
}
int mid=l+r>>1;
if(L<=mid)query(l,mid,o<<1,L,R);
if(mid<R)query(mid+1,r,o<<1|1,L,R);
}
int main(){
bt(1,100000,1);
int t;
inv[1]=1;
for(int i=2;i<=281;++i)
inv[i]=(long long)(mo-mo/i)*inv[mo%i]%mo;//线性推逆元
scanf("%d",&t);
while(t--){
int x;
scanf("%d",&x);
if(x){
int t,k;
scanf("%d%d",&t,&k);
change(1,100000,1,t,k);
}else{
int L,R;
ans.s=1;
ans.p=0;
scanf("%d%d",&L,&R);
query(1,100000,1,L,R);
long long f=ans.s;
for(int i=1;i<=60;++i)
if(ans.p&(1LL<<(i-1))){//计算欧拉函数值
f=f*(prime[i]-1)%mo;
f=f*inv[prime[i]]%mo;
}
printf("%d
",(int)f);
}
}
return 0;
}