洛谷P3987 我永远喜欢珂朵莉～树状数组+vector（暴力）

题目链接：我永远喜欢珂朵莉～

在太阳西斜的这个世界里，置身天上之森。等这场战争结束之后，不归之人与望眼欲穿的众人，人人本着正义之名，长存不灭的过去、逐渐消逝的未来。我回来了，纵使日薄西山，即便看不到未来，此时此刻的光辉，盼君勿忘。

珂朵莉最珂爱了，珂朵莉的题最毒瘤了qwq
不过这道题被我暴力水过去了qwq（以此纪念昨天才补完的末日时小说）
总之先放图qvq

洛咕题解里面好像也有用暴力水过的……不过他们好像用了一些奇怪的操作……比如

vector<vector<int>::iterator> xxx;

并没有看懂……不过从思路来看感觉和我写的差不多……

首先打个表，发现500000500000500000以内数的因数个数最多是200200200个qwq
所以这暗示着珂以暴力水过QWQ
一开始我思路是用一个链表模拟题目的数组，如果有数≤1\le1≤1了就删掉。
但是这样询问仍然是O(n)……O(n)……O(n)……而且发现000的情况很难处理……比较难受……
考虑对于每个iii，用vectorvectorvector记录所有数组中为iii的倍数的下标（从小到大）
假设所有为vvv的倍数的下标都存在mul[v]mul[v]mul[v]里面。
因为是有序的，所以珂以找出一段连续的区间，满足下标都在[l,r][l,r][l,r]区间内。
然后暴力修改即可qwq，如果某个数操作后不是vvv的倍数了，那么把这个数删除。
考虑怎样修改：
因为修改是暴力，所以每次都是单点修改，最后区间查询，所以想到树状数组。
所以修改的时候在树状数组上维护一下即可。

两个细节：

1.1.1.如果从左到右遍历，删除了一个数之后这个数右边的所有元素会向左移……下标也会变……
洛咕上的题解说是把要删掉的元素记下来……最后再删……
其实从右到左遍历就珂以了qwq！{\color{red}\text{其实从右到左遍历就珂以了qwq！}}其实从右到左遍历就珂以了qwq！
因为右到左遍历，即使删除了一个元素，对左边的下标不会有影响qwq

2.2.2.假设一个修改操作是(l,r,v)(l,r,v)(l,r,v)，那么程序只会对mul[v]mul[v]mul[v]进行修改。
假设有一个数能被2v2v2v整除，但它除以vvv后就珂能不能被2v2v2v整除，但这里并没有在mul[2v]mul[2v]mul[2v]中删除这个数qwq
所以暴力修改的时候要判断一下，如果当前数不能被vvv整除，那么continuecontinuecontinue，不做操作。
注意不要把这个数删除……不然会因为常数过大而TLE....TLE....TLE....（珂能是eraseeraseerase常数大？）

毒瘤代码

#include<stdio.h>
#include<cstring>
#include<algorithm>
#include<vector>
#include<math.h>
#define re register int
using namespace std;
typedef long long ll;
int read() {re x=0,f=1;char ch=getchar();while(ch<'0' || ch>'9') {if(ch=='-')  f=-1;ch=getchar();}while(ch>='0' && ch<='9') {x=10*x+ch-'0';ch=getchar();}return x*f;
}
const int Size=500005;
int n,m,a[Size];
vector<int> mul[Size];
void divi(int x,int id) {int sx=sqrt(x);for(re i=1; i<=sx; i++) {if(!(x%i)) {mul[i].push_back(id);if(i*i!=x)   mul[x/i].push_back(id);}}
}
ll tree[Size];
inline int lowbit(int x) {return x&(-x);
}
inline void update(int x,int v) {for(re i=x; i<=n; i+=lowbit(i)) {tree[i]+=v;}
}
inline ll query(int x) {ll ans=0;for(re i=x; i; i-=lowbit(i)) {ans+=tree[i];}return ans;
}
void Update(int l,int r,int x) {if(x==1 || mul[x].empty())    return;int le=lower_bound(mul[x].begin(),mul[x].end(),l)-mul[x].begin();int ri=upper_bound(mul[x].begin(),mul[x].end(),r)-mul[x].begin()-1;for(re i=ri; i>=le; i--) {int now=mul[x][i];if(a[now]%x) {//注意：这时不要erase，否则会TLE.. //mul[x].erase(mul[x].begin()+i);continue;}update(now,a[now]/x-a[now]);a[now]/=x;if(a[now]%x)  mul[x].erase(mul[x].begin()+i);}
}
int main() {n=read();m=read();for(re i=1; i<=n; i++) {a[i]=read();//分解质因数，并统计muldivi(a[i],i);//计入树状数组update(i,a[i]);}while(m--) {int op=read();int l=read();int r=read();if(op==1) {int x=read();Update(l,r,x);} else {printf("%lld\n",query(r)-query(l-1));}}return 0;
}