HDU 4385 Moving Bricks [状态压缩DP]

　　有N个点，每次可以从起点出发到1或者2个点然后回到起点，问走过所有的点需要的最小路程（路程定义为欧氏距离的平方）。

　　显然，到一个点路径为S->i->S,到两个点走过的路径为S->i->j->S。状态压缩然后枚举即可。因为和顺序无关，每次只要枚举最小的没有走的点，然后枚举到这个点返回或者到这个点再去另一个点之后返回的路径长度，相同的状态取最小值即可。

　　

#include <string.h>
#include <stdio.h>
#include <algorithm>
#define INF 0x3f3f3f3f
#define MAXS 1100000
using namespace std;
struct pnt{int x, y;}p[25];
int n, cas, ans[25], anss;
int d[MAXS], pre[MAXS], dis[25][25];

int lowbit(int x){return x&-x;}
bool cmp(const int& x, const int& y) {return lowbit(x) < lowbit(y);}
int sqr(int x){return x * x;}
int getdis(pnt p1, pnt p2){return sqr(p1.x - p2.x) + sqr(p1.y - p2.y);}
int main(){
    //freopen("test.in", "r", stdin);
    scanf("%d", &cas);
    for (int ca = 1; ca <= cas; ca++) {
        scanf("%d%d", &p[0].x, &p[0].y);
        scanf("%d", &n);
        for (int i = 1; i <= n; i++) scanf("%d%d", &p[i].x, &p[i].y);
        for (int i = 0; i <= n; i++) for (int j = 0; j <= n; j++) dis[i][j] = getdis(p[i], p[j]);
        int full = 1<<n;
        for (int i = 1; i < full; i++) d[i] = INF;
        for (int s = 0; s < full; s++) {
            for (int i = 0 ;i < n; i++) if(0 == (s&1<<i)){
                if(d[s|1<<i] > d[s] + 2 * dis[0][i+1]) {
                    d[s|1<<i] = d[s] + 2 * dis[0][i+1];
                    pre[s|1<<i] = s;
                }
                for (int j = i+1, ss = 1<<i|s; j < n; j++) if(0 == (ss&1<<j)) {
                    if (d[ss|1<<j] > d[s] + dis[0][i+1] + dis[i+1][j+1] + dis[j+1][0]) {
                        d[ss|1<<j] = d[s] + dis[0][i+1] + dis[i+1][j+1] + dis[j+1][0];
                        pre[ss|1<<j] = s;
                    }
                }
                break;
            }
        }
        anss = 0;
        for (int x = full -1; x != 0; x = pre[x]) {
            ans[anss++] = x^pre[x];
        }
        sort(ans, ans + anss, cmp);
        printf("Case %d:\n", ca);
        printf("%d\n", d[full - 1]);
        for (int i = 0, x = 0; i < anss; i++) {
            for (int j = 0; j < n; j++){
                if(ans[i]&1<<j)printf("%d%c", j+1, ++x ==n ? '\n' : ' ');
            }

        }
    }
    return 0;
}