P1268 树的重量【构造】

题目描述

树可以用来表示物种之间的进化关系。一棵“进化树”是一个带边权的树，其叶节点表示一个物种，两个叶节点之间的距离表示两个物种的差异。现在，一个重要的问题是，根据物种之间的距离，重构相应的“进化树”。

令N={1..n}，用一个N上的矩阵M来定义树T。其中，矩阵M满足：对于任意的i，j，k，有M[i,j] + M[j,k] >= M[i,k]。树T满足：

1．叶节点属于集合N；

2．边权均为非负整数；

3．dT(i,j)=M[i,j]，其中dT(i,j)表示树上i到j的最短路径长度。

如下图，矩阵M描述了一棵树。

树的重量是指树上所有边权之和。对于任意给出的合法矩阵M，它所能表示树的重量是惟一确定的，不可能找到两棵不同重量的树，它们都符合矩阵M。你的任务就是，根据给出的矩阵M，计算M所表示树的重量。下图是上面给出的矩阵M所能表示的一棵树，这棵树的总重量为15。

输入格式

输入数据包含若干组数据。每组数据的第一行是一个整数n(2<n<30)。其后n-1行，给出的是矩阵M的一个上三角(不包含对角线)，矩阵中所有元素是不超过100的非负整数。输入数据保证合法。

输入数据以n=0结尾。

输出格式

对于每组输入，输出一行，一个整数，表示树的重量。

输入输出样例

输入 #1

5
5 9 12 8
8 11 7
5 1
4
4
15 36 60
31 55
36
0

输出 #1

15
71

和Three Paths on a tree类似地构造一个树上最短路的式子
把整张图看作一张AOE网

通过这张网来建立一棵树，具体实现操作是：
1、把一个为入树的点入树
2、在已经入树的它的邻接点中遍历来更新把这个点入树的最小花费。

显然此时花费是 (dis[i][k] + dis[j][k] - dis[i][j]) / 2;

3、树的总重量加上把这个点入树的最小花费

CODE

#include <bits/stdc++.h>
#define dbg(x) cout << #x << "=" << x << endl
#define eps 1e-8
#define pi acos(-1.0)

using namespace std;
typedef long long LL;

template<class T>inline void read(T &res)
{
    char c;T flag=1;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;res=c-'0';
    while((c=getchar())>='0'&&c<='9')res=res*10+c-'0';res*=flag;
}

namespace _buff {
    const size_t BUFF = 1 << 19;
    char ibuf[BUFF], *ib = ibuf, *ie = ibuf;
    char getc() {
        if (ib == ie) {
            ib = ibuf;
            ie = ibuf + fread(ibuf, 1, BUFF, stdin);
        }
        return ib == ie ? -1 : *ib++;
    }
}

int qread() {
    using namespace _buff;
    int ret = 0;
    bool pos = true;
    char c = getc();
    for (; (c < '0' || c > '9') && c != '-'; c = getc()) {
        assert(~c);
    }
    if (c == '-') {
        pos = false;
        c = getc();
    }
    for (; c >= '0' && c <= '9'; c = getc()) {
        ret = (ret << 3) + (ret << 1) + (c ^ 48);
    }
    return pos ? ret : -ret;
}

const int maxn = 107;

int dis[maxn][maxn];
int n;

int main()
{
    while(scanf("%d",&n) && n) {
        memset(dis, 0, sizeof(dis));
        for ( int i = 1; i <= n; ++i ) {
            for ( int j = i+1; j <= n; ++j ) {
                read(dis[i][j]);
                dis[j][i] = dis[i][j];
            }
        }
        int weight = 0;
        int v = 1;
        for ( int i = 1; i <= n; ++i ) {
            int cost = 0x3f3f3f3f;
            for ( int j = 1; j < i; ++j ) {
                int temp = dis[i][v] + dis[i][j] - dis[v][j];
                cost = min(cost, temp / 2);
                //dbg(cost);
            }
            if(cost != 0x3f3f3f3f)
                weight += cost;
        }
        cout << weight << endl;
    }
    return 0;
}