算法设计与分析书中程序(第07章).docVIP

【程序7-1】 多段图的向前递推算法 templateclass T T GraphT::FMultiGraph(int k, int *p) {//采用程序6-8的邻接表存储图G Tc,*cost=new float[n]; int q, *d=new int[n]; cost[n-1]=0, d[n-1]= -1; //设置向前递推的初值 for (int j=n-2; j=0; j--){ //按n-2, …, 0的次序计算cost和d float min=INFTY; //按式（7-1）计算最小值为cost[j] for (ENodeT *r=a[j]; r; r=r-nextArc) { int v=r-adjVex; if (r-w+cost[v]min) { min=r-w+cost[v]; q=v; } } cost[j]=min; d[j]=q; //q是j在最短子路径上的后继结点 } p[0]=0; p[k-1]=n-1; c=cost[0]; //p[0]和p[n-1]是源点和汇点 for(j=1; j=k-2; j++) p[j]=d[p[j-1]]; //p[i]是最短路径上第i阶段的结点 delete []cost; delete []d; return c; } 【程序7-2】 弗洛伊德算法 templateclass T void MGraphT::Floyd(T** d, int ** path) { int i, j, k; d= new T*[n]; path=new int *[n]; for(i=0; in; i++){ d[i]=new T [n]; path[i]=new int[n]; for (j=0; jn; j++){ //初始化 d[i][j]=a[i][j]; if (i！=j w[i][j]INFTY) path[i][j]=i; else path[i][j]= -1; } } for (k=0; kn; k++) //考察结点k for (i=0; in; i++) for (j=0; jn; j++) if (d[i][k]+d[k][j] d[i][j] ){ d[i][j]=d[i][k]+d[k][j]; path[i][j]=path[k][j]; } } 【程序7-3】 矩阵连乘算法 class MatrixChain { public: MatrixChain(int mSize, int *q); //创建二维数组m和s，一维数组p，并初始化 int MChain(); //一般动态规划法求最优解值 int LookupChain(); //备忘录方法求最优解值（程序7-4） void Traceback(); //构造最优解的公有函数 …… private: void Traceback(int i, int j); //构造最优解的私有递归函数 int LookupChain(int i, int j); //备忘录方法私有递归（程序7-4） int *p, **m, **s, n; }; int MatrixChain::MChain() { //求A[0:n-1]的最优解值 for (int i=0;in; i++) m[i][i]=0; for (int r=2; r=n; r++) for (int i=0; i=n-r; i++) { int j=i+r-1; m[i][j]=m[i+1][j]+p[i]*p[i+1]*p[j+1]; //m[i][j] 的初值 s[i][j]=i; for (int k=i+1; kj; k++) { int t=m[i][k]+m[k+1][j]+p[i]*p[k+1]*p[j+1]; if (tm[i][j]) { m[i][j]=t; s[i][j]=k; } } } return m[0][n-1]; } void MatrixChain::Traceback(int i, int j) { if(i==j) { coutAi; return;} if (is[i][j]) cout(; Traceback(i, s[i][j]); if (is[i][j])cout); if(s[i][j]+1j)cou