40状压DP求助，谢谢

#include<bits/stdc++.h>
using namespace std;
int n;
double x[19],y[19];
double f[1<<19][19];
double d[19][19];
double ans=100000000.0;
int main(){
	cin>>n;
	for(int i=1;i<=n;i++)
		cin>>x[i]>>y[i];
	for(int i=1;i<(1<<19);i++)
		for(int j=0;j<19;j++)
			f[i][j]=100000000.0;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			d[i][j]=(double)sqrt(pow(x[i]-x[j],2)+pow(y[i]-y[j],2));
	for(int i=0;i<n;i++)
		f[1<<i][i+1]=(double)sqrt(pow(x[i+1],2)+pow(y[i+1],2));
	for(int i=1;i<(1<<n);i++)
		for(int j=1;j<=n;j++){
			if((i&(1<<(j-1)))==0) continue;
			for(int k=1;k<=n;k++){
				if((i&(1<<(k-1)))==0||j==k) continue;
				f[i][j]=min(f[i][j],f[i-(1<<(j-1))][k]+d[k][j]);
			}
		}
	for(int i=1;i<=n;i++)
		ans=min(ans,f[(i<<n)-1][i]);
	printf("%.2lf\n",ans);
  	return 0;
}