1. 程式人生 > >#38. 【清華集訓2014】奇數國|線段樹|尤拉函式

#38. 【清華集訓2014】奇數國|線段樹|尤拉函式

最後一個點一直T。。。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define LL long long
const int MAXN=100010;
const int MO=19961993;
const int N=100000;
struct H{LL cnt;int x;}seg[MAXN*4],ans;
int n,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};
void Ins(int now,int l,int r,int loc,int x)
{
	if(l==r)
	{
		seg[now].x=x;
		seg[now].cnt=0;
		for(int i=1;i<=60;i++)
		{
			seg[now].cnt<<=1;
			if(x%prime[i]==0)	
			{
				seg[now].cnt|=1;
				while(x%prime[i]==0) x/=prime[i];
			}
		}	
		return ;
	}
	int mid=(l+r)/2;
	if(loc<=mid) Ins(now*2,l,mid,loc,x);
	else Ins(now*2+1,mid+1,r,loc,x);
	seg[now].cnt=seg[now*2].cnt|seg[now*2+1].cnt;
	seg[now].x=(LL)seg[now*2].x*seg[now*2+1].x%MO;
}
void Q(int now,int l,int r,int s,int t)
{
	if(s<=l && r<=t)
	{
		ans.x=(LL)ans.x*seg[now].x%MO;
		ans.cnt|=seg[now].cnt;
		return ;
	}
	int mid=(l+r)/2;
	if(s<=mid) Q(now*2,l,mid,s,t);
	if(mid+1<=t) Q(now*2+1,mid+1,r,s,t);
}
int F(int x,int p)
{
	if(p==0) return 1;
	int t=F(x,p/2);
	t=(LL)t*t%MO;
	if(p&1) t=(LL)t*x%MO;
	return t;
}
void Build(int now,int l,int r)
{
	if(l==r)
	{
		seg[now].cnt=288230376151711744ll;
		seg[now].x=3;
		return ;
	}
	int mid=(l+r)/2;
	Build(now*2,l,mid);
	Build(now*2+1,mid+1,r);
	seg[now].cnt=seg[now*2].cnt|seg[now*2+1].cnt;
	seg[now].x=(LL)seg[now*2].x*seg[now*2+1].x%MO;
}
int main()
{	
	Build(1,1,N);
	cin >>n;int x,y,opt;
	for(int i=1;i<=n;i++)
	{
		scanf("%d %d %d",&opt,&x,&y);
		if(opt==0)
		{		
			ans.x=1;
			Q(1,1,N,x,y);
			for(int i=60;i>=1;i--)
			{
				if(ans.cnt&1)
					ans.x=(LL)ans.x*F(prime[i],MO-2)%MO*(prime[i]-1)%MO;
				ans.cnt>>=1;
			}
			printf("%d\n",ans.x);
		}
		else
			Ins(1,1,N,x,y);
	}
	return 0;
}