天天写算法之（2-sat+二分）Bomb Game

地址： 点击打开链接

这个题目，如果不是同一个组（设当前点为i，测到的组为j），那么就判断是不是<dis（二分），如果小于，那么就连接i与j^1，这个时候，如果i与j与j^1都行不通，那么就无法实现每组选一个，也就是说当前的dis无法实现（在代码上表现就是构成了强连通分量，想一想i 与  j^1 ,连接了，同时j 与 i^1连接了，而走到j+1的时候，发现i又与j+1^1构边了，同时j+1也和i^1构边了，因此就形成了一个强连通分量，说明不可行），需要再进行二分。

代码:

#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
const int MAXN=210;
const int MAXM=41000;
const double eps=1e-5;
using namespace std;
struct node{
    int x,y;
}p[MAXN];
struct Edge{
    int to,next;
}edge[MAXM<<1];int head[MAXN],tot;
void init(){
    tot=0;
    memset(head,0xff,sizeof(head));
}
void addedge(int u,int v){
    edge[tot].to=v;
    edge[tot].next=head[u];
    head[u]=tot++;
}

int low[MAXN],dfn[MAXN],sccno[MAXN],scc_cnt,dfs_clock;
int Stack[MAXN],top;
bool instack[MAXN];
void dfs(int u){
    int v;
    Stack[top++]=u;
    instack[u]=1;
    low[u]=dfn[u]=++dfs_clock;
    for(int i=head[u];i!=-1;i=edge[i].next){
        v=edge[i].to;
        if(!dfn[v]){
            dfs(v);
            low[u]=min(low[u],low[v]);
        }
        else if(instack[v]&&low[u]>dfn[v]){
            low[u]=dfn[v];
        }
    }
    if(low[u]==dfn[u]){
        ++scc_cnt;
        for(;;){
            v=Stack[--top];
            instack[v]=0;
            sccno[v]=scc_cnt;
            if(v==u) break;
        }
    }
}
bool solveable(int n){
    memset(dfn,0,sizeof(dfn));
    memset(instack,false,sizeof(instack));
    dfs_clock=scc_cnt=top=0;
    for(int i=0;i<n;++i){
        if(!dfn[i]) dfs(i);
    }
    for(int i=0;i<n;i+=2){
        if(sccno[i]==sccno[i^1])
            return false;
    }
    return true;
}
double dist(node a,node b){
    return sqrt((double)(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
void Build_graph(int n,double m){
    init();
    for(int i=0;i<2*n;++i){
        for(int j=0;j<2*n;++j){
            if(i/2!=j/2&&dist(p[i],p[j])<m){
                addedge(i,j^1);
                addedge(j,i^1);
            }
        }
    }
}
int main(){
#ifndef ONLINE_JUDGE
freopen("in.txt","r",stdin);
#endif // ONLINE_JUDGE
    int n;
    while(scanf("%d",&n)!=EOF){
        for(int i=0;i<n;++i){
            scanf("%d%d%d%d",&p[2*i].x,&p[2*i].y,&p[2*i+1].x,&p[2*i+1].y);
        }
        double  l=0.0,r=40000.0;
        while(r-l>=eps){
            double m=(l+r)/2.0;
            Build_graph(n,2*m);
            if(solveable(2*n)) l=m;
            else r=m;
        }
        printf("%.2lf\n",l);
    }
    return 0;
}