斯坦纳树
<strong><span style="font-family:FangSong_GB2312;font-size:14px;">//-------------------- Minimal Steiner Tree-------------------
//G(V,E),A 是V 的一个子集,求至少包含A 中所有点的最小子树
//时间复杂度O( N*M + N*2^A*(2^A+N) )
//INIT: dis[][]距离矩阵; id[]置为集合A 中点的标号;
//CALL: steiner(int n,int m);
const int inf=1000000000;
const int N=55; //总顶点数
const int A=5; //子点数
int id[A]; //包含顶点的序号
int dp[1<<A][N]; //dp[i][j]表示实际图中包含j点的集合i 的MST
bool vis[N];
int dis[N][N]; //距离矩阵
int steiner(int n,int m) {
int i,j,k,top=1<<m;
for (int k = 0; k < n; k++) //最短路算法
for (int i = 0; i < n; i++)
for (int j = 0; j < n; j++)
dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
for(i=0;i<top;i++)
for(j=0;j<n;j++) dp[i][j]=inf;
for(i=0;i<m;i++)
for(j=0;j<n;j++) dp[1<<i][j]=dis[id[i]][j];
for(i=1;i<top;i++) {
if((i&(i-1))==0) continue;
for(k=0;k<n;k++) //init
for(j=(i-1)&i;j>0;j=(j-1)&i)
dp[i][k]=min(dp[i][k],dp[j][k]+dp[i^j][k]);
for(j=0;j<n;j++) //update
for(k=0;k<n;k++)
dp[i][j]=min(dp[i][j],dp[i][k]+dis[k][j]);
}
int ret=inf;
for(i=0;i<n;i++) ret=min(ret,dp[top-1][i]);
return ret;
}
</span></strong>
<strong><span style="font-family:FangSong_GB2312;font-size:14px;">//给定A个点指定两两点联通的斯坦纳树/森林
int steiner(int n,int m) {
int i,j,k,top=1<<m;
for (int k = 0; k < n; k++)
for (int i = 0; i < n; i++)
for (int j = 0; j < n; j++)
dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
for(i=0;i<top;i++)
for(j=0;j<n;j++) dp[i][j]=inf;
for(i=0;i<m;i++)
for(j=0;j<n;j++) dp[1<<i][j]=dis[id[i]][j];
for(i=1;i<top;i++) {
if((i&(i-1))==0) continue;
for(k=0;k<n;k++) //init
for(j=(i-1)&i;j>0;j=(j-1)&i)
dp[i][k]=min(dp[i][k],dp[j][k]+dp[i^j][k]);
for(j=0;j<n;j++) //update
for(k=0;k<n;k++)
dp[i][j]=min(dp[i][j],dp[i][k]+dis[k][j]);
}
int ans[1<<4]; //一共八个点 前四个与后四个对应联通
memset(ans,-1,sizeof(ans));
for (int i=0;i<(1<<4);i++){
int k=i|(i<<4);
int ret=inf;
for(int j=0;j<n;j++) ret=min(ret,dp[k][j]);
ans[i]=ret;
for (int j=i&(i-1);j;j=(j-1)&i){
int kk=j^i;
if (ans[j]!=-1&&ans[kk]!=-1) ans[i]=min(ans[i],ans[j]+ans[kk]);
}
}
return ans[(1<<4)-1];
}</span></strong>
<strong><span style="font-family:FangSong_GB2312;font-size:14px;">//给定A个点前A/2个与后A/2个必须对应联通(但是没有指定如何对应)的斯坦纳树/森林
const int inf=1000000000;
const int N=55; //总顶点数
const int A=11; //子点数
int id[A]; //包含顶点的序号
int dp[1<<A][N]; //dp[i][j]表示实际图中包含j点的集合i 的MST
int d[1<<A];
bool vis[N];
int dis[N][N]; //距离矩阵
int ok(int num,int b){
int hx=0;
for (int i=0;i<b;i++){
if (num&(1<<i)) hx++;
if (num&(1<<(i+b))) hx--;
}
return hx==0;
}
int steiner(int n,int m) {
int i,j,k,top=1<<m;
for (int k = 0; k < n; k++)
for (int i = 0; i < n; i++)
for (int j = 0; j < n; j++)
dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
for(i=0;i<top;i++)
for(j=0;j<n;j++) dp[i][j]=inf;
for(i=0;i<m;i++)
for(j=0;j<n;j++) dp[1<<i][j]=dis[id[i]][j];
for(i=1;i<top;i++) {
d[i]=inf;
if((i&(i-1))==0) continue;
for(k=0;k<n;k++) //init
for(j=(i-1)&i;j>0;j=(j-1)&i)
dp[i][k]=min(dp[i][k],dp[j][k]+dp[i^j][k]);
for(j=0;j<n;j++) //update
for(k=0;k<n;k++)
dp[i][j]=min(dp[i][j],dp[i][k]+dis[k][j]);
int best=0;
for (j=1;j<n;j++) if (dp[i][j]<dp[i][best]) best=j;
d[i]=dp[i][best];
if (ok(i,m/2)) for (int j=(i-1)&i;j;j=(j-1)&i) d[i]=min(d[i],d[i^j]+d[j]);
else d[i]=inf;
}
if (d[top-1]>=inf) return -1;
else return d[top-1];
}</span></strong>
题目:
Poj 3123
Hdu 4085
Zoj 3613
Zoj 3613