【傳送門：BZOJ2663】

簡要題意：

　　給出n個人的坐標和n個寶石的坐標，坐標之間的距離為歐幾裏德距離，只有當一個人與一個寶石的距離<R時，這個人才能控制自己的身軀

　　求出滿足恰好k個人能夠控制自己的身軀的最小的R值和最大的R值（指的是最壞情況下，R的最大值，如果R為正無窮，則輸出+INF）

題解：

　　二分答案+二分圖匹配

　　首先對於最小值就直接二分R，然後如果距離<R的話，就建邊，然後跑二分圖匹配，如果匹配的人數剛好是k就記錄

　　然後是最大值，顯然只有當n=k的時候才會有正無窮的答案，特判一下

　　因為是最壞情況，所以我們當作是n-k個人互相匹配，只有一個人與一個寶石的距離>=R時，才建邊，然後跑完二分圖匹配，如果如果匹配的人數剛好是n-k就記錄

　　這樣就能處理最大值的情況了

參考代碼：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define eps 1e-7
using namespace std;
struct node
{
    double x,y;
}A[51],B[51];
struct edge
{
    int x,y,next;
}a[3100];int len,last[110];
 
void ins(int x,int y)
{
    len++;
    a[len].x=x;a[len].y=y;
    a[len].next=last[x];last[x]=len;
}
int match[110];
int chw[110];
int n,k;
bool findmuniu(int x,int i)
{
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(chw[y]!=i)
        {
            chw[y]=i;
            
 
if(match[y]==0||findmuniu(match[y],i)==true)
            {
                match[y]=x;
                return true;
            }
        }
    }
    return false;
}
int check(double x,int opt)
{
    len=0;memset(last,0,sizeof(last));
    if(opt==1)
    {
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                if(x-sqrt((A[i].x-B[j].x)*(A[i].x-B[j].x)+(A[i].y-B[j].y)*(A[i].y-B[j].y))>eps)
                {
                    ins(i,j+n);
                }
            }
        }
    }
    else
    {
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                if(sqrt((A[i].x-B[j].x)*(A[i].x-B[j].x)+(A[i].y-B[j].y)*(A[i].y-B[j].y))-x>=eps)
                {
                    ins(i,j+n);
                }
            }
        }
    }
    int t=0;
    memset(chw,0,sizeof(chw));
    memset(match,0,sizeof(match));
    for(int i=1;i<=n;i++)
    {
        if(findmuniu(i,i)==true) t++;
    }
    return t;
}
int main()
{
    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++) scanf("%lf%lf",&A[i].x,&A[i].y);
    for(int i=1;i<=n;i++) scanf("%lf%lf",&B[i].x,&B[i].y);
    double l=0.0,r=999999999.9,ans;
    while(r-l>=eps)
    {
        double mid=(l+r)/2;
        int t=check(mid,1);
        if(t>=k)
        {
            if(t==k) ans=mid;
            r=mid-eps;
        }
        else l=mid+eps;
    }
    printf("%.2lf ",ans);
    if(n==k)
    {
        printf("+INF\n");
        return 0;
    }
    l=0.0,r=999999999.9;
    while(r-l>=eps)
    {
        double mid=(l+r)/2;
        int t=check(mid,2);
        if(t==n-k) ans=mid;
        if(t<n-k)
        {
            r=mid-eps;
        }
        else l=mid+eps;
    }
    printf("%.2lf\n",ans);
    return 0;
}

BZOJ2663: [Beijing wc2012]靈魂寶石