HDU3622(二分+2-SAT)

来源：互联网发布：淘宝助理外链图片编辑：程序博客网时间：2024/05/17 03:40

题意不说了，直接讲思路。

首先对半径进行二分，然后再判断炸弹之间的距离是否小于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 26using 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;}

0 0