CF 1114F Please, another Queries on Array?（欧拉函数&线段树）

F. Please, another Queries on Array?（欧拉函数&线段树）

1114F - Please, another Queries on Array?

题目大意

给出一段序列.现可对其进行两种操作.其一是将一段区间l到r之间的所有数乘上某个值x,其二是求一段区间l到r的$\varphi ({\prod }_l^r{a}_i)$其中$\varphi$为求欧拉函数,由于答案可能很大，故对数取$10^9+7$模

解题思路

欧拉函数为积性函数，且有$\varphi (p^k)=p^k(\frac{p-1}{p})$故题目中原所求可以转换成$\varphi ({\prod }_l^r{a}_i)={\prod }_l^r{a}_i\ast \prod \frac{p-1}{p}$其中p为$a_l$到$a_r$之间所有数的所有质因子，通过线段树我们可以简单地维护出区间之积，和该段区间内所有的质因数(这可以通过状态压缩维护出来)

AC代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef pair<int,LL> pil;
const int size=4e5+5;
const int mod=1e9+7;
struct Tree{
	int l,r;
	int sum;
	int tag1;
	LL tag2;
	LL flag; 
}tree[size<<2];
int arr[size];
inline int lson(int x){return x<<1;}
inline int rson(int x){return x<<1|1;}
int prime[] = {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,283,293};
int inv[62];
inline LL quick_pow(LL a,LL b)
{
	LL ans=1;a=a%mod;
	while(b)
	{
		if(b&1) ans=ans*a%mod;
		b>>=1;
		a=a*a%mod;
	}
	return ans;
}
inline void build(int k,int l,int r)
{
	tree[k].l=l,tree[k].r=r;
	tree[k].tag1=1;
	tree[k].tag2=0;
	if(l==r)
	{
		tree[k].sum=arr[l];
		tree[k].flag=0;
		for(int i=0;i<62;i++) if(arr[l]%prime[i]==0) tree[k].flag|=(1LL<<i);
		return ;
	}
	int mid=(l+r)>>1;
	build(lson(k),l,mid);
	build(rson(k),mid+1,r);
	tree[k].sum=(1LL*tree[lson(k)].sum*tree[rson(k)].sum)%mod;
	tree[k].flag=tree[lson(k)].flag|tree[rson(k)].flag;
}
inline void change(int k)
{
	if(tree[k].l!=tree[k].r)
	{
		int ls=lson(k),rs=rson(k);
		tree[ls].sum=1LL*tree[ls].sum*quick_pow(tree[k].tag1,tree[ls].r-tree[ls].l+1)%mod;
		tree[rs].sum=1LL*tree[rs].sum*quick_pow(tree[k].tag1,tree[rs].r-tree[rs].l+1)%mod;
		tree[ls].tag1=1LL*tree[ls].tag1*tree[k].tag1%mod;
		tree[rs].tag1=1LL*tree[rs].tag1*tree[k].tag1%mod;
		tree[ls].flag|=tree[k].tag2;
		tree[rs].flag|=tree[k].tag2;
		tree[ls].tag2|=tree[k].tag2;
		tree[rs].tag2|=tree[k].tag2;
	}
	tree[k].tag2=0;
	tree[k].tag1=1;
}
inline pair<int,LL> query(int k,int l,int r)
{
	if(tree[k].tag1!=1||tree[k].tag2!=0) change(k);
	if(tree[k].l==l&&tree[k].r==r)
	{
		return pil(tree[k].sum,tree[k].flag);
	}
	int mid=(tree[k].l+tree[k].r)>>1;
	if(r<=mid) return query(lson(k),l,r);
	if(l>=mid+1) return query(rson(k),l,r);
	else  {
		pil L=query(lson(k),l,mid),R=query(rson(k),mid+1,r);
		return pil(1LL*L.first*R.first%mod,L.second|R.second);
	}
}
inline void modify(int k,int l,int r,int x,LL flags)
{
	if(tree[k].tag1!=1||tree[k].tag2!=0) change(k);
	if(tree[k].l==l&&tree[k].r==r)
	{
		tree[k].sum=1LL*tree[k].sum*quick_pow(x,tree[k].r-tree[k].l+1)%mod;
		tree[k].tag1=1LL*tree[k].tag1*x%mod;
		tree[k].flag=tree[k].flag|flags;
		tree[k].tag2=tree[k].tag2|flags; 
		return ;
	}
	int mid=(tree[k].l+tree[k].r)>>1;
	if(l>=mid+1) modify(rson(k),l,r,x,flags);
	else if(r<=mid) modify(lson(k),l,r,x,flags);
	else modify(lson(k),l,mid,x,flags),modify(rson(k),mid+1,r,x,flags);
	tree[k].sum=(1LL*tree[lson(k)].sum*tree[rson(k)].sum)%mod;
	tree[k].flag=tree[lson(k)].flag|tree[rson(k)].flag;
}
int main()
{
	for(int i=0;i<62;i++) inv[i]=quick_pow(prime[i],mod-2);
	int n,q;
	string op;
	cin>>n>>q;
	for(int i=1;i<=n;i++) cin>>arr[i];
	build(1,1,n);
	int l,r,x;
	while(q--)
	{
		cin>>op;
		if(op[0]=='M')
		{
			scanf("%d%d%d",&l,&r,&x);
			LL temp=0;
			for(int i=0;i<62;i++) {if(x%prime[i]==0) temp|=(1LL<<i);}
			modify(1,l,r,x,temp);
		}
		else 
		{
			scanf("%d%d",&l,&r);
			pil temp=query(1,l,r);
			LL ans=temp.first;
			for(int i=0;i<62;i++) if((temp.second>>i)&1)ans=1LL*(ans*(prime[i]-1)%mod)*inv[i]%mod;
			printf("%lld\n",ans);
		}
	}
}