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=700;
const int mxn=270000;
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){x=x*10+ch-'0';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=0;
double D[mxn];
//
struct tool{
    int p,r;
}t[mxn];int cnt=1;
int head=1;
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[1].p=read();t[1].r=read();
    n=read();
    int x,y;
    for(i=1;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=0;
    }
    sort(p+1,p+n+1,cmpd);
    for(i=1;i<=n;i+=w){
        L[++tot]=i;R[tot]=min(n,i+w-1);D[tot]=p[R[tot]].dis;
        sort(p+L[tot],p+R[tot]+1,cmpm);
    }
    while(head<=cnt){
        tool now=t[head];
        for(i=1;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=1;
                    }
                break;
            }
            while(L[i]<=R[i] && p[L[i]].m<=now.p){
                if(!p[L[i]].got){
                    p[L[i]].got=1;
                    t[++cnt].p=p[L[i]].p;
                    t[cnt].r=p[L[i]].r;
                }
                L[i]++;
            }
        }
        head++;
    }
    printf("%d\n",cnt-1);
    return 0;
}