给定n个回合,每个回合给你两个点,每个回合只能选择一个点放置炸弹,在n个回合里选出n个点放置炸弹,炸弹的爆炸范围是一个圆形范围,半径可以控制。求满足每个炸弹爆炸的圆形区域不相交的条件下,的所有半径里面的最小半径最大值。
二分枚举两点之间的距离,如果存在两点距离小于等于枚举的距离时,则表示这两个点之间存在矛盾,根据2-sat建图。通过2-sat判断是否满足条件,存在矛盾说明枚举的距离太小,不存在说明没机的距离太大
注意:先按平方算距离 ,最后再开方,会卡精度。
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define clr(x)memset(x,0,sizeof(x))
#define eps 1e-6
int abs(int x){if(x>0)return x;return -x;}
const int maxn=50010;
struct node{
int to;
int next;
}q[50010];
int head[maxn];
int tot;
void add(int s,int u){
q[tot].to=u;
q[tot].next=head[s];
head[s]=tot++;
}
bool ins[maxn];
int color[maxn];
int dfn[maxn],low[maxn],stack[maxn];
int ti,sn,top;
void tarjan(int u){
dfn[u]=low[u]=++ti;
stack[++top]=u;
ins[u]=true;
int i,k;
for(i=head[u];i;i=q[i].next){
k=q[i].to;
if(dfn[k]==0){
tarjan(k);
if(low[k]<low[u])low[u]=low[k];
}
else if(ins[k]&&dfn[k]<low[u])
low[u]=dfn[k];
}
if(dfn[u]==low[u]){
sn++;
do{
k=stack[top--];
ins[k]=false;
color[k]=sn;
}while(k!=u);
}
}
void init(){
clr(low);
clr(dfn);
clr(ins);
clr(head);
tot=1;
top=-1;
sn=ti=0;
}
int n;
struct Point{
int x,y;
}p[500];
double g[500][500];
double dist(Point a,Point b){
return (double)((a.x-b.x)*(a.x-b.x))+(double)((a.y-b.y)*(a.y-b.y));
}
int ok(int mid){
int i,j;
init();
for(i=0;i<n;i++){
for(j=i+1;j<n;j++){
if(g[i][j]<mid+eps){
add(i,j+n);add(j,i+n);
}
if(g[i][j+n]<mid+eps){
add(i,j);add(j+n,i+n);
}
if(g[i+n][j]<mid+eps){
add(i+n,j+n);add(j,i);
}
if(g[i+n][j+n]<mid+eps){
add(i+n,j);add(j+n,i);
}
}
}
for(i=0;i<2*n;i++)
if(dfn[i]==0)tarjan(i);
int flag=0;
for(i=0;i<n;i++)
if(color[i]==color[i+n])flag=1;
if(flag)return 0;
return 1;
}
int main(){
int i,j;
while(scanf("%d",&n)!=EOF){
for(i=0;i<n;i++)
scanf("%d%d%d%d",&p[i].x,&p[i].y,&p[i+n].x,&p[i+n].y);
double maxlen=-1;
for(i=0;i<2*n;i++){
for(j=i;j<2*n;j++){
if(i==j){g[i][j]=0.0;continue;}
g[i][j]=g[j][i]=dist(p[i],p[j]);
if(g[i][j]+eps>maxlen)maxlen=g[i][j];
}
}
double mid,res,l=0.0,r=maxlen;
while(r-l>eps){
mid=(l+r)/2.0;
if(ok(mid)){l=mid;res=mid;}
else r=mid;
}
printf("%.2lf\n",sqrt(res)/2.0);
}
return 0;
}