1433吃奶酪

1433吃奶酪

问至少要跑多少距离，状态压缩DP
一般这种状态压缩DP都和搜索有关系
状态压缩的时间复杂度都是确定的，这一点很nice
首先设f[i][s]表示从i出发遍历集合为s的点的路程最小值，枚举s里其他点进行转移
然后就很简单了，状压好难呜呜呜

#include<cstdio>
#include<cmath>
#include<cstring>
typedef double db;
db x[20],y[20],f[20][35000];
template<class T> T min(T a,T b)
{
 	return a<b?a:b;
}
db dis(int a,int b) 
{
	return sqrt((x[a]-x[b])*(x[a]-x[b])+(y[a]-y[b])*(y[a]-y[b]));//曼哈顿计算 
}
int main()
{
    int n;
	scanf("%d",&n);
    for(int i=1;i<=n;i++) 
		scanf("%lf%lf",&x[i],&y[i]);
    memset(f,127,sizeof(f));
    for(int s=1;s<=(1<<n)-1;s++)
	    for(int i=1;i<=n;i++)
	    {
	        if((s&(1<<(i-1)))==0) continue;//跳到下一个集合
	        if(s==(1<<(i-1))) 
			{
				f[i][s]=0;
				continue;
			}
	        for(int j=1;j<=n;j++)
	        {
	            if((s&(1<<(j-1)))==0||i==j) continue;
	            f[i][s]=min(f[i][s],f[j][s-(1<<(i-1))]+dis(i,j));
	        }
	    }
    db ans=-1;
    for(int i=1;i<=n;i++)
    {
        db s=f[i][(1<<n)-1]+dis(i,0);
        if(ans==-1||ans>s) ans=s;
    }
    printf("%.2lf\n",ans);
    return 0;
}