数学建模- 图与网络模型及方法

-44-第五章图与网络模型及方法§1概论图论起源于18世纪。第一篇图论论文是瑞士数学家欧拉于1736年发表的“哥尼斯堡的七座桥”。1847年，克希霍夫为了给出电网络方程而引进了“树”的概念。1857年，凯莱在计数烷22nnHC的同分异构物时，也发现了“树”。哈密尔顿于1859年提出“周游世界”游戏，用图论的术语，就是如何找出一个连通图中的生成圈，近几十年来，由于计算机技术和科学的飞速发展，大大地促进了图论研究和应用，图论的理论和方法已经渗透到物理、化学、通讯科学、建筑学、生物遗传学、心理学、经济学、社会学等学科中。图论中所谓的“图”是指某类具体事物和这些事物之间的联系。如果我们用点表示这些具体事物，用连接两点的线段（直的或曲的）表示两个事物的特定的联系，就得到了描述这个“图”的几何形象。图论为任何一个包含了一种二元关系的离散系统提供了一个数学模型，借助于图论的概念、理论和方法，可以对该模型求解。哥尼斯堡七桥问题就是一个典型的例子。在哥尼斯堡有七座桥将普莱格尔河中的两个岛及岛与河岸联结起来问题是要从这四块陆地中的任何一块开始通过每一座桥正好一次，再回到起点。当然可以通过试验去尝试解决这个问题，但该城居民的任何尝试均未成功。欧拉为了解决这个问题，采用了建立数学模型的方法。他将每一块陆地用一个点来代替，将每一座桥用连接相应两点的一条线来代替，从而得到一个有四个“点”，七条“线”的“图”。问题成为从任一点出发一笔画出七条线再回到起点。欧拉考察了一般一笔画的结构特点，给出了一笔画的一个判定法则：这个图是连通的，且每个点都与偶数线相关联，将这个判定法则应用于七桥问题，得到了“不可能走通”的结果，不但彻底解决了这个问题，而且开创了图论研究的先河。图与网络是运筹学（OperationsResearch）中的一个经典和重要的分支，所研究的问题涉及经济管理、工业工程、交通运输、计算机科学与信息技术、通讯与网络技术等诸多领域。下面将要讨论的最短路问题、最大流问题、最小费用流问题和匹配问题等都是图与网络的基本问题。我们首先通过一些例子来了解网络优化问题。例1最短路问题（SPP－shortestpathproblem）一名货柜车司机奉命在最短的时间内将一车货物从甲地运往乙地。从甲地到乙地的公路网纵横交错，因此有多种行车路线，这名司机应选择哪条线路呢？假设货柜车的运行速度是恒定的，那么这一问题相当于需要找到一条从甲地到乙地的最短路。例2公路连接问题某一地区有若干个主要城市，现准备修建高速公路把这些城市连接起来，使得从其中任何一个城市都可以经高速公路直接或间接到达另一个城市。假定已经知道了任意两-45-个城市之间修建高速公路的成本，那么应如何决定在哪些城市间修建高速公路，使得总成本最小？例3指派问题（assignmentproblem）一家公司经理准备安排N名员工去完成N项任务，每人一项。由于各员工的特点不同，不同的员工去完成同一项任务时所获得的回报是不同的。如何分配工作方案可以使总回报最大？例4中国邮递员问题（CPP－chinesepostmanproblem）一名邮递员负责投递某个街区的邮件。如何为他（她）设计一条最短的投递路线（从邮局出发，经过投递区内每条街道至少一次，最后返回邮局）？由于这一问题是我国管梅谷教授1960年首先提出的，所以国际上称之为中国邮递员问题。例5旅行商问题（TSP－travelingsalesmanproblem）一名推销员准备前往若干城市推销产品。如何为他（她）设计一条最短的旅行路线（从驻地出发，经过每个城市恰好一次，最后返回驻地）？这一问题的研究历史十分悠久，通常称之为旅行商问题。例6运输问题（transportationproblem）某种原材料有M个产地，现在需要将原材料从产地运往N个使用这些原材料的工厂。假定M个产地的产量和N家工厂的需要量已知，单位产品从任一产地到任一工厂的运费已知，那么如何安排运输方案可以使总运输成本最低？上述问题有两个共同的特点：一是它们的目的都是从若干可能的安排或方案中寻求某种意义下的最优安排或方案，数学上把这种问题称为最优化或优化（optimization）问题；二是它们都易于用图形的形式直观地描述和表达，数学上把这种与图相关的结构称为网络（network）。与图和网络相关的最优化问题就是网络最优化或称网络优化（netwokoptimization）问题。所以上面例子中介绍的问题都是网络优化问题。由于多数网络优化问题是以网络上的流（flow）为研究的对象，因此网络优化又常常被称为网络流（networkflows）或网络流规划等。下面首先简要介绍图与网络的一些基本概念。§2图与网络的基本概念2.1无向图一个无向图(undirectedgraph)G是由一个非空有限集合)(GV和)(GV中某些元素的无序对集合)(GE构成的二元组，记为))(),((GEGVG。其中},,,{)(21nvvvGV称为图G的顶点集（vertexset）或节点集（nodeset），)(GV中的每一个元素),,2,1(nivi称为该图的一个顶点（vertex）或节点（node）；},,,{)(21meeeGE称为图G的边集（edgeset），)(GE中的每一个元素ke(即)(GV中某两个元素jivv,的无序对)记为),(jikvve或ijjikvvvve),,2,1(mk，被称为该图的一条从iv到jv的边（edge）。当边jikvve时，称jivv,为边ke的端点，并称jv与iv相邻（adjacent）；边ke称为与顶点jivv,关联（incident）。如果某两条边至少有一个公共端点，则称这两条边在图G中相邻。边上赋权的无向图称为赋权无向图或无向网络（undirectednetwork）。我们对图和网络不作严格区分，因为任何图总是可以赋权的。一个图称为有限图，如果它的顶点集和边集都有限。图G的顶点数用符号||V或-46-)(G表示，边数用||E或)(G表示。当讨论的图只有一个时，总是用G来表示这个图。从而在图论符号中我们常略去字母G，例如，分别用,,EV和代替)(),(),(GGEGV和)(G。端点重合为一点的边称为环(loop)。一个图称为简单图(simplegraph)，如果它既没有环也没有两条边连接同一对顶点。2.2有向图定义一个有向图（directedgraph或digraph）G是由一个非空有限集合V和V中某些元素的有序对集合A构成的二元组，记为),(AVG。其中},,,{21nvvvV称为图G的顶点集或节点集，V中的每一个元素),,2,1(nivi称为该图的一个顶点或节点；},,,{21maaaA称为图G的弧集（arcset），A中的每一个元素ka(即V中某两个元素jivv,的有序对)记为),(jikvva或),,2,1(nkvvajik，被称为该图的一条从iv到jv的弧（arc）。当弧jikvva时，称iv为ka的尾（tail），jv为ka的头（head），并称弧ka为iv的出弧（outgoingarc），为jv的入弧（incomingarc）。对应于每个有向图D，可以在相同顶点集上作一个图G，使得对于D的每条弧，G有一条有相同端点的边与之相对应。这个图称为D的基础图。反之，给定任意图G，对于它的每个边，给其端点指定一个顺序，从而确定一条弧，由此得到一个有向图，这样的有向图称为G的一个定向图。以下若未指明“有向图”三字，“图”字皆指无向图。2.3完全图、二分图每一对不同的顶点都有一条边相连的简单图称为完全图(completegraph)。n个顶点的完全图记为nK。若YXGV)(，YX，0||||YX（这里||X表示集合X中的元素个数），X中无相邻顶点对，Y中亦然，则称G为二分图(bipartitegraph)；特别地，若YyXx,，则)(GExy，则称G为完全二分图，记成|||,|YXK。2.4子图图H叫做图G的子图（subgraph），记作GH，如果)()(GVHV，)()(GEHE。若H是G的子图，则G称为H的母图。G的支撑子图（spanningsubgraph，又成生成子图）是指满足)()(GVHV的子图H。2.5顶点的度设)(GVv，G中与v关联的边数（每个环算作两条边）称为v的度(degree)，记作)(vd。若)(vd是奇数，称v是奇顶点(oddpoint)；)(vd是偶数，称v是偶顶点(evenpoint)。关于顶点的度，我们有如下结果：(i)Vvvd2)((ii)任意一个图的奇顶点的个数是偶数。2.6图与网络的数据结构网络优化研究的是网络上的各种优化模型与算法．为了在计算机上实现网络优化的-47-算法，首先我们必须有一种方法（即数据结构）在计算机上来描述图与网络。一般来说，算法的好坏与网络的具体表示方法，以及中间结果的操作方案是有关系的。这里我们介绍计算机上用来描述图与网络的5种常用表示方法：邻接矩阵表示法、关联矩阵表示法、弧表表示法、邻接表表示法和星形表示法。在下面数据结构的讨论中，我们首先假设),(AVG是一个简单有向图，mAnV||,||，并假设V中的顶点用自然数n,,2,1表示或编号，A中的弧用自然数m,,2,1表示或编号。对于有多重边或无向网络的情况，我们只是在讨论完简单有向图的表示方法之后，给出一些说明。（i）邻接矩阵表示法邻接矩阵表示法是将图以邻接矩阵（adjacencymatrix）的形式存储在计算机中。图),(AVG的邻接矩阵是如下定义的：C是一个nn的10矩阵，即nnnnijcC}1,0{)(，.),(,0,),(,1AjiAjicij也就是说，如果两节点之间有一条弧，则邻接矩阵中对应的元素为1；否则为0。可以看出，这种表示法非常简单、直接。但是，在邻接矩阵的所有2n个元素中，只有m个为非零元。如果网络比较稀疏，这种表示法浪费大量的存储空间，从而增加了在网络中查找弧的时间。例7对于右图所示的图，可以用邻接矩阵表示为0110010100000100100000110同样，对于网络中的权，也可以用类似邻接矩阵的nn矩阵表示。只是此时一条弧所对应的元素不再是1，而是相应的权而已。如果网络中每条弧赋有多种权，则可以用多个矩阵表示这些权。（ii）关联矩阵表示法关联矩阵表示法是将图以关联矩阵（incidencematrix）的形式存储在计算机中．图),(AVG的关联矩阵B是如下定义的：B是一个mn的矩阵，即mnmnikbB}1,0,1{)(，.,0,),(,,1,),(,,1其它AijkVjAjikVjbik也就是说，在关联矩阵中，每行对应于图的一个节点，每列对应于图的一条弧。如果一个节点是一条弧的起点，则关联矩阵中对应的元素为1；如果一个节点是一条弧的终点，则关联矩阵中对应的元素为1；如果一个节点与一条弧不关联，则关联矩阵中对应的元素为0。对于简单图，关联矩阵每列只含有两个非零元（一个1，一个1）。可以看出，这种表示法也非常简单、直接。但是，在关联矩阵的所有nm个元素中，只有m2个为非零元。如果网络比较稀疏，这种表示法也会浪费大量的存储空间。但由于-48-关联矩阵有许多特别重要的理论性质，因此它在网络优化中是非常重要的概念。例8对于例7所示的图，如果关联矩阵中每列对应弧的顺序为(1,2)，(1,3)，(2,4)，(3,2)，(4,3)，(4,5)，(5,3)和(5,4)，则关联矩阵表示为1110000010110100010110100000110100000011同样，对于网络中的权，也可以通过对关联矩阵的扩展来表示。例如，如果网络中每条弧有一个权，我们可以把关联矩阵增加一行，把每一条弧所对应的权存储在增加的行中。如果网络中每条弧赋有多个权，我们可以把关联矩阵增加相应的行数，把每一条弧所对应的