解题思路:
先二分答案,那么每个点就有了一个匹配范围。
考虑轨道上匹配点的偏转角度,肯定是一个匹配点抵住一个点的匹配边界上,所以一共只有2n种匹配方式,每次暴力网络流是
O(n4logn)
O
(
n
4
l
o
g
n
)
注意到偏转角度不会超过 2π/n 2 π / n ,即最多只会一项每个点一条匹配边的出现,另一条的消失,所以我们可以用扫描线的方法,将操作排序后处理增广和退流操作即可,一次就只会是 O(n) O ( n ) 的复杂度。
由于是二分图,所以退流只会影响三条边,还是很好写的。
时间复杂度 O(n3logn) O ( n 3 l o g n )
#include<bits/stdc++.h>
using namespace std;
int getint()
{
int i=0,f=1;char c;
for(c=getchar();(c!='-')&&(c<'0'||c>'9');c=getchar());
if(c=='-')c=getchar(),f=-1;
for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
return i*f;
}
const int N=405,M=100005,INF=1e9;
const double Pi=acos(-1),eps=1e-9;
int n,S,T,tot,first[N],cur[N],dis[N],nxt[M],to[M],cap[M];
double rad,blk;queue<int>que;
struct point{double x,y;}p[N];
struct data
{
double t;int u,v,op;
data(double _t=0,int _u=0,int _v=0,int _op=0):t(_t),u(_u),v(_v),op(_op){}
inline friend bool operator < (const data &a,const data &b)
{
if(a.t==b.t)return a.op>b.op;
return a.t<b.t;
}
}q[N];
void add(int x,int y,int z)
{
nxt[++tot]=first[x],first[x]=tot,to[tot]=y,cap[tot]=z;
nxt[++tot]=first[y],first[y]=tot,to[tot]=x,cap[tot]=0;
}
bool bfs()
{
for(int i=S;i<=T;i++)dis[i]=-1,cur[i]=first[i];
dis[S]=0,que.push(S);
while(!que.empty())
{
int u=que.front();que.pop();
for(int e=first[u];e;e=nxt[e])
if(dis[to[e]]==-1&&cap[e])dis[to[e]]=dis[u]+1,que.push(to[e]);
}
return dis[T]!=-1;
}
int dfs(int u,int flow)
{
if(u==T)return flow;
int res=0;
for(int &e=cur[u];e;e=nxt[e])
{
int v=to[e];
if(dis[v]==dis[u]+1&&cap[e])
{
int det=dfs(v,min(cap[e],flow-res));
cap[e]-=det,cap[e^1]+=det;
res+=det;if(res==flow)break;
}
}
if(res<flow)dis[u]=-1;
return res;
}
void pop_flow(int u,int v,int &flow)
{
for(int e=first[v],last;e;last=e,e=nxt[e])
if(to[e]==u){e==first[v]?first[v]=nxt[e]:nxt[last]=nxt[e];break;};
int bz=0;
for(int e=first[u],last;e;last=e,e=nxt[e])
{
if(to[e]!=v)continue;
e==first[u]?first[u]=nxt[e]:nxt[last]=nxt[e];
bz=(!cap[e]);break;
}
if(!bz)return;flow--;
for(int e=first[S];e;e=nxt[e])
if(to[e]==u){cap[e]^=1,cap[e^1]^=1;break;}
for(int e=first[T];e;e=nxt[e])
if(to[e]==v){cap[e]^=1,cap[e^1]^=1;break;}
if(bfs())flow+=dfs(S,INF);
}
bool check(double x)
{
memset(first,0,sizeof(first));tot=1;
int res=0,cnt=0;
for(int i=1;i<=n;i++)
{
double dis=sqrt(p[i].x*p[i].x+p[i].y*p[i].y);
if(dis>rad+x||dis<rad-x)return false;
if(dis+rad<=x)for(int j=1;j<=n;j++)add(i,j+n,1);
else
{
double ang=atan2(p[i].x,p[i].y);
double det=acos((rad*rad+dis*dis-x*x)/(2*rad*dis));
double l=ang-det,r=ang+det;
while(l<0)l+=2*Pi;while(r<0)r+=2*Pi;
int L=l/blk,R=r/blk;
q[++cnt]=data(l-L*blk,i,L+1,1),q[++cnt]=data(r-R*blk,i,R+1,0);
L++,R++;
if(l<=r)for(int j=L+1;j<=R;j++)add(i,j+n,1);
else
{
for(int j=L+1;j<=n;j++)add(i,j+n,1);
for(int j=1;j<=R;j++)add(i,j+n,1);
}
}
}
for(int i=1;i<=n;i++)add(S,i,1),add(i+n,T,1);
while(bfs())res+=dfs(S,INF);
if(res>=n)return true;
sort(q+1,q+cnt+1);
for(int i=1;i<=cnt;i++)
{
if(q[i].op)
{
add(q[i].u,q[i].v+n,1);
if(bfs())res+=dfs(S,INF);
if(res>=n)return true;
}
else pop_flow(q[i].u,q[i].v+n,res);
}
return false;
}
int main()
{
//freopen("fleet.in","r",stdin);
//freopen("fleet.out","w",stdout);
n=getint(),rad=getint(),blk=2*Pi/n,S=0,T=2*n+1;
for(int i=1;i<=n;i++)p[i].x=getint(),p[i].y=getint();
double l=0,r=400;
while(r-l>eps)
{
double mid=(l+r)/2;
if(check(mid))r=mid;
else l=mid;
}
printf("%0.8lf\n",l);
return 0;
}