1. 程式人生 > >UVA 12003 分塊 區間計數+單點修改

UVA 12003 分塊 區間計數+單點修改

/**
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;
}