AcWing250 磁力块(分块)
作者:互联网
本题仔细分析发现也是一个二维偏序问题,并且因为强制在线,所以不能用莫队来做,因此考虑分块的做法
这道题,因为题目第一个要求是距离小于半径,因此我们直接算出距离后,通过对距离排序,之后分块
这样的分块就会产生一个性质,在某些块中的所有半径都小于r,而在边界块中,则使用暴力。
在块中我们按质量排序
之后我们只需要给每个块分配一个指针,用于判断质量和吸力的关系,这个指针会一直往后推,不用回溯。
这题还可以替换磁力块,其实就相当于bfs,将已经吸入的都可以当作x,y这个点,只要最后队列空了,说明就结束了
#include<bits/stdc++.h> using namespace std; typedef long long ll; const int N=2e5+10; const int mod=1e9+7; struct node{ ll dis,m; ll r,p; }s[2*N+10],d[510][510]; ll sign[5050]; int e[5050]; int st[5050][5050]; int cnt[5050]; int n,block; ll x,y,x2,y2; bool cmpa(node a,node b){ return a.dis<b.dis; } bool cmpb(node a,node b){ return a.m<b.m; } void init(){ int i; for(i=1;i<=n;i++){ cin>>x2>>y2>>s[i].m>>s[i].p>>s[i].r; s[i].dis=(x2-x)*(x2-x)+(y2-y)*(y2-y); s[i].r*=s[i].r; } block=sqrt(n); sort(s+1,s+1+n,cmpa); for(i=1;i<=n;i++){ d[(i-1)/block+1][(i-1)%block]=s[i]; sign[(i-1)/block+1]=s[i].dis; } for(i=1;i<=(n-1)/block+1;i++){ if(i*block<=n) e[i]=block; else e[i]=n%block; sort(d[i],d[i]+e[i],cmpb); } } bool check(ll a,ll b){ return a<=b; } void bfs(){ queue<node> q; q.push(s[0]); st[0][0]=1; while(q.size()){ node t=q.front(); q.pop(); int i; for(i=1;i<=(n-1)/block+1;i++){ if(sign[i]<=t.r){ while(cnt[i]<e[i]&&check(d[i][cnt[i]].m,t.p)){ if(!st[i][cnt[i]]){ st[i][cnt[i]]=1; q.push(d[i][cnt[i]]); } cnt[i]++; } } else{ for(int j=0;j<e[i];j++){ if(d[i][j].dis<=t.r&&d[i][j].m<=t.p){ if(!st[i][j]){ q.push(d[i][j]); st[i][j]=1; } } } break; } } } } int main(){ ios::sync_with_stdio(0); cin>>x>>y>>s[0].p>>s[0].r>>n; s[0].r*=s[0].r; init(); bfs(); int i; ll ans=0; for(i=1;i<=(n-1)/block+1;i++){ for(int j=0;j<e[i];j++){ if(st[i][j]) ans++; } } cout<<ans<<endl; return 0; }View Code
标签:node,AcWing250,分块,int,ll,x2,磁力,y2,5050 来源: https://www.cnblogs.com/ctyakwf/p/12771699.html