hdu3622:http://acm.hdu.edu.cn/showproblem.php?pid=3622
题意：你有n次，每次你可以在平面上放置一个点，并且每一次都会有两个位置可以选，每一次只能选择其中一个。然后在自己位置上以该点为圆心画圆，这n个圆不能相交，问你最后最小的圆的半径的最大值是多少。
题解：二分+2-sat。对于2-sat的建图始终很迷糊。这里的一次有两个点，并且每一次能选择其中一个，这就相当于把n个点弄成了2*n点。如果i--j距离太近 ，则建边i-->~j,j—>~i.
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
using namespace std;
const int N=;
const int M=;
const double eps=1e-;
struct Edge{
int to,next;
} edge[M];
struct Node{
double x;
double y;
}num1[];
int n,m,cnt,dep,top,atype;
int dfn[N],low[N],vis[N],head[N],st[N],belong[N],in[N],out[N],sum[N];
//sum[i]记录第i个连通图的点的个数,in[i],out[i],表示缩点之后点的入度和初度。
void init(){
cnt=dep=top=atype=;
memset(head,-,sizeof(head));
memset(dfn,,sizeof(dfn));
memset(low,,sizeof(low));
memset(vis,,sizeof(vis));
memset(belong,,sizeof(belong));
memset(in,,sizeof(in));
memset(out,,sizeof(out));
memset(sum,,sizeof(sum));
}
void addedge(int u,int v){
edge[cnt].to=v;
edge[cnt].next=head[u];
head[u]=cnt++;
} void Tarjan(int u){
dfn[u]=low[u]=++dep;
st[top++]=u;
vis[u]=;
for(int i=head[u]; i!=-; i=edge[i].next){
int v=edge[i].to;
if(!dfn[v]){
Tarjan(v);
low[u]=min(low[u],low[v]);
}
else if(vis[v]){
low[u]=min(low[u],dfn[v]);
}
}
int j;
if(dfn[u]==low[u]){
atype++;
do{
j=st[--top];
belong[j]=atype;
sum[atype]++; //记录每个连通分量中点的个数
vis[j]=;
}
while(u!=j);
}
}
bool judge(double mid){
init();
//int base=2*n;
for(int i=;i<=*n;i++){
for(int j=i+;j<=*n;j++){
double dis=sqrt((num1[i].x-num1[j].x)*(num1[i].x-num1[j].x)+(num1[i].y-num1[j].y)*(num1[i].y-num1[j].y));
if(dis<*mid){
if(i<=n&&j<=n){
addedge(i,j+n);
addedge(j,i+n);
}
if(i<=n&&j>n){
addedge(i,j-n);
addedge(j,i+n);
}
if(i>n&&j<=n){
addedge(i,j+n);
addedge(j,i-n);
}
if(i>n&&j>n){
addedge(i,j-n);
addedge(j,i-n);
} }
}
}
for(int i=;i<=*n;i++)
if(!dfn[i])Tarjan(i);
for(int i=;i<=n;i++)
if(belong[i]==belong[i+n])
return false; return true; }
int main(){
while(~scanf("%d",&n)){
for(int i=;i<=n;i++){
scanf("%lf %lf %lf %lf",&num1[i].x,&num1[i].y,&num1[i+n].x,&num1[i+n].y);
}
double l=,r=,ans=;
while(abs(l-r)>eps){
double mid=(l+r)/2.0;
if(judge(mid)){
ans=mid;
l=mid;
}
else{
r=mid;
}
}
printf("%.2f\n",ans);
} }