Save your cats Aizu - 2224

博客讲述了Nicholas的猫咪被邪恶女巫用魔法围栏困住,需用圣水破除围栏拯救猫咪。输入包含魔法桩数量、围栏数量、桩坐标及围栏两端信息,要求输出拯救所有猫咪所需圣水的最小量,还给出了示例输入输出及AC代码。

Save your cats

Aizu - 2224

Problem C:

Nicholas Y. Alford was a cat lover. He had a garden in a village and kept many cats in his garden. The cats were so cute that people in the village also loved them.

One day, an evil witch visited the village. She envied the cats for being loved by everyone. She drove magical piles in his garden and enclosed the cats with magical fences running between the piles. She said “Your cats are shut away in the fences until they become ugly old cats.” like a curse and went away.

Nicholas tried to break the fences with a hummer, but the fences are impregnable against his effort. He went to a church and asked a priest help. The priest looked for how to destroy the magical fences in books and found they could be destroyed by holy water. The Required amount of the holy water to destroy a fence was proportional to the length of the fence. The holy water was, however, fairly expensive. So he decided to buy exactly the minimum amount of the holy water required to save all his cats. How much holy water would be required?

Input

The input has the following format:

N M
x1 y1
.
.
.
xN yN
p1 q1
.
.
.
pM qM

The first line of the input contains two integers N (2 ≤ N ≤ 10000) and M (1 ≤ M). N indicates the number of magical piles and M indicates the number of magical fences. The following N lines describe the coordinates of the piles. Each line contains two integers xi and yi (-10000 ≤ xi, yi ≤ 10000). The following M lines describe the both ends of the fences. Each line contains two integers pj and qj (1 ≤ pj, qjN). It indicates a fence runs between the pj-th pile and the qj-th pile.

You can assume the following:

  • No Piles have the same coordinates.
  • A pile doesn’t lie on the middle of fence.
  • No Fences cross each other.
  • There is at least one cat in each enclosed area.
  • It is impossible to destroy a fence partially.
  • A unit of holy water is required to destroy a unit length of magical fence.

Output

Output a line containing the minimum amount of the holy water required to save all his cats. Your program may output an arbitrary number of digits after the decimal point. However, the absolute error should be 0.001 or less.

Sample Input 1

3 3
0 0
3 0
0 4
1 2
2 3
3 1

Output for the Sample Input 1

3.000

Sample Input 2

4 3
0 0
-100 0
100 0
0 100
1 2
1 3
1 4

Output for the Sample Input 2

0.000

Sample Input 3

6 7
2 0
6 0
8 2
6 3
0 5
1 7
1 2
2 3
3 4
4 1
5 1
5 4
5 6

Output for the Sample Input 3

7.236

Sample Input 4

6 6
0 0
0 1
1 0
30 0
0 40
30 40
1 2
2 3
3 1
4 5
5 6
6 4

Output for the Sample Input 4

31.000

AC代码

//#include<bits/stdc++.h>
#define _CRT_SBCURE_NO_DEPRECATE
#include <set>
#include <map>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <functional>
//#define UP(i,x,y) for(int i=x;i<=y;i++)
//#define DOWN(i,x,y) for(int i=x;i>=y;i--)
#define sdddd(x,y,z,k) scanf("%d%d%d%d", &x, &y, &z, &k)
#define sddd(x,y,z) scanf("%d%d%d", &x, &y, &z)
#define sdd(x,y) scanf("%d%d", &x, &y)
#define sd(x) scanf("%d", &x)
#define mp make_pair
#define pb push_back
#define lson k<<1
#define rson k<<1|1
#define mid (l+r)/2
#define ms(x, y) memset(x, y, sizeof x)
using namespace std;
typedef long long ll;
//typedef unsigned long long ull;
#define MOD 1000000007
const int maxn = 10050;
const int INF = 0x3f3f3f3f;
const ll LINF = 0x3f3f3f3f3f3f3f3f;
int n, m, s, t;
int len;
struct road
{
    int u, v;
    double cost;
};
int arr[maxn][2];
road G[maxn<<2];
int par[maxn];
void init()
{
    for(int i = 0; i < maxn; i++) par[i] = i;
}

int Find(int x)
{
    return par[x]==x? x : par[x] = Find(par[x]);
}
bool same(int x, int y)
{
    int fx = Find(x), fy = Find(y);
    if(fx == fy)return true;
    else        return false;
}
void unite(int x, int y)
{
    int fx = Find(x), fy = Find(y);
    if(fx == fy) return;
    else{
        par[fx] = fy;
        return;
    }
}
bool cmp(road a, road b)
{
    return a.cost>b.cost;
}
double kruskal()
{
    sort(G, G+len, cmp);
    init();
    double res = 0;
    int cnt = 1;
    for(int i = 0; i < len; i++)
    {
        road e = G[i];
        if(!same(e.u, e.v))
        {
            unite(e.u, e.v);
            res+=e.cost;
            cnt++;
        }
    }
    return res;
}
double dist(double x, double y, double u, double v)
{
    return sqrt(abs(x-u)*abs(x-u) + abs(y-v)*abs(y-v));
}
set<pair<int,int> > st;
int main()
{
	//freopen("out.txt", "w", stdout);
    while(~sdd(n,m))
    {
        ms(G,0);
        len = 0;
        int ta, tb;
        for(int i = 1; i <= n; i++)
        {
            sdd(arr[i][0], arr[i][1]);
        }
        double sum = 0;
        for(int i = 1; i <= m; i++)
        {
            sdd(ta, tb);
            if(!st.count(mp(ta, tb)))
            {
                G[len].u = ta;
                G[len].v = tb;
                G[len].cost = dist(arr[ta][0],arr[ta][1],arr[tb][0],arr[tb][1]);
                sum += G[len].cost;
                len++;
                st.insert(mp(ta, tb));
                st.insert(mp(tb, ta));
            }
        }
        double ans = sum - kruskal();
        printf("%.3f\n", ans);
    }
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值