山东师范大学信息科学与工程学院软件工程研究所 徐连诚 2006年11月20日

Slides:



Advertisements
Similar presentations
第 8 章 数组 计算机科学学院 李淮 Tel QQ
Advertisements

1 第 3 章 C++ 中的条件与循环 第 3 次见面! acm.nefu.edu.cn/C++_03.ppt.
While 迴圈 - 不知重複執行次數
数据结构 杨鹏宇 QQ: 版权所有,转载或翻印必究
算法设计与分析 李清勇 教授 办公室:第九教学楼北201.
第十章 依赖于机器的优化 在指令级并行的机器上,程序的运行速度依赖于下面几个因素
第5章 回溯法 欢迎辞.
汽车在( )上行驶.
语文版九年级(下) 多媒体课件.
AI人工智慧報告 黑白棋 班級:資工四乙 學號:498G0009 姓名:盧冠妤.
贪婪算法.
第六章 分支限界法.
国家和我省禽业发展政策 和扶持项目解读 安徽省畜牧兽医局
图论算法 ---最大流问题 长沙市雅礼中学 朱 全 民.
第十章 内部排序 知识点3:快速排序.
第九章 排序 插入排序 交换排序 选择排序 归并排序 基数排序.
資料結構與演算法 課程教學投影片.
第5章 回溯法 “通用的解题法” 欢迎辞.
动态规划 中山纪念中学 宋新波 本模板来源于网络,由第一课件网整理发布,免费分享给大家使用。
第五章 树 东南大学计算机学院 方效林 本课件借鉴了清华大学殷人昆老师 和哈尔滨工业大学张岩老师的课件.
函數(一) 自訂函數、遞迴函數 綠園.
4.1 概述 4.2 类与对象的实现 4.3 对象的初始化和析构 4.4 类的包含 4.5 类模板
计算复杂性和算法分析 计算机科学导论第六讲
第12章 從C到C++語言 12-1 C++語言的基礎 12-2 C++語言的輸出與輸入 12-3 C++語言的動態記憶體配置
書名 Java於資料結構與演算法之實習應用
计算机网络讲义 第5章 批量数据处理—数组 一维数组 排序和查找 二维数组 字符串.
6 使用者函數 6.1 函數定義 宣告函數 呼叫函數 呼叫多個函數 6-6
容斥原理 若干应用 王瑶 张梦微 张雯露 2019/1/11.
C++语言程序设计 第二章 C++简单程序设计.
第二章 插值.
数据结构与数据库 习题课(2) 2016年11月25日.
Chapter12 Graphs Definitions Applications Properties
计算机算法设计与分析(第3版) 王晓东 编著 电子工业出版社.
主讲人: 吕敏 { } Spring 2016 ,USTC 算法基础 主讲人: 吕敏 { } Spring 2016 ,USTC.
教 师:曾晓东 电 话: E_mail: 计算机软件技术基础 教 师:曾晓东 电 话: E_mail:
C++大学基础教程 第3章 C++控制语句 北京科技大学 信息基础科学系.
第四章 第三节 最短路径算法.
Week 2: 程式設計概念與 演算法的效能評估
C语言版 软件与理论教研室 张泽宝 哈尔滨工程大学计算机科学与技术学院.
问题求解 入门.
算法设计复习 内容提要: 算法设计思想回顾(递归和分治、动态规划、贪心算法、回溯法、分支限界法) 经典例子讲解 2019/4/9.
程式結構&語法.
第六章 分支限界法 理解分支限界法的剪枝搜索策略。 掌握分支限界法的算法框架 队列式(FIFO)分支限界法 优先队列式分支限界法.
第5章 回溯法 欢迎辞.
第三章 C++的语句和简单的程序设计 主要内容:
数据结构选讲 北京大学 陈瑜希.
数 据 结 构 与 算 法.
C++语言程序设计 C++语言程序设计 第三章 控制语句 第十一组 C++语言程序设计.
第十二章 施工组织计划技术.
第1章 绪论(二) 教学目标 理解算法的特性及评价标准 掌握算法时间复杂度和空间复杂度的分析方法 1/
程式的時間與空間 Time and Space in Programming
C++语言程序设计教程 第2章 数据类型与表达式 第2章 数据类型与表达式 制作人:杨进才 沈显君.
第2章 认识C语言 教学要点 2. 1 项目二C语言程序识读 2 .2 项目三班级成绩排名 2 .3 知识链接 返回.
二分查找的应用 ——李江贝, 郝琰 2017年9月28日.
資料結構簡介 綠園.
第7章 概率算法 欢迎辞.
第三章 贪心方法 (Greedy Technique)
第八节 算术运算符和算术表达式.
#include <iostream.h>
第六章 贪心算法.
遞迴 Recursion.
第8章 排序 8.1排序基本概念 8.2插入类排序 8.3交换类排序 8.4选择类排序 8.5归并排序 8.6基数排序
第1章 数据结构基础概论 本章主要介绍以下内容 数据结构研究的主要内容 数据结构中涉及的基本概念 算法的概念、描述方法以及评价标准.
主讲人: 吕敏 { } Spring 2016,USTC 算法基础 主讲人: 吕敏 { } Spring 2016,USTC.
第5章 回溯法 欢迎辞.
《数据结构与算法设计》第一部分 面向对象的C++程序设计基础.
近似算法 东南大学计算机学院 方效林.
基本資料型態 變數與常數 運算子 基本的資料處理 授課:ANT 日期:2014/03/03.
第三章 流程控制 程序的运行流程 选择结构语句 循环结构语句 主讲:李祥 时间:2015年10月.
Q1 以下是10個初生嬰兒,首6個月的體重改變紀錄
Presentation transcript:

山东师范大学信息科学与工程学院软件工程研究所 徐连诚 E-Mail:lchxu@163.com 2006年11月20日 算法设计与分析 山东师范大学信息科学与工程学院软件工程研究所 徐连诚 E-Mail:lchxu@163.com 2006年11月20日

第五章 回溯法 学习要点 理解回溯法的深度优先搜索策略 掌握用回溯法解题的算法框架 1)递归回溯最优子结构性质 2)迭代回溯贪心选择性质 3)子集树算法框架 4)排列树算法框架 通过应用范例学习回溯法的设计策略 1)装载问题; 2)批处理作业调度; 3)符号三角形问题 4)n后问题; 5)0-1背包问题; 6)最大团问题; 7)图的m着色问题 8)旅行售货员问题 9)圆排列问题 10)电路板排列问题 11)连续邮资问题

引言 有许多问题,当需要找出它的解集或者要求回答什么解是满足某些约束条件的最佳解时,往往要使用回溯法。 回溯法的基本做法是搜索,或是一种组织得井井有条的、能避免不必要搜索的穷举式搜索法。这种方法适用于解一些组合数相当大的问题。 回溯法在问题的解空间树中,按深度优先策略,从根结点出发搜索解空间树。算法搜索至解空间树的任意一点时,先判断该结点是否包含问题的解:如果肯定不包含,则跳过对该结点为根的子树的搜索,逐层向其祖先结点回溯;否则,进入该子树,继续按深度优先策略搜索。

5.1 回溯法的算法框架 本节介绍回溯法算法框架的有关问题: 一、问题的解空间 二、回溯法的基本思想 三、递归回溯 四、迭代回溯 五、子集树与排列树

一、问题的解空间 应用回溯法解问题时,首先应明确定义问题的解空间。问题的解空间应至少包含问题的一个(最优)解。 问题的解向量:回溯法希望一个问题的解能够表示成一个n元式(x1,x2,…,xn)的形式。 显约束:对分量xi的取值限定。 隐约束:为满足问题的解而对不同分量之间施加的约束。 解空间:对于问题的一个实例,解向量满足显式约束条件的所有多元组,构成了该实例的一个解空间。 注意:同一个问题可以有多种表示,有些表示方法更简单,所需表示的状态空间更小(存储量少,搜索方法简单)。 例如,对于有n种可选物品的0-1背包问题,其解空间由长度为n的0-1向量组成。 n=3时的0-1背包问题用完全二叉树表示的解空间

二、回溯法的基本思想 回溯法的基本步骤: 常用剪枝函数: 关于复杂性: (1)针对所给问题,定义问题的解空间; (2)确定易于搜索的解空间结构; (3)以深度优先方式搜索解空间,并在搜索过程中用剪枝函数避免无效搜索。 常用剪枝函数: 用约束函数在扩展结点处剪去不满足约束的子树; 用限界函数剪去得不到最优解的子树。 关于复杂性: 用回溯法解题的一个显著特征是在搜索过程中动态产生问题的解空间。在任何时刻,算法只保存从根结点到当前扩展结点的路径。如果解空间树中从根结点到叶结点的最长路径的长度为h(n),则回溯法所需的计算空间通常为O(h(n))。而显式地存储整个解空间则需要O(2h(n))或O(h(n)!)内存空间。

生成问题状态的基本方法 扩展结点:一个正在产生儿子的结点称为扩展结点 活结点:一个自身已生成但其儿子还没有全部生成的节点称做活结点 死结点:一个所有儿子已经产生的结点称做死结点 深度优先的问题状态生成法:如果对一个扩展结点R,一旦产生了它的一个儿子C,就把C当做新的扩展结点。在完成对子树C(以C为根的子树)的穷尽搜索之后,将R重新变成扩展结点,继续生成R的下一个儿子(如果存在) 宽度优先的问题状态生成法:在一个扩展结点变成死结点之前,它一直是扩展结点 回溯法:为了避免生成那些不可能产生最佳解的问题状态,要不断地利用限界函数(bounding function)来处死那些实际上不可能产生所需解的活结点,以减少问题的计算量。具有限界函数的深度优先生成法称为回溯法。

示例1 0-1背包问题 n=3, C=30, w={16, 15, 15}, v={45, 25, 25} A Cr=C=30,V=0 B 开始时,Cr=C=30,V=0,A为唯一活结点,也是当前扩展结点 扩展A,先到达B结点 Cr=Cr-w1=14,V=V+v1=45 此时A、B为活结点,B成为当前扩展结点 扩展B,先到达D Cr<w2,D导致一个不可行解,回溯到B 再扩展B到达E E可行,此时A、B、E是活结点,E成为新的扩展结点 扩展E,先到达J Cr<w3,J导致一个不可行解,回溯到E 再次扩展E到达K 由于K是叶结点,即得到一个可行解x=(1,0,0),V=45 K不可扩展,成为死结点,返回到E E没有可扩展结点,成为死结点,返回到B B没有可扩展结点,成为死结点,返回到A A再次成为扩展结点,扩展A到达C Cr=30,V=0,活结点为A、C,C为当前扩展结点 扩展C,先到达F Cr=Cr-w2=15,V=V+v2=25,此时活结点为A、C、F,F成为当前扩展结点 扩展F,先到达L Cr=Cr-w3=0,V=V+v3=50 L是叶结点,且50>45,皆得到一个可行解x=(0,1,1),V=50 L不可扩展,成为死结点,返回到F 再扩展F到达M M是叶结点,且25<50,不是最优解 M不可扩展,成为死结点,返回到F F没有可扩展结点,成为死结点,返回到C 再扩展C到达G Cr=30,V=0,活结点为A、C、G,G为当前扩展结点 扩展G,先到达N,N是叶结点,且25<50,不是最优解,又N不可扩展,返回到G 再扩展G到达O,O是叶结点,且0<50,不是最优解,又O不可扩展,返回到G G没有可扩展结点,成为死结点,返回到C C没有可扩展结点,成为死结点,返回到A A没有可扩展结点,成为死结点,算法结束,最优解X=(0,1,1),最优值V=50 A Cr=C=30,V=0 B w1=16,v1=45 Cr=14,V=45 C Cr=30,V=0 E Cr=14 V=45 D Cr<w2 不可行解 F w2=15,v2=25 Cr=15,V=25 M 25<50 不是 最优 解 J Cr<w3 不可行解 L w3=15,v3=25 Cr=0,V=50 50>45 x=(0,1,1) H I K Cr=14 V=45 x=(1,0,0)

示例2 旅行售货员问题 问题描述:某售货员要到若干城市去推销商品,一直各城市之间的路程,他要选定一条从驻地出发,经过每个城市一遍,最后回到住地的路线,使总的路程最短。 该问题是一个NP完全问题, 有(n-1)!条可选路线 最优解(1,3,2,4,1),最优值25 1 2 3 4 20 6 30 5 10 A B C D E F G H I J K L M N O P Q

三、递归回溯 回溯法对解空间作深度优先搜索,因此,在一般情况下用递归方法实现回溯法。 void backtrack (int t) { if (t>n) output(x); else for (int i=f(n,t);i<=g(n,t);i++) x[t]=h(i); if (constraint(t)&&bound(t)) backtrack(t+1); }

四、迭代回溯 void iterativeBacktrack () { int t=1; while (t>0) 采用树的非递归深度优先遍历算法,可将回溯法表示为一个非递归迭代过程。 void iterativeBacktrack () { int t=1; while (t>0) if (f(n,t)<=g(n,t)) for (int i=f(n,t);i<=g(n,t);i++) { x[t]=h(i); if (constraint(t)&&bound(t)) {if (solution(t)) output(x); else t++;} } else t--;

五、子集树与排列树 遍历子集树需O(2n)计算时间 void backtrack (int t) { if (t>n) output(x); else for (int i=0;i<=1;i++) { x[t]=i; if (legal(t)) backtrack(t+1); } 遍历子集树需O(n!)计算时间 void backtrack (int t) { if (t>n) output(x); else for (int i=t;i<=n;i++) { x[t]=i; if (legal(t)) backtrack(t+1); }

5.2 装载问题 一、问题描述 二、问题分析 三、算法设计 四、算法描述

一、问题描述 有一批共n个集装箱要装上2艘载重量分别为C1和C2的轮船,其中集装箱i的重量为wi,且∑wi≤C1+C2 装载问题要求确定是否有一个合理的装载方案可将这个集装箱装上这2艘轮船。如果有,找出一种装载方案。 容易证明,如果一个给定装载问题有解,则采用下面的策略可得到最优装载方案: (1)首先将第一艘轮船尽可能装满; (2)将剩余的集装箱装上第二艘轮船。

二、问题分析 将第一艘轮船尽可能装满等价于选取全体集装箱的一个子集,使该子集中集装箱重量之和最接近。由此可知,装载问题等价于以下特殊的0-1背包问题。

三、算法设计 解空间:子集树 可行性约束函数(选择当前元素): 上界函数(不选择当前元素): 当前载重量CW+剩余集装箱的重量R ≤当前最优载重量BestW 用回溯法设计解装载问题的O(2n)计算时间算法。在某些情况下该算法优于动态规划算法。

四、算法描述 void backtrack (int i) {// 搜索第i层结点 if (i > n) // 到达叶结点 更新最优解bestx,bestw;return; r -= w[i]; if (cw + w[i] <= c) //搜索左子树 { x[i] = 1; cw += w[i]; backtrack(i + 1); cw -= w[i]; } if (cw + r > bestw) { x[i] = 0; // 搜索右子树 backtrack(i + 1); } r += w[i]; }

5.3 批处理作业调度 一、问题描述 二、实例分析 三、算法描述

一、问题描述 给定n个作业的集合{J1,J2,…,Jn}。每个作业必须先由机器1处理,然后由机器2处理。作业Ji需要机器j的处理时间为tji。对于一个确定的作业调度,设Fji是作业i在机器j上完成处理的时间。所有作业在机器2上完成处理的时间和称为该作业调度的完成时间和。 批处理作业调度问题要求对于给定的n个作业,制定最佳作业调度方案,使其完成时间和达到最小。

二、实例分析 tji 机器1 机器2 作业1 2 1 作业2 3 作业3 这3个作业的6种可能的调度方案是1,2,3;1,3,2;2,1,3;2,3,1;3,1,2;3,2,1; 它们所相应的完成时间和分别是19,18,20,21,19,19。 最佳调度方案是1,3,2,其完成时间和为18。

三、算法描述 解空间:排列树 复杂性:T(n)=O(n!) void Flowshop::Backtrack(int i) { if (i > n) { for (int j = 1; j <= n; j++) bestx[j] = x[j]; bestf = f; } else for (int j = i; j <= n; j++) { f1+=M[x[j]][1]; f2[i]=((f2[i-1]>f1)?f2[i-1]:f1)+M[x[j]][2]; f+=f2[i]; if (f < bestf) { Swap(x[i], x[j]); Backtrack(i+1); f1- =M[x[j]][1]; f- =f2[i]; 解空间:排列树 复杂性:T(n)=O(n!)

5.4 符号三角形问题 一、问题描述 二、问题分析 三、算法描述

一、问题描述 右图是由14个“+”和14个“-”组成的符号三角形。2个同号下面都是“+”,2个异号下面都是“-”。 在一般情况下,符号三角形的第一行有n个符号。 符号三角形问题要求对于给定的n,计算有多少个不同的符号三角形,使其所含的“+”和“-”的个数相同。 + + - + - + + + - - - - + - + + + - - + + - - + - - - +

当前符号三角形所包含的“+”个数与“-”个数 二、问题分析 解向量:用n元组x[1:n]表示符号三角形的第一行。 可行性约束函数: 当前符号三角形所包含的“+”个数与“-”个数 均不超过n*(n+1)/4 无解的判断: n*(n+1)/2为奇数 复杂度分析:计算可行性约束需要O(n)时间,在最坏情况下有 O(2n)个结点需要计算可行性约束,故解符号三角形问题的回溯算法所需的计算时间为 O(n2n)。

三、算法描述 void Triangle::Backtrack(int t) { if ((count>half)||(t*(t-1)/2-count>half)) return; if (t>n) sum++; else for (int i=0;i<2;i++) { p[1][t]=i; count+=i; for (int j=2;j<=t;j++) { p[j][t-j+1]=p[j-1][t-j+1]^p[j-1][t-j+2]; count+=p[j][t-j+1]; } Backtrack(t+1); for (int j=2;j<=t;j++) count-=p[j][t-j+1]; count-=i;

5.5 n后问题 一、问题描述 二、问题分析 三、算法描述

一、问题描述 在n×n格的棋盘上放置彼此不受攻击的n个皇后。按照国际象棋的规则,皇后可以攻击与之处在同一行或同一列或同一斜线上的棋子。n后问题等价于在n×n格的棋盘上放置n个皇后,任何2个皇后不放在同一行或同一列或同一斜线上。 1 Q 2 3 4 5 6 7 8

二、问题分析 解向量:(x1, x2, … , xn) 显约束:xi=1,2, … ,n 隐约束: 1)不同列:xi≠xj  2)不处于同一正、反对角线:|i-j|≠|xi-xj|

三、算法描述 bool Queen::Place(int k) { for (int j=1;j<k;j++) if ((abs(k-j)==abs(x[j]-x[k]))||(x[j]==x[k])) return false; return true; } void Queen::Backtrack(int t) if (t>n) sum++; else for (int i=1;i<=n;i++) { x[t]=i; if (Place(t)) Backtrack(t+1);

5.6 0-1背包问题 解空间:子集树 可行性约束函数: ∑wixi≤C 上界函数: Bound() 算法:略 template<class Typew, class Typep> Typep Knap<Typew, Typep>::Bound(int i) {// 计算上界 Typew cleft = c - cw; // 剩余容量 Typep b = cp; // 以物品单位重量价值递减序装入物品 while (i <= n && w[i] <= cleft) { cleft -= w[i]; b += p[i]; i++; } // 装满背包 if (i <= n) b += p[i]/w[i] * cleft; return b;

5.7 最大团问题 一、基本概念 二、问题分析 三、算法描述 四、进一步改进

一、基本概念 完全子图:给定无向图G=(V,E)。如果UV,且对任意u,vU有(u,v)E,则称U是G的完全子图。 团:G的完全子图U是G的团当且仅当U不包含在G的更大的完全子图中。 G的最大团:是指G中所含顶点数最多的团。 空子图:如果UV且对任意u,vU有(u,v)E,则称U是G的空子图。 独立集:G的空子图U是G的独立集当且仅当U不包含在G的更大的空子图中。 G的最大独立集:是G中所含顶点数最多的独立集。 补图:对于任一无向图G=(V,E)其补图G=(V1,E1)定义为:V1=V,且(u,v)E1当且仅当(u,v)E。 U是G的最大团当且仅当U是G的最大独立集。 1 2 3 4 5

最大团问题的回溯算法backtrack所需的 二、问题分析 解空间:子集树 可行性约束函数: 顶点i到 已选入的顶点集中每一个顶点都有边相连。 上界函数: 有足够多的可选择顶点使得算法有可能 在右子树中找到更大的团。 复杂度分析: 最大团问题的回溯算法backtrack所需的 计算时间显然为O(n2n)。

三、算法描述 void Clique::Backtrack(int i)// 计算最大团 { if (i > n) {// 到达叶结点 for (int j = 1; j <= n; j++) bestx[j] = x[j]; bestn = cn; return;} int OK = 1; // 检查顶点 i 与当前团的连接 for (int j = 1; j < i; j++) if (x[j] && a[i][j] == 0) { // i与j不相连 OK = 0; break;} if (OK) {// 进入左子树 x[i] = 1; cn++; Backtrack(i+1); x[i] = 0; cn--;} if (cn + n - i > bestn) {// 进入右子树 x[i] = 0; Backtrack(i+1);} }

四、进一步改进 选择合适的搜索顺序,可以使得上界函数更有效的发挥作用。例如在搜索之前可以将顶点按度从小到大排序。这在某种意义上相当于给回溯法加入了启发性。 定义Si={vi,vi+1,...,vn},依次求出Sn,Sn-1,...,S1的解。从而得到一个更精确的上界函数,若cn+Si<=max则剪枝。同时注意到:从Si+1到Si,如果找到一个更大的团,那么vi必然属于找到的团,此时有Si=Si+1+1,否则Si=Si+1。因此只要max的值被更新过,就可以确定已经找到最大值,不必再往下搜索了。

5.8 图的m着色问题 一、基本概念 二、问题分析 三、算法描述 四、复杂性分析

一、基本概念 给定无向连通图G和m种不同的颜色。用这些颜色为图G的各顶点着色,每个顶点着一种颜色。是否有一种着色法使G中每条边的2个顶点着不同颜色。这个问题是图的m可着色判定问题。 若一个图最少需要m种颜色才能使图中每条边连接的2个顶点着不同颜色,则称这个数m为该图的色数。 求一个图的色数m的问题称为图的m可着色优化问题。

二、问题分析 解向量:(x1, x2, … , xn)表示顶点i所着颜色x[i] 可行性约束函数:顶点i与已着色的相邻顶点颜色不重复。

三、算法描述 void Color::Backtrack(int t) { if (t>n) { sum++; for (int i=1; i<=n; i++) cout << x[i] << ' '; cout << endl; } else for (int i=1;i<=m;i++) { x[t]=i; if (Ok(t)) Backtrack(t+1); bool Color::Ok(int k) // 检查颜色可用性 { for (int j=1;j<=n;j++) if ((a[k][j]==1)&&(x[j]==x[k])) return false; return true;

∑mi(mn)=nm(mn-1)/(m-1)=O(nmn) 四、复杂性分析 图m可着色问题的解空间树中内结点个数是∑mi(0≤i≤n-1)。 对于每一个内结点,在最坏情况下,用ok检查当前扩展结点的每一个儿子所相应的颜色可用性需耗时O(mn)。因此,回溯法总的时间耗费是 ∑mi(mn)=nm(mn-1)/(m-1)=O(nmn) (0≤i≤n-1)

5.9 旅行售货员问题 解空间:排列树 复杂度分析: void Backtrack(int i) { if (i == n) { if (a[x[n-1]][x[n]] != NoEdge && a[x[n]][1] != NoEdge && (cc + a[x[n-1]][x[n]] + a[x[n]][1] < bestc || bestc == NoEdge)) { for (int j = 1; j <= n; j++) bestx[j] = x[j]; bestc = cc + a[x[n-1]][x[n]] + a[x[n]][1];}} else {// 是否可进入x[j]子树? for (int j = i; j <= n; j++) // 搜索子树 if (a[x[i-1]][x[j]] != NoEdge && (cc+a[x[i-1]][x[i]] < bestc || bestc == NoEdge)) { Swap(x[i], x[j]); cc += a[x[i-1]][x[i]]; Backtrack(i+1); cc -= a[x[i-1]][x[i]]; Swap(x[i], x[j]);}} } 解空间:排列树 复杂度分析: 算法backtrack在最坏情况下可能需要更新当前最优解O((n-1)!)次,每次更新bestx需计算时间O(n),从而整个算法的计算时间复杂性为O(n!)。

5.10 圆排列问题 一、问题描述 二、算法描述 三、算法改进

一、问题描述 给定n个大小不等的圆c1,c2,…,cn,现要将这n个圆排进一个矩形框中,且要求各圆与矩形框的底边相切。圆排列问题要求从n个圆的所有排列中找出有最小长度的圆排列。 例如,当n=3,且所给的3个圆的半径分别为1,1,2时,这3个圆的最小长度的圆排列如图所示。其最小长度为2+4√2。

二、算法描述 void Circle::Backtrack(int t) { if (t>n) Compute(); else for (int j = t; j <= n; j++) { Swap(r[t], r[j]); float centerx=Center(t); if (centerx+r[t]+r[1]<min) {//下界约束 x[t]=centerx; Backtrack(t+1);} Swap(r[t], r[j]);} } float Circle::Center(int t) {// 计算当前所选择圆的圆心横坐标 float temp=0; for (int j=1;j<t;j++) { float valuex=x[j]+2.0*sqrt(r[t]*r[j]); if (valuex>temp) temp=valuex; } return temp; void Circle::Compute(void) {// 计算当前圆排列的长度 float low=0, high=0; for (int i=1;i<=n;i++) { if (x[i]-r[i]<low) low=x[i]-r[i]; if (x[i]+r[i]>high) high=x[i]+r[i]; if (high-low<min) min=high-low;

三、算法改进 上述算法尚有许多改进的余地。例如,象1,2,…,n-1,n和n,n-1, …,2,1这种互为镜像的排列具有相同的圆排列长度,只计算一个就够了,可减少约一半的计算量。 另一方面,如果所给的n个圆中有k个圆有相同的半径,则这k个圆产生的k!个完全相同的圆排列,只计算一个就够了。

5.11 电路板排列问题 一、问题描述 二、算法设计 三、算法描述 四、复杂性分析

一、问题描述 设B={1,2,…,n}是n块电路板的集合。集合L={N1,N2,…,Nm}是这n块电路板的m个连接块。其中每个连接块Ni是B的一个子集,且Ni中的电路板用一根导线连接在一起。不同的连接快可能包含相同的电路板。 设x是n块电路板的一个排列,则电路板的排列密度Density(x)定义为跨越相邻电路板插槽的最大连线数。 电路板排列问题要求度与给定的电路板连接条件(连接块),确定电路板的最佳排列,使其具有最小密度。

二、算法设计 问题空间:排列树

三、算法描述 P152-154

四、复杂性分析 每个结点处Backtrack函数的计算时间:O(m) 结点数:O(n!) 总的计算时间:O(mn!)

5.12 连续邮资问题 本节要点 回顾回溯法的基本思想 连续邮资问题描述 问题分析 算法设计

一、回溯法的基本思想 1、确定问题的解空间 2、找出适当的剪枝函数 3、以深度优先的方式搜索解空间 子集树问题:装载问题、符号三角形问题、0-1背包问题、最大团问题 排列树问题:批处理作业调度、n后问题、旅行售货员问题、圆排列问题、电路板排列问题 其他:图的m着色问题 2、找出适当的剪枝函数 约束函数 限界函数 3、以深度优先的方式搜索解空间 递归回溯 迭代回溯

二、问题描述 假设国家发行了n种不同面值的邮票,并且规定每张信封上最多只允许贴m张邮票。连续邮资问题要求对于给定的n和m的值,给出邮票面值的最佳设计,在1张信封上可贴出从邮资1开始,增量为1的最大连续邮资区间。(NOIP99) 例如,当n=2、m=3时,如果面值分别为1、4,则在l-6之间的每一个邮资值都能得到(当然还有8、9和12);如果面值分别为1、3,则在1-7之间的每一个邮资值都能得到。可以验证当n=2、m=3时,7就是可以得到连续的邮资最大值,面值为l、3。 又如,当n=5和m=4时,面值为(1,3,11,15,32)的5种邮票可以贴出邮资的最大连续邮资区间是1到70。

三、问题分析 基本思路:搜索所有可行解,找出最大连续邮资区间的方案 解向量:用n元组x[1:n]表示n种不同的邮票面值,并约定它们从小到大排列。x[1]=1是唯一的选择。 可行性约束函数:已选定x[1:i-1],最大连续邮资区间是1—r,接下来x[i]的可取值范围是x[i-1]+1—r+1。 如何确定r的值:计算X[1:i]的最大连续邮资区间在本算法中被频繁使用到,因此势必要找到一个高效的方法。考虑到直接递归的求解复杂度太高,我们不妨尝试计算用不超过m张面值为x[1:i]的邮票贴出邮资k所需的最少邮票数y[k]。通过y[k]可以很快推出r的值。事实上,y[k]可以通过递推在O(n)时间内解决: for (int j=0; j<= x[i-2]*(m-1);j++) if (y[j]<m) for (int k=1;k<=m-y[j];k++) if (y[j]+k<y[j+x[i-1]*k]) y[j+x[i-1]*k]=y[j]+k; while (y[r]<maxint) r++;

四、算法设计 P155-156

5.13 回溯法的效率分析 本节要点: 一、影响回溯算法效率的因素 二、重排原理 三、回溯法的效率 四、概率方法

一、影响回溯算法效率的因素 通过前面具体实例的讨论容易看出,回溯算法的效率在很大程度上依赖于以下因素: (1)产生x[k]的时间; (2)满足显约束的x[k]值的个数; (3)计算约束函数constraint的时间; (4)计算上界函数bound的时间; (5)满足约束函数和上界函数约束的所有x[k]的个数。 好的约束函数能显著地减少所生成的结点数。但这样的约束函数往往计算量较大。因此,在选择约束函数时通常存在生成结点数与约束函数计算量之间的折衷。

二、重排原理 对于许多问题而言,在搜索试探时选取x[i]的值顺序是任意的。在其它条件相当的前提下,让可取值最少的x[i]优先。从右图中关于同一问题的2棵不同解空间树,可以体会到这种策略的潜力。 图(a)中,从第1层剪去1棵子树,则从所有应当考虑的分支中一次消去12个分支。对于图(b),虽然同样从第1层剪去1棵子树,却只从应当考虑的分支中消去8个分支。前者的效果明显比后者好。 (a) (b)

三、回溯法的效率 解空间的结构已经选定,影响回溯法效率的前三个因素就可以确定,只剩下生成结点的数目是可变的,它将随问题的具体内容及结点的不同生成方式而变动。即使对同一问题的不同实例,回溯法所产生的结点数也会有很大变化。对于一个实例,回溯法可能只产生O(n)个结点,而对另一个非常相近的实例,回溯法可能就会产生解空间中的所有结点。 如果解空间的结点数是2n或n!,则在最坏情况下,回溯法的时间耗费一般为O(p(n)2n)或O(q(n)n!)。其中p(n)和q(n)均为n的多项式。 n=3,c=30 p=(100,1,1),v=(30,20,20) p=(1,1,100),v=(20,20,30)

四、概率方法 对于一个具体的问题来说,回溯法的有效性往往就体现在当问题实例的规模n较大时,它能够用很少的时间就求出问题的解。而对于一个问题的具体实例我们又很难预测回溯法的算法行为。特别是我们很难估计出回溯法在解这一具体事例时所产生的结点数,这是分析回溯法效率的主要困难。 下面我们介绍一个概率算法,用于克服这一困难: 1、基本思想 2、估算m 3、一个实例:8后问题

1、基本思想 当应用回溯法解某一具体问题的具体实例I时,可用蒙特卡罗方法来估算回溯法将要产生的结点数目。 该方法的基本思想是:在解空间树上产生一条随机的路径,然后沿此路径估算满足约束条件的结点总数m: 设x是所产生的随机路径上的一个结点,且位于解空间树的第i层上; 对于x的所有儿子结点用约束函数检测出满足条件的结点数目mi; 路径上下一个结点是从x的这mi个孩子中随机选取的; 这条路径一直延伸到一个叶结点或没有满足条件的子结点为止; 通过这些mi的值,就可以估算出解空间树中满足约束条件的结点总数m。 在回溯法求解问题的所有解时,这个数特别有用,因为在这种情况下,解空间中所有满足条件的结点都必须生成;若只要求用回溯法找出问题的一个解,则所生成的结点数一般只是这m个结点中的一小部分,此时用m来估计回溯法生成的结点数就过于保守。

2、估算m 估算m的公式:

3、实例:8后问题 P159

课后习题 习题 5-14,5-16,5-18,5-27,5-29