CF920F SUM and REPLACE 题解

题目链接

本题是很典型的一道 线段树/分块与暴力结合的题
与本题类似的还有如下几题：
SP2713 GSS4
P4145 上帝造题的七分钟2 / 花神游历各国
CF438D The Child and Sequence

个人本题是这几题中包装的最好的，也是唯一能交题解的一道

本题的修改操作非常不好处理，也不容易合并，比较难用传统的线段树方法结局。但是这里的修改操作有一个特点，就是有效的修改次数非常少

显然地，一个数 $n$ 最多有 $2\sqrt{n}$ 个因数，也就是说题目中的 $D(i)\le 2\sqrt{i}$
那么，如果每次让 $a_i=D(a_i)$，就相当于让 $a_i$ 开了一个根号
这样的话，$a_i$ 会快速地减小。假如 $a_i$ 等于 $10^6$（题目中的最大值），那么在不超过 $5$ 次操作后，$a_i$ 就变成 $2$ 了，并且以后无论修改多少次，$a_i$ 永远是 $2$。

那么，每次修改时我们可以暴力地在线段树中递归，如果当前递归到的区间的最大值 $\le 2$，就停止递归，否则就递归到边界进行单点修改。时间复杂度 $\mathcal{O}(n\log n)$

为了美观以及方便理解，这里只贴出未开long long的代码

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
const int Maxn=3e5+10;
const int Maxm=Maxn<<2;
int sum[Maxm],maxv[Maxm];
int a[Maxn],f[1000010];
int n,m;
inline int read()
{
	int s=0,w=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
	while(ch>='0' && ch<='9')s=(s<<3)+(s<<1)+(ch^48),ch=getchar();
	return s*w;
}
inline void push_up(int k)
{
	sum[k]=sum[k<<1]+sum[k<<1|1];
	maxv[k]=max(maxv[k<<1],maxv[k<<1|1]);
}
void modify(int k,int l,int r,int x,int y)
{
	if(maxv[k]<=2)return;
	if(l==r)
	{
		a[l]=f[a[l]];
		sum[k]=maxv[k]=a[l];
		return;
	}
	int mid=(l+r)>>1;
	if(x<=mid)modify(k<<1,l,mid,x,y);
	if(mid<y)modify(k<<1|1,mid+1,r,x,y);
	push_up(k);
}
int query(int k,int l,int r,int x,int y)
{
	if(x<=l && r<=y)return sum[k];
	int mid=(l+r)>>1,ret=0;
	if(x<=mid)ret=query(k<<1,l,mid,x,y);
	if(mid<y)ret+=query(k<<1|1,mid+1,r,x,y);
	return ret;
}
void build(int k,int l,int r)
{
	if(l==r)
	{
		sum[k]=maxv[k]=a[l];
		return;
	}
	int mid=(l+r)>>1;
	build(k<<1,l,mid);
	build(k<<1|1,mid+1,r);
	push_up(k);
}
int main()
{
	// freopen("in.txt","r",stdin);
	n=read(),m=read();
	for(int i=1;i<=n;++i)
	a[i]=read();
	build(1,1,n);
	int cnt=1e6;
	for(int i=2;i<=cnt;++i)
	{
		for(int j=i;j<=cnt;j+=i)
		++f[j];
		++f[i];
	}
	while(m--)
	{
		int opt=read(),x=read(),y=read();
		if(opt==1)modify(1,1,n,x,y);
		else printf("%d\n",query(1,1,n,x,y));
	}
	return 0;
}