版权声明:存在错误或者不清楚的地方还望指出 https://blog.csdn.net/hypHuangYanPing/article/details/82947151
/**
UVA 12003
链接:https://vjudge.net/problem/UVA-12003
区间计数+单点修改
*/
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=5e5+7;
int n,q,blo,pos[maxn];
ll s[maxn];
vector<ll>vec[800];
int main(){
int n,m,u;scanf("%d %d %d",&n,&m,&u);
blo=(int)sqrt(n*1.0);
for(int i=1;i<=n;i++) scanf("%lld",&s[i]);
for(int i=1;i<=n;i++) pos[i]=(i-1)/blo+1;
for(int i=1;i<=n;i++) vec[pos[i]].push_back(s[i]);
for(int i=1;i<=pos[n];i++) sort(vec[i].begin(),vec[i].end());
while(m--){
int l,r,v,p;scanf("%d %d %d %d",&l,&r,&v,&p);
int ret=0;
for(int i=pos[l]+1;i<=pos[r]-1;i++) ret+=lower_bound(vec[i].begin(),vec[i].end(),v)-vec[i].begin();
if(pos[l]!=pos[r]){
for(int i=l;i<=blo*pos[l];i++) if(s[i]<v) ret++;
for(int i=(pos[r]-1)*blo+1;i<=r;i++) if(s[i]<v) ret++;
}
else{
for(int i=l;i<=r;i++) if(s[i]<v) ret++;
}
ll val=1ll*u*ret/(r-l+1);
vec[pos[p]].clear();
s[p]=val;
for(int i=(pos[p]-1)*blo+1;i<=min(n,pos[p]*blo);i++) vec[pos[p]].push_back(s[i]);
sort(vec[pos[p]].begin(),vec[pos[p]].end());
}
for(int i=1;i<=n;i++)
printf("%lld\n",s[i]);
return 0;
}