bzoj1601 [Usaco2008 Oct]灌水

http://www.elijahqi.win/archives/1006
Description

Farmer John已经决定把水灌到他的n(1<=n<=300)块农田，农田被数字1到n标记。把一块土地进行灌水有两种方法，从其他农田饮水，或者这块土地建造水库。 建造一个水库需要花费wi(1<=wi<=100000),连接两块土地需要花费Pij(1<=pij<=100000,pij=pji,pii=0). 计算Farmer John所需的最少代价。

Input

*第一行：一个数n

*第二行到第n+1行：第i+1行含有一个数wi

*第n+2行到第2n+1行：第n+1+i行有n个被空格分开的数，第j个数代表pij。

Output

*第一行：一个单独的数代表最小代价.

Sample Input

4
5
4
4
3
0 2 2 2
2 0 3 3
2 3 0 4
2 3 4 0
Sample Output

9

输出详解：

Farmer John在第四块土地上建立水库，然后把其他的都连向那一个，这样就要花费3+2+2+2=9

水库建立在哪里我们不好处理 不妨我们建立一个超级大水库 然后和每个点都有w[i]的权值做边 这样生成的树我们可以认为生成树的时候就可以方便的统计是否有水库的情况了
再次熟练一下prim （蒟蒻的微笑）

#include<cstdio>
#include<queue>
#include<vector>
#include<cstring>
#define N 330
using namespace std;
#define pa pair<int,int>
inline int read(){
    int x=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x; 
}
priority_queue<pa,vector<pa>,greater<pa> >q;
int cost[N][N],f[N],ans,n;bool visit[N];
void prim(){
    memset(f,0x7f,sizeof(f));f[0]=0;q.push(make_pair(0,0));
    while (!q.empty()){
        int x=q.top().second;q.pop();
        if (visit[x]==true) continue;visit[x]=true;ans+=f[x];
        for (int y=1;y<=n;++y){
            if (visit[y]) continue;
            if (cost[x][y]&&cost[x][y]<f[y]) {
                f[y]=cost[x][y];q.push(make_pair(f[y],y));
            }
        }
    }
}
int main(){
    freopen("bzoj1601.in","r",stdin);
    n=read();for (int i=1;i<=n;++i) cost[0][i]=read();
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n;++j) cost[i][j]=read();
    prim();
    printf("%d",ans);
    return 0;
}