Newcoder 132 E.简单数据结构1(BIT+数论)

186 篇文章 0 订阅
54 篇文章 0 订阅

Description

给一个长为 n n n的序列, m m m次操作,每次操作:

1.区间 [ l , r ] [l,r] [l,r] x x x

2.对于区间 [ l , r ] [l,r] [l,r],查询 a l a l + 1 a l + 2 . . .   m o d   p a_l^{a_{l+1}^{a_{l+2}...}}\ mod\ p alal+1al+2... mod p,一直到 a r a_r ar

请注意每次的模数不同。

Input

第一行两个整数 n , m n,m n,m表示序列长度和操作数

接下来一行, n n n个整数,表示这个序列

接下来 m m m行,可能是以下两种操作之一:

操作 1 1 1:区间 [ l , r ] [l,r] [l,r]加上 x x x

操作 2 2 2:查询区间 [ l , r ] [l,r] [l,r]的那个式子 m o d   p mod\ p mod p的值

( n , m ≤ 5 ⋅ 1 0 5 , x ≤ 2 ⋅ 1 0 9 , p ≤ 2 ⋅ 1 0 7 ) (n,m\le 5\cdot 10^5,x\le 2\cdot 10^9,p\le 2\cdot 10^7) (n,m5105,x2109,p2107)

Output

对于每个询问,输出一个数表示答案

Sample Input

6 4
1 2 3 4 5 6
2 1 2 10000007
2 2 3 5
1 1 4 1
2 2 4 10

Sample Output

1
3
1

Solution

树状数组维护区间修改,由指数循环定理 A B % C = A B % φ ( C ) + φ ( C ) % C , B ≥ φ ( C ) A^B\%C=A^{B\%\varphi(C)+\varphi(C)}\%C,B\ge \varphi(C) AB%C=AB%φ(C)+φ(C)%C,Bφ(C)和欧拉函数的快速收敛,直接递归求即可,注意有时候指数可能不用多加 φ ( C ) \varphi(C) φ(C),故需要判断其幂指数上的值和 φ ( C ) \varphi(C) φ(C)的大小关系

Code

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
#define maxn 20000007
int n,m,prime[maxn],res,phi[maxn];
void get_euler(int n=2e7)
{
    res=0;
    for(int i=2;i<=n;i++)
    {
        if(!phi[i])prime[res++]=i,phi[i]=i-1;
        for(int j=0;j<res&&prime[j]*i<=n;j++)
        {
            if(i%prime[j]==0) 
            {
            	phi[i*prime[j]]=phi[i]*prime[j];
            	break;
            }
            phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
} 
struct BIT 
{
	#define lowbit(x) (x&(-x))
	ll b[500005];
	int n;
	void init(int _n)
	{
		n=_n;
		for(int i=1;i<=n;i++)b[i]=0;
	}
	void update(int x,int v)
	{
		if(x>n)return ;
		while(x<=n)
		{
			b[x]+=v;
			x+=lowbit(x);
		}
	}
	ll query(int x)
	{
		ll ans=0;
		while(x)
		{
			ans+=b[x];
			x-=lowbit(x);
		}
		return ans;
	}
}bit;
P Pow(ll a,int b,int p)
{
	int ans=1,flag=0;
	if(a==1)flag=(p==1);
	else
	{
		ll temp=1;
		for(int i=1;i<=b;i++)
		{
			temp=temp*a;
			if(temp>=p)
			{
				flag=1;
				break;
			}
		}
	}
	a%=p;
	while(b)
	{
		if(b&1)ans=(ll)ans*a%p;
		a=(ll)a*a%p;
		b>>=1;
	}
	return P(ans,flag);
}
P Solve(int l,int r,int p)
{
	ll a=bit.query(l);
	if(a==1||l==r||p==1)return P(a%p,a>=p?1:0);
	P b=Solve(l+1,r,phi[p]);
	if(b.second)b.first+=phi[p];
	return Pow(a,b.first,p);
}
int main()
{
	get_euler();
	scanf("%d%d",&n,&m);
	bit.init(n);
	for(int i=1;i<=n;i++)
	{
		int a;
		scanf("%d",&a);
		bit.update(i,a);
		bit.update(i+1,-a);
	}
	while(m--)
	{
		int op,l,r,p,x;
		scanf("%d%d%d%d",&op,&l,&r,&p);
		if(op==1)bit.update(l,p),bit.update(r+1,-p);
		else printf("%d\n",Solve(l,r,p).first);
	}
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值