小花梨的数组

ProblemE

小花梨的数组

时间限制：1000ms

空间限制：512MB

Description

小花梨得到了一个长度为n的数组a，现在要对它进行三种操作：

⚫1 l  r对所有的i∈[l,r],a[i]=a[i]∗minprime(a[i])

⚫2 l r对所有的i∈[l,r],a[i]=a[i]/minprime(a[i])

⚫3x求a[x]的值

minprime(x)= 1                      (x=1) 

                      x的最小素因子(x≥2)

现在给出初始数组a，对其进行m次操作，对于第三种操作输出结果。

Input

第一行输入两个正整数n,m，表示数组的长度以及操作次数(1≤n,m≤100000)

第二行输入n个正整数表示数组푎(1≤ai≤1000000)

接下来m行表示m次操作，每行输入格式为"1  l   r"或者"2  l  r"，或者"3x"，对应上述三种操作。1≤l,r,x≤n,l≤r

Output

对于第三种操作输出答案即可，答案对109+7进行取模。

先要求出所有的素因子，然后用线段树维护的方法求出除法的次数和乘法的次数

对于每一次除法，我们先判断这个数之前有没有被乘过 如果被乘过，那么就抵消一次乘法。

详细见代码

AC代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int MOD=1e9+7;
const int maxn=1e5+5;
ll k,n,m,mul[maxn<<2],divn[maxn<<2],a[maxn];
vector<ll>f[maxn];
long long  noprime[10000005];//表示某一个数是不是素数
long long prime[10000005];//表示第几个素数是...
void lln()
{
	noprime[0]=1;//0不是素数
	noprime[1]=1;//1不是素数
	k=0;
	for(ll i=2; i<=1000015; i++)
	{
		if(noprime[i]==0)//如果i是素数
			prime[k++]=i;//保存在prime里
		for(ll j=0; j<k&&i*prime[j]<=1000015; j++)//从第一个素数到已知的最后一个素数
		{
			noprime[i*prime[j]]=1;//合数一定不是素数
			if(i%prime[j]==0) break;//即比一个合数大的质数和该合数的乘积可用一个更大的合数和比其小的质数相乘得到
		}
	}
}
void pushdown(ll l,ll r,ll res)
{
	if(mul[res]||divn[res])
	{
		if(mul[res<<1]>=divn[res]) mul[res<<1]-=divn[res];//如果左儿子的乘法次数大于父亲的除法次数就直接减掉 
		else divn[res<<1]+=(divn[res]-mul[res<<1]),mul[res<<1]=0;//不然的话 就在左儿子的除法次数里加上剩下的除法次数 乘法次数就变成了0 
		if(mul[res<<1|1]>=divn[res]) mul[res<<1|1]-=divn[res];//右儿子同上 
		else divn[res<<1|1]+=(divn[res]-mul[res<<1|1]),mul[res<<1|1]=0;
		mul[res<<1|1]+=mul[res];//除法处理完了就处理乘法 乘法直接加上即可 
		mul[res<<1]+=mul[res];
		mul[res]=0;
		divn[res]=0;
	}
}
void update(ll l,ll r,ll res,ll L,ll R,ll p)
{
	if(L<=l&&R>=r)
	{
		if(p==1) mul[res]++;//如果是乘法就直接加上 
		else//如果是除法 
		{
			if(mul[res]!=0) mul[res]--;//如果乘法次数有多 ，那么直接减掉一次乘法即可 
			else divn[res]++;//不然就加上一次除法 
		}
		return;
	}
	pushdown(l,r,res);
	ll mid=(l+r)/2;
	if(L<=mid) update(l,mid,res<<1,L,R,p);
	if(R>mid) update(mid+1,r,res<<1|1,L,R,p);
}
void query(ll l,ll r,ll res,ll& mull,ll& divv,ll x)
{
	if(l==r)
	{
		mull=mul[res];
		divv=divn[res];
		return ;
	}
	pushdown(l,r,res);
	ll mid=(l+r)/2;
	if(x<=mid) query(l,mid,res<<1,mull,divv,x);
	else query(mid+1,r,res<<1|1,mull,divv,x);
}
ll pow(ll a, ll b, ll m)
{
	ll ans = 1;
	while(b)
	{
		if(b & 1)ans = (a * ans) % m;
		b >>= 1;
		a = (a * a) % m;
	}
	return ans;
}
int  main()
{
	lln();//先打一个素数表 
	scanf("%lld%lld",&n,&m);
	for(ll i=1; i<=n; i++)
	{
		scanf("%lld",&a[i]);
		ll x=a[i];
		for(ll j=0; j<k&&prime[j]*prime[j]<=x; j++)
		{
			while(x%prime[j]==0)
			{
				x/=prime[j];
				f[i].push_back(prime[j]);//把a[i]所有能整除的素数输入到vector中 
			}
		}
		if(x!=1) f[i].push_back(x);
	}
	while(m--)
	{
		ll op,x,l,r;
		scanf("%lld",&op);
		if(op==1)
		{
			scanf("%lld%lld",&l,&r);
			update(1,n,1,l,r,1);//1代表乘法操作 
		}
		else if(op==2)
		{
			scanf("%lld%lld",&l,&r);
			update(1,n,1,l,r,2);//2代表除法操作 
		}
		else
		{
			scanf("%lld",&x);
			ll mull=0,divv=0;
			query(1,n,1,mull,divv,x);//求出这个点被除了几次 乘了几次 
			if(divv>=f[x].size()) printf("1\n");//如果被除的次数大于他的所有素因子数 就可以直接输出1 
			else
			{
				ll ans=1,temp=f[x][divv];
				for(ll i=divv; i<f[x].size(); i++)
				{
					ans*=f[x][i];//后面的这些数是没有被影响的，所以可以直接乘上去 
				}
				ans=ans*pow(temp,mull,MOD)%MOD;//在用快速幂 求出乘的次数 乘上ans 可以得出最后的答案 
				printf("%lld\n",ans);
			}
		}
	}
	return 0;
}