CDQ分治&整体二分练习2——P3527 [POI2011]MET-Meteors

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/Love_mona/article/details/81283629

蒟蒻的垂死挣扎

这题就是整体二分的裸题了，把修改和询问离线下来丢进去二分就好了。

注意一下： 1.这题是个环     2.这题多了一个hack数据是中途会暴long long的

第一个特判一下就好了，第二个的话因为询问有值域的存在所以超出值域的直接break出来即可。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<queue>
#include<map>
#include<set>
#include<vector>
#define ll long long
#define MOD 1000000007
#define N 300010
#define RG register

using namespace std;

inline ll read(){
    RG ll x=0,o=1; RG char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
    if(ch=='-') o=-1,ch=getchar();
    while(ch>='0'&&ch<='9') x=((x<<3)+(x<<1))+ch-'0',ch=getchar();
    return x*o;
}

vector<int> V[N];
int n,m,k,ans[N];
struct mona { int l,r; ll w; } s[N];
struct joker{ ll p; int id; } q[N],q1[N],q2[N];
struct Tree { ll tr[N];
    inline void Modify(int x,ll val){
        for(RG int i=x;i<=m;i+=(i&-i)) tr[i]+=val;
    }
    inline ll Query(int x){ ll sum=0;
        for(RG int i=x;i;i-=(i&-i)) sum+=tr[i]; return sum;
    }
} T;
inline void Fall(int i,int op){
    if(s[i].l>s[i].r) T.Modify(1,s[i].w*op);
    T.Modify(s[i].l,s[i].w*op),T.Modify(s[i].r+1,-s[i].w*op);
}
inline void Divide(int l,int r,int ql,int qr){
    if(ql>qr) return ;
    if(l==r) { for(RG int i=ql;i<=qr;++i) ans[q[i].id]=l; return ; }
    RG int mid=(l+r)>>1; for(RG int i=l;i<=mid;++i) Fall(i,1);
    RG int t1=0,t2=0;
    for(RG int i=ql;i<=qr;++i){
        RG ll sum=0; 
        for(RG int j=0,l=V[q[i].id].size();j<l;++j){
            sum+=T.Query(V[q[i].id][j]);
            if(sum>1e9) break ;
        }
        if(sum>=q[i].p)  q1[++t1]=q[i];
        else q[i].p-=sum,q2[++t2]=q[i];
    }   for(RG int i=l;i<=mid;++i) Fall(i,-1);
    for(RG int i=1;i<=t1;++i) q[i+ql-1]=q1[i];
    for(RG int i=1;i<=t2;++i) q[i+ql+t1-1]=q2[i];
    Divide(l,mid,ql,ql+t1-1),Divide(mid+1,r,ql+t1,qr);
}

int main(){
    n=read(),m=read();
    for(RG int i=1;i<=m;++i) V[read()].push_back(i);
    for(RG int i=1;i<=n;++i) q[i]=(joker) { read(),i } ;
    k=read();
    for(RG int i=1;i<=k;++i) s[i]=(mona) { read(),read(),read() };
    s[++k]=(mona) { 1,m,2e9 };
    Divide(1,k,1,n);
    for(RG int i=1;i<=n;++i) if(ans[i]<k) printf("%d\n",ans[i]); else puts("NIE");
}