洛谷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-l行，给出的是矩阵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

像这样的构造题目一般的技巧就是把一些点的信息合并然后减小规模，

最后得到答案。

考虑如何合并两个叶子，对于两个有相同父亲的叶子A，B，

设他们的父亲为C，设A到C的距离为x，B到C的距离为y。

那么对于所有其他的叶子D，我们有M[D,A] - M[D,B] = x - y,

且x + y = M[A,B]。

我们枚举这样的A，B,检查他们是否满足条件，如果满足，

那么我们可以解出x，y，然后把A,B合并起来。

合并A，B就是把A到C的距离调整为M[A,B]，

然后修正各个点到A的距离，即M[D,A] = M[D,A] + y。

最后合并到只剩下两个点时他们之间的距离即为答案。

附代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#define MAXN 35
#define MAX 999999999
using namespace std;
int n,ans,a[MAXN][MAXN];
inline int read(){
	int date=0,w=1;char c=0;
	while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();}
	while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();}
	return date*w;
}
int main(){
    while(1){
        n=read();
        if(n==0)break;
        memset(a,0,sizeof(a));
        for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++)
        a[j][i]=a[i][j]=read();
        ans=a[1][2];
        for(int i=3;i<=n;i++){
            int minn=MAX;
            for(int j=2;j<=i-1;j++)minn=min(minn,(a[1][i]+a[i][j]-a[1][j])/2);
            ans+=minn;
        }
        printf("%d\n",ans);
    }
	return 0;
}