bzoj3813: 奇数国

9 篇文章 0 订阅
9 篇文章 0 订阅

该村子在财产上的求和运算等同于我们的乘法运算

该村子在财产上的求和运算等同于我们的乘法运算

该村子在财产上的求和运算等同于我们的乘法运算

一开始直接把它忽略了。。。导致题意都没看懂,直接看题解一句话题意,听说是最简单的一道。。。

考虑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;
}


  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值