bzoj3813 奇数国 线段树+欧拉函数

Description

---

在一片美丽的大陆上有100000个国家，记为1到100000。这里经济发达，有数不尽的账房，并且每个国家有一个银行。某大公司的领袖在这100000个银行开户时都存了3大洋，他惜财如命，因此会不时地派小弟GFS清点一些银行的存款或者让GFS改变某个银行的存款。该村子在财产上的求和运算等同于我们的乘法运算，也就是说领袖开户时的存款总和为3100000。这里发行的软妹面额是最小的60个素数（p1=2,p2=3,…,p60=281），任何人的财产都只能由这60个基本面额表示，即设某个人的财产为fortune（正整数），则fortune=p1^k1*p2k2*…p60^K60。
领袖习惯将一段编号连续的银行里的存款拿到一个账房去清点，为了避免GFS串通账房叛变，所以他不会每次都选择同一个账房。GFS跟随领袖多年已经摸清了门路,知道领袖选择账房的方式。如果领袖选择清点编号在[a,b]内的银行财产，他会先对[a,b]的财产求和（计为product），然后在编号属于[1,product]的账房中选择一个去清点存款，检验自己计算是否正确同时也检验账房与GFS是否有勾结。GFS发现如果某个账房的编号number与product相冲，领袖绝对不会选择这个账房。怎样才算与product不相冲呢？若存在整数x,y使得numberx+producty=1，那么我们称number与product不相冲，即该账房有可能被领袖相中。当领袖又赚大钱了的时候，他会在某个银行改变存款，这样一来相同区间的银行在不同的时候算出来的product可能是不一样的，而且领袖不会在某个银行的存款总数超过1000000。
现在GFS预先知道了领袖的清点存款与变动存款的计划，想请你告诉他，每次清点存款时领袖有多少个账房可以供他选择，当然这个值可能非常大，GFS只想知道对19961993取模后的答案。

x≤100000，当ai=0时0≤ci−bi≤100000

Solution

---

考虑那个柿子的含义，有解当且仅当gcd(prod,num)=1
因此我们实际上要求区间积的φ，质数只有60个，然后就十分显然了

跑得好慢啊，预处理质数的逆元和次幂才卡过去QUQ

Code

---

#include <stdio.h>
#include <string.h>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)

typedef long long LL;
const int MOD=19961993;
const int N=100005;

struct rec {
	int r[61];

	rec operator +(rec b) const {
		rec a=*this,c;
		rep(i,1,60) c.r[i]=a.r[i]+b.r[i];
		return c;
	}
} sum[N<<2];

bool not_prime[N];

int prime[N];

LL ny[N],pw[61][N];

int read() {
	int x=0,v=1; char ch=getchar();
	for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
	for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
	return x*v;
}

LL ksm(LL x,LL dep) {
	LL ret=1;
	for (;dep;dep>>=1) {
		(dep&1)?(ret=ret*x%MOD):0;
		x=x*x%MOD;
	}
	return ret;
}

void pre_work() {
	for (int i=2;prime[0]<60;++i) {
		if (!not_prime[i]) prime[++prime[0]]=i;
		for (int j=1;i*prime[j]<N&&j<=prime[0];++j) {
			not_prime[i*prime[j]]=1;
			if (i%prime[j]==0) break;
		}
	}
}

void modify(int now,int tl,int tr,int x,rec v) {
	if (tl==tr) return (void) (sum[now]=v);
	int mid=(tl+tr)>>1;
	if (x<=mid) modify(now<<1,tl,mid,x,v);
	else modify(now<<1|1,mid+1,tr,x,v);
	sum[now]=sum[now<<1]+sum[now<<1|1];
}

rec query(int now,int tl,int tr,int l,int r) {
	if (tl==l&&tr==r) return sum[now];
	int mid=(tl+tr)>>1;
	if (r<=mid) return query(now<<1,tl,mid,l,r);
	if (l>mid) return query(now<<1|1,mid+1,tr,l,r);
	rec qx=query(now<<1,tl,mid,l,mid);
	rec qy=query(now<<1|1,mid+1,tr,mid+1,r);
	return qx+qy;
}

void build(int now,int tl,int tr) {
	if (tl==tr) {
		sum[now].r[2]=1;
		return ;
	}
	int mid=(tl+tr)>>1;
	build(now<<1,tl,mid); build(now<<1|1,mid+1,tr);
	sum[now]=sum[now<<1]+sum[now<<1|1];
}

int main(void) {
	build(1,1,100000);
	pre_work();
	rep(i,1,prime[0]) {
		ny[i]=ksm(prime[i],MOD-2);
		pw[i][0]=1;
		rep(j,1,100000) pw[i][j]=pw[i][j-1]*prime[i]%MOD;
	}
	for (int T=read();T--;) {
		int opt=read(),l=read(),r=read();
		if (opt==0) {
			LL ans=1;
			rec ret=query(1,1,100000,l,r);
			rep(i,1,60) if (ret.r[i]) {
				ans=ans*(prime[i]-1)%MOD*ny[i]%MOD*pw[i][ret.r[i]]%MOD;
			}
			printf("%lld\n", ans);
		} else { rec v;
			rep(i,1,60) {
				v.r[i]=0;
				if (r%prime[i]==0) {
					int cnt=0;
					for (;r%prime[i]==0;) r/=prime[i],cnt++;
					v.r[i]=cnt;
				}
			}
			modify(1,1,100000,l,v);
		}
	}
	return 0;
}