交巡警服務(wù)平臺的設(shè)置與調(diào)度11年B題

1、全國大學(xué)生數(shù)學(xué)建模競賽承 諾 書我們仔細(xì)閱讀了中國大學(xué)生數(shù)學(xué)建模競賽的競賽規(guī)則.我們完全明白，在競賽開始后參賽隊(duì)員不能以任何方式（包括電話、電子郵件、網(wǎng)上咨詢等）與隊(duì)外的任何人（包括指導(dǎo)教師）研究、討論與賽題有關(guān)的問題。我們知道，抄襲別人的成果是違反競賽規(guī)則的, 如果引用別人的成果或其他公開的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻(xiàn)的表述方式在正文引用處和參考文獻(xiàn)中明確列出。我們鄭重承諾，嚴(yán)格遵守競賽規(guī)則，以保證競賽的公正、公平性。如有違反競賽規(guī)則的行為，我們將受到嚴(yán)肅處理。我們參賽選擇的題號是（從A/B/C/D中選擇一項(xiàng)填寫）： B 我們的參賽報(bào)名號為（如果賽區(qū)設(shè)置報(bào)名號的話）：

2、所屬學(xué)校（請?zhí)顚懲暾娜?西北大學(xué) 參賽隊(duì)員 (打印并簽名) ：1. 張舒岱 2. 劉羽 3. 張成悟 指導(dǎo)教師或指導(dǎo)教師組負(fù)責(zé)人 (打印并簽名)： 日期： 2014 年 8 月 10日全國大學(xué)生數(shù)學(xué)建模競賽編 號 專 用 頁賽區(qū)評閱編號（由賽區(qū)組委會評閱前進(jìn)行編號）：賽區(qū)評閱記錄（可供賽區(qū)評閱時使用）：評閱人評分備注全國統(tǒng)一編號（由賽區(qū)組委會送交全國前編號）：全國評閱編號（由全國組委會評閱前進(jìn)行編號）：交巡警服務(wù)平臺的設(shè)置與調(diào)度摘要交巡警服務(wù)平臺位置的選取以及劃分交巡警服務(wù)平臺的管轄范圍對于處理突發(fā)事件有非常大的影響?，F(xiàn)階段，一般依據(jù)經(jīng)驗(yàn)選取服務(wù)平臺位置及劃分管轄區(qū)域。所以如何科學(xué)合理

3、處理的交巡警服務(wù)平臺的設(shè)置與調(diào)度問題具有十分重要的現(xiàn)實(shí)意義。本文研究了交巡警服務(wù)平臺的設(shè)置與調(diào)度問題。具體討論了在給定的區(qū)域A內(nèi)，如何合理的設(shè)置交巡警服務(wù)平臺的管轄區(qū)域；發(fā)生特殊事件時應(yīng)如何調(diào)動服務(wù)平臺警力以快速封鎖區(qū)域A；應(yīng)該增加多少數(shù)量交巡警服務(wù)平臺以及在哪個位置增加。本文建立最短路模型、0-1整數(shù)規(guī)劃模型，利用軟件解決了分配各平臺管轄范圍、調(diào)度警務(wù)資源以及合理設(shè)置交巡警服務(wù)平臺這三個方面的問題。在解決分配各平臺管轄范圍問題時，本文建立了最短路模型。通過求解各個路口到交巡警平臺的距離是否滿足最低時間限制，解決交巡警服務(wù)平臺分配管轄范圍的問題。本文在軟件上運(yùn)用算法進(jìn)行求解，給出了中心城區(qū)A的

4、20個服務(wù)平臺的管轄范圍，并求得到達(dá)最近的交巡警服務(wù)平臺的時間超過3分鐘的6個路口。在解決調(diào)度警務(wù)資源快速封鎖城區(qū)的問題時，本文建立了0-1整數(shù)規(guī)劃模型。以封鎖城區(qū)所用時間最少為限制條件，利用lingo軟件編程求解，給出了該區(qū)交巡警服務(wù)平臺警力合理的調(diào)度方案，并求得對13個交通要道實(shí)現(xiàn)全封鎖最短需要8.01分鐘。在解決交巡警服務(wù)平臺的選址問題時，本文建立了雙目標(biāo)0-1整數(shù)規(guī)劃模型?？紤]到建設(shè)新的服務(wù)平臺需要投入更多的成本和警務(wù)資源，還需平衡各個服務(wù)平臺的工作量。因此，以增加服務(wù)平臺數(shù)最小和服務(wù)平臺工作量方差最小為目標(biāo)，建立了雙目標(biāo)0-1整數(shù)規(guī)劃模型。解出增加的服務(wù)平臺數(shù)為4個，新增的服務(wù)平臺具

5、體位置為A29,A39,A48,A88。 本文所提供的模型考慮到均衡各個交巡警服務(wù)平臺的工作量和新建服務(wù)臺的成本，使結(jié)果更加合理符合需求，可以推廣到任何一個市區(qū)甚至更廣范圍內(nèi)的交巡警服務(wù)平臺的設(shè)置與調(diào)度問題的解決中。也可以廣泛應(yīng)用于社區(qū)衛(wèi)生室、公共衛(wèi)生間、消防救火中心等社會服務(wù)部門的選址問題，對實(shí)際有指導(dǎo)意義。關(guān)鍵詞：Dijkstra算法 雙目標(biāo)0-1整數(shù)規(guī)劃模型 Lingo編程一、問題重述“有困難找警察”，是家喻戶曉的一句流行語。警察肩負(fù)著刑事執(zhí)法、治安管理、交通管理、服務(wù)群眾四大職能。為了更有效地貫徹實(shí)施這些職能，需要在市區(qū)的一些交通要道和重要部位設(shè)置交巡警服務(wù)平臺。每個交巡警服務(wù)平臺的職

6、能和警力配備基本相同。由于警務(wù)資源是有限的，如何根據(jù)城市的實(shí)際情況與需求合理地設(shè)置交巡警服務(wù)平臺、分配各平臺的管轄范圍、調(diào)度警務(wù)資源是警務(wù)部門面臨的一個實(shí)際課題。試就某市設(shè)置交巡警服務(wù)平臺的相關(guān)情況，建立數(shù)學(xué)模型分析研究下面的問題：附件1中的附圖1給出了該市中心城區(qū)A的交通網(wǎng)絡(luò)和現(xiàn)有的20個交巡警服務(wù)平臺的設(shè)置情況示意圖，相關(guān)的數(shù)據(jù)信息見附件2。請為各交巡警服務(wù)平臺分配管轄范圍，使其在所管轄的范圍內(nèi)出現(xiàn)突發(fā)事件時，盡量能在3分鐘內(nèi)有交巡警（警車的時速為60km/h）到達(dá)事發(fā)地。對于重大突發(fā)事件，需要調(diào)度全區(qū)20個交巡警服務(wù)平臺的警力資源，對進(jìn)出該區(qū)的13條交通要道實(shí)現(xiàn)快速全封鎖。實(shí)際中一個平臺

7、的警力最多封鎖一個路口，請給出該區(qū)交巡警服務(wù)平臺警力合理的調(diào)度方案。根據(jù)現(xiàn)有交巡警服務(wù)平臺的工作量不均衡和有些地方出警時間過長的實(shí)際情況，擬在該區(qū)內(nèi)再增加2至5個平臺，請確定需要增加平臺的具體個數(shù)和位置。2、 模型假設(shè)（1）每個交巡警服務(wù)平臺的職能和警力配備基本相同；（2）警車的行駛速度恒定，不考慮實(shí)際交通狀況的影響； （3）交巡警服務(wù)平臺接到報(bào)警后能立即出警，中間沒有延誤；（4）每個節(jié)點(diǎn)只能被一個服務(wù)平臺管轄；（5）一個平臺的警力最多封鎖一個路口。三、符號說明賦權(quán)連通圖頂點(diǎn)i到頂點(diǎn)j的距離頂點(diǎn)i到頂點(diǎn)j的距離0-1整型變量第i個交巡警服務(wù)平臺調(diào)度到第j個交通路口的情況；封鎖要道所需時間最長的

8、服務(wù)平臺的出警時間中心城區(qū)A的第個交叉路口節(jié)點(diǎn)的發(fā)案率；第i個交巡警服務(wù)平臺到第j個交通路口的最短距離距離的偏差限四、問題分析交巡警服務(wù)平臺位置的選取以及劃分交巡警服務(wù)平臺的管轄范圍是一件非常重要的事情。由于種種原因，現(xiàn)在交巡警服務(wù)平臺的選址及管轄區(qū)域劃分大多根據(jù)經(jīng)驗(yàn)進(jìn)行。缺乏科學(xué)系統(tǒng)的位置選取與管轄區(qū)劃分，造成了警務(wù)資源浪費(fèi)、突發(fā)事件處理不及時等多種問題。另外，警務(wù)資源是有限的，設(shè)置交巡警服務(wù)平臺也要耗費(fèi)大量的資源。所以如何科學(xué)的規(guī)劃管轄區(qū)，合理的設(shè)立新的交巡警服務(wù)平臺具有重大的意義。本文意在根據(jù)現(xiàn)有的城區(qū)道路圖與交巡警服務(wù)平臺位置，根據(jù)到達(dá)突發(fā)事件地點(diǎn)使用時間最少原則進(jìn)行交巡警服務(wù)平臺管轄

9、范圍劃分，根據(jù)使用資金盡量少及各交巡警服務(wù)平臺的工作量盡量一致的原則設(shè)立新的交巡警服務(wù)平臺。在問題一中有三個子問題需要解決。（1） 要對20個交巡警服務(wù)平臺分配管轄范圍，使其在所管轄的范圍內(nèi)出現(xiàn)突發(fā)事件時，盡量能在3分鐘內(nèi)有交巡警（警車的時速為60km/h）到達(dá)事發(fā)地。即計(jì)算各交巡警服務(wù)平臺與各個地點(diǎn)的距離。將巡警在3分鐘內(nèi)到達(dá)事發(fā)地轉(zhuǎn)化為交巡警服務(wù)平臺距離事發(fā)地距離不超過3km。這是典型的最短路模型。對于這個問題，我們采用算法。（2） 當(dāng)重大突發(fā)事件發(fā)生后，要對中心城區(qū)A的20個交巡警服務(wù)平臺的警力資源進(jìn)行調(diào)度，對進(jìn)出該區(qū)的13條交通要道實(shí)現(xiàn)快速全封鎖，其關(guān)鍵在于合理調(diào)度警務(wù)資源使得封鎖全部

10、要道所需的總時間達(dá)到最小，也就是使得出警時間最長的服務(wù)平臺所需的時間盡可能的小。實(shí)際中一個平臺的警力最多封鎖一個路口，給出該區(qū)交巡警服務(wù)平臺警力合理的調(diào)度方案，我們采用0-1模型進(jìn)行求解，給出該區(qū)交巡警服務(wù)平臺警力合理的調(diào)度方案。（3） 針對現(xiàn)有的中心城區(qū)A的20個交巡警服務(wù)平臺進(jìn)行分析后，需要新增加25個服務(wù)平臺以解決工作量不平衡和部分路口節(jié)點(diǎn)出警時間過長的問題。這屬于交巡警服務(wù)平臺選址問題。一方面考慮采用集合覆蓋模型，目的是在滿足所有節(jié)點(diǎn)3分鐘內(nèi)都有警方到達(dá)的條件下，使新增設(shè)的服務(wù)平臺數(shù)目盡可能得小，從而降低了建設(shè)成本。另一方面也要考慮新增設(shè)服務(wù)平臺后，能夠解決服務(wù)平臺工作量不平衡的問題，

11、所以把盡可能均衡各個服務(wù)平臺工作量作為第二個目標(biāo) 。因此考慮需要建立一個兩目標(biāo)的0-1整數(shù)規(guī)劃模型。五、模型的建立和求解問題1.1A區(qū)交巡警服務(wù)平臺管轄范圍的分配為各交巡警服務(wù)平臺分配管轄范圍，使其在所管轄的范圍內(nèi)出現(xiàn)突發(fā)事件時，盡量能在3分鐘內(nèi)有交巡警（警車的時速為60km/h）到達(dá)事發(fā)地。即計(jì)算各交巡警服務(wù)平臺與各個地點(diǎn)的距離。將巡警在3分鐘內(nèi)到達(dá)事發(fā)地轉(zhuǎn)化為交巡警服務(wù)平臺距離事發(fā)地距離不超過3km。我們建立最短路模型：，即求出離每個地點(diǎn)距離最短的交巡警服務(wù)平臺，將這個地點(diǎn)分配給距離最短的交巡警平臺管轄。如何求每個地點(diǎn)之間的距離，我們運(yùn)用了軟件中的最短路函數(shù)“（）”函數(shù)。將題目附件中各節(jié)點(diǎn)

12、坐標(biāo)轉(zhuǎn)化為矩陣，進(jìn)而建立稀疏矩陣，然后利用最短路函數(shù)提取出符合3分鐘路程要求的矩陣，最后進(jìn)行整理，即可得所求。模型的求解：先用算法求解出各交巡警服務(wù)平臺到各個路口節(jié)點(diǎn)的最短距離，利用軟件進(jìn)行運(yùn)算（運(yùn)算程序見附件一），代入數(shù)據(jù),求得管轄范圍如下：表一 A區(qū)各服務(wù)平臺的管轄范圍服務(wù)平臺節(jié)點(diǎn)編號管轄的節(jié)點(diǎn)11、42、44、79、8022、66、74、75、76、7833、43、54、55、68、7044、57、58、60、62、6355、49、5366、48、50、51、52、56、5977、30、3488、35、36、37、47、99、3210101111、25、26、2712121313、21、

13、22、23、2414141515、311616、33、45、461717、40、411818、71、72、77、84、85、87、88、89、911919、64、65、67、69、73、82、832020、81、86（原始表格見附件二）把出警時間不超過3分鐘，轉(zhuǎn)化為服務(wù)平臺距離所管轄的路口距離不超過3千米。由此檢驗(yàn)得六個路口（ 28，29，38，39，61，92）不滿足出警時間要求。到達(dá)最近的交巡警服務(wù)平臺的時間超過3分鐘的6個路口如下表：路口標(biāo)號282938396192歸屬標(biāo)號1515162720距離47.58150.00534.05936.82241.90236.013問題1.2A區(qū)交巡警

14、服務(wù)平臺警力調(diào)度方案對于重大突發(fā)事件，需要調(diào)度全區(qū)20個交巡警服務(wù)平臺的警力資源，對進(jìn)出該區(qū)的13條交通要道實(shí)現(xiàn)快速全封鎖。要求得該如何分配警力資源封鎖路口我們使用0-1整數(shù)規(guī)劃模型。設(shè)表示第個交巡警服務(wù)平臺調(diào)度到第j個交通路口的情況，即：其中根據(jù)對問題的分析，要實(shí)現(xiàn)對要道的快速全封鎖，所以模型的目標(biāo)是使封鎖所有要道的總時間最短，其關(guān)鍵在于控制封鎖要道所需時間最長的服務(wù)平臺的出警時間，使之達(dá)到最小值。設(shè)表示封鎖要道所需時間最長的服務(wù)平臺的出警時間Aij表示20個交警平臺到13個交通路口的距離封鎖要道所需時間最長的服務(wù)平臺的出警時間最短：min ；在13個交通路口上，每個路口都必須有一個交巡警：

15、 ；每個交巡警服務(wù)平臺至多只能去一個路口：；每個交巡警到達(dá)路口的距離均小于最后一個到達(dá)路口的交警平臺與該路口的距離： 。建立模型如下： 目標(biāo)函數(shù)： min ;：模型的求解：我們利用lingo軟件進(jìn)行目標(biāo)函數(shù)的求解。主程序如下：（具體程序見附件三）由程序結(jié)果可得如下表格：封鎖方案表格距離（km）節(jié)點(diǎn)走向時間（m）距離（km）節(jié)點(diǎn)走向時間（m）2.476482.472.608162.603.105303.105.409105.408.01.7298.017.7010227.704.7015284.703.8011243.803.982383.986.7416146.740.354620.353.2

16、614213.26012120下圖是所求得各巡邏點(diǎn)到個路口的時間表：取其中的最大值為8.01546.故控制封鎖要道所需時間最短約為8.01分鐘。問題1.3擬增A區(qū)交巡警服務(wù)平臺對于擬增A區(qū)交巡警服務(wù)平臺，首先分析現(xiàn)有的交巡警服務(wù)平臺的分布，發(fā)現(xiàn)存在交巡警服務(wù)平臺工作量不均衡和部分交通路口出警時間過長的問題。這就要求新增加幾個服務(wù)平臺后，使得各個交巡警服務(wù)平臺的工作量盡可能相同以及使各個交通路口出警時間都被控制在3分鐘內(nèi)。新建服務(wù)平臺需要成本，所以需要合理確定服務(wù)平臺的選址，使需要建立的服務(wù)平臺的數(shù)目最小。由此，我們參考集合覆蓋模型（在一定的區(qū)域內(nèi)，設(shè)置最小數(shù)量的設(shè)施來覆蓋其中所有的點(diǎn)），建立了

17、一個兩目標(biāo)0-1整數(shù)規(guī)劃模型。（見附件四）目標(biāo)函數(shù)： 約束條件：(i=1,2,92 =1,2,92)且 每個交通路口都有一個平臺管轄： 所有平臺到其所管轄的交通路口最短距離中的最大值不超過距離偏差,即各個交通路口出警時間都被控制在3分鐘內(nèi)： (i=1,2,92 j=1,2,92)模型的求解：第三個子問題所建立的是雙目標(biāo)的0-1整數(shù)規(guī)劃模型，第一目標(biāo)為增加的服務(wù)平臺最少，第二目標(biāo)為各個服務(wù)平臺每天的服務(wù)強(qiáng)度方差最小。為了求解模型，先只考慮第一目標(biāo)，然后再在第一目標(biāo)最優(yōu)的情況下給出第二目標(biāo)最優(yōu)的解。在lingo軟件中運(yùn)行程序，從24至30范圍內(nèi)取出若干個偏差限與所對應(yīng)的目標(biāo)函數(shù)值。由結(jié)果，增加5個

18、平臺，標(biāo)號與坐標(biāo)分別為：29（246,337）;39（371,333）;48（315,374）;51（348.5,380.5）;88（444.5,383）由于該模型發(fā)案率的均衡性不是很好，增加的5個平臺不一定是劃分A區(qū)的最優(yōu)解。因此要進(jìn)行進(jìn)一步的分析處理。將25個平臺代入1.1的模型中，與原A區(qū)的劃分結(jié)果相比較（以dij>30的節(jié)點(diǎn)個數(shù)的多少為標(biāo)準(zhǔn)），逐個去掉平臺個數(shù)，再重復(fù)以上操作，最終得出結(jié)果。使用和求A區(qū)偏差限相同的方法分別確定增加5個平臺、4個平臺、3個平臺時的偏差限，即最好均衡性分別為1.9、1.85、1.85，并設(shè)計(jì)表格進(jìn)行對比，表格如下：增加的點(diǎn)數(shù)均衡性最優(yōu)解時間超過3分鐘

19、的路口數(shù)增加的平臺51.9886.2554529 39 48 51 8841.85954.6909429 39 48 8831.851170.171429 39 48權(quán)衡均衡性、最優(yōu)解及增加點(diǎn)個數(shù)所需花費(fèi)的資源，增加四個點(diǎn)數(shù)是最節(jié)省資源且效果最好的。故：交巡警服務(wù)平臺增加四個，它們分別為：A29 A39 A48 A88六、模型的評價(jià)6.1模型的優(yōu)點(diǎn)：（1） 對于問題一的第一個子問題：題目規(guī)定以最短時間為目標(biāo)，可轉(zhuǎn)化為最短距離不超過3km。可是對于有一些節(jié)點(diǎn)，無論怎樣安排都不能達(dá)到警方到達(dá)事故現(xiàn)場為3分鐘之內(nèi)。我們對于這些節(jié)點(diǎn)單獨(dú)處理，并根據(jù)已知條件計(jì)算出這幾個節(jié)點(diǎn)應(yīng)該屬于哪個交巡警服務(wù)平臺。（

20、2） 對于問題一的第三個子問題：該問題要求確定增加2-5個交巡警服務(wù)平臺的方案。一方面，根據(jù)限制條件，巡警到達(dá)其管轄節(jié)點(diǎn)應(yīng)在3分鐘之內(nèi)。另一方面，應(yīng)平均合理安排各個交巡警服務(wù)平臺的工作量，不至于出現(xiàn)某個服務(wù)平臺工作量過大而某個平臺工作量過小的情況。6.2模型的缺點(diǎn)：對于問題一的第一個子問題，建立的模型目標(biāo)單一，并沒有考慮到各個服務(wù)平臺交巡警工作梁的不均衡性。在分配結(jié)果中可以看到，部分交巡警的工作量很小而部分交巡警的工作量很大。這顯然是不合理的。事實(shí)上，可以建立一個多目標(biāo)規(guī)劃問題，以交巡警到達(dá)事故現(xiàn)場和工作量的方差為目標(biāo)，使得分配管轄范圍的方案更加合理。七、模型的推廣本文所提供的模型可以推廣到任

21、何一個市區(qū)甚至更廣范圍內(nèi)的交巡警服務(wù)平臺的設(shè)置與調(diào)度問題的解決中。也可以廣泛應(yīng)用于社區(qū)衛(wèi)生室、公共衛(wèi)生間、消防救火中心等社會服務(wù)部門的選址問題。在實(shí)際運(yùn)用中有很高的實(shí)用性及通用性。八、參考文獻(xiàn)1殷代君，廣義最大覆蓋模型在應(yīng)急設(shè)施選址中的應(yīng)用研究，中外企業(yè)家，2010年第3期：169-170,2010。2卓金武，MATLAB在數(shù)學(xué)建模中的應(yīng)用，北京：北京航空航天大學(xué)出版社，2011。3韓中庚，數(shù)學(xué)建模方法及其應(yīng)用，北京：高等教育出版社，2009。4謝金星，薛毅，優(yōu)化建模與lingo軟件，北京：清華大學(xué)出版社，2005。附件一：function list,AdjMatrix,distance =

22、getList()edge = csvread('edge.csv');range = size(edge);%刪除非A區(qū)的邊f(xié)or i =range:-1:1 if (edge(i,1)>92)|(edge(i,2)>92) edge(i,:) = ; endend node = csvread('node.csv');S = edge(:,1);%起點(diǎn)向量E = edge(:,2);%終點(diǎn)向量 range = size(S);W = zeros(range,1);%權(quán)向量for i=1:1:range a = node(S(i,1),1),nod

23、e(S(i,1),2); b = node(E(i,1),1),node(E(i,1),2); distance = norm(a-b); W(i,1) = distance;end%建立稀疏矩陣AdjMatrix = sparse(S;E,E;S,W;W,92,92);%求最短路徑矩陣distance = graphallshortestpaths(AdjMatrix,'directed',false);%提取出符合3分鐘路程要求的矩陣A = distance(1:20,21:92);rangeI,rangeJ = size(A);for i=1:1:rangeI for j

24、=1:1:rangeJ if A(i,j)>30 A(i,j) = NaN; end endend %整理list = zeros(72,1);for i=1:1:72 Max = 0; R = NaN; for j=1:1:20 if A(j,i)>Max Max = A(j,i); R = j; end end list(i,1) = R;endT = linspace(1,20,20)'list = T;list;end附件三：title 0-1; sets: p/1.13/:a; l/1.20/:b; link(l,p):c,d; endsets data: d=2

25、22.36,160.28,92.87,192.93,210.96,225.02,228.93,190.01,195.16,120.83,58.81,118.50,48.85, 204.64,141.30,73.88,173.95,191.97,206.03,211.21,172.29,177.44,103.11,39.82,103.10,60.35, 183.52,127.67,60.26,160.32,178.35,192.41,190.09,151.17,156.32,82.00,60.94,81.98,43.93, 219.97,150.09,82.67,182.73,200.76,21

26、4.82,226.54,162.27,155.35,81.03,48.61,73.96,3.50, 176.28,129.70,62.28,162.35,177.50,191.55,182.85,113.07,106.15,31.83,94.21,24.76,52.55, 176.59,130.00,62.59,162.65,177.80,191.86,183.16,113.37,106.46,32.14,94.52,25.06,53.37, 149.15,109.01,41.60,141.66,150.36,164.42,155.72,85.70,80.15,5.83,73.53,12.90

27、,79.92, 140.93,94.34,26.92,126.99,142.14,156.19,147.50,102.28,104.93,30.61,58.85,30.99,86.77, 130.11,82.74,15.33,115.39,131.32,145.38,136.68,97.76,107.24,34.92,47.26,41.99,93.37, 75.87,127.76,69.57,95.11,77.08,91.13,82.44,141.95,151.44,79.11,101.50,86.19,147.61, 37.91,83.37,113.95,50.72,32.70,46.75,

28、38.05,186.33,195.82,123.50,145.88,130.57,191.99, 0.00,119.50,145.43,86.85,68.83,64.77,35.92,217.81,227.30,154.98,177.36,162.05,223.47, 59.77,59.73,127.15,27.08,9.06,5.00,23.85,228.08,237.57,165.25,161.21,172.32,213.32, 119.50,0.00,67.42,32.65,50.68,64.73,83.59,180.50,189.17,114.84,101.48,121.91,153.59, 170.30,132.98,65.56,165.63,171.51,185.56,176.87,47.52,57.01,44.01,97.50,51.09,118.10, 145.43,67.42,0.00,100.07,118.09,132.15,151.00,113.08,121.75,47.43,3