bzoj3813: 奇數國

該村子在財產上的求和運算等同於我們的乘法運算

該村子在財產上的求和運算等同於我們的乘法運算

該村子在財產上的求和運算等同於我們的乘法運算

一開始直接把它忽略了。。。導致題意都沒看懂，直接看題解一句話題意，聽說是最簡單的一道。。。

考慮phi(x)=x*∏(pi-1)/pi，維護區間和，區間每個質數是否存在(long long)，即可得到答案。

#include<iostream>
#include<cstdio>
#define Hz 19961993
#define N 100005
#define n 100000
#define ll long long
using namespace std;
int pr[60]={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};
ll inv[N],sum[N<<2],bit[N<<2],Sum,Bit,T,p,x,y;
void add(int k,int l,int r,int x,int y,ll p)
{
	if (l==r){sum[k]=y;bit[k]=p;return;}
	int mid=(l+r)>>1,L=k<<1,R=L|1;
	if (x<=mid) add(L,l,mid,x,y,p);
	else add(R,mid+1,r,x,y,p);
	sum[k]=sum[L]*sum[R]%Hz;bit[k]=bit[L]|bit[R];
}
void qry(int k,int l,int r,int x,int y)
{
	if (x<=l&&r<=y)
		{Sum=Sum*sum[k]%Hz;Bit|=bit[k];return;}
	int mid=(l+r)>>1;
	if (x<=mid) qry(k<<1,l,mid,x,y);
	if (y>mid) qry(k<<1|1,mid+1,r,x,y);
}
ll part(ll x)
{
	ll ans=0;
	for (ll i=0,j=1;i<60;i++,j<<=1)
		if (x%pr[i]==0) ans|=j;
	return ans;
}
ll Get(ll x,ll y)
{
	for (ll i=0,j=1;i<60;i++,j<<=1)
		if (y&j) x=x*inv[pr[i]]%Hz*(pr[i]-1)%Hz;
	return x;
}
int main()
{
	inv[1]=1;
	for (int i=2;i<=n;i++)
		inv[i]=(Hz-Hz/i)*inv[Hz%i]%Hz;
	p=part(3);
	for (int i=1;i<=n;i++)
		add(1,1,n,i,3,p);
	
	scanf("%lld",&T);
	while(T--)
	{
		scanf("%lld%lld%lld",&p,&x,&y);
		if (p==1)
		{
			p=part(y);
			add(1,1,n,x,y,p);
		}
		else
		{
			Sum=1;Bit=0;
			qry(1,1,n,x,y);
			printf("%lld\n",Get(Sum,Bit));
		}
	}
	return 0;
}