2021icpc网络赛第二场L题

感谢Akari Kito大佬的博客让我写出了这道题。

分析题目提到区间乘和区间询问，一眼看上去就是线段树，主要的难点就是在于怎么去区间维护欧拉函数。（本蒟蒻就死在了欧拉函数上）。

写本题根据欧拉函数的一个性质，

//orz是我储存欧拉函数的数组

if(p是质数)
{
    if(i%p==0)orz[i*p]=orz[i]*p;
    else orz[i*p]=orz[i]*(p-1);
}

仔细观察这个性质，我们现在要区间乘一个数字，根据乘法交换律，我们可以把区间乘的数字拆分成一个个质数，然后再利用上面的这个性质，就可以做到直接在线段树上修改区间维护的欧拉函数。

那么现在的问题来了，我们怎么判断这个区间该直接乘以p还是直接乘以（p-1）呢，或者是这段区间里面既有需要乘以p的又有需要乘以（p-1）的呢。这里有要涉及到势能线段树了，要是你不懂的话可以去找大佬的博客学一学然后再回来写这个题（因为本蒟蒻解释不清楚势能线段树，所以不在这里解释）。

对于区间乘p还是乘（p-1）我们在每个维护区间中加一个bitset<30>tg（100以内的质数只有25个，bitsset的每一个1表示当前这个数是100内第i个质数的倍数，0表示不是第i个质数的倍数），

这样子我们通过当前这段区间的bitset就可以知道该乘p还是该乘（p-1）还是该继续往下递归了；

其实吧bitset换成bool【30】也是可以的，但是bitset比bool快而且一个&操作就行，bool还需要一个for循环，所以我用的是bitset；

我们再建立一个cnt【105】【30】，用来记录1-100分解之后的质因数有哪几个；

下面贴上ac代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+10;
const ll mod=998244353;
bitset<30>st[105]; //st[i][j]表示数字i能否被第j个质数整除
struct node
{
	int l,r;
	bitset<30>tg;
	ll sum,lazy;	
} t[N<<2];
int n,m;
int a[N];
int pri[30];//100以内的质数
int cnt[105][30];//cnt[i][j]表示数字i质因数分解之后含有几个pri[j]
int orz[105];//欧拉函数

int phi(int x)
{
	int ans=x;
	for(int i=2;i*i<=x;i++)
	{
		if(x%i==0)
		{
			ans=ans/i*(i-1);
			while(x%i==0)x/=i;
		}
	}
	if(x>1)ans=ans/x*(x-1);
	return ans;
}

void init()
{
	int tot=0;
	for(int i=2;i<=100;i++)
	{
		bool flag=0;
		for(int j=2;j<=sqrt(i);j++)
		{
			if(i%j==0)
			{
				flag=1;
				break;
			}
		}
		if(!flag)pri[++tot]=i;	
	}
	for(int i=1;i<=100;i++)
	{
		for(int j=1;j<=25;j++)
		{
			int temp=i;
			while(temp%pri[j]==0)
			{
				cnt[i][j]++;
				temp/=pri[j];
			}
			st[i][j]=cnt[i][j];	
		}	
	}	
} 

void pushup(int k)
{
	t[k].sum=(t[k<<1].sum+t[k<<1|1].sum)%mod;
	t[k].tg=t[k<<1].tg&t[k<<1|1].tg;
}

void f(int k,ll v)
{
	t[k].sum=(t[k].sum*v)%mod;
	t[k].lazy=(t[k].lazy*v)%mod;
}

void pushdown(int k)
{
	f(k<<1,t[k].lazy);
	f(k<<1|1,t[k].lazy);
	t[k].lazy=1;
}

void build(int k,int l,int r)
{
	t[k].l=l;t[k].r=r;t[k].lazy=1;
	if(l==r)
	{
		t[k].sum=1ll*orz[a[l]];
		t[k].tg=st[a[l]];
		return ;
	}
	int mid=(l+r)>>1;
	build(k<<1,l,mid);
	build(k<<1|1,mid+1,r);
	pushup(k);
}

void modify(int k,int l,int r,int pos,int cnt)  //pos表示的是100内的第几个质数；
{
	if(l<=t[k].l&&t[k].r<=r&&t[k].tg[pos])
	{
		for(int i=1;i<=cnt;i++)
		{
			t[k].sum=t[k].sum*pri[pos]%mod;
			t[k].lazy=t[k].lazy*pri[pos]%mod;
		}
		return ;
	}
	if(t[k].l==t[k].r)
	{
		t[k].sum=t[k].sum*(pri[pos]-1)%mod;
		t[k].tg[pos]=1;
		for(int i=1;i<cnt;i++)t[k].sum=t[k].sum*pri[pos]%mod;
		return ;
	}
	if(t[k].lazy>1)pushdown(k);
	int mid=(t[k].l+t[k].r)>>1;
	if(l<=mid)modify(k<<1,l,r,pos,cnt);
	if(r>mid)modify(k<<1|1,l,r,pos,cnt);
	pushup(k);
}

ll query(int k,int l,int r)
{
	if(l<=t[k].l&&t[k].r<=r)return t[k].sum%mod;
	if(t[k].lazy>1)pushdown(k);
	ll ans=0;
	int mid=(t[k].l+t[k].r)>>1;
	if(l<=mid)ans=(ans+query(k<<1,l,r))%mod;
	if(r>mid)ans=(ans+query(k<<1|1,l,r))%mod;
	return ans;
}
int main()
{
	init();
	orz[1]=1;
	for(int i=1;i<=100;i++)orz[i]=phi(i);
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)scanf("%d",&a[i]);
	build(1,1,n);
	int pos,l,r,c;
	while(m--)
	{
		scanf("%d%d%d",&pos,&l,&r);
		if(pos==1)
		{
			printf("%lld\n",query(1,l,r));
		}
		else
		{
			scanf("%d",&c);
			for(int i=1;i<=25;i++)
			{
				if(cnt[c][i])modify(1,l,r,i,cnt[c][i]);
			}
		}
	}
	return 0;
}