Contest Hunter #46 T1 磁力块 [分块]

描述

在一片广袤无垠的原野上，散落着N块磁石。每个磁石的性质可以用一个五元组(x,y,m,p,r)描述，其中x,y表示其坐标，m是磁石的质量，p是磁力，r是吸引半径。若磁石A与磁石B的距离不大于磁石A的吸引半径，并且磁石B的质量不大于磁石A的磁力，那么A可以吸引B。
小取酒带着一块自己的磁石L来到了这篇原野的(x0,y0)处，我们可以视为磁石L的坐标为(x0,y0)。小取酒手持磁石L并保持原地不动，所有可以被L吸引的磁石将会被吸引过来。在每个时刻，他可以选择更换任意一块自己已经获得的磁石（当然也可以是自己最初携带的L磁石）在(x0,y0)处吸引更多的磁石。小取酒想知道，他最多能获得多少块磁石呢？

输入格式

第一行五个整数x0,y0,pL,rL,N，表示小取酒所在的位置，磁石L磁力、吸引半径和原野上散落磁石的个数。
接下来N行每行五个整数x,y,m,p,r，描述一块磁石的性质。

输出格式

输出一个整数，表示最多可以获得的散落磁石个数（不包含最初携带的磁石L）。

样例输入

0 0 5 10 5
5 4 7 11 5
-7 1 4 7 8
0 2 13 5 6
2 -3 9 3 4
13 5 1 9 9

样例输出

3

数据范围与约定

• 对于30%的数据，1<=N<=1000。
• 对于100%的数据，1<=N<=250000，-10^9<=x,y<=10^9，1<=m,p,r<=10^9。

将石块按距离分块，每块内按质量排序。挨个尝试用已有的石头吸引剩下的石头即可。

一直觉得分块就是稍有些技巧的暴力……

事实上这题的标解分块还真不如我的花式暴力跑得快23333

花式暴力：http://www.cnblogs.com/SilverNebula/p/5929674.html

 /*by SilverN*/
 #include<iostream>
 #include<algorithm>
 #include<cstring>
 #include<cstdio>
 #include<cmath>
 using namespace std;
 const int w=;
 const int mxn=;
 int read(){
     int x=,f=;char ch=getchar();
     while(ch<'' || ch>''){if(ch=='-')f=-;ch=getchar();}
     while(ch>='' && ch<=''){x=x*+ch-'';ch=getchar();}
     return x*f;
 }
 struct node{
     double dis;
     int m,p,r;
     bool got;
 }p[mxn];
 int hd[mxn];
 int cmpd(const node a,const node b){
     return a.dis<b.dis;
 }
 int cmpm(const node a,const node b){
     return a.m<b.m;
 }
 //
 int L[mxn],R[mxn],tot=;
 double D[mxn];
 //
 struct tool{
     int p,r;
 }t[mxn];int cnt=;
 int head=;
 int n;
 int x0,yy0;
 double dist(int x1,int y1,int x2,int y2){
     return sqrt( (double)(x1-x2)*(x1-x2)+(double)(y1-y2)*(y1-y2) );
 }
 int main(){
     x0=read();yy0=read();
     int i,j;
     t[].p=read();t[].r=read();
     n=read();
     int x,y;
     for(i=;i<=n;i++){
         x=read();y=read();
         p[i].dis=dist(x,y,x0,yy0);
         p[i].m=read();    p[i].p=read();    p[i].r=read();
         p[i].got=;
     }
     sort(p+,p+n+,cmpd);
     for(i=;i<=n;i+=w){
         L[++tot]=i;R[tot]=min(n,i+w-);D[tot]=p[R[tot]].dis;
         sort(p+L[tot],p+R[tot]+,cmpm);
     }
     while(head<=cnt){
         tool now=t[head];
         for(i=;i<=tot;i++){
             if(D[i]>now.r){
                 for(j=L[i];j<=R[i];j++)
                     if(!p[j].got && p[j].dis<=now.r && p[j].m<=now.p){
                         t[++cnt].p=p[j].p;
                         t[cnt].r=p[j].r;
                         p[j].got=;
                     }
                 break;
             }
             while(L[i]<=R[i] && p[L[i]].m<=now.p){
                 if(!p[L[i]].got){
                     p[L[i]].got=;
                     t[++cnt].p=p[L[i]].p;
                     t[cnt].r=p[L[i]].r;
                 }
                 L[i]++;
             }
         }
         head++;
     }
     printf("%d\n",cnt-);
     return ;
 }