HDU 3920 Clear All of Them I 状态压缩DP 2011 Multi-University Training Contest 9 - Host by BJTU

/*
参考：http://blog.csdn.net/racebug2010/article/details/6683963
第一道状态压缩DP
state的二进制表示（如110011）表示状态第1，2，5，6的点已经被消灭
DP[state]表示为state状态时的最小距离和
DP[state]=min（DP[state],DP[state^(1<<i)^(1<<j)]+dis(p[i],p[j])+dis(my,p[i])）
其中i表示state的第i位(这一位必须为1)
所以只要枚举两个state为1的位i,j进行状态转移（纯粹这样会超时）
由于110011可以由状态000011，010001，010010，100001，100010，110000这些状态而来
但000011和110011转换成110011的实质是一样的（都是1，2配对 5，6配对）
状态可以由六个减少到000011，010001，010010 三个
所以不需要枚举i,j,只需固定i，枚举j就行了
*/
#include<iostream>
#include <stdlib.h>
#include <algorithm>
#include <stdio.h>
#include<vector>
#include<cmath>;
using namespace std;
typedef __int64 i64;
const int M=1000000007;
struct Point
{
    double x,y;
}my,p[22];
int n;
double dp[1<<21];
double dis(Point a,Point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
bool cmp(Point a,Point b)
{
    return dis(my,a)<dis(my,b);
}
//优化前***********************************************************************************
// double DP(int state)
// {
    // if(dp[state]>=0)return dp[state];
    // dp[state]=1<<30;
    // for(int i=0;i<n;i++)
    // {
        // if(state&(1<<i))
        // {
            // for(int j=i+1;j<n;j++)
            // if(state&(1<<j))
            // dp[state]=min(dp[state],DP(state^(1<<i)^(1<<j))+dis(p[i],p[j])+dis(my,p[i]));
        // }
    // }
    // if(dp[state]==(1<<30))
    // dp[state]=0;
    // return dp[state];
// }

//优化后***********************************************************************************
double DP(int state)
{
    if(dp[state]>=0)return dp[state];
    dp[state]=1<<30;
    int pos=-1;
    for(int i=0;i<n;i++)
    {
        if(state&(1<<i))
        {
            if(pos==-1)pos=i;
            else 
            {
                dp[state]=min(dp[state],DP(state^(1<<i)^(1<<pos))+dis(my,p[pos])+dis(p[pos],p[i]));
            }
        }
    }
    if(dp[state]==(1<<30))
    dp[state]=0;
    return dp[state];
}
int main()
{
    int Case;
    scanf("%d",&Case);
    for(int k=1;k<=Case;k++)
    {
        scanf("%lf%lf",&my.x,&my.y);
        scanf("%d",&n);
        n=2*n;
        for(int i=0;i<n;i++)
        scanf("%lf%lf",&p[i].x,&p[i].y);
        sort(p,p+n,cmp);
        for(int i=0;i<(1<<n);i++)
        dp[i]=-1;
        printf("Case #%d: %.2lf\n",k,DP((1<<n)-1));
    }
    return 0;
}