bzoj 3110 K大数查询

第一道整体二分，因为只需要知道每个询问区间中比mid大的数有多少个，就可以直接用线段树区间加，区间求和了。

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cstring>
  4 #include<algorithm>
  5 #define int long long
  6 #define N 50005
  7 #define ls x*2,l,mid
  8 #define rs x*2+1,mid+1,r
  9 using namespace std;
 10 int n,m;
 11 struct node
 12 {
 13     int sign,x,y,c;
 14 }op[N];
 15 int q[N],ans[N];
 16 int a[N*8],lazy[N*8];
 17 void push_down(int x,int l,int mid,int r)
 18 {
 19     if(lazy[x]!=0)
 20     {
 21         lazy[x*2]+=lazy[x];
 22         lazy[x*2+1]+=lazy[x];
 23         a[x*2]+=lazy[x]*(mid-l+1);
 24         a[x*2+1]+=lazy[x]*(r-mid);
 25         lazy[x]=0;
 26     }
 27     return ;
 28 }
 29 int qur(int x,int l,int r,int ll,int rr)
 30 {
 31     if(l>=ll&&r<=rr)return a[x];
 32     int mid=(l+r)>>1;
 33     push_down(x,l,mid,r);
 34     if(rr<=mid)return qur(ls,ll,rr);
 35     if(ll>mid)return qur(rs,ll,rr);
 36     return qur(ls,ll,rr)+qur(rs,ll,rr);
 37 }
 38 void gai(int x,int l,int r,int ll,int rr,int d)
 39 {
 40     if(l>=ll&&r<=rr)
 41     {
 42         lazy[x]+=d;
 43         a[x]+=d*(r-l+1);
 44         return ;
 45     }
 46     int mid=(l+r)>>1;
 47     push_down(x,l,mid,r);
 48     if(ll<=mid)gai(ls,ll,rr,d);
 49     if(rr>mid)gai(rs,ll,rr,d);
 50     a[x]=a[x*2]+a[x*2+1];
 51 }
 52 int tmp[2][N];
 53 void solve(int L,int R,int l,int r)
 54 {
 55     if(L>R)return ;
 56     if(l==r)
 57     {
 58         for(int i=L;i<=R;i++)if(op[q[i]].sign==2)ans[q[i]]=l;
 59         return ;
 60     }
 61     int mid=(l+r)>>1;
 62     int cnt1=0,cnt2=0;
 63     for(int i=L;i<=R;i++)
 64     {
 65         if(op[q[i]].sign==1)
 66         {
 67             if(op[q[i]].c>mid)
 68             {
 69                 gai(1,1,n,op[q[i]].x,op[q[i]].y,1);
 70                 tmp[1][++cnt2]=q[i];
 71             }
 72             else tmp[0][++cnt1]=q[i];
 73         }
 74         else
 75         {
 76             int y=qur(1,1,n,op[q[i]].x,op[q[i]].y);
 77             if(y>=op[q[i]].c)
 78             {
 79                 tmp[1][++cnt2]=q[i];
 80             }
 81             else op[q[i]].c-=y,tmp[0][++cnt1]=q[i];
 82         }
 83     }
 84     for(int i=L;i<=R;i++)
 85     {
 86         if(op[q[i]].sign==1)
 87         {
 88             if(op[q[i]].c>mid)gai(1,1,n,op[q[i]].x,op[q[i]].y,-1);
 89         }
 90     }
 91     int l1=L+cnt1-1;
 92     for(int i=1;i<=cnt1;i++)q[L+i-1]=tmp[0][i];
 93     for(int i=1;i<=cnt2;i++)q[l1+i]=tmp[1][i];
 94     solve(L,l1,l,mid);solve(l1+1,R,mid+1,r);
 95 }
 96 signed main()
 97 {
 98     scanf("%lld%lld",&n,&m);
 99     for(int i=1;i<=m;i++)
100     {
101         scanf("%lld%lld%lld%lld",&op[i].sign,&op[i].x,&op[i].y,&op[i].c);
102         q[i]=i;
103     }
104     solve(1,m,1,n);
105     for(int i=1;i<=m;i++)
106     {
107         if(op[i].sign==2)printf("%lld\n",ans[i]);
108     }
109     return 0;
110 }