算法设计与分析——电路布线(动态规划)

一、问题描述

在一块电路板的上下两端分别有n个接线柱。根据电路设计,要求用导线 (i,π(i)),将上端接线柱 i 与下端接线柱 π(i) 相连,
如图,其中 π(i),1<=i<=n,是(1,2……,n)的一个排列。导线(i,π(i))称为该电路板上的第i条连线。对于任何 1<=i<s<=n,第i条连线和第s条连线相交的充分且必要条件是 π(i) > π(s)。

在制作电路板时,要求将这n条线分布到若干个绝缘层上,在同一层上的连线不能相交。电路布线问题要确定将哪些连线安排在第一层上,使得该层上有尽可能多的连线。换句话说,该问题要求确定导线集Nets = { (i,π(i)),1<=i<=n }的最大不相交子集。

 二、算法思路

1、最优子结构性质

N(i,j)表示上面节点i与下面节点j连线的左侧区域内(包括i j连线)的连线集合,MNS(i,j)表示连线左侧区域内最大不相交连线子集,Size(i,j)表示MNS(i,j)集合中连线的个数。

在这里注意N(i,j)和MNS(i,j)表示的都是集合!!内存储的是连线,Size(i,j)存储的才是最大不相交连线的个数!!

2、递归计算最优值

当i=1的时候很好理解。

当i>1时,我们还是看上面的那个连线图。

当j<π(i)时,t[8][9]=t[7][9]。

当j>=π(i)时,t[7][9]=[6][8]+1。琢磨一下是不是符合表达式?

void MNS(int C[],int n,int **size)  
{  
    for(int j=0;j<C[1];j++)  
    {  
        size[1][j]=0;  
    }  
  
    for(int j=C[1]; j<=n; j++)  
    {  
        size[1][j]=1;  
    }  
  
    for(int i=2; i<n; i++)  
    {  
        for(int j=0; j<C[i]; j++)  
        {  
            size[i][j]=size[i-1][j];//当i<c[i]的情形  
        }  
        for(int j=C[i]; j<=n; j++)  
        {  
            //当j>=c[i]时,考虑(i,c[i])是否属于MNS(i,j)的两种情况  
            size[i][j]=max(size[i-1][j],size[i-1][C[i]-1]+1);  
        }  
    }  
    size[n][n]=max(size[n-1][n],size[n-1][C[n]-1]+1);  
}  

3、构造最优解

 

 红色标明的就是算法选择的路径(向上优先,也可以用向左优先,答案都是四个,但值会有一点不同)。在斜角值改变时可以取得所求的子集。即 9->10,7->9, 5->5, 3->4

void Traceback(int C[],int **size,int n,int Net[],int& m)  
{  
    int j=n;  
    m=0;  
    for(int i=n;i>1;i--)  
    {  
        if(size[i][j]!=size[i-1][j])//此时,(i,c[i])是最大不相交子集的一条边  
        {  
            Net[m++]=i;  
            j=C[i]-1;//更新扩展连线柱区间  
        }  
    }  
    if(j>=C[1])//处理i=1的情形  
    {  
        Net[m++]=1;  
    }  
}  

4、计算复杂性

5、演示代码

#include <iostream>
using namespace std;

const int N = 10;

void MNS(int C[],int n,int **size);
void Traceback(int C[],int **size,int n,int Net[],int& m);

int main()
{
    int c[] = {0,8,7,4,2,5,1,9,3,10,6};//下标从1开始
    int **size = new int *[N+1];

    for(int i=0; i<=N; i++)
    {
        size[i] = new int[N+1];
    }

    MNS(c,N,size);

    cout<<"电路布线最大不相交连线数目为:"<<size[N][N]<<endl;

    int Net[N],m;
    Traceback(c,size,N,Net,m);

    cout<<"最大不相交连线分别为:"<<endl;
    for(int i=m-1; i>=0; i--)
    {
        cout<<"("<<Net[i]<<","<<c[Net[i]]<<") ";
    }
    cout<<endl;
    return 0;
}

void MNS(int C[],int n,int **size)
{
    for(int j=0;j<C[1];j++)
    {
        size[1][j]=0;
    }

    for(int j=C[1]; j<=n; j++)
    {
        size[1][j]=1;
    }

    for(int i=2; i<n; i++)
    {
        for(int j=0; j<C[i]; j++)
        {
            size[i][j]=size[i-1][j];//当i<c[i]的情形
        }
        for(int j=C[i]; j<=n; j++)
        {
            //当j>=c[i]时,考虑(i,c[i])是否属于MNS(i,j)的两种情况
            size[i][j]=max(size[i-1][j],size[i-1][C[i]-1]+1);
        }
    }
    size[n][n]=max(size[n-1][n],size[n-1][C[n]-1]+1);
}

void Traceback(int C[],int **size,int n,int Net[],int& m)
{
    int j=n;
    m=0;
    for(int i=n;i>1;i--)
    {
        if(size[i][j]!=size[i-1][j])//此时,(i,c[i])是最大不相交子集的一条边
        {
            Net[m++]=i;
            j=C[i]-1;//更新扩展连线柱区间
        }
    }
    if(j>=C[1])//处理i=1的情形
    {
        Net[m++]=1;
    }
}

 

posted @ 2019-10-12 11:54  王陸  阅读(3007)  评论(0编辑  收藏