HDU3362+状态压缩

dp[ i ]表示该状态下得所需花费。

 /*
 状态压缩dp
 dp[i] = min( dp[ i-j ]+cost[ j ] );
 由i-j的状态转到i的状态
 */
 #include<stdio.h>
 #include<string.h>
 #include<stdlib.h>
 #include<iostream>
 #include<math.h>
 using namespace std;
 const int maxn = ;
 const double inf = 99999999.0;
 const int maxm = (<<);
 const double eps = 1e-;
 struct node{
     double x,y;
 }dot[ maxn ];
 double dp[ maxm ];
 double dis[ maxn ][ maxn ];
 double dist( int i,int j ){
     return sqrt( ( dot[i].x-dot[j].x )*( dot[i].x-dot[j].x )+( dot[i].y-dot[j].y )*( dot[i].y-dot[j].y ) );
 }
 int main(){
     int n;
     while( scanf("%d",&n),n ){
         int sum,cnt,flag;
         cnt = ;
         sum = ;
         for( int i=;i<n;i++ ){
             scanf("%lf%lf%d",&dot[i].x,&dot[i].y,&flag);
             if( flag== ){
                 sum += (<<i);
                 cnt++;
             }
         }
         if(( n>&&cnt< )||( n==&&cnt== )){
             printf("No Solution\n");
             continue;
         }
         for( int i=;i<n;i++ )
             for( int j=;j<n;j++ )
                 dis[ i ][ j ] = dist( i,j );
         int N = (<<n);
         for( int i=;i<N;i++ )
             dp[ i ] = inf;
         dp[ sum ] = 0.0;
         for( int i=sum;i<N;i++ ){
             if( dp[i]>inf ) continue;
             for( int j=;j<n;j++ ){
                 if( i&(<<j) )//j是固定的
                     continue;
                 double min1 = inf;
                 double min2 = inf;
                 for( int k=;k<n;k++ ){//找出min1,min2
                     if( i&(<<k) ){
                         if( min1>dis[ j ][ k ] ){
                             min2 = min1;
                             min1 = dis[ j ][ k ];
                         }
                         else if( min2>dis[ j ][ k ] ){
                             min2 = dis[ j ][ k ];
                         }
                     }
                 }
                 dp[ i|(<<j) ] = min( dp[ i|(<<j) ],dp[i]+min1+min2 );
             }
         }
         printf("%.6lf\n",dp[ N- ]);
     }
     return ;
 }