P3987 我永远喜欢珂朵莉~

最新推荐文章于 2022-01-13 15:02:30 发布

crazy morning

最新推荐文章于 2022-01-13 15:02:30 发布

阅读量95

点赞数

分类专栏：线段树

本文链接：https://blog.csdn.net/weixin_49843717/article/details/117590398

版权

线段树专栏收录该内容

10 篇文章 1 订阅

订阅专栏

这篇博客介绍了如何使用线段树和因数分解的方法解决一个区间修改的问题。作者首先指出暴力求解无法通过，然后提出了一个优化方案，即为每个数的因数创建一个vector来存储可以整除的下标，并在修改时从这些向量中截取并更新。同时，考虑到某些除法可能改变数的因数性质，需要额外维护。最后，博主给出了线段树的实现代码，并展示了如何结合因数向量进行区间修改和查询操作。

摘要由CSDN通过智能技术生成

题目

思路

~~Ynoi的题居然不是分块？？~~
本题是大学的弱化版~~但鄙人的代码被大学卡掉了~~
但是这题还是过得去的

显然暴力过不去~~这不废话嘛~~
所以我们考虑用另一种做法来表示下标。
考虑对所有数的因数分别开一个vector（平衡树也可以，但是……长而且难打，所以这题算是C党福利？）存a中可以被整除的下标，然后每次修改就从这里截取一段出来修改，然后除完还可以被整除的再放回去。
但也有一种情况：16/4=4，这里除完就不能被8整除了，所以我们维护时还要再考虑现在能不能除，能除再考虑。
然后用线段树维护区间和就好。

线段树是远古码风，请结合更新函数内部代码，查询函数内部代码理解
code:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<vector>
using namespace std;
inline int read()
{
	int x=0,f=1;
	char c;
	c=getchar();
	while (c!='-'&&(c>'9'||c<'0')) c=getchar();
	if (c=='-')
	{
		f=-1;
		c=getchar();
	}
	while (c<='9'&&c>='0')
	{
		x=x*10+c-'0';
		c=getchar();
	}
	return x*f;
}
int n,m,a[100005];
vector<int> s[500005];
long long tree[2000008],lazy[2000008];
int f,L,R;
void jianshu(int l,int r,int id)
{
	if (l==r)
	{
		tree[id]=a[l];
		return;
	}
	jianshu(l,(l+r)>>1,id*2);
	jianshu(((l+r)>>1)+1,r,id*2+1);
	tree[id]=tree[id<<1]+tree[(id<<1)+1];
	return; 
}
void down(int l,int r,int id)
{
	if (lazy[id])
	{
		int mid=(l+r)>>1;
		lazy[id*2]+=lazy[id],lazy[id*2+1]+=lazy[id];
		tree[id*2]+=lazy[id]*(mid-l+1);
		tree[id*2+1]+=lazy[id]*(r-mid); 
		lazy[id]=0;
	}
	return;
}
void qjgx(int l,int r,int id,int x)
{
	if (l>=L&&R>=r)
	{
		tree[id]+=x*(r-l+1);
		lazy[id]+=x;
		return;
	}
	down(l,r,id);
	int mid=(l+r)>>1;
	if (L<=mid) qjgx(l,mid,id*2,x);
	if (R>mid) qjgx(mid+1,r,id*2+1,x);
	tree[id]=tree[id*2]+tree[id*2+1];
	return;
}
long long qjcx(int l,int r,int id)
{
	if (l>=L&&R>=r)
	{
		return tree[id];
	}
	down(l,r,id);
	int mid=(l+r)>>1;
	long long ans=0;
	if (L<=mid) ans+=qjcx(l,mid,id*2);
	if (R>mid) ans+=qjcx(mid+1,r,id*2+1);
	return ans;
}
vector<int>::iterator x,y,it;
vector<vector<int>::iterator> d;
int main()
{
	n=read(),m=read();
	for (int i=1;i<=n;i++)
	{
		a[i]=read();
		for (int j=1;j*j<=a[i];j++)
		{
			if (a[i]%j) continue;
			s[j].push_back(i);
			if (a[i]/j!=j) s[a[i]/j].push_back(i);
		}
	}
	jianshu(1,n,1);
	while (m--)
	{
		f=read(),L=read(),R=read();
		if (f==1)
		{
			f=read();
			if (f==1||s[f].empty()) continue;
			x=lower_bound(s[f].begin(),s[f].end(),L);
			if (x==s[f].end()) continue;
			y=upper_bound(s[f].begin(),s[f].end(),R);
			d.clear();
			for (it=x;it!=y;it++)
			{
				if (a[*it]%f==0)
				{
					int wj=a[*it]/f;
					L=*it,R=*it;
					qjgx(1,n,1,-a[*it]+wj);
					a[*it]=wj;
					if (a[*it]%f) d.push_back(it);
				}
			}
			for (int i=d.size()-1;i>0;i--)
			{
				s[f].erase(d[i]);
			}
		}
		else
		{
			printf("%lld\n",qjcx(1,n,1));
		}
	}
	return 0;
}