UVALive 4999 Electricity Connection

最小斯坦纳树：求图G的一个真子集里面的所有点连通，且总边权值和最小的生成树称为MinimalSteinerTree（最小斯坦纳树）。

模板：
dp[u][state]表示以u为根，指定集合中的点的连通状态为state的生成树的最小总权值。
1.通过连通状态的子集进行转移。
dp[u][state]=min{ dp[u][sub1]+dp[u][sub2] | sub1^sub2 == state }

2.在当前枚举的连通状态下，对该连通状态进行松弛操作。
dp[u][state]=min{ dp[v][state]+e[u][v] }

代码：

#include<bits/stdc++.h>using namespace std;#define CLR(A, X) memset(A, X, sizeof(A))const int INF = 0x3f3f3f3f;bool inq[70];char mp[10][10];int w[70], dir[4][2] = {1,0, -1,0, 0,-1, 0,1};int hou[1<<8], dp[70][1<<8];queue<int> Q;void spfa(int S) {    for(int i = 0; i < 64; i++) {        inq[i] = 0;        if(dp[i][S] != INF) {            Q.push(i);            inq[i] = 1;        }    }    while(!Q.empty()) {        int u = Q.front(); Q.pop();        inq[u] = 0;        for(int i = 0; i < 4; i++) {            int x = u/8+dir[i][0], y = u%8+dir[i][1];            if(x<0||y<0||x>=8||y>=8) continue;            int v = x*8+y;            int t = dp[u][S]+1+w[v];            if(dp[v][S] > t) {                dp[v][S] = t;                if(!inq[v]) {                    inq[v] = 1;                    Q.push(v);                }            }        }    }}int main() {    int T;    scanf("%d", &T);    while(T--) {        CLR(dp, INF); CLR(hou, -1);        int pl, pw, r, k = 0;        scanf("%d%d", &pl, &pw);        for(int i = 0; i < 8; i++) {            scanf("%s", mp[i]);            for(int j = 0; j < 8; j++) {                int u = i*8+j;                if(mp[i][j] == 'G') r = u, w[u] = 0;                else if(mp[i][j] == '.') w[u] = pl;                else if(mp[i][j] == 'W') w[u] = pw;                else { w[u] = 0; hou[1<<k] = u; k++; }            }        }        int all = (1<<k)-1;        for(int S = 0; S <= all; S++) {            if(~hou[S]) dp[hou[S]][S] = 0;            for(int i = 0; i < 64; i++) {                for(int S0 = (S-1)&S; S0; S0 = (S0-1)&S) {                    dp[i][S] = min(dp[i][S], dp[i][S0]+dp[i][S^S0]-w[i]);                }            }            spfa(S);        }        static int cas = 0;        printf("Case %d: %d\n", ++cas, dp[r][all]);    }    return 0;}