整体二分——BZOJ2527/Luogu3527 [POI2011]MET-Meteors

题面：洛谷3527 BZOJ2527
这题首先可以想到二分答案但是呢有很多点需要二分
同时每个点的二分状态都差不多
那么我们就一起二分（整体二分）
每次二分出答案以后对于每个国家计算有几颗流星
这个树状数组和线段树都可以维护（我用树状数组调了一个上午QAQ）
然后超过希望总数的点放到左边去，没有的放到右边去
然后分别继续二分
对然后统计答案即可

#include<bits/stdc++.h>using namespace std;typedef long long lll;int n,m,k,nedge=0,p[300001],nex[300001],head[300001];int v[300001],L[300001],R[300001],t[300001],ans[300001];lll p1[300001],p2[300001],rp=0,f[300001];bool jzq[300001];inline int lowbit(int x){return(x&-x);}inline void add(int x,lll s){for(;x<=m;x+=lowbit(x))f[x]+=s;}inline lll sum(int x){lll ans=0;for(;x>0;x-=lowbit(x))ans+=f[x];return ans;}inline void addedge(int a,int b){p[++nedge]=b;nex[nedge]=head[a];head[a]=nedge;}void jia(int x,int v){    int l=L[x],r=R[x],vv=v*t[x];    if(r>=l)add(r+1,-vv),add(l,vv);    else add(1,vv),add(r+1,-vv),add(l,vv);}void erfen(int l,int r,int ll,int rr){    if(rr<ll)return;    if(l==r){        for(int i=ll;i<=rr;i++)ans[p1[i]]=l;        return;    }    int mid=l+r>>1;    while(rp<mid)jia(++rp,1);while(rp>mid)jia(rp--,-1);    lll s;    for(int i=ll;i<=rr;i++){        p2[i]=p1[i];s=0;        for(int k=head[p1[i]];k;k=nex[k]){            s+=sum(p[k]);            if(s>=v[p1[i]])break;        }        jzq[p1[i]]=(s>=v[p1[i]])?1:0;    }    s=ll-1;int ss=0;    for(int i=ll;i<=rr;i++)if(jzq[p2[i]])p1[++s]=p2[i];    for(int i=ll;i<=rr;i++)if(!jzq[p2[i]])p1[++ss+s]=p2[i];    erfen(l,mid,ll,s);    erfen(mid+1,r,s+1,rr);}int main(){    scanf("%d%d",&n,&m);    for(int i=1;i<=m;i++){        int x;scanf("%d",&x);        addedge(x,i);    }    for(int i=1;i<=n;i++)scanf("%d",&v[i]);    scanf("%d",&k);    for(int i=1;i<=k;i++)scanf("%d%d%d",&L[i],&R[i],&t[i]);    L[++k]=1;R[k]=m;t[k]=1e9;    for(int i=1;i<=n;i++)p1[i]=i;    erfen(1,k,1,n);    for(int i=1;i<=n;i++)if(ans[i]!=k)printf("%d\n",ans[i]);    else puts("NIE");    return 0;}