[bsoj1922] 机器蛇

---

题目描述

　　在未来的某次战争中，我军计划了一次军事行动，目的是劫持敌人的航母。由于这个计划高度保密，你只知道你所负责的一部分：机器蛇的通信网络。计划中要将数百条机器蛇投放到航母的各个角落里。由于航母内部舱室、管线错综复杂，且大部分由金属构成，因此屏蔽效应十分强烈，况且还要考虑敌人的大强度电子干扰，如何保持机器蛇间的联系，成了一大难题。每条机器蛇的战斗位置由作战计划部门制定，将会及时通知你。每条机器蛇上都带有接收、发射系统，可以同时与多条机器蛇通讯。由于整个系统承载的数据量庞大，需要一个固定的通讯网络。情报部门提供了极其详尽的敌方航母图纸，使你对什么地方有屏蔽了如指掌。
　　请你设计一个程序，根据以上信息构造通讯网络，要求信息可以在任意两条机器蛇间传递，同时为了避免干扰，通讯网络的总长度要尽可能的短。

---

输入格式

　　第一行是一个整数n（n≤200）表示参战的机器蛇总数。
　　以下n行，每行两个整数xi，yi，为第i支机器蛇的战斗位置。
　　接下来一行是一个整数m（m≤100）表示航母内部可能产生屏蔽的位置。
　　最后m行，每行四个整数ai，bi，ci，di，表示线段(ai，bi)-(ci，di)处可能有屏蔽，也就是说通讯网络不能跨越这条线段。

---

输出格式

仅一个实数，表示建立的通讯网的最短长度，保留3位小数。如果不能成功建立通讯网，请输出-1.000。

---

样例数据

样例输入

3
1 3
3 1
5 5
1
0 0 3 3

样例输出

8.944

---

题目分析

最小生成树
建边时用计算几何判断线段是否相交

---

源代码

#include<algorithm>
#include<iostream>
#include<iomanip>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
inline const int Get_Int() {
    int num=0,bj=1;
    char x=getchar();
    while(x<'0'||x>'9') {
        if(x=='-')bj=-1;
        x=getchar();
    }
    while(x>='0'&&x<='9') {
        num=num*10+x-'0';
        x=getchar();
    }
    return num*bj;
}
struct Vector;
struct Point {
    double x,y;
};
double Dist(Point a,Point b) {
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
struct Vector {
    double x,y;
};
Vector operator - (Point a,Point b) {
    Vector c;
    c.x=b.x-a.x;
    c.y=b.y-a.y;
    return c;
}
double Cross(Vector a,Vector b) { //叉积
    return a.x*b.y-b.x*a.y;
}
bool Segment_Intersect(Point A,Point B,Point C,Point D) { //判断两线段是否相交
    return (Cross((C-A),(D-A))*Cross((C-B),(D-B))<0)&&(Cross((A-C),(B-C))*Cross((A-D),(B-D))<0);
}

const int maxn=50005;
struct Edge {
    int x,y;
    double v;
    bool operator < (Edge b) const {
        return v<b.v;
    }
} a[maxn];
struct Kruskal {
    int n,m,cnt;
    double ans;
    int father[maxn];
    void init(int n,int m) {
        this->n=n;
        this->m=m;
        ans=cnt=0;
    }
    void AddEdge(int from,int to,double dist) {
        a[++cnt].x=from;
        a[cnt].y=to;
        a[cnt].v=dist;
    }
    int Get_Father(int x) {
        if(father[x]==x)return x;
        return father[x]=Get_Father(father[x]);
    }
    bool main() {
        int cnt=0;
        sort(a+1,a+m+1);
        for(int i=1; i<=n; i++)father[i]=i;
        for(int i=1; i<=m; i++) {
            int f1=Get_Father(a[i].x),f2=Get_Father(a[i].y);
            if(f1!=f2) {
                father[f1]=f2;
                ans+=a[i].v;
                cnt++;
                if(cnt==n-1)break;
            }
        }
        if(cnt!=n-1)return false;
        return true;
    }
};

Point P[205],Start[205],End[205];
Kruskal krus;
int n,m,sum=0;
int main() {
    n=Get_Int();
    for(int i=1; i<=n; i++) {
        P[i].x=Get_Int();
        P[i].y=Get_Int();
    }
    m=Get_Int();
    for(int i=1; i<=m; i++) {
        Start[i].x=Get_Int();
        Start[i].y=Get_Int();
        End[i].x=Get_Int();
        End[i].y=Get_Int();
    }
    for(int i=1; i<=n; i++)
        for(int j=i+1; j<=n; j++) {
            int bj=1;
            for(int k=1; k<=m; k++)
                if(Segment_Intersect(P[i],P[j],Start[k],End[k])) {
                    bj=0;
                    break;
                }
            if(bj==0)continue;
            sum++;
            krus.AddEdge(i,j,Dist(P[i],P[j]));
        }
    krus.init(n,sum);
    if(!krus.main())puts("-1.000");
    else printf("%0.3lf\n",krus.ans);
    return 0;
}

---