Uva10911 最优匹配问题

题目链接

集合上的动态规划，设d(i,S)表示把前i个点中，位于集合S中的元素两两配对的最小距离和,状态转移方程：d(i,S)=min{|P_iP_j|+d(i-1,S-{i}-{j}) | j∈S}, |P_iP_j|表示P_i和P_j间的距离，而上面的方程可以进一步简化，阶段i不必要保存，因为i已经包含在集合S中了，而且S表示的是前i个点的集合，所以i是集合中最大的元素，则d(S)表示“把S中的元素两两配对的的最小距离和”转移方程为 d(S)=min{|P_iP_j|+d(S-{i}-{j})|j∈S,i=max(S)}

#include <iostream>
#include <cmath>
#include <vector>
#include <algorithm>
using namespace std;

const int maxn = 25;
const int INF = 1<<30;
int f[maxn][2];
double d[100000];
char str[maxn];
int n;

double dist(int i, int j)
{
    return sqrt((f[i][0] - f[j][0])*(f[i][0] - f[j][0]) +
        (f[i][1] - f[j][1])*(f[i][1] - f[j][1]));
}

int main()
{
    int Case = 0;
    while (scanf("%d", &n) && n != 0)
    {
        getchar();
        for (int i = 0; i < 2 *n ; i++)
        {
            scanf("%s %d %d", str, &f[i][0], &f[i][1]);
        }

        n = 2 * n;
        d[0] = 0;
        int i, j, S;

        for (S = 1; S < (1 << n);S++)//子集枚举用二进制表示
        {           
            d[S] = INF;
            for (i = 0; i < n; i++)
            {
                if (S & (1 << i))
                    break;
            }
            for (j = i + 1; j < n; j++)
            {
                if (S & (1 << j))
                    d[S] = min(d[S], dist(i, j) + d[S ^ (1 << i) ^ (1 << j)]);//d[S ^ (1 << i) ^ (1 << j)]表示集合S去掉Pi和Pj后剩下的集合
            }
        }
        printf("Case %d: %.2f\n", ++Case, d[(1 << n) - 1]);
    }
    return 0;
}