bzoj5316: [Jsoi2018]绝地反击【计算几何+二分+扫描线+网络流】

传送门

解题思路：

先二分答案，那么每个点就有了一个匹配范围。
考虑轨道上匹配点的偏转角度，肯定是一个匹配点抵住一个点的匹配边界上，所以一共只有2n种匹配方式，每次暴力网络流是$O(n^4logn)$

注意到偏转角度不会超过$2\pi/n$，即最多只会一项每个点一条匹配边的出现，另一条的消失，所以我们可以用扫描线的方法，将操作排序后处理增广和退流操作即可，一次就只会是$O(n)$的复杂度。

由于是二分图，所以退流只会影响三条边，还是很好写的。

时间复杂度$O(n^3logn)$

#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;
}