bzoj3276磁力 两种要求下的最大值：分块or线段树+拓扑

进阶指南上的做法是分块的。。

但是线段树搞起来也挺快，将磁石按照距离排序，建立线段树，结点维护区间质量最小值的下标

进行拓扑，每次在可行的范围内在线段树中找到质量最小的下标取出，取出后再将线段树对应的点设置成0

查询时找区间不为0最小值的下标即可

#include<cstdio>
#include<algorithm>
#define N 250010
typedef long long ll;
int n,i,x0,y0,nowp,x,y,r,c,v[N<<],tmp,h=,t,q[N];ll nowr;
struct P{int m,p;ll d,r;}a[N];
inline bool cmp(P x,P y){return x.d<y.d;}
inline void read(int&a){
  char c;bool f=;a=;
  while(!((((c=getchar())>='')&&(c<=''))||(c=='-')));
  if(c!='-')a=c-'';else f=;
  while(((c=getchar())>='')&&(c<=''))(a*=)+=c-'';
  if(f)a=-a;
}
inline ll sqr(ll x){return x*x;}
inline int lower(){
  int l=,r=n,t=,mid;
  while(l<=r)if(a[mid=(l+r)>>].d<=nowr)l=(t=mid)+;else r=mid-;
  return t;
}
inline int merge(int x,int y){
  if(!x)return y;
  if(!y)return x;
  return a[x].m<a[y].m?x:y;
}
inline void up(int x){v[x]=merge(v[x<<],v[x<<|]);}
void build(int x,int a,int b){
  if(a==b){v[x]=a;return;}
  int mid=(a+b)>>;
  build(x<<,a,mid),build(x<<|,mid+,b),up(x);
}
void change(int x,int a,int b,int c){
  if(a==b){v[x]=;return;}
  int mid=(a+b)>>;
  c<=mid?change(x<<,a,mid,c):change(x<<|,mid+,b,c);
  up(x);
}
void ask(int x,int a,int b){
  if(b<=c){tmp=merge(tmp,v[x]);return;}
  int mid=(a+b)>>;
  ask(x<<,a,mid);
  if(c>mid)ask(x<<|,mid+,b);
}
int main(){
  read(x0),read(y0),read(nowp),read(r),read(n),nowr=sqr(r);
  for(i=;i<=n;i++){
    read(x),read(y),read(a[i].m),read(a[i].p),read(r);
    a[i].d=sqr(x-x0)+sqr(y-y0),a[i].r=sqr(r);
  }
  std::sort(a+,a+n+,cmp),build(,,n);
  if(c=lower())while(){
    tmp=,ask(,,n);
    if(!tmp||a[tmp].m>nowp)break;
    change(,,n,q[++t]=tmp);
  }
  while(h<=t){
    nowp=a[q[h]].p,nowr=a[q[h++]].r;
    if(c=lower())while(){
      tmp=,ask(,,n);
      if(!tmp||a[tmp].m>nowp)break;
      change(,,n,q[++t]=tmp);
    }
  }
  return printf("%d",t),;
}