题意不说了,直接讲思路。
首先对半径进行二分,然后再判断炸弹之间的距离是否小于2*半径,如果是,那么就连接i->j^1和j->i^1,然后用强连通判断可行性。
#include <stdio.h>
#include <string.h>
#include <string>
#include <iostream>
#include <algorithm>
#include <vector>
#include <math.h>
#include <map>
#include <queue>
#include <stack>
#include <set>
#define M 205
#define LL long long
#define Ld __int64
#define eps 0.00001
#define INF 999999999
#define MOD 112233
#define MAX 26
using namespace std;
struct Node
{
double x,y;
};
struct Node node[M];
vector<int> G[M];
//dfn数组表示dfs时到达i点的时间,indx表示时间
int dfn[M],low[M],sccno[M],scc_cnt;
int indx;
int num[M];
stack<int> s;
void Tarjan(int u)
{
indx++;
dfn[u]=low[u]=indx; //为结点u设定次序编号和low初值
s.push(u); //将结点u压入栈中
for(int i=0;i<G[u].size();i++) //枚举每条边
{
int v=G[u][i];
if(!dfn[v]) //若结点v未被访问过
{
Tarjan(v); //继续往下找
low[u]=min(low[u],low[v]);
}
else if(!sccno[v]) //若结点v在栈中
{
low[u]=min(low[u],dfn[v]);
}
}
if(low[u]==dfn[u]) //如果结点u是强连通分量的根
{
scc_cnt++;
for(;;)
{
int x=s.top(); //讲v退栈,为该强连通分量中一个顶点
s.pop();
sccno[x]=scc_cnt;
num[scc_cnt]++; //记录每个强连通分量的结点数
if(x==u)
break;
}
}
}
void find_scc(int n)
{
indx=scc_cnt=0;
memset(num,0,sizeof(num));
memset(sccno,0,sizeof(sccno));
memset(dfn,0,sizeof(dfn));
for(int i=0;i<n*2;i++)
if(!dfn[i])
Tarjan(i);
}
double distant(Node a,Node b)
{
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double maxn(double a,double b)
{
if(a<b)
return b;
return a;
}
int main()
{
int n;
while(~scanf("%d",&n))
{
for(int i=0;i<n*2;i+=2)
{
scanf("%lf%lf%lf%lf",&node[i].x,&node[i].y,&node[i^1].x,&node[i^1].y);
}
double l=0.0,r=40000.0;
double ans=-1;
while(r-l>eps)
{
double mid=(l+r)/2;
for(int i=0;i<n*2;i++)
G[i].clear();
for(int i=0;i<n*2;i+=2)
{
for(int j=i+2;j<n*2;j++)
{
if(distant(node[i],node[j])<2*mid)
{
G[i].push_back(j^1);
G[j].push_back(i^1);
}
}
}
for(int i=1;i<n*2;i+=2)
{
for(int j=i+1;j<n*2;j++)
{
if(distant(node[i],node[j])<2*mid)
{
G[i].push_back(j^1);
G[j].push_back(i^1);
}
}
}
find_scc(n);
int flag=0;
for(int i=0;i<n*2;i+=2)
{
if(sccno[i]==sccno[i+1])
{
flag=1;
break;
}
}
if(flag)
{
r=mid;
}
else
{
l=mid;
ans=mid;
}
}
printf("%.2lf\n",ans);
}
return 0;
}