斯坦纳树

斯坦纳树

<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