AcWing 250 磁力快(分块)
在一片广袤无垠的原野上,散落着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)。
数据范围
\(1≤N≤250000, −10^9≤x,y≤10^9, 1≤m,p,r≤10^9\)
输入样例:
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
分析
将磁石按照质量排序,分成t块,然后对于一个磁石一定有一个k,使得前k-1块所有的磁石的质量都小于它的磁力,k+1块之后的磁石质量都比它磁力大,而对块内按照距离排序之后,只需要挨个扫就可以了,扫过的标记、加入队列,调整区间。但对于第k块,因为不是所有的质量都小于当前磁块的磁力,所以要全部扫一遍,直到距离大于磁力半径,该块不能被调整区间
#include <bits/stdc++.h>
using namespace std;
const int N = 250010;
struct node{
int x,y,m,p,r;
double dis;
}a[N];
int x0,y,n;
int L[N],R[N],v[N];
double S(int x){return 1.0*x * x;}
bool cmp(node a,node b){
return a.m < b.m;
}
bool cmp2(node a,node b){
return a.dis < b.dis;
}
int main(){
scanf("%d%d%d%d%d",&x0,&y,&a[0].p,&a[0].r,&n);
for(int i=1;i<=n;i++){
scanf("%d%d%d%d%d",&a[i].x,&a[i].y,&a[i].m,&a[i].p,&a[i].r);
a[i].dis = sqrt(S(a[i].x-x0)+S(a[i].y-y));
}
int t = sqrt(n);
for(int i=1;i<=t;i++){
L[i] = (i-1) * t + 1;
R[i] = i * t;
}
if(R[t] < n){
t ++;L[t] = R[t-1] + 1;R[t] =n;
}
sort(a+1,a+n+1,cmp);
for(int i=1;i<=t;i++){
sort(a+L[i],a+R[i]+1,cmp2);
}
int res = 0;
queue<int> q;
q.push(0);
while(!q.empty()){
int u = q.front();q.pop();
res++;
int r = a[u].r,p = a[u].p;
bool flag = true;//标记是否为前k-1块
for(int i=1;i<=t;i++){
if(!flag)break;
int j = L[i];
for(;j<=R[i];j++){
if(v[j])continue;
if(p >= a[j].m){
if(r >= a[j].dis){
v[j] = 1;q.push(j);
}else{
break;
}
}else{
flag = false;continue;//是第k块
}
}
if(flag)L[i] = j;//不是第k块就调整左坐标点
}
}
printf("%d\n",res-1);//减去原始的磁石
return 0;
}
注:转载请注明出处