185. [USACO Oct08] 挖水井

185. [USACO Oct08] 挖水井(点击转到COGS)

输入文件：water.in   输出文件：water.out  
时间限制：1 s   内存限制：128 MB

描述

农夫约翰决定给他的N（1<=N<=300）个牧场浇水，这些牧场被自然的命名为1..N。他可以给一个牧场引入水通过在这个牧场挖一口井或者修一条管道使这个牧场和一个已经有水的牧场连接。

在牧场i挖一口井的花费是w_i（1<=w_i<=100000）。修建一条水管连接牧场i和牧场j的花费是p_ij（1<=p_ij<=100000;p_ij=p_ji;p_ii=0）。

请确定农夫约翰为了完成浇灌所有的牧场所需的最小的总花费。

题目名称：water

输入格式：

第1行：一个单独的整数n。

第2..n+1行：第i+1行包含一个单独的整数w_i。

第n+2..2n+1行：第n+1+i行包含n个用空可分开的整数；其中第j个数是p_ij。

输入样例（file water.in）：

4

5

4

4

3

0 2 2 2

2 0 3 3

2 3 0 4

2 3 4 0

输入说明：

这里有4个牧场，修井和修管道的代价如图。

输出格式：

•     第1行：一个单独的整数，表示花费。

输出样例（file water.out）：

9

输出说明：

农夫约翰可以在第4个牧场修井，并且将每个牧场和第一个连接起来，这样，花费是3+2+2+2=9。

思路：

将挖水井看成边(连向自己的)，利用并查集，求出各最小生成树

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #define maxn 123456789
 using namespace std;
 struct edge
 {
     int fr,to,w;
 }E[];//size > 300*300
 int n,x,tot,num;
 int w[],father[];
 int find(int x)
 {
     if(father[x]!=x) father[x]=find(father[x]);
     return father[x];
 }
 void unionn(int x,int y)
 {
     father[x]=y;
 }
 void add(int x,int y,int z)
 {
     num++;
     E[num].fr=x;
     E[num].to=y;
     E[num].w=z;
 }
 bool cmp(const edge &a,const edge &b)
 {
     return a.w<b.w;
 }
 int main()
 {
     freopen("water.in","r",stdin);
     freopen("water.out","w",stdout);
     scanf("%d",&n);
     for(int i=;i<=n;i++)
       father[i]=i;
     for(int i=;i<=n;i++)
       scanf("%d",&w[i]);
     for(int i=;i<=n;i++)
       for(int j=;j<=n;j++)
       {
           scanf("%d",&x);
         if(x)
           add(i,j,x);
       }
     for(int i=;i<=n;i++)
     {
         add(i,n+,w[i]);//E2:从i向n+1(任一点)挖井的费用(把边存进去)
         add(n+,i,w[i]);//E3:从n+1向i挖井的费用 E2的反图
     }
     sort(E+,E+num+,cmp);
     int k=;
     for(int i=;i<=num;i++)
     {
         int r1=find(E[i].fr);
         int r2=find(E[i].to);
         if(r1!=r2)
         {
             unionn(r1,r2);
             k++;
             tot+=E[i].w;
         }
         if(k==n)//最小生成树的n-1条边+至少一个水井
           break;
     }
     printf("%d",tot);
     fclose(stdin);fclose(stdout);
     return ;
 }