2011年数学建模B题国家一等奖

小时

2011高教社杯全国大学生数学建模竞赛
+ \7 n6 r( r8 S9 D+ Y
承  诺  书
* d: |; h( L- [3 g! @7 w
我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
0 h- m' r  O' X$ U+ }- a我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题.
* w* W5 W( b( M/ _2 R" }我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出.
" B$ b1 e  T+ {我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性.如有违反竞赛规则的行为,我们将受到严肃处理.

' `$ h+ }& L6 R" P0 N. Q我们参赛选择的题号是(从A/B/C/D中选择一项填写): B                     
- L1 K2 u6 C1 B( A        我们的参赛报名号为(如果赛区设置报名号的话):                           
3 i- d6 y/ Y* E" _9 @& W+ \8 }所属学校(请填写完整的全名):            
  t/ S- m  [+ t" D+ x3 d参赛队员 (打印并签名) :1.                                       
                       2.                                       
. b" p' K3 D* ]/ J                       3.                                         
指导教师或指导教师组负责人  (打印并签名):               
, w1 _1 _! z$ c' H, C% M; ]
                                              日期:  2011  年  9 月 11 日

$ h& C- v/ m% [8 s# [& \2 v
+ }" f& N* c) ^& b, v
1 |" w& Z# a6 e8 U3 {

" t+ }" M9 E$ d9 l8 r" `, u- C5 m赛区评阅编号(由赛区组委会评阅前进行编号):

7 i$ K/ r" \/ o8 ^( x" ]2011高教社杯全国大学生数学建模竞赛
, w( w4 X# L) X2 o8 A
3 G, [; S6 Z- {; {编 号 专 用 页
4 X$ [* r) b) m: M& n+ O' D- }

  d9 ?+ X7 S7 R) e5 ]9 U
9 G1 d5 {- V: a6 k赛区评阅编号(由赛区组委会评阅前进行编号):


+ i2 T. D4 |# q) E
/ O4 l/ e, j8 \赛区评阅记录(可供赛区评阅时使用):
2 g! l9 ?9 L$ z( F4 G5 H
0 E" J1 s% k8 |/ N; o+ g评
  t6 j" G; ?0 V: v阅
人
; e2 h/ R0 P" A$ B9 N  s                                                                               
  p9 A4 V! U; f0 w0 o0 L

4 \" A- p; w% Q, o" Z评
分
8 `5 T1 u1 o5 S; E3 U+ c                                                                               
" o8 [9 @9 U! {
备
) Y. |" r3 e3 r+ S+ K9 V* x" m9 i注
                                                                               

$ A, I8 O5 m/ q2 v) ~


全国统一编号(由赛区组委会送交全国前编号):
: r! r. J, S, k6 }, u8 I& P
. o. g0 s/ W) @  @. S( Q) z9 @



1 s8 {4 X, j8 l; E) w全国评阅编号(由全国组委会评阅前进行编号):


0 d! Z' |" B! G) Y; M" q+ D
/ q& d% F+ q# c6 p" t
2 ]* c. \( [8 Y' L$ B
& y; J' F/ H+ T( |; z/ _) f交巡警服务平台的设置与调度
8 d" r  E/ ^9 F摘   要
( I7 X8 k* u  b由于警务资源是有限的,所以根据城市的实际情况与需求,合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是有关部门面临的一个实际课题.本文着力于通过所给资料,寻找最优化的交巡台设置与调度方案. 按照设置交巡警服务平台的原则和任务,我们首先对问题1用Floyd算法,提出最佳的交巡警服务平台管辖区域划分方案,缩短了出警时间,平衡了工作量,然后采用回溯法,给出了应对突发事件的警力比较合理调度方案；对于问题2,我们将其归结为全局的配置问题,首先用优化后的Floyd算法对该市现有六城区的交巡警服务平台设置进行改进,其次以时间最短、围堵区域最小为原则,提出了应对重大刑事案件的最佳围堵方案.
对于问题1,本文将最短时间问题转化为单向最短路径问题.我们没有运用经典的求最短距的Dijkstra算法,采取时间复杂度更简便的Floyd算法,应用Matlab编程,以出警时间最短为原则,将72个交通节点分配给20个交巡警服务平台；对于出现突发事件,本文采用回溯法,以最节省警力、实现全区封锁联动时间（即封锁路口最长时间）最短为目标,成功的实现了应对突发事件时警力的合理调度；对于某些交巡警服务平台工作量大、出警时间过长等问题,本文利用Mathematica对附表2中的数据进行分析,整理分析A区各节点事故发生率后,利用图论的相关知识,提出应增设4个服务平台,基本实现警力的最优配置.最后,借助于Matlab和Mathematica软件,对附件中所提供的数据进行了筛选,去除异常数据,对残缺数据进行适当补充,并从中随机抽取了3组数据(每组8个采样)对理论结果进行了数据模拟,结果显示,理论结果与数据模拟结果吻合良好.
而对于问题2,我们对附件中所提供的A,B,C,D,E,F六城区的数据进行了整合与分析,并做出了直观的图表.遵循警情主导警务原则、快速出警原则、方便与安全原则,并结合辖区地域特征、人口分布、交通状况、治安状况和未来城市发展规划等实际情况,在充分考虑现有警力和财力并确保安全的条件下,科学分析现有平台的数量和具体位置的合理性.数据显示C区和F区的事故发生率较高、交巡警服务平台工作量高于全市平均水平、交巡警服务平台平均每天出警时间过长,针对以上问题我们再次利用均衡二分法,并考虑区域边界处的设点拥挤问题,提出了在C区增设5个交巡平台、F区增设1个交巡平台.对于该市地点P(第32个节点)处发生了重大刑事案件的围堵问题,本文将其归结为资源调配问题.本文合理假设了犯罪嫌疑人的车行驶速度(分三种情况考虑:等于警车速度,警车速度的二倍,警车速度的一半),确定三分钟后犯罪嫌疑人逃逸的可能覆盖范围,从而利用回溯法的思想采用Matlab编程确定犯罪嫌疑人的车的所有可能位置.以时间最短、围堵区域最小为原则,采用改进的穷举算法,快速地形成围堵区域,并实现了围堵区域最小的目的.实现了资源调配问题的优化决策.
1 L, {! W' S+ [$ ^考虑到该城市未来发展规划,只需对本文所建模型进行适当改进即可,在此不进行详细解答.
( }3 P3 l$ G1 l

2 Q2 w0 |2 ^, ^8 w) ]+ J; n关键词    最短路径  Floyd算法  回溯法  穷举法  优化决策

8 R* z4 h" r! N& m; l
3 Q0 b% ], g/ ]9 w/ p
" k  }8 k! H. E- S8 e7 ]目   录
: r6 r" v! Z2 f) }9 o* N交巡警服务平台的设置与调度        1
0 N2 ^) u  t" \; q3 j" L( |" V5 X摘   要        1
- y  Y" j% w0 b: J1.问题重述        1
* W5 j/ e* i( X) v" U* ^2.问题分析        1
$ P3 q' S. B. z2.1对于问题一的分析        1
/ O8 ?  K8 c& a) v$ h2.2对问题二的分析        1
3.模型假设        2
4.定义与符号说明        2
5.模型的建立与求解        2
5.1 问题一的模型        2
0 ]8 o7 H: Y; j. c( s5.1.1 模型建立        2
7 J. a! Y! `! I/ j7 C  |5.1.2 模型求解        3
: _9 W% J9 k) u5.2 问题二的模型        8
5.2.1 模型建立        8
5.2.2 模型求解        9
) X6 u0 s8 j; c8 g) h" q7.模型的评价与推广        10
8. 附件        10
0 \$ [% |* D+ {$ \) m附件1:用Floyd算法分配个服务平台管辖区域        10
5 i! r3 s' d; C. P) ^) N( Z# g! P附件2:邻接矩阵的matlab实现程序        22
附件3:围堵方案的java实现程序        29
) `  N; b8 J( D; B% e& V附件4:全区的交巡警平台有效覆盖范围(有效代表三分钟内可以到达)        30
5 Q, g9 q' F* x' J附件5:用Mathmatica求数据均值与方差        30
附件6:输入任意两点的坐标,输出两点间距离        30
附件7:A区各线路距离        31


4 b# z% Y/ m* X( r


8 ^6 ?" e" h6 A; O( B: s
1 v- b0 W6 m3 }
1 q) z, ]$ o, ~& `( e2 }1.问题重述
6 L! f6 l1 H9 v0 A2 G; r5 ~“有困难找警察”,是家喻户晓的一句流行语.警察肩负着刑事执法、治安管理、交通管理、服务群众四大职责.为了更有效地贯彻实施这些职能,需要在市区的一些交通要道、人员密集区和重要部位设置交巡警服务平台.每个交巡警服务平台的职能和警力配备基本相同.由于警务资源的有限性,根据城市的实际情况与需求,合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题.本文着力于寻找最优化的设置与调度方案.
3 G- o( ^; t, `/ v5 ~6 o* n问题1要求合理分配交巡警服务平台的管辖范围,使其在所管辖的范围内出现突发事件时,尽量能在3分钟内有交巡警(警车的时速为60km/h)到达事发地；对于重大突发事件,给出该区交巡警服务平台警力合理的调度方案,尽快封锁道路；拟在该区内再增加2至5个平台,以减少出警时间、平均工作量,确定需要增加平台的具体个数和位置.
问题2要求分析研究该市现有交巡警服务平台设置方案的合理性并给出解决方案；如果该市地点P(第32个节点)处发生了重大刑事案件,在案发3分钟后接到报警,犯罪嫌疑人已驾车逃跑.为了快速搜捕嫌疑犯,给出调度全市交巡警服务平台警力资源的最佳围堵方案.
2.问题分析
本题所要解决的是A区以及全市的安巡警服务平台设置与调度问题,根据现实生活状况,我们首先要考虑的是警力资源的限制,即要使得所布置的警力尽可能的少.其次是在交巡台数量最少的情况下,力求警员到达现场的时间在3分钟以内,解决突发状况.
2.1对于问题一的分析
该市中心城区A的交通网络有92个节点和20个交巡警服务平台,要求当突发事件发生时,尽量能在3分钟内有交巡警到达事发地,已知警车的时速为V=60km/h,我们将最短时间转化为最短路问题,应用Floyd算法,求解出A区距离每一节点最近的交巡台,即将该节点分配给该交巡台.
对于重大突发事件,要实现对进出该区的13条交通要道进行快速封锁,即需调度交巡台尽快到达13个节点,重复Floyd算法,找出最近交巡台,即可找出调配方案.但需注意的是,有的出入口本来就有交巡台,但为了达最优化,需进行重新分配,故应用回溯法,找到调度方案.
* a+ c4 b6 _6 _' B8 }0 @现有交巡台工作量不均衡和有些地方出警时间过长,统计A区各个交巡台案发率,计算均值与方差,在案发率较高地带增设交巡台,平衡工作量,尽量缩短出警时间.
2.2对问题二的分析
! H' m, Y& _& d对于问题二,是对问题一的进一步改进与推广,在遵循警情主导警务原则,快速出警原则与方便与安全原则,结合辖区地域特征、人口分布和治安状况等实际情况,充分考虑现有警力和财力并确保安全的条件下,设置交巡平台,重复上一问的做法,评估交巡平台的合理性.对于改进方案,应考虑城区内部工作量,城区之间的联系以及城市边界的警力调度.
2 F- f9 ^$ O' B3 |* u对于突发状况的围堵方案,应在最短时间内对可能逃逸区域进行合围,最小范围内缩小包围圈.
3.模型假设
/ U1 X1 ]* Y% q' l" q- @1.假设题中所给数据均真实可靠.
2.出警时道路恒畅通(无交通事故、交通堵塞等发生),警车行驶正常,警车及肇事车辆行驶时均以60km/h匀速行驶,转弯处不需要花费时间.
9 y1 u2 G  Y) _' w4 k8 P) \6 h3.事故均发生在路口节点,两节点连线上认为没有事故发生.
" y2 k) O' w5 U8 m9 T1 K9 d4.每条线路行驶都是双向的.
) `: U+ O0 p6 s7 B, R5.考虑肇事车辆在P点向各个方向逃逸的概率相等.
6.在整个行驶中,车辆只在主要干道行驶.
; D3 c1 U1 a3 {0 X9 m+ Z7.发生事故时,忽略反应调度时间.
3 {' ]  H5 F0 A* a4 T: |4.定义与符号说明           
$ u0 I" C. U& Y5 X4 y+ ~6 w
任意两个标志点 与 之间的距离


0 P' c1 c/ k3 g( u5 r& D1 a标志点间的距离组成的距离矩阵
2 z1 a0 M  ]- C9 H+ @3 h+ h
# V. K+ z7 ^, L3 ~: l标志点的邻接矩阵
. Y* @; K) `0 T  L) X
3 v/ x1 C2 {$ M; c# H, W邻接矩阵的元素

+ C% }$ U* R; L) m2 |5 y. |9 P相邻标志点间的距离矩阵
8 @! z8 a* Y) y9 Q# ]
相邻标志点 与 间的距离

* P* X0 O3 z( q, ~; ]& I  J) }0 X6 {
% v# H2 h5 O- t. Q9 ~标志点的权值矩阵
" A' q! v/ G" A# a6 L) W  
/ J$ L% y0 D4 J1 Z" A5 Q标志点间的最短距离矩阵

标志点 与 之间的最短距离
: s; q% x! p" |4 X1 M
6 N) ^* g0 J- F8 S% v, J3 u7 Q
1 B- n& v4 O9 f( @* ~2 @肇事车辆逃逸速度

) }6 z. U1 n8 Z# t, J( S8 V' S/ ^: _  
. Z/ ^/ U! H5 s) Y5 G5.模型的建立与求解
5.1 问题一的模型
5 @; g' B! }$ I  A' z: G% k5.1.1 模型建立
此问是关于最短路径的模型分析及MATLAB的实现A区道路状况及交巡台的设置如图1所示.本文应用Floyd算法,通过构造距离矩阵,依次找出距离每一节点最近的交巡台,使得有事故发生时,交巡警在最短时间内到达事故现场,以此为依据分配管辖区域.如果道路不通时,认为两端节点的距离为无穷.
                   图1  A区各节点及服务平台示意图
  E* R. `: R3 f, Q7 x% n2 w当有重大突发事件时,要对进出该区的13条交通要道进行快速封锁,固定13个出入口,应用回溯法,找到距离节点最近的交巡平台.封锁时间决定于最后到达节点的时间,由于一个平台的警力最多封锁一个路口,至少需调动13个平台的警力.
为达到工作量的均衡和出警时间尽可能的短,需进行优化决策.考虑每一节点案发率的不同,在A区增设2到5个平台,使得每一平台的工作量均衡,平均出警时间大体相同.
9 }* i% s- h9 ?5 y$ m) B! h5.1.2 模型求解
3 u# w; _# v; L) ]首先我们可以根据题中所给的各个标志点的坐标,用matlab计算出任意两点之间的直线距离,得到92*92的距离矩阵:
6 a" R2 b& p9 f5 y( C7 d' T  _
根据题中的分布图,我们可以得到各标志点的邻接矩

" G6 ^, l% R/ J6 w# c5 `即如果两个点相邻,则邻接矩阵中相对应的元素的值为1,否则为0；例如:3和44这两个点相邻,那么  .
$ Z9 m( w: ^4 W1 ]# H6 O: X    根据Floyd算法,我们是要求出任意两节点之间的距离,所以我们需要得到相邻两个结点的直线距离.我们可以利用距离矩阵的元素 与 的点乘积得到相邻标志点间的距离矩阵:

对于D中不相邻点间距离0改为无穷大(Inf)从而得到节点与节点间的权值矩阵:
: ^9 n5 R  Q5 i9 L8 `
* ~; Y% z  I$ J, c3 s3 @1 P  ?即如果15和10之间不相邻,也即不能直接到达,那么D中的 和 都将变成 和 等于无穷大(Inf),否则则等于D中相应元素的数据.
7 j' {0 x8 s2 }运用Floyd算法求出任意两点间最短距离,得到最短距离矩阵 :
& J  U/ J8 K2 k$ P$ @" Z
1 w' J8 {; m1 i2 G, p由Floyd算法,运行MATLAB程序,可统计出距离每一节点最近的交巡台的位置,MATLAB运行结果如表1所示.带括号的节点为发生事故时任意交巡台都不能在三分钟内赶到节点.
交巡台—节点        距离        交巡台—节点        距离
13—21        27.0831         4—57        18.6815
& b% {7 Z' m" n13—22        9.0554        6—58        23.8414
/ M% v/ [, w) t4 W1 R! q13—23        5.0000         6—59        16.0312
) S, v4 A! Y- F) K) Y$ ~! t2 C13—24        23.8537         4—60        17.9240
12—25        17.8885         4—(61)        52.1055
3 h- q1 @7 D0 O; _' H11—26        9.0000         4—62        3.5000
. S3 b6 K. g" V. R; S3 P11—27        16.4330         4—63        10.3087
9 x$ r4 T# d& M  N& F: L15—(28)        47.5184         4—64        9.3632
4 `# H! y: Y3 P. P, U0 j/ a& [15—(29)        57.0052        3—65        15.2398
7—30        5.8310         3—66        18.4012
6 u4 i, f/ B8 f. U$ ]; o9—31        20.5572        1—67        14.9158
7—32        11.4018        1—75—68        10.7927
8—33        8.2765        1—69        5.0000
9—34        5.0249         2—70        8.6023
' }0 Y2 \7 h0 x& i/ v; o9—35        4.2426         1—74—71        11.2650
2 v; T. r7 W- _16—36        6.0828         2—72        16.4031
16—37        11.1818         18—73        19.7231
16—(38)        34.0588         1—74        6.2650
* n1 N6 D4 {1 h1 |7 \9 @5 h2—(39)        36.8219        1—75        6.2650
: I* s& V& w0 P- Y" q3 r2—40        19.1442         1—76        9.8005
; d6 e: T: f  p: T7 g, ]: u17—41        8.5000         19—77        9.8489
17—42        9.8489        1—78        6.4031
/ d' V, P, S/ }! e+ c2—43        8.0000        19—79        4.4721
2—44        9.8468        18—80        8.0623
9—45        10.9508         18—81        6.7082
8—46        9.3005        18—82        10.7935
7—47        12.8062        18—83        5.3852
7—48        12.9021        20—84        11.7522
5—49        5.0000        20—85        4.4721
5—50        8.4853        20—86        3.6050
5 P& o# T8 o  x3 u: i! _& ^5—51        12.8932        20—87        14.6511
5—52        17.1944        20—88        12.9464
5—53        11.7082        20—89        14.7522
3—54        22.7089        18—90        19.5256
; @+ W7 u5 _+ R3—55        12.6590        20—91        16.0060
5—56        21.4370        20—(92)        36.0060
# l% i# q% x$ |表1  该市A区指定节点到交巡警服务平台最短距离
由上表可初步确定A区20个交巡台的管辖范围,如表2所示.带括号的节点为发生事故时任意交巡台都不能在三分钟内赶到节点.
交巡台序号        辖区内节点        辖区内案发率        交巡台序号        辖区内节点        辖区内案发率
2 L& X; S8 i% [6 t3 i& L1        67 68 69 71 74 75 76 78         9.4        2         40 43 44 70 72 39        9.7
3        54 55 65 66        5.6        4        57 60 62 63 64        6.6
5        49 50 51 52 53 56        7.7        6        58 59        4.5
7        30 32 47 48 61        9        8        33 46        5
9        31 34 35 45        8.2        10                 1.6
11        26 27        4.6        12         25        4
13        21 22 23 24        8.5        14                 2.5
$ A0 T$ j+ @3 C, w' i5 m2 w15        (28) (29)        4.8        16        36 37 (38)         5
17        41 42        5.3        18        73 80 81 82 83        7
! p- T0 i/ N4 p0 E' J4 `19        77 79        3.4        20        84 85 86 87 88 89 91 90 (92)        11.5
# C# t5 j9 `( L0 n& s表2  该市A区交巡警服务平台所管辖交叉路口清单


4 R7 ~2 n2 @" r0 ^图2  A区各交巡台管辖区域示意图
    需要说明的是,同一条路整体归一个交巡台管理.
当有重大突发事件时,固定13个进出A区的节点,运用回溯法,结合上表,找到距离节点最近的交巡台,以此来达到总体时间的最短,我们一共可以得到四个方案,在这个过程中可以发现,有些交巡台要避免去最近的节点封锁而去较远的节点,以此来节省警力.具体封锁方案如表3、表4所示.最短调度时间均为8.0155.
方案一:
交巡台        过程        出入口节点号
2        路径        40 39
38
        时间        3.9822min       
/ L- b( T5 @0 z4 `5 z4        路径                62
. @3 ]/ M! d- n5 Y2 a        时间        0.3500min       
, j) c- T3 i3 B# ]* f6        路径        47 48        30
# S; w7 J* Y7 @( u* h, R( A        时间        3.1829min       
& |) Z7 T3 N# P7        路径        30        29
3 J7 ]' i6 E( t& `  U3 y        时间        8.0154min       
. {/ c# W- f$ f' D6 x7 P8        路径        47        48
" E" D3 Q7 W  j        时间        3.0995min       
% i5 [! v' d6 q, ^' `) p$ e' V9        路径        35 36        16
/ R4 O0 k; h# |$ [( U" o        时间        1.5083min       
10        路径        26 27        12
4 q8 ~  M0 S" e2 @2 @        时间        7.5863min       
* h$ V' z8 M8 h' u11        路径                22
) ]5 D$ T5 J6 h- c$ @4 ~5 A3 e+ r/ w        时间        3.2696min       
  G2 L% W( P7 A12        路径        25        24
5 z! I; c/ {  W( y7 H" S& Z# H        时间        3.5916min       
13        路径                23
        时间        0.5000min       
14        路径                21
        时间        3.2649min       
15        路径                28
6 h8 I3 U: ]8 I6 l. J# g        时间        4.7518min       
16        路径                14
        时间        6.7417min       
表3  A区突发事件封锁方案一
方案二        方案三        方案四
路口标号        平台号        路口标号        平台号        路口标号        平台号
12        13        12        13        12        10
14        16        14        23        14        16
16        6        16        9        16        6
9 T' a5 W1 y+ ~! r, Q8 x+ d21        14        21        11        21        14
22        10        22        10        22        12
23        11        23        14        23        13
, h/ {. ^# b$ m- u* x$ ?24        12        24        12        24        11
' H' ^% C3 H+ K% w2 C2 P28        15        28        15        28        15
% p* y+ b5 U' k# f2 t2 ]' F29        7        29        7        29        7
30        8        30        6        30        9
2 g, f! l% V- ~8 B2 v* d: C1 W38        19        38        17        38        1
8 [4 {, A4 K! v: A. ]5 ~8 I48        5        48        9        48        8
62        20        62        20        62        2
表4  A区突发事件封锁方案二、三、四
% W: [0 z7 U" J) B  ^5 b. Z在对交巡台均衡工作量,加快出警时间方面,综合各节点的案发率、交巡台到其辖区内任一节点的路程进行综合评估,做出优化决策.
# B, d. N) U$ \/ p4 m% ~7 Z在案发率较高地带增设交巡台,以缓解周围交巡台的工作压力,为达均衡工作量的目的,将32号节点从7号交巡台归到8号交巡台,44号节点从2号交巡台归到3号交巡台,39号节点从2号交巡台归到16号交巡台,47号节点从6号交巡台归到7号交巡台,61号节点从7号交巡台归到4号交巡台.这样,A区每交巡台平均每天处理案件数从6.1950件、方差6.8289降到每天处理5.1917件、方差2.2182,极大的协调了工作量.对于个别节点的重新划分,会增加出警时间,但在总体上平均每天的出警时间大大缩短了.
0 ]$ f: E# P$ S: q1 G5 x& B0 p综上考虑,共增设4个交巡台,重新分配的结果如表4所示.
交巡台序号        负责区域内的节点        管辖区域内的案发率        平均每天出警时间
1        69 71 74 75 78        6.6        3.5069
2        40 43 70 72         7.2        6.6736
3 O2 y( s2 ?7 H9 Y4 x3        54 55 44        5.2        4.9314
4        57 60 61 62 63        6.4        7.7015
. H+ Y" n$ @  p2 `$ m) X5        49 50 51 52 53 56        7.7        5.9456
6        58 59 47        6.1        6.3949
7        30 48        5.9        3.5258
8        33 46           5.0             2.2748
9        35 45        4.9        2.1271
% k9 U( m$ S: T4 k! l1 _& Q& }! r10                1.6       
  U* M' E4 H8 S: t11        26 27        4.6        2.3946
6 |: |& J+ U4 f; ], t12        25        4.0        2.8622
13        23 24        5.7        3.8239
14                2.5       
( M4 C7 Y- y, a9 A15        (28) (29)        4.8        14.1580
16        36 37 (38) (39)        6.4        10.0562
0 |. n' Z) T% _+ l9 \9 `17        41 42        5.3        2.5689
18        73 80 81 83        5.9        3.8438
19        77 79        3.4        1.1457
20        85 86 87 (92)        6.4        5.5354
: ]  W. b' h  x, m% n! U0 n21        22        2.8        2.5239
8 g9 y1 t: W1 _" D! K% o! G31        32 34        4.9        4.3962
2 @9 n: h7 n# B6 a# T% m$ f66        64 65 67 68 76         5.1        2.6655
6 V7 d* |' A0 P0 U, l$ M90        82 84 88 89 91        6.2        3.2171
表5  优化后的A区交巡台管辖区域示意图

" I$ s$ o& q% m4 V图3  A区增设平台示意图
    图中方块所示节点即为增设平台处.
5.2 问题二的模型
5.2.1 模型建立
* u3 U- E. p( J5 f对于问题2,对附件中所提供的A,B,C,D,E,F六城区的数据进行整合,做出直观的图表.遵循警情主导警务原则、快速出警原则、方便与安全原则,结合辖区地域特征、人口分布、交通状况、治安状况和未来城市发展规划等实际情况,充分考虑现有警力和财力并确保安全,科学分析现有平台的数量和具体位置的合理性.
$ i( k; Q) y6 L对于该市地点P(第32个节点)处发生了重大刑事案件的围堵问题,本文将其归结为资源调配问题.本文合理假设了犯罪嫌疑人的车行驶速度(分三种情况考虑:等于警车速度,警车速度的二倍,警车速度的一半),并确定三分钟后犯罪嫌疑人的车行驶的最远距离,从而利用回溯法的思想采用Matlab编程确定犯罪嫌疑人的车的所有可能位置.以时间最短、围堵区域最小为原则,采用改进的双层Floyd算法,快速地形成围堵区域,并使围堵区域尽可能的小.
4 ]5 U9 D5 w% [4 j5.2.2 模型求解
全市整体状况如表5所示,数据显示C区和F区的事故发生率较高、交巡警服务平台工作量高于全市平均水平且交巡警服务平台平均每天出警时间过长,针对以上问题本文再次利用问题1的Floyd算法,并考虑区域边界处的设点拥挤问题,本文提出了在C区增加5个服务平台、在F区增加1个服务平台.
全市六个城区        城区面积        城区人口        平台数        平均人口        全区案发率        各区平台案发率均值
" ]9 C* r& _! YA        22        60        20        2.727        124.5        6.625
; }: ]. c4 `9 q6 Z+ ?. CB        103        21        8        0.204        66.4        8.3
3 s( b3 k- o/ V" a) nC        221        49        17        0.223        187.2        11.012
D        383        73        9        0.191        67.8        7.533
E        432        76        15        0.176        119.4        7.96
+ d  q; X" I5 F* cF        274        53        11        0.193        109.2        9.927
均值                53.3333        13.3333        0.619
& m5 I) a4 L- t* H  I% h9 ^112.4167
8.5595
. Z' I1 F. p% ]表6  全市整体状况

图4  全市增设交巡台位置示意图(方块所示区域)
& H2 b7 d0 ], N: ~6 e" z! \对于P点发生重大刑事案件,动用全市警力进行围堵,我们希望使得包围圈尽可能的小,由于犯罪嫌疑人的车速度未知,我们分以下三种情况进行考虑:
1)当犯罪嫌疑人的车速与警车速度同,即 .
6 s6 U' K+ W& `; w1 B1 P0 e运用穷举法,对肇事车辆可能的逃逸路线进行分析,以3分钟路程为半径,找到肇事车辆逃逸的覆盖范围,如图5所示,其中实线表示可能路径,在此范围内有8、9、10、15号共4个交巡平台,保证这4个平台警力不动,组成第一范围包围圈.

) s/ |5 E! X* F% A5 D: d图5 肇事车逃逸3分钟内覆盖区域示意图
. i% d! I: Y0 ~* c1 H. u    进一步分析可能的逃逸路线,调度16号交巡台到36号节点,2号交巡台到3号节点,3号交巡台到55号节点,6号交巡台到47号节点,组成第二组半包围,保证对A区的封锁.若肇事车辆经36号节点逃往16号节点,则会与16号交巡台在途中相遇.
对于从32号节点经7号节点逃逸到30号和47号节点,存在从A区逃往其他城区的可能,需调动其他城区交巡台的支援.将C区119号交巡台调度到237号节点,将D区320号交巡台调度到371号节点,321号交巡台经368号、369号节点到370号节点,至此,在全市范围内实现全面封锁.
2)当犯罪嫌疑人的车速比警车车速小,即 ,我们令
8 W) g8 V8 E5 _9 |" X    方法同1),寻找分钟逃逸范围内所覆盖的全部交巡台,经过整合分析,保持7号、8号、9号、15号共4个交巡台原地封锁,10号交巡台到34号节点封锁,6号交巡台到47号节点封锁,16号交巡台到36号节点封锁,3号交巡台经55号节点到46号节点进行封锁,2号交巡台经3号节点到45号节点封锁,在此过程中,10号、2号和3号交巡台会在途中与肇事车辆相遇.
3)当犯罪嫌疑人的车速比警车车速大,即 ,我们令
由于肇事车辆逃逸速度较快,可能会逃逸到C区和F区,故需调动C区和F区警力进行围堵.
   A区将20号交巡台调到62号节点,16号交巡台调到36号节点,2号交巡台经40号到39号节点,17号交巡台调到41号节点,15号、10号、4号、3号、5号、7号、8号、9号交巡台原地封锁,其余交巡台向其邻近的路口节点进行增援.经过分析,肇事车辆可能由28号、48号、30号进入C区及A、D两区的交汇地带,或由16号节点逃逸到F区,在此,对C区、D区、F区交巡台进行如下调配,实现全市封锁:
   C区:240号交巡台调到239号节点,170号交巡台调到225节点,167号交巡台调到259节点.
   D区:320号交巡台调度到371号节点,321号交巡台经368号、369号节点到370号节点.
" K! M  o7 i- ]   F区:477号交巡台调度到501号节点,518号交巡台调到521号节点,478号节点调到527号节点,484号节点到571号节点.
7.模型的评价与推广
本文避免了时间复杂度较复杂的Dijkstra算法,选用Floyd算法,在求最短路径上提高了效率,代码编写简单.
) C. R0 y) H0 q& \' M模型的建立思路清晰,遵循可操作性、科学性、可比性原则,该模型建立出了在较理想状态下交巡警平台的最优设置,减少出警时间,均衡工作量,提高工作效率,在遇突发事件时,可尽快实现道路封锁,给生活中交巡警平台的设立予参考,具有一定的实际应用价值,也可以应用于其他适用区域.模型的运算由矩阵、向量的运算组成,易于用数学软件求解和验证.
" g! ]4 r. {' A/ p, x本模型较好的解决了交巡警平台的最优选址问题,当事故发生时,交巡警可以第一时间到达事发地点,有效的改善了交巡警在执行任务中的效率,在经济迅猛发展的今天,城市加速扩张,人口迅速增长,交巡警平台的设置是平安城市的最好保障.该模型也可运用到其他最优选址问题中去,比如关于消防救援工作最优路径问题、重大生产安全事故应急救援问题、公共交通的最优路径问题等. 同时也可利用该模型算法拓展模型在其他领域的适用范围.
该模型也有一定的局限性,如现实中不能时刻都保证道路的畅通性.既不能保证出警的时间总是维持在3分钟之内.忽略了实际地形对于车速的影响以及实际生活中存在的不定因素.
! |/ A% q9 `5 Y' i* g: q2 M
) }+ b) ~) [) K, h参考文献
& \' e) A2 }) i* @9 ][1]         [徐孝凯,王凤禄],《数据结构简明教程》第二版,北京:清华大学出版社,2005年4月1日
[2]         [李建中,骆吉洲],《华章数学译丛》第二版,北京:机械工业出版社,2002年6月
[3]         [陈庆华等],《组合最优化技术及其应用》第1版,北京:国防科技大学出版社,1989年8月
[4]         [W.T.Tutte],《Graph Theory》,英国:Cambridge University Press,2001年3月1日
! [3 }$ O9 x: j+ U1 U% x! v5 M8.附件
附件1:用Floyd算法分配个服务平台管辖区域
/ n1 L8 w# m8 _1 c5 B9 |[area_a _x,area_a _y]=find (location_all _daolu<=92);
* e/ R7 n# P3 H0 ` road_index _a=[area_a _x,area_a _y];
road_index _a
' ^0 L; ]7 P: V) e: n# H+ s# n$ ]1 B
1 p* F! `8 G4 G( Y, F& w" xa1=find (road_index _a (:,2)==1);
a2=find (road_index _a (:,2)==2);
% L4 r' j3 r' {; mA=road_index _a (a1,1);B=road_index _a (a2,1);[c,ia,ib] = intersect (A,B);
c
$ Z7 x( `9 q& o. O' I
" F2 u9 r$ h; G% A' I. dsize (c)
9 I: _3 D. s) [! e
for i=1:140
    tt=c (i);
    uu=location_all _daolu (tt,;
0 W5 H8 w* D/ s! b7 |    uu1=uu (1);uu2=uu (2);
    vv1=location_a _zuobiao (uu1,;
    vv2=location_a _zuobiao (uu2,;
    ww1=[vv1 (1),vv2 (1)];
    ww2=[vv1 (2),vv2 (2)];
+ V* [7 X6 O, s0 l: D    line (ww1,ww2)
end
) x( K; G9 X9 K6 v for i=1:140
    tt=c (i);
    uu=location_all _daolu (tt,;
& P/ U* Q) O9 o% ?( W5 y2 I1 K    uu1=uu (1);uu2=uu (2);
4 L( U+ |  v2 G    vv1=location_a _zuobiao (uu1,;
    vv2=location_a _zuobiao (uu2,;
    ww1=[vv1 (1),vv2 (1)];
3 _4 p% j0 d1 w    ww2=[vv1 (2),vv2 (2)];
    line (ww1,ww2,'k')
( D  \$ F9 [- q: `4 \4 C4 o1 r2 {end
! C# `+ p1 G5 J3 Y0 O. B5 t% ??? Error using ==> line
% String argument is an unknown option.
2 o' M5 W' z: h# H& `' l
8 ?1 Y6 H) i& I( F3 qfor i=1:140
2 d; ^7 B- }) w+ G+ c- P7 E& m    tt=c (i);
    uu=location_all _daolu (tt,;
    uu1=uu (1);uu2=uu (2);
2 f0 j5 l; W7 A4 T    vv1=location_a _zuobiao (uu1,;
    vv2=location_a _zuobiao (uu2,;
; v5 r+ v) F! Q: X5 o" q, @    ww1=[vv1 (1),vv2 (1)];
    ww2=[vv1 (2),vv2 (2)];
    line (ww1,ww2,'Color',[.8 .8 .8])
end
save data_b _problem
2 u9 w$ D6 X. B# H) t; `
- M4 i- t% y' U/ p# F8 ^" |% A区节点间的邻接矩阵
load data_b _problem;
matric_lingjie=zeros (92,92);
/ y* |, v. g& \+ b% |# w  ]) P[xx,yy]=find (matric_lingjie==0);matric_lingjie (xx,yy)=inf;
for i=1:92
5 e' s% R6 {; n' p& X7 ]        matric_lingjie (i,i)=0;
7 w  V5 u9 g8 _9 ]end

* R( e6 ?  p) V" s7 A5 L. T! }% ^
for i=1:140
    tt=c (i);
    uu=location_all _daolu (tt,;
  D' T7 M- ]  s; z$ o( T" e    uu1=uu (1);uu2=uu (2);                       % 端点序号
    vv1=location_a _zuobiao (uu1,;           % 第一个端点坐标
6 `8 g& K& `: J: `    vv2=location_a _zuobiao (uu2,;           % 第二个端点坐标
+ F& M7 ?2 M" T" k% g    % 计算端点间距离
( T" l, {! h/ x6 j: x    distance=sqrt ((vv1 (1)-vv2 (1))^2+(vv1 (2)-vv2 (2))^2);
* s) n# w% N% o6 s6 T4 m' ?( ]0 {  [    matric_lingjie (uu1,uu2)=distance;
    matric_lingjie (uu2,uu1)=distance;        % 赋值给邻接矩阵
1 [- V) [4 U9 F1 @6 z+ m9 t4 g' Hend

! s" i3 Y& ?/ R[D,R]=floyd (matric_lingjie);
4 P, T8 s  k" P! C
) C- a; Z  X6 ~' ]4 ^8 omatric_fenkuai=D (1:20,;
: p& i# t2 l* ^* L
) b, I( ^' I7 y1 V, ?+ ]3 o" V" N- Wfor i=1:20
    for j=1:92
9 {5 Q9 |. U1 S& f) j        if matric_fenkuai (i,j)>30
: }: \4 l- H" R; S            matric_fenkuai (i,j)=0;
# ]5 W: R/ k9 c) [& @        end
    end
end

' S0 @0 V7 {/ ?) X' {6 R( o8 sti=zeros (1,92);
5 Y: }3 K& I- E% F, z" fti (1)=text (location_a _zuobiao (1,1),location_a _zuobiao (1,2)+1.5,'1');
ti (2)=text (location_a _zuobiao (2,1),location_a _zuobiao (2,2)+1.5,'2');
ti (3)=text (location_a _zuobiao (3,1),location_a _zuobiao (3,2)+1.5,'3');
ti (4)=text (location_a _zuobiao (4,1),location_a _zuobiao (4,2)+1.5,'4');
$ y, W. ~9 H# q8 s' Oti (5)=text (location_a _zuobiao (5,1),location_a _zuobiao (5,2)+1.5,'5');
" `* b1 _6 X8 Z# K- Ati (6)=text (location_a _zuobiao (6,1),location_a _zuobiao (6,2)+1.5,'6');
ti (7)=text (location_a _zuobiao (7,1),location_a _zuobiao (7,2)+1.5,'7');
ti (8)=text (location_a _zuobiao (8,1),location_a _zuobiao (8,2)+1.5,'8');
4 r5 C  V, p( ?! P) U/ ]ti (9)=text (location_a _zuobiao (9,1),location_a _zuobiao (9,2)+1.5,'9');
2 O3 _& T: o# X* [, Nti (10)=text (location_a _zuobiao (10,1),location_a _zuobiao (10,2)+1.5,'10');
ti (11)=text (location_a _zuobiao (11,1),location_a _zuobiao (11,2)+1.5,'11');
5 y8 B! u4 {/ W3 a5 T1 C- k. jti (12)=text (location_a _zuobiao (12,1),location_a _zuobiao (12,2)+1.5,'12');
ti (13)=text (location_a _zuobiao (13,1),location_a _zuobiao (13,2)+1.5,'13');
, a6 t& O1 b" Vti (14)=text (location_a _zuobiao (14,1),location_a _zuobiao (14,2)+1.5,'14');
7 o* |4 t1 q: B0 j" W. ]9 s  u. lti (15)=text (location_a _zuobiao (15,1),location_a _zuobiao (15,2)+1.5,'15');
8 D* u8 Z. G5 f1 U5 ]* Qti (16)=text (location_a _zuobiao (16,1),location_a _zuobiao (16,2)+1.5,'16');
% U+ v+ B2 Q7 Y% }ti (17)=text (location_a _zuobiao (17,1),location_a _zuobiao (17,2)+1.5,'17');
# M9 `! U% ~& c2 A; z; y7 Qti (18)=text (location_a _zuobiao (18,1),location_a _zuobiao (18,2)+1.5,'18');
  f. W4 v5 [, l6 A5 t; Wti (19)=text (location_a _zuobiao (19,1),location_a _zuobiao (19,2)+1.5,'19');
ti (20)=text (location_a _zuobiao (20,1),location_a _zuobiao (20,2)+1.5,'20');
* c2 f9 v; G. F0 [4 f% [5 Dti (21)=text (location_a _zuobiao (21,1),location_a _zuobiao (21,2)+1.5,'21');
6 C" ~8 W6 e6 E+ G2 j% V2 Nti (22)=text (location_a _zuobiao (22,1),location_a _zuobiao (22,2)+1.5,'22');
ti (23)=text (location_a _zuobiao (23,1),location_a _zuobiao (23,2)+1.5,'23');
ti (24)=text (location_a _zuobiao (24,1),location_a _zuobiao (24,2)+1.5,'24');
ti (25)=text (location_a _zuobiao (25,1),location_a _zuobiao (25,2)+1.5,'25');
5 O' ], X9 c+ ~: d' Xti (26)=text (location_a _zuobiao (26,1),location_a _zuobiao (26,2)+1.5,'26');
1 K0 @' X) Q- |$ B) ]6 Cti (27)=text (location_a _zuobiao (27,1),location_a _zuobiao (27,2)+1.5,'27');
ti (28)=text (location_a _zuobiao (28,1),location_a _zuobiao (28,2)+1.5,'28');
/ t  w! ~2 U7 D1 @& oti (29)=text (location_a _zuobiao (29,1),location_a _zuobiao (29,2)+1.5,'29');
ti (30)=text (location_a _zuobiao (30,1),location_a _zuobiao (30,2)+1.5,'30');
5 b1 Z: W" `4 @' q0 i' c  Dti (31)=text (location_a _zuobiao (31,1),location_a _zuobiao (31,2)+1.5,'31');
ti (32)=text (location_a _zuobiao (32,1),location_a _zuobiao (32,2)+1.5,'32');
: o3 c* R" a: A$ d) lti (33)=text (location_a _zuobiao (33,1),location_a _zuobiao (33,2)+1.5,'33');
) A8 _2 v* [3 u$ a3 @) Yti (34)=text (location_a _zuobiao (34,1),location_a _zuobiao (34,2)+1.5,'34');
+ b- ^/ Q1 n$ n% a0 d( \4 g9 tti (35)=text (location_a _zuobiao (35,1),location_a _zuobiao (35,2)+1.5,'35');
  t6 X: a5 g5 ^1 {6 E( S" Ati (36)=text (location_a _zuobiao (36,1),location_a _zuobiao (36,2)+1.5,'36');
1 I4 m4 e- w1 k- l% I  h6 K2 t0 j$ G) ?ti (37)=text (location_a _zuobiao (37,1),location_a _zuobiao (37,2)+1.5,'37');
ti (38)=text (location_a _zuobiao (38,1),location_a _zuobiao (38,2)+1.5,'38');
1 ?# n& R$ n) O8 Nti (39)=text (location_a _zuobiao (39,1),location_a _zuobiao (39,2)+1.5,'39');
ti (40)=text (location_a _zuobiao (40,1),location_a _zuobiao (40,2)+1.5,'40');
, i8 O3 X7 M* Z( [4 Rti (41)=text (location_a _zuobiao (41,1),location_a _zuobiao (41,2)+1.5,'41');
ti (42)=text (location_a _zuobiao (42,1),location_a _zuobiao (42,2)+1.5,'42');
ti (43)=text (location_a _zuobiao (43,1),location_a _zuobiao (43,2)+1.5,'43');
( g0 W- X# z7 Z! N8 h" @ti (44)=text (location_a _zuobiao (44,1),location_a _zuobiao (44,2)+1.5,'44');
9 R, ~! q% I2 z* [- W. d% y4 ^ti (45)=text (location_a _zuobiao (45,1),location_a _zuobiao (45,2)+1.5,'45');
ti (46)=text (location_a _zuobiao (46,1),location_a _zuobiao (46,2)+1.5,'46');
ti (47)=text (location_a _zuobiao (47,1),location_a _zuobiao (47,2)+1.5,'47');
! u/ ]8 j/ M5 a9 U6 b& Sti (48)=text (location_a _zuobiao (48,1),location_a _zuobiao (48,2)+1.5,'48');
ti (49)=text (location_a _zuobiao (49,1),location_a _zuobiao (49,2)+1.5,'49');
, p# f: l# [* `+ Y" x$ ?ti (50)=text (location_a _zuobiao (50,1),location_a _zuobiao (50,2)+1.5,'50');
ti (51)=text (location_a _zuobiao (51,1),location_a _zuobiao (51,2)+1.5,'51');
: \8 f3 S. C2 s4 O* Cti (52)=text (location_a _zuobiao (52,1),location_a _zuobiao (52,2)+1.5,'52');
ti (53)=text (location_a _zuobiao (53,1),location_a _zuobiao (53,2)+1.5,'53');
- T/ f" T" f: K' W/ G/ eti (54)=text (location_a _zuobiao (54,1),location_a _zuobiao (54,2)+1.5,'54');
9 c6 z, @2 {) J+ u0 z7 q% ~ti (55)=text (location_a _zuobiao (55,1),location_a _zuobiao (55,2)+1.5,'55');
- L+ P5 W0 ?3 x' k& x; v" kti (56)=text (location_a _zuobiao (56,1),location_a _zuobiao (56,2)+1.5,'56');
7 R+ n' ?  P9 S) qti (57)=text (location_a _zuobiao (57,1),location_a _zuobiao (57,2)+1.5,'57');
ti (58)=text (location_a _zuobiao (58,1),location_a _zuobiao (58,2)+1.5,'58');
: o( T) p7 b2 s. I; hti (59)=text (location_a _zuobiao (59,1),location_a _zuobiao (59,2)+1.5,'59');
& w2 n; A5 R( x  A8 wti (60)=text (location_a _zuobiao (60,1),location_a _zuobiao (60,2)+1.5,'60');
- G$ F' S( x- q- }( Q- O" Rti (61)=text (location_a _zuobiao (61,1),location_a _zuobiao (61,2)+1.5,'61');
ti (62)=text (location_a _zuobiao (62,1),location_a _zuobiao (62,2)+1.5,'62');
ti (63)=text (location_a _zuobiao (63,1),location_a _zuobiao (63,2)+1.5,'63');
! r7 _; I- q$ H, E* pti (64)=text (location_a _zuobiao (64,1),location_a _zuobiao (64,2)+1.5,'64');
& {1 |  R% b7 F6 i) Y# x3 `ti (65)=text (location_a _zuobiao (65,1),location_a _zuobiao (65,2)+1.5,'65');
ti (66)=text (location_a _zuobiao (66,1),location_a _zuobiao (66,2)+1.5,'66');
9 r2 C$ m- m8 Z9 s2 B3 mti (67)=text (location_a _zuobiao (67,1),location_a _zuobiao (67,2)+1.5,'67');
4 V& S* A5 O3 n9 Tti (68)=text (location_a _zuobiao (68,1),location_a _zuobiao (68,2)+1.5,'68');
ti (69)=text (location_a _zuobiao (69,1),location_a _zuobiao (69,2)+1.5,'69');
* v. z/ x1 W6 z$ q. o# C$ U( Oti (70)=text (location_a _zuobiao (70,1),location_a _zuobiao (70,2)+1.5,'70');
! M5 j( a' a( N( _, Fti (71)=text (location_a _zuobiao (71,1),location_a _zuobiao (71,2)+1.5,'71');
ti (72)=text (location_a _zuobiao (72,1),location_a _zuobiao (72,2)+1.5,'72');
ti (73)=text (location_a _zuobiao (73,1),location_a _zuobiao (73,2)+1.5,'73');
$ j7 M  O) C+ {ti (74)=text (location_a _zuobiao (74,1),location_a _zuobiao (74,2)+1.5,'74');
# C$ g( ^( y8 E+ ~: F: Xti (75)=text (location_a _zuobiao (75,1),location_a _zuobiao (75,2)+1.5,'75');
ti (76)=text (location_a _zuobiao (76,1),location_a _zuobiao (76,2)+1.5,'76');
" T$ x+ u. A4 d. n/ Cti (77)=text (location_a _zuobiao (77,1),location_a _zuobiao (77,2)+1.5,'77');
ti (78)=text (location_a _zuobiao (78,1),location_a _zuobiao (78,2)+1.5,'78');
ti (79)=text (location_a _zuobiao (79,1),location_a _zuobiao (79,2)+1.5,'79');
ti (80)=text (location_a _zuobiao (80,1),location_a _zuobiao (80,2)+1.5,'80');
ti (81)=text (location_a _zuobiao (81,1),location_a _zuobiao (81,2)+1.5,'81');
ti (82)=text (location_a _zuobiao (82,1),location_a _zuobiao (82,2)+1.5,'82');
ti (83)=text (location_a _zuobiao (83,1),location_a _zuobiao (83,2)+1.5,'83');
ti (84)=text (location_a _zuobiao (84,1),location_a _zuobiao (84,2)+1.5,'84');
ti (85)=text (location_a _zuobiao (85,1),location_a _zuobiao (85,2)+1.5,'85');
5 f$ |9 a, I4 P3 l2 W, Q6 V! Fti (86)=text (location_a _zuobiao (86,1),location_a _zuobiao (86,2)+1.5,'86');
3 }- I0 @9 }8 @# a* V8 Hti (87)=text (location_a _zuobiao (87,1),location_a _zuobiao (87,2)+1.5,'87');
+ S- t  o0 Y) S5 s& B0 f; q2 {ti (88)=text (location_a _zuobiao (88,1),location_a _zuobiao (88,2)+1.5,'88');
ti (89)=text (location_a _zuobiao (89,1),location_a _zuobiao (89,2)+1.5,'89');
. T5 S  \/ X7 P6 Iti (90)=text (location_a _zuobiao (90,1),location_a _zuobiao (90,2)+1.5,'90');
ti (91)=text (location_a _zuobiao (91,1),location_a _zuobiao (91,2)+1.5,'91');
+ r- J- I; i/ n7 B' ?ti (92)=text (location_a _zuobiao (92,1),location_a _zuobiao (92,2)+1.5,'92');


) l/ A$ a& Z7 L6 Mfenpei=[1        1
  c. [+ a# ~4 c1        67
1        68
1 D# Z2 F; Y' `, u1        69
1        71
5 y( Z' R/ `! d1        73
0 Z( B8 b- p* o! c! F& P2 {1        74
) s& V7 ~* O: [  `: F1        75
1        76
1        78
2        2
2        39
/ i# c; n& v; D) C( e2        40
/ }3 e3 `1 y& c1 t5 I& s2        43
2        44
$ x) _1 Z: {' _4 L, h5 E2        70
+ g% O- j' Z8 T' C1 i! p4 b0 P. l( ~7 S2        72
4 t' N* w0 c' j  X/ b3        3
3 a: `! `( m4 T* D( @9 P- `7 d# E' r3        54
3        55
6 m- j3 C6 b$ M* x3        65
3        66
' z* u5 u: P$ H9 Z4        4
: |* a$ [) X; B6 H0 W4        57
4        60
+ A7 S' e" }7 |2 B: `4        62
4        63
4        64
0 s: T! X) B' f3 M4 {  M- q5        5
. T6 M6 h: @9 H" Y1 T' q) z9 ?5        49
" n  C# T2 \* O# |& {! L* s- ]+ ~- Z5        50
& I% g* m' B5 }4 h" j2 I5        51
6 `: r4 S) L2 k) d$ B+ J5        52
: ^6 V$ a% Y. j, [5        53
5        56
5        58
5        59
8 [' ^' l* U' P& y- c6        6
' t5 g0 S7 m4 j- }8 R7        7
5 P# P' J1 S( P2 R4 Q7        30
7        32
7        47
0 X# m. b" d7 j) B7        48
3 ]% ~+ D7 a3 p  v7        61
8        8
2 H$ `; s7 M* J- X/ Y. ]8        33
8        46
1 f* ~* P5 }1 x  O& W  t9        9
; @! B5 _3 f7 A# F9        31
9        34
1 Q( u& R9 y3 Q1 E6 Q9        35
- H) S) y. i; o% d9        45
' `0 W4 h% ?! Y& s9 e10        10
& B4 D! S6 r) b1 t9 S. P' l11        11
11        26
11        27
12        12
& y7 h7 W# d, l$ s" P3 F: l1 q12        25
13        13
4 Z' j$ V% T8 t2 E/ v/ W13        21
13        22
13        23
13        24
14        14
15        15
# D  B8 z0 k$ x/ R3 P% `15        28
" k! g, J  p* L2 i* D/ h15        29
; v) L0 x9 \. t' Q8 A16        16
16        36
16        37
/ V' I( d6 L0 {1 }- r! j# g16        38
; p7 b9 A- V% R+ L& r17        17
17        41
17        42
7 m* V: A3 b9 V7 t& {' d) j, J7 j18        18
0 e  O1 L3 ~7 V5 ~18        80
18        81
. k8 ~/ a6 O" D( U% W* F18        82
18        83
9 J/ M* k' J0 s% g$ ]( n# y" J19        19
; d# n# B7 T# |7 D9 T( K19        77
19        79
20        20
20        84
20        85
20        86
20        87
& X7 h, n1 n( r, d. X4 ~20        88
20        89
20        90
- t2 H3 K' a! r$ ?) B& h& \  F20        91
20        92];
( ^, }* O! {' j3 [3 Z  L
( p1 i1 c7 |6 l2 z1 H$ Pfor i=1:92
    if fenpei (i,1)==1
* w" Y! f$ p  W( A# [; H        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'b');hold on;
    elseif fenpei (i,1)==2
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'g');hold on;
    elseif fenpei (i,1)==3
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'r');hold on;
4 D% l. g) n  ?, p; y    elseif fenpei (i,1)==4
: s1 K0 ?( H. ]( m9 Z! u+ g        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'c');hold on;
* o; X! m6 a. S& g; V$ x! k    elseif fenpei (i,1)==5
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'m');hold on;
    elseif fenpei (i,1)==6
& E$ _- W% T* S# z( T& [/ T        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'y');hold on;
- S9 W* p' L0 ^1 ?5 a2 Y7 y' I9 o    elseif fenpei (i,1)==7
2 d) f6 `5 W1 |        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'k');hold on;
    elseif fenpei (i,1)==8
" v( O% `7 z' ^3 J- x/ y        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'b+');hold on;
5 O/ M, K+ ]; u' V; t( q5 Z/ o, [    elseif fenpei (i,1)==9
; @- }9 @5 Z+ t% X0 v        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bo');hold on;
' C$ Z! N0 x  v    elseif fenpei (i,1)==10
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bs');hold on;
    elseif fenpei (i,1)==11
( ?3 ^0 q3 m% x; v9 ?7 M2 \* n  ]        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'g+');hold on;
    elseif fenpei (i,1)==12
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'go');hold on;
0 ~: i! _  a" C    elseif fenpei (i,1)==13
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'gs');hold on;
$ s7 d: `, b0 e3 p4 @! u- I& b    elseif fenpei (i,1)==14
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'r+');hold on;
    elseif fenpei (i,1)==15
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'ro');hold on;
. A! F. Z# [$ V& r9 w* ^* |" l    elseif fenpei (i,1)==16
$ e3 c. i" T- d0 J        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'rs');hold on;
    elseif fenpei (i,1)==17
: t9 m) Z) o1 F, U        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'c+');hold on;
9 V, w$ p6 `8 ~9 \& q' ~; R1 D    elseif fenpei (i,1)==18
% Z0 {5 U1 o! b+ i1 Y' f        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'co');hold on;
    elseif fenpei (i,1)==19
) L  F' j* Z+ X6 q        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'cs');hold on;
. s. v5 ^* p- S" o0 d( Q7 D    else
, W6 j; }- K6 ]- i+ o7 z5 X        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'m+');hold on;
    end
end
( J2 e1 q9 ^! m# z: B
/ X! g1 x. l/ Q, q4 n7 hfor i=1:92
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'b');hold on;
end
' [6 T7 @  p6 B  |7 qaxis ([190,460,250,410]);
for i=1:140
# x3 k0 f) n8 R4 _: o    tt=c (i);
    uu=location_all _daolu (tt,;
$ q& Q  p0 O6 `6 y/ ]) p    uu1=uu (1);uu2=uu (2);
    vv1=location_a _zuobiao (uu1,;
    vv2=location_a _zuobiao (uu2,;
3 I4 _& x* }% i" C0 V, P    ww1=[vv1 (1),vv2 (1)];
    ww2=[vv1 (2),vv2 (2)];
8 Y5 ?7 C" D( \# J' {    line (ww1,ww2,'Color',[.5 .5 .5])
- \' E2 R( [) e) b& r+ K4 Xend
* {& Z. S- `9 `5 _. A& {: P) l
5 p- {  T! P# oti=zeros (1,92);
ti (1)=text (location_a _zuobiao (1,1),location_a _zuobiao (1,2)+1.5,'1');
ti (2)=text (location_a _zuobiao (2,1),location_a _zuobiao (2,2)+1.5,'2');
7 {8 s0 {3 u+ ^+ ]6 S+ q( f- kti (3)=text (location_a _zuobiao (3,1),location_a _zuobiao (3,2)+1.5,'3');
; @2 B; R9 L( x3 t/ U; ^ti (4)=text (location_a _zuobiao (4,1),location_a _zuobiao (4,2)+1.5,'4');
ti (5)=text (location_a _zuobiao (5,1),location_a _zuobiao (5,2)+1.5,'5');
ti (6)=text (location_a _zuobiao (6,1),location_a _zuobiao (6,2)+1.5,'6');
5 D4 H1 I- k& U; lti (7)=text (location_a _zuobiao (7,1),location_a _zuobiao (7,2)+1.5,'7');
( m* F" y9 _; F. ]8 Bti (8)=text (location_a _zuobiao (8,1),location_a _zuobiao (8,2)+1.5,'8');
  u1 K# d. l/ r% m$ jti (9)=text (location_a _zuobiao (9,1),location_a _zuobiao (9,2)+1.5,'9');
5 ?/ F7 N2 m/ t1 ~$ J4 a! @ti (10)=text (location_a _zuobiao (10,1),location_a _zuobiao (10,2)+1.5,'10');
, e& C+ |$ W% D5 ?ti (11)=text (location_a _zuobiao (11,1),location_a _zuobiao (11,2)+1.5,'11');
2 f& y7 b+ M" P/ [; U6 a+ fti (12)=text (location_a _zuobiao (12,1),location_a _zuobiao (12,2)+1.5,'12');
ti (13)=text (location_a _zuobiao (13,1),location_a _zuobiao (13,2)+1.5,'13');
/ A" Q! ^1 n$ z# F% ?" k1 u% dti (14)=text (location_a _zuobiao (14,1),location_a _zuobiao (14,2)+1.5,'14');
& ?% \$ e5 s6 M( }, q# i+ H# Zti (15)=text (location_a _zuobiao (15,1),location_a _zuobiao (15,2)+1.5,'15');
$ U1 q7 h6 F% V& f4 K5 {  Rti (16)=text (location_a _zuobiao (16,1),location_a _zuobiao (16,2)+1.5,'16');
ti (17)=text (location_a _zuobiao (17,1),location_a _zuobiao (17,2)+1.5,'17');
ti (18)=text (location_a _zuobiao (18,1),location_a _zuobiao (18,2)+1.5,'18');
ti (19)=text (location_a _zuobiao (19,1),location_a _zuobiao (19,2)+1.5,'19');
ti (20)=text (location_a _zuobiao (20,1),location_a _zuobiao (20,2)+1.5,'20');
  h4 d( {( H* B; Q8 ^ti (21)=text (location_a _zuobiao (21,1),location_a _zuobiao (21,2)+1.5,'21');
" _& {# G7 _9 D) b( P/ hti (22)=text (location_a _zuobiao (22,1),location_a _zuobiao (22,2)+1.5,'22');
8 G; _' [. |  jti (23)=text (location_a _zuobiao (23,1),location_a _zuobiao (23,2)+1.5,'23');
) @9 i% h5 k; S, x0 @6 m0 iti (24)=text (location_a _zuobiao (24,1),location_a _zuobiao (24,2)+1.5,'24');
ti (25)=text (location_a _zuobiao (25,1),location_a _zuobiao (25,2)+1.5,'25');
  A. m- b9 Q! o; s' b+ b* Z1 {ti (26)=text (location_a _zuobiao (26,1),location_a _zuobiao (26,2)+1.5,'26');
ti (27)=text (location_a _zuobiao (27,1),location_a _zuobiao (27,2)+1.5,'27');
1 t/ d/ c4 q# r( n) Sti (28)=text (location_a _zuobiao (28,1),location_a _zuobiao (28,2)+1.5,'28');
9 _. E/ @$ a) d# Uti (29)=text (location_a _zuobiao (29,1),location_a _zuobiao (29,2)+1.5,'29');
1 S# R! B8 a/ C: j+ Q8 H0 S$ ?ti (30)=text (location_a _zuobiao (30,1),location_a _zuobiao (30,2)+1.5,'30');
9 }5 _8 G& A+ B) Z3 Yti (31)=text (location_a _zuobiao (31,1),location_a _zuobiao (31,2)+1.5,'31');
ti (32)=text (location_a _zuobiao (32,1),location_a _zuobiao (32,2)+1.5,'32');
) A/ g8 \! @$ N! ]; R" Q7 R1 `6 Gti (33)=text (location_a _zuobiao (33,1),location_a _zuobiao (33,2)+1.5,'33');
ti (34)=text (location_a _zuobiao (34,1),location_a _zuobiao (34,2)+1.5,'34');
  h6 P# L4 x9 q) {! B9 @ti (35)=text (location_a _zuobiao (35,1),location_a _zuobiao (35,2)+1.5,'35');
ti (36)=text (location_a _zuobiao (36,1),location_a _zuobiao (36,2)+1.5,'36');
ti (37)=text (location_a _zuobiao (37,1),location_a _zuobiao (37,2)+1.5,'37');
ti (38)=text (location_a _zuobiao (38,1),location_a _zuobiao (38,2)+1.5,'38');
- _- e, b* s6 M+ cti (39)=text (location_a _zuobiao (39,1),location_a _zuobiao (39,2)+1.5,'39');
ti (40)=text (location_a _zuobiao (40,1),location_a _zuobiao (40,2)+1.5,'40');
$ k, Y( A/ o4 C( qti (41)=text (location_a _zuobiao (41,1),location_a _zuobiao (41,2)+1.5,'41');
9 M3 a! @$ b! Jti (42)=text (location_a _zuobiao (42,1),location_a _zuobiao (42,2)+1.5,'42');
7 s' t: c* v3 T& Hti (43)=text (location_a _zuobiao (43,1),location_a _zuobiao (43,2)+1.5,'43');
# }, @; k# z/ D) S! ~8 w* pti (44)=text (location_a _zuobiao (44,1),location_a _zuobiao (44,2)+1.5,'44');
ti (45)=text (location_a _zuobiao (45,1),location_a _zuobiao (45,2)+1.5,'45');
ti (46)=text (location_a _zuobiao (46,1),location_a _zuobiao (46,2)+1.5,'46');
ti (47)=text (location_a _zuobiao (47,1),location_a _zuobiao (47,2)+1.5,'47');
$ G: r# ?# S3 V! Zti (48)=text (location_a _zuobiao (48,1),location_a _zuobiao (48,2)+1.5,'48');
% z3 z' V8 ?: O+ H6 s# Hti (49)=text (location_a _zuobiao (49,1),location_a _zuobiao (49,2)+1.5,'49');
ti (50)=text (location_a _zuobiao (50,1),location_a _zuobiao (50,2)+1.5,'50');
: b2 w. c! ]! rti (51)=text (location_a _zuobiao (51,1),location_a _zuobiao (51,2)+1.5,'51');
ti (52)=text (location_a _zuobiao (52,1),location_a _zuobiao (52,2)+1.5,'52');
ti (53)=text (location_a _zuobiao (53,1),location_a _zuobiao (53,2)+1.5,'53');
! \7 W) }2 m( S5 ?5 e3 bti (54)=text (location_a _zuobiao (54,1),location_a _zuobiao (54,2)+1.5,'54');
ti (55)=text (location_a _zuobiao (55,1),location_a _zuobiao (55,2)+1.5,'55');
ti (56)=text (location_a _zuobiao (56,1),location_a _zuobiao (56,2)+1.5,'56');
ti (57)=text (location_a _zuobiao (57,1),location_a _zuobiao (57,2)+1.5,'57');
ti (58)=text (location_a _zuobiao (58,1),location_a _zuobiao (58,2)+1.5,'58');
2 H/ |7 v1 k. P# [) V0 |) Tti (59)=text (location_a _zuobiao (59,1),location_a _zuobiao (59,2)+1.5,'59');
ti (60)=text (location_a _zuobiao (60,1),location_a _zuobiao (60,2)+1.5,'60');
ti (61)=text (location_a _zuobiao (61,1),location_a _zuobiao (61,2)+1.5,'61');
ti (62)=text (location_a _zuobiao (62,1),location_a _zuobiao (62,2)+1.5,'62');
ti (63)=text (location_a _zuobiao (63,1),location_a _zuobiao (63,2)+1.5,'63');
ti (64)=text (location_a _zuobiao (64,1),location_a _zuobiao (64,2)+1.5,'64');
- f6 {  c7 m# _& F( y8 }ti (65)=text (location_a _zuobiao (65,1),location_a _zuobiao (65,2)+1.5,'65');
2 M* ~; }9 n7 b5 R' x* jti (66)=text (location_a _zuobiao (66,1),location_a _zuobiao (66,2)+1.5,'66');
9 j' U5 {, d# @$ y* N: mti (67)=text (location_a _zuobiao (67,1),location_a _zuobiao (67,2)+1.5,'67');
7 d! _  C+ b  O, U- Fti (68)=text (location_a _zuobiao (68,1),location_a _zuobiao (68,2)+1.5,'68');
ti (69)=text (location_a _zuobiao (69,1),location_a _zuobiao (69,2)+1.5,'69');
3 Y8 t; L" `7 s. `0 e' M: q3 I1 oti (70)=text (location_a _zuobiao (70,1),location_a _zuobiao (70,2)+1.5,'70');
ti (71)=text (location_a _zuobiao (71,1),location_a _zuobiao (71,2)+1.5,'71');
ti (72)=text (location_a _zuobiao (72,1),location_a _zuobiao (72,2)+1.5,'72');
. c3 }* N* }* d% k- L3 T0 Z6 rti (73)=text (location_a _zuobiao (73,1),location_a _zuobiao (73,2)+1.5,'73');
ti (74)=text (location_a _zuobiao (74,1),location_a _zuobiao (74,2)+1.5,'74');
7 O# x) n3 i) n3 F, E% D7 T0 Vti (75)=text (location_a _zuobiao (75,1),location_a _zuobiao (75,2)+1.5,'75');
ti (76)=text (location_a _zuobiao (76,1),location_a _zuobiao (76,2)+1.5,'76');
ti (77)=text (location_a _zuobiao (77,1),location_a _zuobiao (77,2)+1.5,'77');
7 R+ W' |; K5 u" W8 }* B0 lti (78)=text (location_a _zuobiao (78,1),location_a _zuobiao (78,2)+1.5,'78');
ti (79)=text (location_a _zuobiao (79,1),location_a _zuobiao (79,2)+1.5,'79');
ti (80)=text (location_a _zuobiao (80,1),location_a _zuobiao (80,2)+1.5,'80');
) H1 D( r2 o! U5 s( [% `ti (81)=text (location_a _zuobiao (81,1),location_a _zuobiao (81,2)+1.5,'81');
ti (82)=text (location_a _zuobiao (82,1),location_a _zuobiao (82,2)+1.5,'82');
5 @" m5 M% u* V* |ti (83)=text (location_a _zuobiao (83,1),location_a _zuobiao (83,2)+1.5,'83');
ti (84)=text (location_a _zuobiao (84,1),location_a _zuobiao (84,2)+1.5,'84');
* D% |  K- b: X* pti (85)=text (location_a _zuobiao (85,1),location_a _zuobiao (85,2)+1.5,'85');
ti (86)=text (location_a _zuobiao (86,1),location_a _zuobiao (86,2)+1.5,'86');
ti (87)=text (location_a _zuobiao (87,1),location_a _zuobiao (87,2)+1.5,'87');
ti (88)=text (location_a _zuobiao (88,1),location_a _zuobiao (88,2)+1.5,'88');
; V. |' c- V$ |- {ti (89)=text (location_a _zuobiao (89,1),location_a _zuobiao (89,2)+1.5,'89');
ti (90)=text (location_a _zuobiao (90,1),location_a _zuobiao (90,2)+1.5,'90');
ti (91)=text (location_a _zuobiao (91,1),location_a _zuobiao (91,2)+1.5,'91');
ti (92)=text (location_a _zuobiao (92,1),location_a _zuobiao (92,2)+1.5,'92');
& ^3 X& S$ k! h6 ^4 X   
' q8 Q0 U- x/ a9 h: u* Shelp voronoi
2 ^3 a7 A( ~, |# c. l[vx,vy] = voronoi (location_a _zuobiao _x,location_a _zuobiao _y);

; G0 N. ^9 M8 g# u3 `" Y- Afor i=1:92
    if fenpei (i,1)==1
1 ?! i/ E% t, g; ?7 ~        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bh');hold on;
    elseif fenpei (i,1)==2
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'gh');hold on;
    elseif fenpei (i,1)==3
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'rh');hold on;
* u( y' X! l! `' Y# p4 k    elseif fenpei (i,1)==4
5 Y( S4 ~9 g5 W( W, ]        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'ch');hold on;
    elseif fenpei (i,1)==5
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'mh');hold on;
. i  j, r4 W! q0 A+ w, c    elseif fenpei (i,1)==6
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'yh');hold on;
    elseif fenpei (i,1)==7
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'kh');hold on;
. n: B- W& c. ~9 M# g' @& s    elseif fenpei (i,1)==8
3 h/ a2 f! S4 A+ t4 t/ J0 ^        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'b+');hold on;
' G) r1 ]2 [  @    elseif fenpei (i,1)==9
, i) |" G+ s, w1 [+ Q        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bo');hold on;
+ Y% R5 e, y* {  l7 f, Q    elseif fenpei (i,1)==10
) L, s% `2 k  ]3 q% |        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bs');hold on;
    elseif fenpei (i,1)==11
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'g+');hold on;
    elseif fenpei (i,1)==12
  a3 P  s! O& Y& n! I  H        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'go');hold on;
    elseif fenpei (i,1)==13
$ u, J& l) u5 I: x  J        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'gs');hold on;
    elseif fenpei (i,1)==14
! w- t1 r" o  s. {) F; H' v& p0 v        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'r+');hold on;
    elseif fenpei (i,1)==15
) G7 M5 ?) U0 q% v! [        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'ro');hold on;
' W! Q/ e0 b4 t$ |2 u0 X7 ]/ Q    elseif fenpei (i,1)==16
  u  }% ?" m) @. Y1 u: z        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'rs');hold on;
    elseif fenpei (i,1)==17
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'c+');hold on;
0 M  y0 g8 n- B1 A. [4 T    elseif fenpei (i,1)==18
* Z% Y: X9 `, J+ C    plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'co');hold on;
% f1 z0 e8 B& G    elseif fenpei (i,1)==19
" r# i/ C0 c/ P* \; R9 L, I        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'cs');hold on;
    else
& f- I# |& x: t1 L; N        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'m+');hold on;
    end
end
1 \- f) S6 @4 r# y; F* ^2 W! f附件2:邻接矩阵的matlab实现程序
data=[
" w2 m1 Z9 `8 z9 D, j# y. W1   75
1   78
9 F! c, d$ i, l4 R1 |2   44
3   45
* r8 y  f" I! L& |3   65
4   39
4   63
5   49
# _" p& ?4 v5 l5   50
6 C" O$ s- P2 \( L1 q0 z. L6   59
% H6 `- b% ^% U/ i" U7   32
: e- W7 {0 l- G' h) b/ K+ ~  a7   47
8   9
" o& r# H0 h0 |7 t* q" |" {. R8   47
9   35
0 F, K% j& _( F+ w/ K10  34
) G( ?' @9 W4 N& V/ E11  22
11  26
12  25
12  471
8 {% H! i+ x9 J, y/ Y14  21
15  7
15  31
, i' {- c2 j+ p$ l; K4 c16  14
2 _) `: [8 x  ?/ o# P) w16  38
17  40
17  42
17  81
* d- `1 N' q& k8 w18  81
  k9 J  n1 }) \, p4 b18  83
6 f* a" E7 z: Z% j: |1 v8 Z19  79
20  86
21  22
22  372
22  13
" ?  V- D) P/ M& L9 _% E( r  L23  13
2 t. w+ [5 h' h) Z2 f23  383
* V3 J9 n! k$ F/ Z) R& T24  13
24  25
  P; X+ ^0 f5 b25  11
26  27
' u9 n! T2 V( A, R. E26  10
27  12
28  29
28  15
9 H" b1 w; q/ ~: X9 ^% W1 l0 a% I29  30
30  7
3 k+ s6 I6 O3 L/ k30  48
; j8 u( V- a& I$ c  J+ D31  32
. _& m1 |: z, K# I3 X& ]6 Z31  34
* h! g& }- }1 J( M6 h32  33
33  34
33  8
4 f2 u2 w$ q0 D. y2 M34  9
35  45
  M' @8 l5 Z# f, x( i/ Z' l36  35
. _4 L1 j* [: D. F; o6 @$ Q& n8 G9 o36  37
36  16
8 `& @% i% Y9 [' [/ \. N36  39
4 f; y8 j8 c; x; f3 M4 s37  7
38  39
38  41
39  40
" e0 ~6 v2 `  J. Z+ i0 w40  2
41  17
41  92
! I7 A- e. R/ O! M6 j42  43
* Q$ c! n3 p/ ]) W; h& J43  2
' _9 |/ k1 ]% r, U. m/ O: a; e43  72
44  3
45  46
46  8
46  55
7 N; t0 N/ i' _" U47  48
) ]; u# h  X2 m( @47  6
. h( o. [- M( O; J7 j; s47  5
8 X, }7 C% L# ]$ z5 }48  61
49  50
; |/ S* Z8 D5 k- f49  53
& k% ]) I  U4 r50  51
4 N. M( U3 T2 i51  52
5 p& x* G- U- ^9 [, S. n51  59
52  56
$ Z& O9 Z  h7 S5 s53  52
$ n: x: \2 j) b9 v+ F+ F53  54
54  55
1 ~0 x1 s6 J* e% O- ?4 ]2 S. g$ S54  63
2 `& o: S9 U9 @55  3
56  57
3 Y9 v+ d9 ]/ c5 Q' Z1 ]57  58
57  60
8 i, B; N/ E3 \5 \9 i% r* ]8 b57  4
58  59
60  62
& \) ?% f2 Q/ O61  60
5 F% z# k1 ^! u62  4
7 q2 ]7 u8 W8 V" n  S" g62  85
63  64
5 A5 c' O8 r7 h6 D0 r7 J; f64  65
64  76
1 v! h4 i- Y# a3 s6 f65  66
66  67
66  76
8 z! V- B% G% r" a" C% W6 Y' D8 }; D/ ?67  44
% q1 ?3 {# N8 P/ @67  68
7 K0 D2 y2 o3 D% P68  69
68  75
6 Q  L; g7 n& A" R" _, R; H. N& y7 l69  70
69  71
  l# R" Y% `2 x! A( e' @: B69  1
: u4 P/ O1 m4 o% m- Y70  2
70  43
71  72
71  74
" K# h* z  r% ?0 ^( i# @7 D72  73
73  74
73  18
% _' l4 S2 V, D1 z0 @  x6 n74  1
74  80
75  76
76  77
77  78
7 Y: T3 H  M  M! A- o77  19
78  79
79  80
* q' a; ^% @# k- J9 S; [; K80  18
) ?8 C2 ^  M* ?81  82
82  83
82  90
83  84
84  85
( T" w! {+ T: o5 I5 q85  20
8 h, A8 q& D" s+ @( X86  87
2 {' L, Q0 v4 ?3 z86  88
87  88
5 ^" v) R' p5 T) C0 R87  92
; B3 n" w* c2 [88  89
88  91
6 j# m9 Q1 C% X' P0 c3 _' K89  20
89  84
. }5 U; [( ~3 h7 d/ I6 Z# }0 P89  90
90  91
91  92
];
x=data(:,1);
" P' \9 {8 Y' ?! ty=data(:,2);
z=[6.264982043
6.403124237
9.486832981
( G, c; U  x* e1 u# i7 L% c42.46469122
15.23975065
45.60975773
10.30776406
% ], G; V% Y: D, ^6 ^5
8.485281374
16.03121954
11.40175425
12.80624847
11.5974135
+ s1 ~2 x2 q  j20.79663434
4.242640687
% J" x9 s& n% l49.21635907
% ]# g$ C/ Z7 `1 W: {; H/ k, n32.69556545
+ q. b% @# g4 U, Q; l3 X0 T9
17.88854382
$ r# P6 O* y! W* A7 Y384.4697647
, d1 {# x4 r2 k9 H: ~32.64965543
38.18376618
40
67.41661516
: [# V, n) C3 A! p34.05877273
7 ^' @* U0 V5 J. ?, T/ S% Q26.87936011
9.848857802
40.22437072
6.708203932
/ Y# l. H2 h9 @5.385164807
2 i$ Y& a# ?4 M- a( I1 z4.472135955
3.605551275
18.02775638
358.0460864
9.055385138
5
347.6348659
23.85372088
18.02775638
- [7 ^# q1 K& T- F9 v! h4 q* ]" k20.02498439
7.433034374
/ E3 V+ G! q: M/ `* v35.38361203
33.04920574
9.486832981
3 \" ~4 B& }7 D3 G: v; J2 R# K/ G47.51841748
. g  r) I- R$ ]9 s" m$ S74.3236167
* _3 {2 m4 a5 |" s$ G2 C5.830951895
7.071067812
11.70469991
1 C% [! Q+ V- @. }+ I$ E, Z15.53222457
+ B' o" D  A1 {0 m3 {/ H5.099019514
7.566372975
8.276472679
% M3 u0 q7 @0 ?5 w3 k. p5.024937811
) Z' U& P, ^1 _% P$ y6.708203932
5
2 c; D6 t+ k! q- h+ e6 _9 o5.099019514
$ \; i" ]# F! {  F6.08276253
- g; z" ^% @# s35.0142828
30.41381265
1 g# C% B0 i8 f6 Y8 h3
40.07804885
17.67766953
% Q3 a8 D& ?* {19.14418972
8.5
46.31684359
8.062257748
8
8.062257748
! l9 T3 t* Z& t11.62970335
) [" ^& n. j* M$ ]5 m2 w- A6
+ }' c( G  z0 y" _9.300537619
29.42787794
: `! U. O2 W& e1 |4 p10.19803903
14.56021978
56.26944108
29
6 W" k0 T) g1 x9 ]8 Z2 q10.44030651
2 Q! o5 \6 g1 x. {6.708203932
3.807886553
2 O$ N- k8 h$ \# G9 S" u* r4.301162634
5 [& ~: E- H; z" |' @. N4 X* {2.915475947
5 A& h3 }* t5 P* X' Q" Y4.242640687
, V% A: M/ I7 K# }$ K/ E8.544003745
22.8035085
10.04987562
24.18677324
; s3 S/ w  d( X6 ~12.6589889
8 v% Z2 X  K7 a+ Q. r0 a$ n12.3794184
) Y8 k  r: I$ [4 Y& c7.5
6 A( F0 t- U, @1 Q, u8.139410298
18.68154169
7.810249676
7 B& m+ }% z/ A3 B( ]13.89244399
* r4 |" L  p; _34.71310992
3.5
. i+ y  L& a+ c6 d60.01666435
9.055385138
5.830951895
0 E# M) U+ ?2 p. B13.15294644
6 q7 Q/ d  Z' l2 J) U3.16227766
4.242640687
9.219544457
14.76482306
4.123105626
7 W& E  o- \1 q" z8 d) m7.071067812
6 q9 v. i- X( D" X2 U4.527692569
5.385164807
6.403124237
6 {& k; L. \$ ~, w8 @/ u5
8.602325267
! y1 P- K. {) T8 V; y, Y1 C; `7.615773106
5
6.103277808
8.062257748
4.031128874
# e! y" S% `$ N- c8 A, |19.72308292
8 H+ e" X% m% c0 {& t  A6.264982043
( q- J2 @) M+ H9 h2 ^: d+ [16.91892432
& n4 D& h; ?+ H* x3.535533906
4.472135955
! J5 P$ F% p; @" S- w10
9.848857802
* C- U. x, U. U+ s" g" v2 d6.708203932
4.472135955
8.062257748
5.024937811
/ s1 M7 V. A* q+ @5 W- r2 {5.408326913
8.732124598
0 [8 e  a9 W! s9.848857802
7.280109889
+ B2 E7 P$ k& s* u0 l7 C9 ]9 a4.472135955
, ]# v" P1 u5 C* H11.04536102
* L. w1 W6 z3 S/ l5 k6 c, ~: C% B6 V9.340770846
4.031128874
3 T5 @6 a5 G9 g9 ~+ m. c$ b7 i21.37755833
4.031128874
# i* a) Q* }7 E3.041381265
9.486832981
0 `% R$ w0 |% _( N4 o; x3
( c- i3 m' j, T/ h1 w8 k% a3.535533906
4.74341649
; o' p$ j* S8 q/ W/ J20.02498439
];
7 q+ o/ b0 Q/ }, [2 Yxx=zeros(92);
' i7 Y$ \! Y% ~3 ?: c* R# qfor i=1:143
    xx(x(i),y(i))=z(i);
    xx(y(i),x(i))=z(i);
end
for i=1:92
1 @" t8 A% S: A    for j=1:92
    if xx(i,j)==0
        xx(i,j)=10000;
" R' J! ~: _' V+ t    end
/ f0 k2 q$ c# n+ A& x& F" h) k    end
$ j5 R9 |3 _" `  t8 e/ aend
for i=1:92
3 c; t# q9 J) t( \    xx(i,i)=0;
end
! e, R3 e9 x' Z- U, ^9 i
附件3:围堵方案的java实现程序
public void SearchMax(ArrayList<SerachNode>  searchnode ,int[][]quantu,ArrayList<Node> allNode)
{
while(searchnode.size()>0)
2 W3 O1 K! d' E0 d& f{
SerachNode tmp=searchnode.get(0);
searchnode.remove(0);

for(int j=0;j<quantu.length;j++) //扩展当前结点
{
int i=tmp.getBianhao();
if(quantu[i][j]>0)
$ K9 T$ c# d. T9 }{

SerachNode newNode=new SerachNode(allNode.get                                                  (j),tmp.getCurrentQuanzhi()+quantu[i][j]) ;
6 x3 i+ ?0 n8 t: }+ [' o
0 D% ]% n% ^# k3 rif(!bounderMAX(newNode))
/ H0 e7 r5 G9 ~{
MaxresultSet.add(newNode);
4 x3 j& `+ t* J. b3 Mcontinue;
/ X+ i& [8 A" i5 x# U0 [+ Z5 u
- U& p6 V7 D( x# }& s( l9 p- ~}
. T2 y. U/ v# O; @
- ?2 G( \2 ^8 u* T  J
4 Z/ [( |# W5 S6 Q3 ?$ k( pif(!this.contain(searchnode, newNode))
{
searchnode.add(newNode);
}
}}}
, r& E: S, w( ~8 X. e! T
2 ?" g* r2 r5 s. o, A附件4:全区的交巡警平台有效覆盖范围(有效代表三分钟内可以到达)
( y2 S8 K/ a3 ?9 k

, \3 t( a, G" R. Y附件5:用Mathmatica求数据均值与方差
3 F0 {# ?7 ?5 ~" G- v) Wdata={9.1,5.6,7.6,9,8.2,4.6,8.5,4.8,5.3,3.4,8.3,7.2,4.5,5,1.6,4,2.5,6.4,7.9,10.6};
2 [" r1 u( c7 YMean[data]
$ U1 w7 f+ E8 J1 ^5 A$ x7 [Variance[data]
9 r; s" T; `, z4 @; N
9 x9 B# b5 g1 J( k3 U, D) x% [附件6:输入任意两点的坐标,输出两点间距离
程序中以节点10和34间距离为例（Mathematica）
2 o8 b" }* {1 D$ n4 D# qIn[1]:
Clear["Glibal'*"]
x1=328;
" N' t! w7 ]8 C  Tx2=282;
( b, Z" ~6 }# T, o; L+ ~) j; ^" dy1=342.5;
+ V7 U3 }. f$ Gy2=325;
d=Sqrt[(x1-x2)^2+(y1-y2)^2]
1 w! S. Z5 F0 L/ jOut[1]:49.2164
4 m1 j+ E  @% U7 m2 `
$ R, Z( y+ B& l1 N0 r, y# \6 c' l- I$ Y  C附件7:A区各线路距离
% ^! t5 @% D7 O  Y7 S+ ]路线起点
% \5 ~8 l0 r2 s# T8 M标号        路线终点
标号        起点
横坐标        起点
纵坐标        终点
横坐标        终点
' Q" x$ z; V/ B0 l6 q纵坐标        距离
& a3 b, S+ e# b# L& g7 Y1        75        413        359        418.5        356        6.2649820431
1        78        413        359        417        364        6.4031242374
2        44        403        343        394        346        9.4868329805
3        45        383.5        351        342        342        42.464691215
3        65        383.5        351        395        361        15.239750654
5 e& E* h  ?# v1 }& _4        39        381        377.5        371        333        45.609757728
4        63        381        377.5        391        375        10.307764064
$ s# c5 k0 Y. [4 J' B5        49        339        376        342        372        5
5 f/ P: V- e0 Y+ G9 ~5        50        339        376        345        382        8.4852813742
% o. e1 a: B2 q8 N6        59        335        383        351        382        16.031219542
+ P$ i: a. c& d3 l1 Y3 i) j7        32        317        362        326        355        11.401754251
7        47        317        362        325        372        12.806248475
8        9        334.5        353.5        333        342        11.597413505
0 d2 x" J0 k; P- Q8        47        334.5        353.5        325        372        20.796634343
9        35        333        342        336        339        4.2426406871
10        34        282        325        328        342.5        49.216359069
11        22        247        301        234        271        32.695565449
11        26        247        301        256        301        9
; O" ]- q) D9 [. F12        25        219        316        227        300        17.88854382
12        471        219        316        155        316        64
14        21        280        292        251        277        32.649655435
15        7        290        335        317        362        38.183766184
15        31        290        335        314        367        40
16        14        337        328        280        292        67.416615163
16        38        337        328        371        330        34.058772732
17        40        415        335        388.5        330.5        26.879360111
17        42        415        335        419        344        9.8488578018
17        81        415        335        438        368        40.224370722
18        81        432        371        438        368        6.7082039325
18        83        432        371        434        376        5.3851648071
  K- o7 K2 s: g19        79        418        374        420        370        4.472135955
6 t+ }# q; H, @' I) m: w20        86        444        394        447        392        3.6055512755
# Z* o' G# p4 c3 K1 z4 _6 Q& R21        22        251        277        234        271        18.027756377
22        372        234        271        232.5        264        7.1589105316
! C# e' @. i$ s2 f: ~4 A* H22        13        234        271        225        270        9.0553851381
23        13        225        265        225        270        5
, ]* |! m8 v1 h9 v) N8 w; r23        383        225        265        192        264        33.015148038
24        13        212        290        225        270        23.853720884
24        25        212        290        227        300        18.027756377
- |! R& y7 p1 Y! ^25        11        227        300        247        301        20.024984395
7 Z  w  O( {$ `) ]" _' V# b4 n. E26        27        256        301        250.5        306        7.4330343737
$ W: j/ j7 G4 D* a# W( j: j1 {26        10        256        301        282        325        35.383612026
3 T7 `3 W* t. a5 B% O$ h1 _/ A27        12        250.5        306        219        316        33.049205739
& ]# i4 H; d# I/ [" h& S+ `. L: b28        29        243        328        246        337        9.4868329805
28        15        243        328        290        335        47.518417482
29        30        246        337        314        367        74.323616704
" _. h! M# b6 y30        7        314        367        317        362        5.8309518948
6 m; a: w( n; K7 X# R$ J/ ^) P, k+ N! W30        48        314        367        315        374        7.0710678119
& \' d& {3 N. q1 Q31        32        315        351        326        355        11.704699911
31        34        315        351        328        342.5        15.532224567
/ X5 y0 z0 O( z" V" e* ~32        33        326        355        327        350        5.0990195136
4 R+ W4 }. m8 z3 f" I33        34        327        350        328        342.5        7.5663729752
0 b$ l7 Y6 g* N/ H- ^33        8        327        350        334.5        353.5        8.2764726786
34        9        328        342.5        333        342        5.0249378106
35        45        336        339        342        342        6.7082039325
36        35        336        334        336        339        5
36        37        336        334        331        335        5.0990195136
8 J1 z9 T/ z) W% I8 {36        16        336        334        337        328        6.0827625303
36        39        336        334        371        333        35.0142828
/ y, u# g, {% a1 D37        7        331        335        317        362        30.413812651
! \) U5 d5 q3 `" l9 m38        39        371        330        371        333        3
38        41        371        330        411        327.5        40.078048855
2 ~6 B$ d  B/ r% Q( a2 q39        40        371        333        388.5        330.5        17.67766953
! j0 c- x9 g3 R* V# t; ^9 B+ G40        2        388.5        330.5        403        343        19.144189719
41        17        411        327.5        415        335        8.5
41        92        411        327.5        444        360        46.316843588
  _1 X+ t/ @- t7 T& ^4 d42        43        419        344        411        343        8.0622577483
* i9 p5 U5 |' {) c43        2        411        343        403        343        8
# O) o: ~  U# `' K- S7 [43        72        411        343        418        347        8.0622577483
44        3        394        346        383.5        351        11.62970335
45        46        342        342        342        348        6
46        8        342        348        334.5        353.5        9.3005376189
% R& X: M; V0 A, f' U" ^) {% |46        55        342        348        371        353        29.427877939
' E' d- t! V1 x' T47        48        325        372        315        374        10.198039027
  u' R) D6 g0 y: A# D4 F47        6        325        372        339        376        14.560219779
47        5        325        372        381        377.5        56.269441085
48        61        315        374        335        395        29
7 z. q0 J3 a1 \2 i- Q: Q3 V# b: t49        50        342        372        345        382        10.440306509
; X/ `; K( o3 k6 j# x5 k7 N49        53        342        372        348        369        6.7082039325
50        51        345        382        348.5        380.5        3.8078865529
51        52        348.5        380.5        351        377        4.3011626335
8 |# I8 M/ E* k: S- I2 b7 m* x1 h51        59        348.5        380.5        351        382        2.9154759474
. `& j: f+ v5 ]* v# P52        56        351        377        354        374        4.2426406871
53        52        348        369        351        377        8.5440037453
53        54        348        369        370        363        22.803508502
  ]$ U) s7 L% f; t; F54        55        370        363        371        353        10.049875621
3 C: i4 l- [* L. J: W54        63        370        363        391        375        24.186773245
55        3        371        353        383.5        351        12.658988901
6 I9 X& l0 z( d# I$ l' h3 V56        57        354        374        363        382.5        12.379418403
8 n1 h, ^- L0 r' n- m57        58        363        382.5        357        387        7.5
57        60        363        382.5        369        388        8.139410298
8 X3 Z- B" A! g57        4        363        382.5        381        377.5        18.681541692
1 s8 J6 A  q. g58        59        357        387        351        382        7.8102496759
60        62        369        388        381        381        13.892443989
% r/ V8 A, a( }% a. \5 V61        60        335        395        369        388        34.713109915
62        4        381        381        381        377.5        3.5
62        85        381        381        440        392        60.016664352
+ v8 w5 U4 y, g" r2 |63        64        391        375        392        366        9.0553851381
64        65        392        366        395        361        5.8309518948
64        76        392        366        405        368        13.152946438
65        66        395        361        398        362        3.1622776602
0 B) z5 t8 E  K2 o* w" J0 ]66        67        398        362        401        359        4.2426406871
66        76        398        362        405        368        9.2195444573
! z8 b' @& x" Z% g8 M67        44        401        359        394        346        14.76482306
67        68        401        359        405        360        4.1231056256
( e1 g7 r! P2 H68        69        405        360        410        355        7.0710678119
) r2 ~# `1 M% C9 y# \68        75        405        360        405.5        364.5        4.5276925691
69        70        410        355        408        350        5.3851648071
69        71        410        355        415        351        6.4031242374
69        1        410        355        413        359        5
70        2        408        350        403        343        8.602325267
  j$ n8 k+ A3 L# o2 P3 {' }; f% P70        43        408        350        411        343        7.6157731059
4 Q( |! h/ @) ]71        72        415        351        418        347        5
5 Q' e2 e* W! m2 a- h7 h2 [2 b  x4 S71        74        415        351        418.5        356        6.1032778079
6 W) V* N1 Z. g& F% g& h+ W72        73        418        347        422        354        8.0622577483
73        74        422        354        418.5        356        4.0311288741
5 B1 H6 Y; C" \" ^8 t( T4 K73        18        422        354        432        371        19.723082923
74        1        418.5        356        413        359        6.2649820431
74        80        418.5        356        424        372        16.918924316
7 o( j& F" l4 w7 G5 N9 A. d* L  j8 @75        76        405.5        364.5        405        368        3.5355339059
76        77        405        368        409        370        4.472135955
77        78        409        370        417        364        10
77        19        409        370        418        374        9.8488578018
78        79        417        364        420        370        6.7082039325
79        80        420        370        424        372        4.472135955
80        18        424        372        432        371        8.0622577483
, \9 a" U( z3 u4 A$ Z: {81        82        438        368        438.5        373        5.0249378106
82        83        438.5        373        434        376        5.4083269132
2 W4 D" e" V6 V  T2 @7 y82        90        438.5        373        440.5        381.5        8.7321245983
83        84        434        376        438        385        9.8488578018
84        85        438        385        440        392        7.2801098893
7 _; e( ~( o9 u. J85        20        440        392        444        394        4.472135955
86        87        447        392        448        381        11.045361017
86        88        447        392        444.5        383        9.3407708461
' w& |7 u! @0 x) c7 I87        88        448        381        444.5        383        4.0311288741
3 V" p) d& L6 ^1 d) K87        92        448        381        444        360        21.377558326
88        89        444.5        383        441        385        4.0311288741
; n$ X* I( `& h3 |  `6 W' l88        91        444.5        383        445        380        3.0413812651
% n! R6 J4 J5 ?9 a2 o- W89        20        441        385        444        394        9.4868329805
89        84        441        385        438        385        3
  ~/ @# n, X) t89        90        441        385        440.5        381.5        3.5355339059
( F" o2 |( }3 A+ n# C0 z90        91        440.5        381.5        445        380        4.7434164903
; c# u* p1 |8 J1 M0 K' f91        92        445        380        444        360        20.024984395

Dou1

太详细了，谢谢楼主

莫西莫西

灰常好，就是看着太费劲

李梦龙33

谢谢楼主分享，正在做。。。

霜之咏叹调

haobaniyignla  

往东就亮

谢谢楼主分享！！！！

qqwhw2012

好东西
111

467857726de

不错哦  谢谢楼主

hbdkfk2

谢谢楼主！！！！！

liupeng723911

楼主你太好了.........