二分+2-sat——hdu3062

hdu3622升级版

注意要保留两位小数

/*
给定n对圆心(x,y),要求从每对里找到一个点画圆，不可相交
使得最小半径最大
二分答案，设最小半径为r
然后两两配对一次进行判断，在2-sat上连边即可
*/
#include<bits/stdc++.h>
using namespace std;
#define esp 1e-6
#define N 405
#define M 405*405
struct Node{double x,y;}p[N][];
struct Edge{int to,nxt;}e[M<<];
int n,m,head[N],tot;
void add(int u,int v){
    e[tot].to=v;e[tot].nxt=head[u];head[u]=tot++;
}

int dfn[N],low[N],cnt,id[N],ind,stk[N],top,ins[N];
void tarjan(int x){
    low[x]=dfn[x]=++ind;
    stk[++top]=x;ins[x]=;
    for(int i=head[x];i!=-;i=e[i].nxt){
        int v=e[i].to;
        if(!dfn[v]){
            tarjan(v);
            low[x]=min(low[x],low[v]);
        }
        else if(ins[v])
            low[x]=min(low[x],low[v]);
    }
    if(low[x]==dfn[x]){
        int y;cnt++;
        do{
            y=stk[top--];
            ins[y]=;
            id[y]=cnt;
        }while(x!=y);
    }
}

void init(){
    cnt=tot=ind=top=;
    memset(head,-,sizeof head);
    memset(low,,sizeof low);
    memset(dfn,,sizeof dfn);
    memset(ins,,sizeof ins);
}

int check(Node a,Node b,double r){
    double dis=sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
    if(dis>=*r)return ;
    return ;
}
int judge(double r){//判断半径r是否可行
    init();
    for(int i=;i<n-;i++)
        for(int j=i+;j<n;j++){
            if(!check(p[i][],p[j][],r))
                add(i*,j*+),add(j*,i*+);
            if(!check(p[i][],p[j][],r))
                add(i*,j*),add(j*+,i*+);
            if(!check(p[i][],p[j][],r))
                add(i*+,j*+),add(j*,i*);
            if(!check(p[i][],p[j][],r))
                add(i*+,j*),add(j*+,i*);
        }
    for(int i=;i<*n;i++)
        if(!dfn[i])tarjan(i);
    for(int i=;i<n;i++)
        if(id[i*]==id[i*+])return ;
    return ;
}

int main(){
    while(cin>>n){
        for(int i=;i<n;i++){
            cin>>p[i][].x>>p[i][].y;
            cin>>p[i][].x>>p[i][].y;
        }
        double l=,r=,mid,ans=;
        while(l+esp<=r){
            mid=(l+r)/;
            if(judge(mid))
                ans=mid,l=mid;
            else r=mid;
        }
        printf("%.2lf\n",ans);
    }
}