洛谷P4140 奇数国

Description

在一片美丽的大陆上有100000个国家,记为1到100000。这里经济发达,有数不尽的账房,并且每个国家有一个银行。某大公司的领袖在这100000个银行开户时都存了3大洋,他惜财如命,因此会不时地派小弟GFS清点一些银行的存款或者让GFS改变某个银行的存款。该村子在财产上的求和运算等同于我们的乘法运算,也就是说领袖开户时的存款总和为3100000。这里发行的软妹面额是最小的60个素数(p1=2,p2=3,…,p60=281),任何人的财产都只能由这60个基本面额表示,即设某个人的财产为fortune(正整数),则fortune=p1k1*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取模后的答案。


Input

第一行一个整数x表示领袖清点和变动存款的总次数。

接下来x行,每行3个整数ai,bi,ci。ai为0时表示该条记录是清点计划,领袖会清点bi到ci的银行存款,你需要对该条记录计算出GFS想要的答案。ai为1时表示该条记录是存款变动,你要把银行bi的存款改为ci,不需要对该记录进行计算。


Output

输出若干行,每行一个数,表示那些年的答案。


Solution

ax+by =1 那么gcd(a,b)=1
再想想小于等于x且与x互质的数个数 那不就是欧拉函数嘛
题目六十个质数已经明示了

欧拉函数是不完全积性函数

ϕ ( p ) = p − 1 \phi (p)=p-1 ϕ(p)=p1 其中p是质数
其中对于 ϕ ( x ) \phi (x) ϕ(x) ϕ ( p ) \phi (p) ϕ(p)
若x%p==0 则 ϕ ( x ∗ p ) = ϕ ( x ) ∗ p \phi (x*p)=\phi (x)*p ϕ(xp)=ϕ(x)p
否则 ϕ ( x ∗ p ) = ϕ ( x ) ∗ ϕ ( p ) \phi (x*p)=\phi (x)*\phi(p) ϕ(xp)=ϕ(x)ϕ(p)

那么直接60棵线段树维护每个银行存款质因数分解情况 查询直接快速幂乘就行了。


Code

#include<bits/stdc++.h>
using namespace std;
const int mod=19961993;
struct data{
	int a[61];
}tree[400010];
int p[101],cnt,n;
bool vis[1010];
void build(int id,int l,int r){
	if(l==r){
		tree[id].a[2]=1;
		return;
	}
	int mid=(l+r)>>1;
	build(id<<1,l,mid);
	build(id<<1|1,mid+1,r);
	for(int i=1;i<=60;i++)
	tree[id].a[i]=tree[id<<1].a[i]+tree[id<<1|1].a[i];
}
void update(int id,int l,int r,int x,int v){
	//cout<<id<<" "<<l<<" "<<r<<" "<<x<<" "<<v<<endl;
	if(l==r){
		int s=v,tot=0;
		for(int i=1;i<=60;i++){
			tot=0;
			while(s%p[i]==0)
			s/=p[i],tot++;
			tree[id].a[i]=tot;
		}
		return;
	}
	int mid=(l+r)>>1;
	if(x<=mid) update(id<<1,l,mid,x,v);
	else update(id<<1|1,mid+1,r,x,v);
	for(int i=1;i<=60;i++)
	tree[id].a[i]=tree[id<<1].a[i]+tree[id<<1|1].a[i];
}
data query(int id,int l,int r,int ql,int qr){
	if(ql<=l&&r<=qr) return tree[id];
	data ans,u;
	ans.a[0]=1;
	for(int i=1;i<=60;i++)
	ans.a[i]=0;
	int mid=(l+r)>>1;
	if(ql<=mid) ans=query(id<<1,l,mid,ql,qr);
	if(qr>mid){
		u=query(id<<1|1,mid+1,r,ql,qr);
		if(ans.a[0]==1) ans=u;
		else for(int i=1;i<=60;i++)
		ans.a[i]+=u.a[i];
	}
	return ans;
}
long long fpow(int x,int y){
	long long s=1,base=x;
	while(y){
		if(y&1) s=(s*base)%mod;
		base=(base*base)%mod;
		y>>=1;
	}
	return s;
}
void work(int n){
	for(int i=2;i<n;i++){
		if(!vis[i]) p[++cnt]=i;
		if(cnt==60) break;
		for(int j=1;j<=cnt&&p[j]*i<n;j++){
			vis[p[j]*i]=1;
			if(i%p[j]==0) break;
		}
	}
}
void qry(int x,int y){
	data u=query(1,1,1e5,x,y);
	long long sum=1;
	for(int i=1;i<=60;i++)
	if(u.a[i]) sum=((((p[i]-1)%mod*fpow(p[i],u.a[i]-1))%mod)*sum)%mod;//,cout<<i<<" "<<p[i]<<" "<<u.a[i]<<endl;
	printf("%lld\n",sum%mod);
}
int main(){
	int x,y,z;
	work(1000);
	scanf("%d",&n);
	build(1,1,1e5);
	for(int i=1;i<=n;i++){
		scanf("%d%d%d",&x,&y,&z);
		if(x==0) qry(y,z);
		else update(1,1,1e5,y,z);
	}
}
/*
6
0 1 3
1 1 5
0 1 3
1 1 7
0 1 3
0 2 3
*/
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值