hdu 4385 Moving Bricks

题目大意：

一个人在搬砖，初始位置在(x0,y0)，目的地也是(x0,y0)，给你n（n<=20）块砖的位置(xi,yi)，规定两点之间时间耗费为欧几里德距离的平方，此人一次可以搬不超过两块砖，问你最少的时间花费并给出方案。

题解：

状态DP。f[i]表示i集合，然后f[i]=min(f[i^(1<<(j-1))]+2*d[j][0],f[i^(1<<(j-1))^(1<<(k-1))]+d[k][0]+d[j][0]+d[k][j]),其中：

((i>>(j-1))&1)==1，((i>>(k-1))&1)==1。

在开一个way[i]记录方案

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int inf=1<<29;
int f[1500000];
int way[1500000];
int d[30][30],x[30],y[30];
int sec,n,num;
void dp(int i)
{
    if(f[i]<inf)return;
    int x=i,j=1,k;
    while(x!=0)
    {
        if(((i>>(j-1))&1)==1)
        {
            dp(i^(1<<(j-1)));
            if(f[i^(1<<(j-1))]+2*d[j][0]<f[i])
            {
            f[i]=f[i^(1<<(j-1))]+2*d[j][0];
            way[i]=i^(1<<(j-1));
            }
            k=j+1;
            int y=x>>1;
            while(y!=0)
            {
            if(((i>>(k-1))&1)==1 && d[j][k]<d[j][0]+d[k][0])//cut
            {
                dp(i^(1<<(j-1))^(1<<(k-1)));
                if(f[i^(1<<(j-1))^(1<<(k-1))]+d[k][0]+d[j][0]+d[k][j]<f[i])
                {
                f[i]=f[i^(1<<(j-1))^(1<<(k-1))]+d[k][0]+d[j][0]+d[k][j];
                way[i]=i^(1<<(j-1))^(1<<(k-1));
                }
            }
            k++;
            y=y>>1;
            }
        }
        x=x>>1;
        j++;
    }
}
void out(int x)
{
    int t=x^way[x],i=1;
    while(t!=0)
    {
        if((t&1)==1)
        {
        num++;
        if(num<n)printf("%d ",i);
        else printf("%d\n",i);
        }
        i++;
        t=t>>1;
    }
    if(way[x]>0)out(way[x]);
}
int main()
{
    scanf("%d",&sec);
    for(int z=1;z<=sec;z++)
    {
        scanf("%d%d",&x[0],&y[0]);
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        scanf("%d%d",&x[i],&y[i]);
        for(int i=0;i<=n;i++)
         for(int j=0;j<=n;j++)
         d[i][j]=(x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]);
        for(int i=0;i<=(1<<n);i++)
        f[i]=inf;
        f[0]=0;
        dp((1<<n)-1);
        printf("Case %d:\n",z);
        printf("%d\n",f[(1<<n)-1]);
        num=0;
        out((1<<n)-1);
    }
    return 0;
}