点集配对问题

最新推荐文章于 2024-04-16 09:50:49 发布

zzti_xiaowei

最新推荐文章于 2024-04-16 09:50:49 发布

阅读量2.2k

点赞数 1

分类专栏： Acm--动态规划状压dp

本文链接：https://blog.csdn.net/zzti_xiaowei/article/details/82700450

版权

Acm--动态规划同时被 2 个专栏收录

28 篇文章 0 订阅

订阅专栏

状压dp

5 篇文章 1 订阅

订阅专栏

<摘自刘汝佳训练指南P61>

题意： 空间里n个点（n为偶数），把它们配成n/2对点，使得每个点恰好在一个点对中，要求所有点对中，两点距离之和最小值，其中n<=20。

思路：
dp(s)：表示集合s配对后的最小距离和。
状态转移方程： dp(s) = min(dp{S - {i} - {j}+ |Pi Pj| | j∈S, j > i, i = min{S}}

很巧妙的一点就是i取的是min{s}！！！如果枚举所有可能，会有很多重复取法，复杂度会*n！！！

假如有12345678共八个坐标，其中的某一取法为[1,2]一组，[3,4]一组,[5,6]一组,[7,8]一组。
那么其它取法有：
(1,2)->(3,4)->(5,6)->(7,8)
(5,6)->(7,8)->(1,2)->(3,4)
(1,2)->(3,4)->(7,8)->(5,6)
(3,4)->(5,6)->(7,8)->(1,2)
(7,8)->(1,2)->(3,4)->(5,6)
etc…

它们的结果是一样的，只是取法顺序不同，所以没有必要枚举所有结果！！！

代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
const int inf=0x3f3f3f3f;
const int Max_n=20;

int n;
struct node{
    int x,y;
}no[Max_n];
double dp[1<<Max_n];

double sq(int i,int j){
    return sqrt((no[i].x-no[j].x)*(no[i].x-no[j].x)+(no[i].y-no[j].y)*(no[i].y-no[j].y));
}

int main()
{
    scanf("%d",&n);
    for(int i=0;i<n;i++)
        scanf("%d%d",&no[i].x,&no[i].y);
    for(int i=0;i<(1<<n);i++)
        dp[i]=inf;
    dp[0]=0;
    for(int i=1;i<(1<<n);i++){
        int k=0;
        while(!(i&1<<k))k++;
        for(int j=k+1;j<n;j++){
            if(i&1<<j)
//              dp[i]=min(dp[i],dp[i^(1<<k)^(1<<j)]+sq(k,j));
                dp[i]=min(dp[i],dp[i-(1<<k)-(1<<j)]+sq(k,j));
        }
    }
    printf("%.3f\n",dp[(1<<n)-1]);
    return 0;
}