[清华集训2014]奇数国

Posted Mrsrz

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[清华集训2014]奇数国相关的知识,希望对你有一定的参考价值。

题目:BZOJ3813、UOJ#38。

题目大意:某国家有100000个银行,初始每个银行只有3块钱。有两个操作:①修改某个银行存的钱;②求某段区间内所有钱的“总和”的欧拉函数值。

这些钱数的质因子只会包含前60个质数,这个国家的加法相当于我们的乘法(什么奇怪的国家( ̄_, ̄ ))。

解题思路:原题对第②种操作的描述是,总和为product,让你求编号1~product内,多少数number满足$number·x+product·y=1$。

我们将其转化为$number·x\equiv 1(mod\ product)$,要使这个同余方程有解,当且仅当number和product互质。

那么这个问题就转化为[1,product]里与product互质的数有多少个,即$\varphi(product)$。

然后利用公式$\varphi(n)=n*\prod_{i}(1-\frac{1}{p_i})=n*\prod_{i}\frac{p_i-1}{p_i}$求解即可。

区间维护,显然用线段树,那么还要保存每段出现的质因子怎么办?由于最多出现60个质数,我们用一个long long的每一位表示一个质数,然后用或运算(or,|)即可实现相乘操作。

然后它又要模一个19961993,所以求欧拉函数时要用到这60个质数的逆元,而19961993是个质数,我们线性推一波逆元即可。

概括算法为:线段树+位运算+数论。

所以总时间复杂度为$O(t\log (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\n",(int)f);
		}
	}
	return 0;
}

 

以上是关于[清华集训2014]奇数国的主要内容,如果未能解决你的问题,请参考以下文章

[清华集训2014]奇数国

UOJ#38清华集训2014奇数国

[清华集训2014]奇数国

清华集训2014 做题记录

奇数国[清华集训2014 Day1]

清华集训 2014--奇数国(线段树&欧拉函数&乘法逆元&状态压缩)