codeforces 8C - Looking for Order 状态DP

int  dist[maxn][maxn] ;
int  dp[1<<maxn]  , father[1<<maxn] ;

pair<int , int> pt[maxn]  , p ;

int  gtd(int i , int j){
return   (pt[i].first - pt[j].first) * (pt[i].first - pt[j].first)
+ (pt[i].second - pt[j].second) * (pt[i].second - pt[j].second) ;
}

int  main(){

int  n  ,  x , y ;

while(cin>>p.first>>p.second){
scanf("%d" , &n) ;
for(int i = 0 ; i < n ; i++)  scanf("%d%d" , &pt[i].first , &pt[i].second) ;
pt
= p ;

for(int i = 0 ; i <= n ; i++)
for(int j = i+1 ; j <= n ; j++) dist[i][j] = dist[j][i] = gtd(i , j) ;

memset(dp , -1 , sizeof(dp)) ;
dp[0] = 0  ;

for(int i = 0 ; i < (1<<n) ; i++){
if(dp[i] == -1) continue ;

for(int j = 0 ; j < n ; j++){
if((i & (1<<j)) == 0){
int t = i ^ (1<<j) ;
int dis = dp[i] + 2 * dist[j]
;
if(dp[t] == -1 || dp[t] > dis){
dp[t] = dis ;
father[t] = i ;
}
for(int k = 0 ; k < n ; k++){
if((t & (1<<k)) == 0){
int t2 = t ^ (1<<k) ;
int dis = dp[i] + dist[j]
+ dist[j][k] + dist[k]
;
if(dp[t2] == -1 || dp[t2] > dis){
dp[t2] = dis ;
father[t2] = i ;
}
}
}
break ;
}
}
}

cout<< dp[(1<<n) - 1]<<endl ;
int u = (1<<n) - 1 ;
vector<int> ans ; ans.clear() ;
while(u){
int f = father[u] ;
int c = f ^ u ;
ans.push_back(0) ;
vector<int> v ; v.clear() ;
for(int i = 0 ; i < n ; i++){
if(c & (1 << i)) v.push_back(i+1) ;
}
for(int i = v.size()-1 ; i >= 0 ; i--) ans.push_back(v[i]) ;
u = f ;
}
ans.push_back(0)  ;

printf("%d" , ans[ans.size() - 1]) ;
for(int i = ans.size() - 2 ; i >= 0 ; i--) printf(" %d" , ans[i]) ;
puts("") ;
}
return 0 ;
}