日韩av黄I国产麻豆传媒I国产91av视频在线观看I日韩一区二区三区在线看I美女国产在线I麻豆视频国产在线观看I成人黄色短片

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 >

Python小白的数学建模课-07.选址问题

發布時間:2025/3/15 36 豆豆
生活随笔 收集整理的這篇文章主要介紹了 Python小白的数学建模课-07.选址问题 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

選址問題是要選擇設施位置使目標達到最優,是數模競賽中的常見題型。

小白不一定要掌握所有的選址問題,但要能判斷是哪一類問題,用哪個模型。

進一步學習 PuLP工具包中處理復雜問題的字典格式快捷建模方法。

歡迎關注『Python小白的數學建模課 @ Youcans』系列,每周持續更新



1. 選址問題

選址問題是指在某個區域內選擇設施的位置使所需的目標達到最優。選址問題也是一種互斥的計劃問題。

例如投資場所的選址:企業要在 m 個候選位置選擇若干個建廠,已知建廠費用、運輸費及 n 個地區的產品需求量,應如何進行選址。

選址問題是運籌學中經典的問題之一,選址問題在生產生活、物流、甚至軍事中都有著非常廣泛的應用,如工廠、倉庫、急救中心、消防站、垃圾處理中心、物流中心、導彈倉庫的選址等。更重要的,選址問題也是數模競賽的熱點問題。

選址是重要的長期決策,選址的好壞直接影響到服務方式、服務質量、服務效率、服務成本等,從而影響到利潤和市場競爭力,選址問題的研究有著重大的經濟、社會和軍事意義。

選址問題有四個基本要素:設施、區域、距離和優化目標。

歡迎關注『Python小白的數學建模課 @ Youcans』系列,每周持續更新
Python小白的數學建模課-01.新手必讀
Python小白的數學建模課-02.數據導入
Python小白的數學建模課-03.線性規劃
Python小白的數學建模課-04.整數規劃
Python小白的數學建模課-05.0-1規劃
Python小白的數學建模課-06.固定費用問題
Python小白的數學建模課-07.選址問題
Python小白的數學建模課-09.微分方程模型
Python小白的數學建模課-10.微分方程邊值問題
Python小白的數學建模課-12.非線性規劃
Python小白的數學建模課-15.圖論的基本概念
Python小白的數學建模課-16.最短路徑算法
Python小白的數學建模課-17.條件最短路徑算法


1.1 設施

選址問題加粗樣式中所說的設施,在具體題目中可以是工廠、倉庫、服務站等形式。

1.2 區域

選址問題中所說的區域,在具體題目中可以是工廠、車間的內部布局,也可以是給定的某個地區、甚至空間范圍。

按照規劃區域的特征,可以分為連續選址問題和離散選址問題。連續選址問題,設施可以布局在區域內的任意位置,就要求出最優選址的坐標;離散選址問題,只能從若干候選位置中進行選擇,運籌學中的選址問題通常是這類離散選址問題。

1.3 距離

選址問題中所說的距離,是指設施到服務對象之間的距離,在具體題目中也可以是某個選址位置的服務時間、成本、覆蓋范圍。如果用圖論方法求解,通常就是連接頂點的邊的權值。

當問題所關注的是設施到服務對象之間的距離時,如果問題給出的不是頂點之間的距離,而是設施的位置坐標,要注意不是只有歐式距離,對于不同問題也可能是球面距離、曼哈頓距離、切比雪夫距離。

1.4 優化目標

選址問題要求選擇最好的選址位置,但選址位置只是決策變量,選擇的最終目的通常是實現加權距離最短、費用最小、利潤最大、時間最短,這才是優化問題的目標函數。

按照目標函數的特點,可以分為:中位問題,要求總成本最小;中心問題,服務于每個客戶的最大成本最小;反中心問題:服務于每個客戶的最小成本最大。



2. 常見選址問題及建模

2.1 P-中位問題(P-median problem)

P-中位問題,假設有 N 個候選服務站和 M 個需求點,要從 N 個候選服務站中選擇 P 個,使所有需求點到最近的服務站的加權距離 dijd_{ij}dij? 的總和最小。需求點 i 的權值,通常是指該需求點的需求量。

這是一個 MinSum 問題,定義決策變量 xjx_jxj? 為選中的服務站,yijy_{ij}yij? 將各需求點匹配到最近的服務站:
xj={1,服務站j被選中0,服務站j未被選中x_j = \begin{cases} 1,& 服務站\ j\ 被選中\\ 0,& 服務站\ j\ 未被選中 \end{cases} xj?={10??j??j??
yij={1,需求點i由服務站j服務0,需求點i不由服務站j服務y_{ij} = \begin{cases} 1,& 需求點\ i\ 由服務站\ j\ 服務\\ 0,& 需求點\ i\ 不由服務站\ j\ 服務 \end{cases} yij?={10??i??j??i??j??
可以建立數學模型如下:

min∑i∈M∑j∈Nwidijyijs.t.:{∑j∈Nxj=P∑j∈Nyij=1,?iyij?xj≤0,?i,jxj∈{0,1},yij∈{0,1}min\;\sum_{i \in M}\sum_{j \in N} w_i d_{ij}y_{ij}\\ s.t.:\;\begin{cases} \sum_{j \in N} x_{j} = P\\ \sum_{j \in N} y_{ij} = 1,\forall i\\ y_{ij} - x_j \leq 0,\forall i,j\\ x_{j} \in \{0,1\}, \;y_{ij} \in \{0,1\} \end{cases} miniM?jN?wi?dij?yij?s.t.:??????????jN?xj?=PjN?yij?=1?iyij??xj?0?i,jxj?{0,1},yij?{0,1}?

其中:j 為服務站,i 為需求點,wiw_iwi? 為需求點 i 的需求量, dijd_{ij}dij? 為需求點 i 到服務站 j 的距離。


2.2 P-中心問題

P-中心問題,假設有 N 個候選服務站和 M 個需求點,要從 N個候選服務站中選擇 P個,使任一需求點到最近的服務站的最大距離最小。

這是一個 MinMax 問題,需要最小化任何需求點與其鄰近設施點的最大距離。P-中位問題追求總和最小,可以理解為發展經濟總量優先;P-中心問題關注最差個體的最好結果,可以理解為優先進行扶貧。

定義決策變量 xjx_jxj? 為選中的服務站,yijy_{ij}yij? 將各需求點匹配到最近的服務站:
xj={1,服務站j被選中0,服務站j未被選中x_j = \begin{cases} 1,& 服務站\ j\ 被選中\\ 0,& 服務站\ j\ 未被選中\end{cases} xj?={10??j??j??

yij={1,需求點i由服務站j服務0,需求點i不由服務站j服務y_{ij} = \begin{cases} 1,& 需求點\ i\ 由服務站\ j\ 服務\\ 0,& 需求點\ i\ 不由服務站\ j\ 服務\end{cases} yij?={10??i??j??i??j??

可以建立數學模型如下:

minDs.t.:{∑j∈Nwidijyij≤D,?i∑j∈Nxj=P∑j∈Nyij=1,?iyij?xj≤0,?i,jxj∈{0,1},yij∈{0,1}min\; D\\ s.t.:\;\begin{cases} \sum_{j \in N} w_i d_{ij} y_{ij} \leq D, \forall i\\ \sum_{j \in N} x_{j} = P\\ \sum_{j \in N} y_{ij} = 1, \forall i\\ y_{ij} - x_j \leq 0, \forall i,j\\ x_{j} \in \{0,1\}, \;y_{ij} \in \{0,1\} \end{cases} minDs.t.:????????????????jN?wi?dij?yij?D,?ijN?xj?=PjN?yij?=1,?iyij??xj?0,?i,jxj?{0,1},yij?{0,1}?

其中:j 為服務站,i 為需求點, dijd_{ij}dij? 為需求點 i 到服務站 j 的距離。如果只求需求點到最近的服務站的最大距離,則 wi=1w_i = 1wi?=1 ;如果要求任一需求點到最近的服務站的最大運費,則 wiw_iwi? 為需求點 i 的需求量,即加權最大距離。


2.3 集合覆蓋問題

覆蓋模型適用于一些特殊場景,例如消防中心、救護車、巡邏車等應急設施的區位選址問題。覆蓋問題分為集合覆蓋問題(Set covering problem)和最大覆蓋問題(Maximal covering problem)。

集合覆蓋問題研究滿足覆蓋所有需求點顧客的前提下,服務站的最少個數或建設費用最小的問題。假設有 N 個候選服務站和 M 個需求點,已知每個服務站的服務范圍(或服務容量),要從 N個候選服務站中選擇若干個,使所有需求點得到服務(到所屬服務站的距離或時間小于給定的臨界值),服務站的個數最少或成本最小。

定義參數 aija_{ij}aij? 為每個服務站的覆蓋范圍:

aij={1,服務站j可以覆蓋需求點i0,服務站j不能覆蓋需求點ia_{ij} = \begin{cases} 1,& 服務站\ j\ 可以覆蓋需求點\ i\\ 0,& 服務站\ j\ 不能覆蓋需求點\ i \end{cases} aij?={10??j??i?j??i?

定義決策變量 xjx_jxj? 為選中的服務站:
xj={1,服務站j被選中0,服務站j未被選中x_j = \begin{cases} 1,& 服務站\ j\ 被選中\\ 0,& 服務站\ j\ 未被選中\end{cases} xj?={10??j??j??

可以建立數學模型如下:

min∑j∈Ncjxjs.t.:{∑j∈Nixj≥1,?i∈Mxj∈{0,1}min\; \sum_{j \in N} c_j x_{j}\\ s.t.:\;\begin{cases} \sum_{j \in N_i} x_j \geq 1, \forall i \in M\\ x_{j} \in \{0,1\} \end{cases} minjN?cj?xj?s.t.:{jNi??xj?1,?iMxj?{0,1}?

其中:j 為服務站,i 為需求點,cjc_jcj? 為服務站 j 的建設費用(最少個數問題中不需要考慮),Ni={j:aij=1}N_i=\{j:a_{ij}=1\}Ni?={j:aij?=1} 是覆蓋需求點 i 的候選服務站的集合。


2.4 最大覆蓋問題

最大覆蓋問題研究在已知服務站的數目和服務半徑的條件下,如何設立 P個服務站使得可接受的服務需求最大的問題。

定義決策變量 xjx_jxj? 為選中的服務站:
xj={1,服務站j被選中0,服務站j未被選中x_j = \begin{cases} 1,& 服務站\ j\ 被選中\\ 0,& 服務站\ j\ 未被選中 \end{cases} xj?={10??j??j??
KaTeX parse error: Undefined control sequence: \ at position 57: …,& 需求點\ i\ 未被覆蓋\? ?\end{cases}

可以建立數學模型如下:

max∑i∈Niwizis.t.:{zi≤∑j∈Nixj,?i∈M∑j∈Mxj=pxj∈{0,1},zi∈{0,1}max\; \sum_{i \in N_i} w_i z_i\\ s.t.:\; \begin{cases} z_i \leq \sum_{j \in N_i} x_j , \forall i \in M\\ \sum_{j \in M} x_j = p\\ x_{j} \in \{0,1\}, z_{i} \in \{0,1\} \end{cases} maxiNi??wi?zi?s.t.:??????zi?jNi??xj?,?iMjM?xj?=pxj?{0,1},zi?{0,1}?

其中:j 為服務站,i 為需求點,wiw_iwi? 為需求點 i 的需求量。


![在這里插入圖片描述](https://img-blog.csdnimg.cn/e808c4a986574475b00d925506c28cdb.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3lvdWNhbnM=,size_16,color_FFFFFF,t_70#pic_center)

2.5 其它選址問題

其它選址問題,在數學建模中應用相對較少,限于篇幅不能逐一介紹其數學模型。在此將各模型的特點簡要介紹,以便判斷問題的類型。

帶固定費用和容量限制的選址問題

服務站建站的固定費用和服務站的容量(能力)限制這兩個因素具有很強的實際意義,經常作為基本選址問題的深化研究課題。
無容量限制的固定費用下的選址問題,就是去掉服務站個數的約束,并將固定建站費用加到 P-中位問題的目標函數上。

選址分配問題

選址分配問題類似于 P-中位問題,有 m 個服務站需要選址,n 個已知位置的顧客分配給不同的設施,已知每個服務站的能力和每個顧客的需求,要求服務站的選址和顧客對服務站的分配,使顧客與所分配服務站的距離總和最小。

隨機選址問題
服務站的運行時間、建設成本、需求點位置、需求數量等全部或部分參數是不確定的,但服從某種隨機分布。

動態選址問題
研究未來若干時間段內服務站的最優選址問題,在不同時間段內動態選址模型的參數是變化的,但在某一時間段內模型參數是確定的。

競爭選址問題
研究考慮市場上存在兩個以上同類產品或服務的提供者,或服務站提供多個產品或服務。



2.6 選址問題的求解算法與編程實現

設施選址問題通常是是 NP 問題,不存在多項式時間算法。常用的近似解法有:

線性規劃舍入算法,相當于整數規劃問題的求解算法。首先給出原問題的整數規劃模型,然后求解相應的線性規劃松弛問題得到分數最優解,根據可行要求對分數最優解進行改造,構造原問題的整數可行解。

原始對偶算法,首先找到對偶問題的一個可行解,再根據該對偶可行解構造原始問題的整數可行解,不斷調整對偶問題的可行解,直到找到最優解為止。

局部搜索算法:給定初始可行解,定義適當的鄰域,通過引入恰當的調整策略,在鄰域中得到改進的可行解,依次迭代,直到調整策略不能改進為止

啟發式算法或隨機優化算法。

本節作為線性規劃問題系列的一篇,仍然選擇 PuLP工具包求解選址問題。很多選址問題適合用圖論方法描述和求解,這將在后續課程中進行介紹。



3. 案例 1:PuLP求解指派問題

說明:本案例是指派問題,不是選址問題。因指派問題未單獨成文,因此將該案例放在本文中。

另外,本案例給出了 PuLP 工具包使用字典方式快捷編程的使用方法,這在選址問題中是非常方便的。

3.1 游泳接力賽的指派問題

游泳隊中 A、B、C、D 四名運動員組成 4x100米混合泳接力隊,運動員各種泳姿的成績如下表所示 (單位:秒):

隊員\項目自由泳蛙泳蝶泳仰泳
A56746163
B63696571
C57776367
D55766262

如何安排 A、B、C、D 四名運動員的泳姿,才最有可能取得好成績?


3.2 指派問題建模分析

引入 0-1 變量 xijx_{ij}xij?

xi,j={0,第i人不游第j種姿勢1,第i人游第j種姿勢,i,j=1,...4x_{i,j} = \begin{cases} 0,第\;i\;人不游第\;j\;種姿勢\\ 1,第\;i\;人游第\;j\;種姿勢,i,j=1,...4 \end{cases} xi,j?={0ij姿1ij姿i,j=1,...4?

指派問題的數學模型就可以描述為:
minf(x)=∑i=1n∑j=1n(cijxij)s.t.:{∑j=1nxij=1,i=1,...,n∑i=1nxij=1,j=1,...,nxij=0,1,i,j=1,...,nmin\;f(x) = \sum_{i=1} ^n \sum_{j=1} ^n (c_{ij} x_{ij})\\ s.t.:\;\begin{cases} \sum_{j=1} ^n x_{ij} = 1,i=1,...,n\\ \sum_{i=1} ^n x_{ij} = 1,j=1,...,n\\ x_{ij} = 0,1,i,j=1,...,n \end{cases} minf(x)=i=1n?j=1n?(cij?xij?)s.t.:??????j=1n?xij?=1i=1,...,ni=1n?xij?=1j=1,...,nxij?=0,1i,j=1,...,n?

其中:

ci,j=(56746163636965715777636755766262)c_{i,j}=\left( \begin{matrix} 56 & 74 & 61 & 63 \\ 63 & 69 & 65 & 71 \\ 57 & 77 & 63 & 67 \\ 55 & 76 & 62 & 62 \end{matrix} \right) ci,j?=?????56635755?74697776?61656362?63716762??????

3.3 指派問題模型求解的編程

模型求解,用標準模型的優化算法對模型求解,得到優化結果。模型求解的編程步驟如下:

(0)導入 PuLP庫函數

import pulp

(1)定義一個規劃問題

AssignLP = pulp.LpProblem("Assignment_problem_for_swimming_relay_race", sense=pulp.LpMinimize)

pulp.LpProblem 用來定義問題的構造函數。參數 sense 指定問題求目標函數的最小值/最大值 。

(2)定義決策變量

rows = cols = range(0, 4)x = pulp.LpVariable.dicts("x", (rows, cols), cat="Binary")

pulp.LpVariable 用來定義決策變量的函數。參數 cat 設定變量類型,’ Binary ’ 表示 0/1 變量。

注意,指派問題、選址問題中都涉及 N*M 維矩陣變量,變量個數很多,如果逐一定義非常冗長,而且容易出錯、不便修改。本例使用 pulp.LpVariable.dicts 提供的字典格式定義了 4*4 個變量 xijx_{ij}xij?,使程序大為簡化。

(3)添加目標函數

scoreM = [[56,74,61,63],[63,69,65,71],[57,77,63,67],[55,76,62,62]]AssignLP += pulp.lpSum([[x[row][col]*scoreM[row][col] for row in rows] for col in cols])

本例程在語句內使用兩重 for 循環遍歷列表實現所有變量的線性組合 ,使程序大為簡化。

(4)添加約束條件

for row in rows:AssignLP += pulp.lpSum([x[row][col] for col in cols]) == 1 # sum(x(i,j),j=1,4)=1, i=1,4for col in cols:AssignLP += pulp.lpSum([x[row][col] for row in rows]) == 1 # sum(x(i,j),i=1,4)=1, j=1,4

快捷方法對于約束條件的定義與對目標函數的定義相似,使用字典定義參數,使用循環定義約束條件,使程序簡單、結構清楚。

(5)求解和結果輸出

AssignLP.solve() # youcansprint(AssignLP.name)member = ["隊員A","隊員B","隊員C","隊員D"]style = ["自由泳","蛙泳","蝶泳","仰泳"]if pulp.LpStatus[AssignLP.status] == "Optimal": # 獲得最優解xValue = [v.varValue for v in AssignLP.variables()]# [0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0]xOpt = np.array(xValue).reshape((4, 4)) # 將 xValue 格式轉換為 4x4 矩陣print("最佳分配:" )for row in rows:print("{}\t{} 參加項目:{}".format(xOpt[row],member[row],style[np.argmax(xOpt[row])]))print("預測最好成績為:{}".format(pulp.value(AssignLP.objective)))

xValue 獲得的是列表變量,通過 numpy 的 reshape() 函數轉換為 4*4 矩陣,便于格式化輸出。


3.4 指派問題 Python 例程

# mathmodel08_v1.py # Demo08 of mathematical modeling algorithm # Solving assignment problem with PuLP. # Copyright 2021 Youcans, XUPT # Crated:2021-06-02 # Python小白的數學建模課 @ Youcansimport pulp # 導入 pulp 庫 import numpy as np# 主程序 def main():# 問題建模:"""決策變量:x(i,j) = 0, 第 i 個人不游第 j 種姿勢x(i,j) = 1, 第 i 個人游第 j 種姿勢i=1,4, j=1,4目標函數:min time = sum(sum(c(i,j)*x(i,j))), i=1,4, j=1,4約束條件:sum(x(i,j),j=1,4)=1, i=1,4sum(x(i,j),i=1,4)=1, j=1,4變量取值范圍:x(i,j) = 0,1 """# 游泳比賽的指派問題 (assignment problem)# 1.建立優化問題 AssignLP: 求最小值(LpMinimize)AssignLP = pulp.LpProblem("Assignment_problem_for_swimming_relay_race", sense=pulp.LpMinimize) # 定義問題,求最小值# 2. 建立變量rows = cols = range(0, 4)x = pulp.LpVariable.dicts("x", (rows, cols), cat="Binary")# 3. 設置目標函數scoreM = [[56,74,61,63],[63,69,65,71],[57,77,63,67],[55,76,62,62]]AssignLP += pulp.lpSum([[x[row][col]*scoreM[row][col] for row in rows] for col in cols])# 4. 施加約束for row in rows:AssignLP += pulp.lpSum([x[row][col] for col in cols]) == 1 # sum(x(i,j),j=1,4)=1, i=1,4for col in cols:AssignLP += pulp.lpSum([x[row][col] for row in rows]) == 1 # sum(x(i,j),i=1,4)=1, j=1,4# 5. 求解AssignLP.solve() # youcans# 6. 打印結果print(AssignLP.name)member = ["隊員A","隊員B","隊員C","隊員D"]style = ["自由泳","蛙泳","蝶泳","仰泳"]if pulp.LpStatus[AssignLP.status] == "Optimal": # 獲得最優解xValue = [v.varValue for v in AssignLP.variables()]# [0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0]xOpt = np.array(xValue).reshape((4, 4)) # 將 xValue 格式轉換為 4x4 矩陣print("最佳分配:" )for row in rows:print("{}\t{} 參加項目:{}".format(xOpt[row],member[row],style[np.argmax(xOpt[row])]))print("預測最好成績為:{}".format(pulp.value(AssignLP.objective)))returnif __name__ == '__main__': # Copyright 2021 YouCans, XUPTmain() # Python小白的數學建模課 @ Youcans

3.5 Python 例程運行結果

Welcome to the CBC MILP Solver Version: 2.9.0 Build Date: Feb 12 2015 Result - Optimal solution foundAssignment_problem_for_swimming_relay_race 最佳分配: [0. 0. 1. 0.] 隊員A 參加項目:蝶泳 [0. 1. 0. 0.] 隊員B 參加項目:蛙泳 [1. 0. 0. 0.] 隊員C 參加項目:自由泳 [0. 0. 0. 1.] 隊員D 參加項目:仰泳 預測最好成績為:249.0

4. 案例 2:PuLP求解選址問題

4.1 消防站的選址問題

例題 2:某城市有 8 個區,每個區最多建一個消防站,擬建設消防站到各區的最長時間如下表所示?,F要求任何區域發生火警時,消防車能在 10分鐘內趕到。在此條件下盡量減少消防站數量,應該在哪幾個區建設消防站?

區域12345678
1712182024262528
214581516181818
319941410221613
41415151018151418
5201812209251412
6182120162061015
722182015161559
830221520141886

4.2 選址問題建模分析

首先判斷這是一個集合覆蓋問題,要求從 8 個候選消防站中選擇若干個,在所有需求點得到服務的時間都小于臨界值 10分鐘的條件下,選擇消防站的數量最少。本問題不考慮各候選站點建設費用的差異,即不帶權重。

定義參數 RijR_{ij}Rij? 為每個消防站的覆蓋范圍:
Rij={1,消防站j可以覆蓋區域i0,消防站j不能覆蓋區域iR_{ij} = \begin{cases} 1,& 消防站\ j\ 可以覆蓋區域\ i\\ 0,& 消防站\ j\ 不能覆蓋區域\ i \end{cases} Rij?={10??j??i?j??i?

由擬建消防站到各區的最長時間表可以得到參數 RijR_{ij}Rij? 如下表:

區域12345678
110000000
201100000
301101000
400010000
500001000
600000110
700000011
800000011

定義決策變量 xjx_jxj? 為選中的服務站:

xj={1,消防站j被選中0,消防站j未被選中x_j = \begin{cases} 1,& 消防站\ j\ 被選中\\ 0,& 消防站\ j\ 未被選中\end{cases} xj?={10??j??j??

可以建立數學模型如下:

minf(x)=∑j=18xjs.t.:{∑j=18xjRij≥1,i=1,..8xj∈{0,1},j=1,..8min\; f(x) = \sum_{j=1}^8 x_{j}\\ s.t.:\;\begin{cases} \sum_{j=1}^8 x_j R_{ij} \geq 1, &i=1,..8\\ x_{j} \in \{0,1\}, &j=1,..8 \end{cases} minf(x)=j=18?xj?s.t.:{j=18?xj?Rij?1,xj?{0,1},?i=1,..8j=1,..8?

選址問題的模型求解,用標準模型的優化算法對模型求解,得到優化結果。

模型求解的編程步驟與指派問題是一致的,且在例程中給出了詳細的注釋,就不再進行逐項解釋了。

需要注意的是,選址問題的決策變量、參數、約束條件的數量較大(N*M),如果對變量、約束條件逐個進行定義,編程過程將是非常冗長和痛苦的,因此需要使用列表、字典等快捷方式進行定義。對于更大規模的問題,模型中的數據要通過讀取數據文件獲得,就更需要采用這種方式來編程。


4.3 選址問題 Python 例程

# mathmodel09_v1.py # Demo08 of mathematical modeling algorithm # Solving set covering problem with PuLP. # Copyright 2021 Youcans, XUPT # Crated:2021-06-06 # Python小白的數學建模課 @ Youcansimport pulp # 導入 pulp 庫# 主程序 def main():# 問題建模:"""決策變量:x(j) = 0, 不選擇第 j 個消防站x(j) = 1, 選擇第 j 個消防站, j=1,8目標函數:min fx = sum(x(j)), j=1,8約束條件:sum(x(j)*R(i,j),j=1,8) >=1, i=1,8變量取值范圍:x(j) = 0,1"""# 消防站的選址問題 (set covering problem, site selection of fire station)# 1.建立優化問題 SetCoverLP: 求最小值(LpMinimize)SetCoverLP = pulp.LpProblem("SetCover_problem_for_fire_station", sense=pulp.LpMinimize) # 定義問題,求最小值# 2. 建立變量zones = list(range(8)) # 定義各區域 youcansx = pulp.LpVariable.dicts("zone", zones, cat="Binary") # 定義 0/1 變量,是否在該區域設消防站# 3. 設置目標函數SetCoverLP += pulp.lpSum([x[j] for j in range(8)]) # 設置消防站的個數# 4. 施加約束reachable = [[1, 0, 0, 0, 0, 0, 0, 0],[0, 1, 1, 0, 0, 0, 0, 0],[0, 1, 1, 0, 1, 0, 0, 0],[0, 0, 0, 1, 0, 0, 0, 0],[0, 0, 0, 0, 1, 0, 0, 0],[0, 0, 0, 0, 0, 1, 1, 0],[0, 0, 0, 0, 0, 0, 1, 1],[0, 0, 0, 0, 0, 0, 1, 1]] # 參數矩陣,第 i 消防站能否在 10分鐘內到達第 j 區域for i in range(8):SetCoverLP += pulp.lpSum([x[j]*reachable[j][i] for j in range(8)]) >= 1# 5. 求解SetCoverLP.solve()# 6. 打印結果print(SetCoverLP.name)temple = "區域 %(zone)d 的決策是:%(status)s" # 格式化輸出if pulp.LpStatus[SetCoverLP.status] == "Optimal": # 獲得最優解for i in range(8):output = {'zone': i+1, # 與問題中區域 1~8 一致'status': '建站' if x[i].varValue else '--'}print(temple % output)print("需要建立 {} 個消防站。".format(pulp.value(SetCoverLP.objective)))returnif __name__ == '__main__': # Copyright 2021 YouCans, XUPTmain() # Python小白的數學建模課 @ Youcans

4.4 Python 例程運行結果

Welcome to the CBC MILP Solver Version: 2.9.0 Build Date: Feb 12 2015 Result - Optimal solution foundSetCover_problem_for_fire_station 區域 1 的決策是:建站 區域 2 的決策是:-- 區域 3 的決策是:建站 區域 4 的決策是:建站 區域 5 的決策是:-- 區域 6 的決策是:建站 區域 7 的決策是:建站 區域 8 的決策是:-- 需要建立 5.0 個消防站

5. 小結

  • 關于規劃問題,我們從線性規劃、整數規劃、0-1規劃到一些特殊類型問題,用 5節課進行了介紹,到這里就暫告一段落了。后面根據需要,可能還會講非線性規劃,實際上主要是非線性優化問題了。
  • 雖然各種規劃問題的求解算法差別很大,但我們所用的編程實現方法都是基于 PuLP工具包,編程步驟都是一致的。
  • 本系列集中體現了與其它課程的區別,沒有展開講算法的實現步驟,而是重點講編程方法的選擇、建立模型方程的過程和編程實現的步驟,這主要是為了便于小白學習和掌握。
  • 【本節完】


    版權聲明:

    歡迎關注『Python小白的數學建模課 @ Youcans』 原創作品

    原創作品,轉載必須標注原文鏈接:(https://blog.csdn.net/youcans/article/details/117650843)。

    Copyright 2021 Youcans, XUPT

    Crated:2021-06-06


    歡迎關注 『Python小白的數學建模課 @ Youcans』 系列,持續更新
    Python小白的數學建模課-01.新手必讀
    Python小白的數學建模課-02.數據導入
    Python小白的數學建模課-03.線性規劃
    Python小白的數學建模課-04.整數規劃
    Python小白的數學建模課-05.0-1規劃
    Python小白的數學建模課-06.固定費用問題
    Python小白的數學建模課-07.選址問題
    Python小白的數學建模課-09.微分方程模型
    Python小白的數學建模課-10.微分方程邊值問題
    Python小白的數學建模課-12.非線性規劃
    Python小白的數學建模課-15.圖論的基本概念
    Python小白的數學建模課-16.最短路徑算法
    Python小白的數學建模課-17.條件最短路徑算法
    Python小白的數學建模課-18.最小生成樹問題
    Python小白的數學建模課-19.網絡流優化問題
    Python小白的數學建模課-20.網絡流優化案例
    Python小白的數學建模課-21.關鍵路徑法
    Python小白的數學建模課-22.插值方法
    Python小白的數學建模課-23.數據擬合全集
    Python小白的數學建模課-A1.國賽賽題類型分析
    Python小白的數學建模課-A2.2021年數維杯C題探討
    Python小白的數學建模課-A3.12個新冠疫情數模競賽賽題及短評
    Python小白的數學建模課-B2. 新冠疫情 SI模型
    Python小白的數學建模課-B3. 新冠疫情 SIS模型
    Python小白的數學建模課-B4. 新冠疫情 SIR模型
    Python小白的數學建模課-B5. 新冠疫情 SEIR模型
    Python小白的數學建模課-B6. 新冠疫情 SEIR改進模型
    Python數模筆記-PuLP庫
    Python數模筆記-StatsModels統計回歸
    Python數模筆記-Sklearn
    Python數模筆記-NetworkX
    Python數模筆記-模擬退火算法


    總結

    以上是生活随笔為你收集整理的Python小白的数学建模课-07.选址问题的全部內容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。

    日韩一区二区三区在线观看 | 亚洲天天看 | 91精品久久久久久 | 97香蕉久久国产在线观看 | 蜜臀av麻豆| 亚洲精品在线国产 | 日韩一级片大全 | 国产精品久久电影网 | 视频在线观看入口黄最新永久免费国产 | 亚洲码国产日韩欧美高潮在线播放 | 欧美激情一区不卡 | 日日操天天操狠狠操 | 狂野欧美激情性xxxx | 日韩免费成人 | 又粗又长又大又爽又黄少妇毛片 | 国产午夜精品av一区二区 | 美女搞黄国产视频网站 | 亚洲精品乱码白浆高清久久久久久 | 91中文字幕 | 国产日韩视频在线观看 | 国产性xxxx | 五月激情综合婷婷 | 91久久精| 激情欧美xxxx | www黄com| 9色在线视频 | 麻豆视频在线免费 | 成人毛片在线观看 | 国产婷婷一区二区 | 国产成人一区二区三区影院在线 | 亚洲天天摸日日摸天天欢 | 欧美日韩国产一区二区三区在线观看 | 日韩在线视频一区二区三区 | 超碰在线中文字幕 | 亚洲男男gaygay无套同网址 | 久久免费播放 | 97超碰在线资源 | 伊人一级 | 国产视频在线免费观看 | 国产视频1区2区3区 久久夜视频 | 日本午夜免费福利视频 | 韩日av一区二区 | 日韩影视大全 | 99久久精品午夜一区二区小说 | 国产乱码精品一区二区三区介绍 | 国产精品黑丝在线观看 | 国产偷v国产偷∨精品视频 在线草 | 久久99国产精品久久99 | 日日夜夜噜噜噜 | 81国产精品久久久久久久久久 | 日韩精品中字 | 天天操天天射天天爱 | 人人干人人搞 | 91porny九色91啦中文 | 国产一级高清视频 | 欧美日韩免费一区二区三区 | 操久久网| 久久这里只有精品首页 | 蜜臀av.com | 最新中文字幕在线观看视频 | 在线精品在线 | 91看片在线播放 | 国产美腿白丝袜足在线av | 国产精品九九九 | 日韩中字在线 | 国产精品1区2区3区 久久免费视频7 | 国产色久 | 国产精品毛片久久久久久 | 91成人免费在线视频 | 久久免费视频这里只有精品 | 韩国精品在线 | 免费视频资源 | 99c视频高清免费观看 | 91麻豆精品国产91久久久无需广告 | 九九热在线精品视频 | 日韩欧美电影在线 | 亚洲高清在线视频 | 五月综合激情婷婷 | 精品视频久久久久久 | 99久久99久久免费精品蜜臀 | 69av国产 | 狠狠色丁香婷婷综合久久片 | 在线观看国产v片 | 精品一区二区三区四区在线 | 色婷婷色 | 亚洲成人av在线播放 | 国产精品成人免费一区久久羞羞 | 国产精品入口a级 | 中文字幕 二区 | 在线日韩视频 | 18国产精品福利片久久婷 | 国产中文在线播放 | a午夜在线 | 毛片久久久 | 开心色激情网 | 涩涩资源网 | 国产精品免费大片视频 | 热久久99这里有精品 | 成人av在线看 | 天天操天天舔天天干 | 久久免费精品一区二区三区 | 国产亚洲成av片在线观看 | 久草在线视频国产 | 亚洲成人av免费 | 日日夜夜免费精品 | 天天插天天| 国产成人精品网站 | 亚洲视频 在线观看 | 日韩大片在线免费观看 | 九九视频精品在线 | 日韩欧美一区二区三区黑寡妇 | 在线播放精品一区二区三区 | 国产精品毛片一区二区 | 激情av资源网 | 亚洲中字幕 | 国产精品福利在线 | 欧美日韩在线精品一区二区 | 狠狠色丁香婷婷综合久小说久 | 密桃av在线 | 亚洲在线a | 激情五月在线 | av一区二区三区在线 | av亚洲产国偷v产偷v自拍小说 | 国产成人精品国内自产拍免费看 | 婷婷综合激情 | 久久综合色8888 | 国产成本人视频在线观看 | 国产一级二级在线 | 国产一级片久久 | 成人黄色大片 | 精品久久99| 久久综合偷偷噜噜噜色 | 久久久91精品国产一区二区精品 | 亚洲人成人在线 | 99精品成人 | 亚洲无吗天堂 | 91日韩在线播放 | 亚洲精区二区三区四区麻豆 | 激情综合亚洲精品 | 99国产精品一区 | 久久精品视频在线看 | 日日干美女| 懂色av一区二区三区蜜臀 | 五月天亚洲综合 | 操操日日| 五月婷在线观看 | 伊人五月在线 | 亚洲免费在线观看视频 | 欧美高清视频不卡网 | 狠狠干夜夜操天天爽 | 综合在线亚洲 | 国产精品一区二区久久 | 久久精品黄 | 国产精品欧美激情在线观看 | 欧美日韩国产一区二区三区 | 91麻豆精品国产91久久久久久久久 | 黄色1级毛片 | 天天av综合网| 在线综合 亚洲 欧美在线视频 | 国产 日韩 中文字幕 | 国内一级片在线观看 | 国产精品久久久久免费观看 | 91成人免费视频 | 麻豆播放 | 亚洲精品视频网站在线观看 | 免费视频资源 | 色婷婷骚婷婷 | 少妇bbb好爽| 亚洲一级片免费观看 | 色婷婷成人 | 精品99免费视频 | 亚洲天堂精品视频在线观看 | 久久人人爽av| 国产精品1区2区 | 成年人网站免费在线观看 | 国产视频中文字幕 | 国产精品久久久久久久久久久久午夜 | 免费a一级 | 国产二区精品 | 国产一级黄色片免费看 | 免费黄色网止 | 久久夜av| 国产精品理论在线观看 | 人人草人人草 | 亚洲一区二区三区毛片 | 涩涩网站在线看 | 久久www免费人成看片高清 | 国产91精品一区二区麻豆亚洲 | 国产在线精品国自产拍影院 | av网址在线播放 | 亚洲色综合 | 狠狠色丁香婷综合久久 | 欧美视频在线二区 | av在线小说 | 九九精品无码 | 91免费观看 | 国产综合婷婷 | 在线免费观看视频a | 天天操综合网站 | 在线免费观看黄色 | 美女免费视频观看网站 | 欧美日韩免费观看一区=区三区 | 天天操网站 | 成人一级片免费看 | 婷婷久久综合网 | 成人影音av| 免费中文字幕 | 久久综合婷婷国产二区高清 | 干干干操操操 | 亚洲理论电影 | 精品一区久久 | 大荫蒂欧美视频另类xxxx | 免费av观看 | 国产福利91精品一区二区三区 | av大片免费在线观看 | 日韩高清一二区 | 日韩视频免费看 | 中文字幕视频一区 | 国产精品久久久久久久毛片 | 色丁香久久 | 美女久久久久久久久久久 | 国产美女视频一区 | 国内精品免费 | 国产日韩在线播放 | 免费观看丰满少妇做爰 | 一级片免费在线 | 97国产超碰| 国产一级精品视频 | 99精品久久只有精品 | 久久av中文字幕片 | 涩涩资源网 | 久久草网站 | av免费观看网站 | 天天综合色 | 黄色精品免费 | 国产成人免费精品 | 国产裸体bbb视频 | 国内精品久久久久久久久久清纯 | 亚洲精品99久久久久中文字幕 | 国产免费片 | 开心激情网五月天 | 国产午夜精品免费一区二区三区视频 | 国产黄在线观看 | 韩日电影在线观看 | 久草免费在线视频 | 永久av免费在线观看 | 久久成人综合视频 | 日日夜夜人人精品 | 欧美另类重口 | 欧美性生活一级片 | 国产aa免费视频 | 欧美综合在线观看 | 丁香六月婷婷激情 | 欧美日韩视频在线一区 | 亚洲韩国一区二区三区 | 美女网站黄免费 | 国产成人久久 | 99av国产精品欲麻豆 | 久草视频在线资源 | 日韩一区二区免费视频 | 91精品办公室少妇高潮对白 | 手机av网站| 久久国产精品小视频 | 超碰在线公开免费 | 久久久久久久久久久久久久电影 | 天天干,狠狠干 | 激情久久婷婷 | 国产一区私人高清影院 | 国产人成一区二区三区影院 | 99久久精品免费一区 | 国产精在线 | 久久综合免费 | 国产精品久久久久久久av电影 | 91亚色视频在线观看 | 日韩精品一区二区三区高清免费 | 成人污视频在线观看 | 久久色视频 | 免费在线成人 | 色狠狠操 | 91九色精品女同系列 | 黄色福利视频网站 | 国产黄色成人av | 久久综合久久综合这里只有精品 | 国产美女免费 | 六月丁香婷 | 天天天天天天天操 | 久久国产a | 在线91播放 | 91av手机在线 | 亚洲国产成人久久综合 | 99精品视频免费看 | 中文区中文字幕免费看 | 日本高清免费中文字幕 | 日韩av看片 | 国产精品乱码久久久 | 亚洲 欧美 变态 国产 另类 | www.久久色 | 日韩色av色资源 | 91桃色国产在线播放 | 久久久久影视 | 国产剧情一区 | 黄色a一级片 | 欧美日韩视频在线观看免费 | 婷婷视频在线观看 | 香蕉网址 | 亚洲美女免费视频 | 国产精品96久久久久久吹潮 | 在线天堂中文在线资源网 | 国产精品免费一区二区三区在线观看 | av免费看网站 | 久久尤物电影视频在线观看 | 奇米影视在线99精品 | 中文字幕在线视频第一页 | 综合婷婷丁香 | 国产精品 中文在线 | 久久午夜鲁丝片 | 中文字幕在线播出 | 日韩视频免费 | 色操插 | 国产成人一级 | 久久在线看 | 一区二区三区免费在线观看视频 | 日韩三级免费 | 国产日韩精品视频 | 欧美性色黄| 午夜.dj高清免费观看视频 | 很黄很色很污的网站 | 99久热在线精品视频观看 | 丝袜av一区| 免费av黄色 | 久草精品电影 | 一级黄色片在线 | 在线成人国产 | 九九电影在线 | www.看片网站 | 国产成人三级 | 国产午夜三级一区二区三桃花影视 | 香蕉久草在线 | 最近中文字幕在线播放 | 久久综合狠狠综合久久激情 | 在线不卡的av | 超碰电影在线观看 | 亚洲日韩精品欧美一区二区 | 亚洲精品白浆高清久久久久久 | 婷婷综合亚洲 | 夜夜爽天天爽 | 免费精品国产va自在自线 | 国产福利久久 | 成人免费网站在线观看 | 在线电影 你懂得 | 日韩欧美视频在线 | 国产精品 中文在线 | 成人超碰97 | 欧美国产日韩在线观看 | 五月婷婷六月丁香 | 久久精品在线免费观看 | 五月婷婷国产 | 久久黄页 | www久| 久久精品中文字幕少妇 | 手机av电影在线观看 | 91亚洲精品久久久蜜桃 | 天天色天天射天天综合网 | 狠狠狠狠狠狠 | 久久久九色精品国产一区二区三区 | 国内一级片在线观看 | 精品91| 国产成人精品亚洲精品 | 久久免费国产视频 | 美国三级黄色大片 | 国产中文视频 | 日韩免费av网址 | 国产精品久久久久永久免费观看 | 久久久一本精品99久久精品66 | 国内成人精品2018免费看 | 国产1级视频 | 久久综合狠狠综合久久综合88 | aaa日本高清在线播放免费观看 | 青草草在线 | 在线观看国产亚洲 | 91九色最新地址 | 日本久久久久久久久久 | 日韩美精品视频 | 麻豆手机在线 | 免费三级网 | 国产一区成人在线 | 九九精品在线观看 | 丁香导航 | 国内精品久久久久久中文字幕 | 国产91精品一区二区麻豆亚洲 | 色99视频 | 免费美女久久99 | 91福利国产在线观看 | 蜜臀av性久久久久蜜臀aⅴ涩爱 | 视色网站| 免费在线播放 | 91刺激视频 | 亚洲高清激情 | 免费看国产曰批40分钟 | 久草视频免费观 | 高清免费av在线 | 久久国产精品99久久久久久丝袜 | 伊人婷婷色| 日本黄色免费看 | 91精品在线看 | 国产成人久久精品亚洲 | 亚洲黄色区 | 国产日韩三级 | 婷婷精品进入 | 在线a视频免费观看 | 天天插夜夜操 | 婷婷在线观看视频 | 国产丝袜网站 | 国产成人久久久久 | 国产一在线精品一区在线观看 | 色在线国产 | 激情婷婷欧美 | 正在播放一区二区 | 欧美午夜久久久 | 婷婷av电影 | 国产精品欧美 | 91在线在线观看 | 亚洲精品国产片 | 成人在线一区二区三区 | 精品电影一区 | 欧美三级高清 | 人人澡人人舔 | 97操碰 | 国产精品九九久久久久久久 | 久久6精品| 国产剧情av在线播放 | 日韩欧美精品在线 | 亚洲人成在线观看 | 麻豆视频在线免费 | 欧美日韩一区二区免费在线观看 | www.99久久.com| 成人久久网 | 精品国产一区二区三区久久久蜜臀 | 中文字幕视频一区二区 | 欧美在线观看视频一区二区三区 | 人人干网站| 色资源中文字幕 | 免费国产一区二区 | 国产精品自产拍 | 精品高清美女精品国产区 | 五月天网页| 一二三精品视频 | 久99久精品视频免费观看 | 三级av免费| 少妇搡bbb | 天天射射天天 | 欧美最猛性xxxxx亚洲精品 | 欧美一级性 | 丁香高清视频在线看看 | 国产精品久久久久影院 | 久久天天躁狠狠躁夜夜不卡公司 | 亚洲老妇xxxxxx | 黄色精品免费 | 亚洲精品国偷自产在线91正片 | 99久久久国产精品免费99 | 岛国片在线 | 日韩在线视频在线观看 | 亚洲国产成人av网 | 99在线免费视频 | 亚洲成av人影院 | 伊人婷婷久久 | 婷婷六月丁 | 国产破处精品 | 中文字幕日韩有码 | 国产精品免费一区二区 | 播五月综合 | 婷婷日韩 | 国内外成人在线 | 亚洲国产成人在线观看 | 成人蜜桃视频 | 亚洲精品999 | 国产黄色高清 | 免费日韩 精品中文字幕视频在线 | 丁香六月久久综合狠狠色 | 久久www免费人成看片高清 | 69绿帽绿奴3pvideos | 天天操偷偷干 | av噜噜噜在线播放 | 91麻豆网| 久草国产视频 | 99视频+国产日韩欧美 | 免费视频久久久久久久 | 国产在线一区观看 | 在线看片日韩 | 涩涩在线| 日韩视频区 | 免费成人在线网站 | 国产精品久久久久久久电影 | av亚洲产国偷v产偷v自拍小说 | 日韩中文字幕视频在线 | 一区二区久久 | 精品亚洲免费视频 | 国产香蕉久久精品综合网 | 免费日韩 精品中文字幕视频在线 | 日韩免费在线网站 | 国产不卡精品视频 | 国产一区欧美日韩 | 五月婷丁香网 | a午夜在线| 天天视频色版 | 草久视频在线观看 | 中日韩欧美精彩视频 | 夜夜摸夜夜爽 | 国产美女精品人人做人人爽 | 精品99在线视频 | 2023国产精品自产拍在线观看 | 久久久久久久久久久网站 | 国精产品999国精产品视频 | 天天av在线播放 | 婷婷激情小说网 | 免费成人在线观看 | 国产精品久久久久永久免费 | 国产手机精品视频 | 天堂资源在线观看视频 | aaa亚洲精品一二三区 | 一级黄色片在线观看 | 99亚洲国产精品 | 久久精品亚洲一区二区三区观看模式 | 成人国产网站 | 久久久九色精品国产一区二区三区 | 色综合天天做天天爱 | 免费视频成人 | 亚洲一区二区麻豆 | 国产天天爽 | 色综合天天色综合 | 日韩成片 | 最近中文字幕在线中文高清版 | 免费日韩高清 | 天天伊人狠狠 | 精品在线观看一区二区 | 色亚洲网| 亚洲性xxxx | 在线成人欧美 | 麻豆你懂的 | 婷婷国产一区二区三区 | 97色资源 | 午夜.dj高清免费观看视频 | 国产片免费在线观看视频 | 1024手机看片国产 | 中文字幕在线观看免费 | 国产精品麻豆果冻传媒在线播放 | 天天综合网 天天综合色 | 91c网站色版视频 | 91桃色在线观看视频 | 日产乱码一二三区别免费 | 久久一区二 | 国产精品成人免费 | 一二区av | 国产 日韩 欧美 自拍 | 亚洲精品成人在线 | 国内精品久久久久久久久久 | 99操视频| 久久神马影院 | 色窝资源 | 九九日九九操 | 超碰国产在线观看 | 国产成人av福利 | 久久成人一区二区 | 国产香蕉视频在线观看 | 久久亚洲人 | 九九久久精品 | 玖玖综合网| 国产二区免费视频 | 天堂av在线免费观看 | 国产精品a久久 | 五月激情电影 | 五月婷婷在线播放 | 性色av一区二区三区在线观看 | 香蕉视频一级 | 国产视频精品免费 | 久久成人一区二区 | 亚洲国产视频a | 免费观看www7722午夜电影 | 日韩区欠美精品av视频 | 国产日本亚洲高清 | 四虎国产精品永久在线国在线 | 激情综合亚洲精品 | 国内视频在线观看 | 欧美色图88 | 免费在线播放黄色 | 亚洲精选视频免费看 | 美女国产精品 | 亚洲伊人成综合网 | 99精品影视 | 三级黄色欧美 | 欧美日本日韩aⅴ在线视频 插插插色综合 | 精品影院一区二区久久久 | 欧美日韩免费网站 | 色综合色综合久久综合频道88 | 在线视频区 | 精品久久久久久国产91 | 久在线| 青青视频一区 | 欧美日韩一区二区三区视频 | 中文久草| 亚洲国产操 | 国产美女视频黄a视频免费 久久综合九色欧美综合狠狠 | 亚洲精品免费播放 | 中文字幕av最新 | 黄色小说免费在线观看 | 久久久久久国产精品久久 | 免费看成人a | 99视频精品全部免费 在线 | 国产精品久久久久久999 | 免费在线观看一区二区三区 | 韩国av一区二区三区在线观看 | 国产一区二区三区久久久 | 国产精品综合av一区二区国产馆 | 国产精品视频app | 国产一区成人在线 | 欧美一级片 | 日韩精品中文字幕久久臀 | 一级黄色免费网站 | 婷婷丁香花五月天 | 天天操天天射天天操 | 欧美一区二区三区在线视频观看 | 欧美激情精品久久久久 | 手机av网站 | 伊人久久一区 | 美女国产精品 | 激情av在线资源 | 色免费在线 | 我爱av激情网 | 黄色资源网站 | 色综合中文综合网 | 91av视频免费在线观看 | 久久精彩免费视频 | 亚洲精品国产精品乱码不99热 | 成人av影视| 欧美成人精品欧美一级乱 | 日本性久久 | 日韩欧美国产精品 | 欧美aaaxxxx做受视频 | 奇米影视在线99精品 | 99久久精品视频免费 | 在线视频一区二区 | 91系列在线 | 在线免费色视频 | 免费h精品视频在线播放 | 国产精品久久久久久久久久久久午夜 | 久久最新视频 | 色视频网页 | 成人av在线影院 | 天天插视频 | 免费黄色a网站 | 91亚洲精品在线 | 99爱视频在线观看 | 精品久久久99 | 一区精品在线 | 国产1级毛片 | 日韩欧美区 | 91自拍成人 | 九色琪琪久久综合网天天 | 韩国在线一区二区 | 人人舔人人插 | 亚洲婷婷在线视频 | 91最新在线视频 | 日本精品视频在线观看 | 激情久久久久久久久久久久久久久久 | 欧美色图p | av超碰在线观看 | 又色又爽的网站 | 欧美精品在线观看免费 | 综合国产在线 | 91丨九色丨91啦蝌蚪老版 | 永久免费毛片在线观看 | 91在线视频免费播放 | 欧美精品亚州精品 | 久久国产精品久久国产精品 | 久久图| 91av原创| 免费国产黄线在线观看视频 | 国产精品一区二区在线观看 | av大全在线观看 | 久久66热这里只有精品 | 国产只有精品 | 激情av一区二区 | 欧美一区二区三区在线视频观看 | 日韩免费电影网 | 亚洲干视频在线观看 | 成人三级网址 | 国产精品一区二区久久精品爱微奶 | 天天久久综合 | 一区二区三区免费播放 | 久久久久久久久久久久电影 | 97热视频 | a√资源在线 | 天天色视频 | 久久久久久久久免费视频 | 99视频免费看 | 久久无码av一区二区三区电影网 | 亚洲最大色 | 91成人精品一区在线播放69 | 亚洲伊人色 | 大胆欧美gogo免费视频一二区 | 中文一区在线观看 | 美女免费视频观看网站 | 欧洲亚洲激情 | 在线a亚洲视频播放在线观看 | 天天舔天天射天天操 | 又污又黄的网站 | 在线免费av网 | 成人xxxx| 91免费视频黄 | av电影免费在线看 | 永久中文字幕 | 国产精品不卡av | 丁香婷婷综合色啪 | 国产精品99视频 | 国产色区 | 亚洲无线视频 | 高清av网站 | 综合久久网站 | 婷婷在线免费观看 | 91视频在线免费看 | 免费看一级片 | 在线国产高清 | 亚洲一区在线看 | 在线欧美小视频 | 夜夜操天天操 | 激情丁香 | 日韩一区二区三区高清免费看看 | 精品一区电影国产 | 美女精品 | 粉嫩高清一区二区三区 | 一区二区欧美激情 | 奇米影视四色8888 | 日韩欧美在线综合网 | 亚洲精品视频在线观看免费 | 亚洲天堂色婷婷 | 深爱激情综合网 | 在线观看日韩精品 | 亚洲免费永久精品国产 | 免费视频一级片 | 欧美一二三四在线 | 成人免费看视频 | 这里只有精品视频在线 | 911久久香蕉国产线看观看 | 亚洲黄色在线观看 | 天堂网av 在线 | 亚洲精品免费观看视频 | 中文电影网 | 波多野结衣最新 | 97超在线 | 国产97视频在线 | 婷婷精品国产欧美精品亚洲人人爽 | 国产精品18久久久久久久网站 | 久久久久久久久久免费视频 | 免费在线观看av不卡 | 精品在线免费观看 | 精品免费观看 | 亚洲性视频 | 粉嫩av一区二区三区四区在线观看 | 日本公妇色中文字幕 | 成人久久视频 | 欧美作爱视频 | 91在线视频观看免费 | 婷婷视频在线观看 | 日韩欧美一区二区在线观看 | 国产自制av | 亚洲国产99| 久久综合久色欧美综合狠狠 | 免费观看久久 | 97超碰免费 | 操久久免费视频 | 一区二区精 | 97超碰国产精品 | 91热爆视频| 久久久国产影视 | 中文字幕4 | 日韩视频一区二区在线 | 欧美另类v| 狠狠地日 | 国产精品每日更新 | 久久99电影| 精品久久久久久亚洲综合网 | 欧美日韩国产区 | 色伊人网 | 在线视频手机国产 | 96av麻豆蜜桃一区二区 | 色姑娘综合网 | 不卡的av中文字幕 | 国产精品18久久久久久首页狼 | 国产精品免费观看视频 | 99精品国产一区二区三区麻豆 | 一区二区三区在线视频111 | 久久午夜鲁丝片 | 毛片基地黄久久久久久天堂 | 香蕉久久久久久av成人 | 在线观看免费一区 | 成人亚洲精品国产www | 在线看中文字幕 | 欧美综合干 | 国产在线一区二区 | 中文字幕频道 | 精品在线看 | 欧美日韩xx | 一级片黄色片网站 | 三级小视频在线观看 | 欧美另类性 | 国产看片网站 | 婷婷丁香九月 | 女人18毛片90分钟 | 久久er99热精品一区二区 | 亚洲欧美日韩一二三区 | 91精品国产99久久久久 | 日本久久久久 | 国产色在线视频 | 最近中文字幕高清字幕在线视频 | 天天拍天天操 | 色噜噜噜噜 | 婷婷午夜天 | 99视频在线免费观看 | 欧美精品久久久久久久久老牛影院 | 最新中文字幕在线观看视频 | 黄色三几片 | 日韩精品在线免费播放 | 最近更新好看的中文字幕 | 欧美肥妇free| 日韩精品中文字幕久久臀 | 精品免费国产一区二区三区四区 | 亚洲区另类春色综合小说校园片 | 色橹橹欧美在线观看视频高清 | 日韩在线免费高清视频 | 日韩高清免费在线观看 | a级国产片 | 免费看国产精品 | 免费在线观看不卡av | 国内精品久久久久影院男同志 | 奇米影视在线99精品 | 亚洲mv大片欧洲mv大片免费 | 中文字幕电影高清在线观看 | 日日成人网 | 6080yy精品一区二区三区 | 欧美久久久久久久久久 | 久久艹在线观看 | 亚洲三级毛片 | 人人澡人人爽 | 99在线播放 | 视频在线播放国产 | 天天射天天干天天操 | 欧美精品中文字幕亚洲专区 | 日韩精品一区二区在线观看 | 日韩成人看片 | 久久精品激情 | 欧美日韩国产综合一区二区 | 午夜在线资源 | 三级在线视频播放 | 日韩免费在线看 | 91精品国产综合久久福利 | 亚洲永久精品一区 | 国产精品美女久久久久久2018 | 综合久久久久久久 | 国产专区视频在线观看 | 中文字幕在线视频国产 | 亚洲国产成人精品在线观看 | av久久在线 | 久久艹免费 | 在线成人免费电影 | 国产一区视频在线观看免费 | 99视频在线精品 | 激情综合一区 | 亚洲日日射 | 最近最新中文字幕 | 三级在线国产 | 一区二区视频在线看 | 日本特黄一级片 | 丁香影院在线 | 午夜精品一区二区三区四区 | 不卡av免费在线观看 | 2019中文最近的2019中文在线 | 99久久www| 免费又黄又爽视频 | 在线观看av网 | 精品久久久久久亚洲综合网 | 天天射天天爱天天干 | 日韩精品一区二区在线视频 | 超级碰碰免费视频 | 96视频在线| 国产欧美三级 | 丁香花中文在线免费观看 | 日韩三级av| 久久狠狠一本精品综合网 | 日韩伦理一区二区三区av在线 | 天天色天 | 精品久久九九 | 午夜精品久久久久久久99热影院 | 久草免费在线观看 | 日日干av | 麻豆av一区二区三区在线观看 | 成人在线免费视频 | 国产这里只有精品 | 久久夜夜操 | 在线日本v二区不卡 | 国产日产av | av在线不卡观看 | 在线综合色 | 欧美日韩综合在线 | 亚洲精品一区二区三区在线观看 | 久久久久99精品国产片 | 国产精品h在线观看 | 国产成人在线免费观看 | 正在播放国产精品 | 人人精品久久 | 精品国产一区二区三区久久久久久 | 欧美老女人xx | 草久草久 | 国产高清视频在线观看 | 在线观看韩日电影免费 | 天天综合网 天天综合色 | 欧美在线观看禁18 | 国产黄色免费在线观看 | 777xxx欧美 | 成人黄色在线观看视频 | 国产网站色 | 欧美另类交在线观看 | av大全在线观看 | 亚洲 av网站 | 久久 亚洲视频 | 国产资源| 婷婷在线免费观看 | 久久天堂亚洲 | adn—256中文在线观看 | 天天操天天干天天 | 91在线最新| 日韩一二三区不卡 | 黄色字幕网 | 伊人久久精品久久亚洲一区 | 国产韩国精品一区二区三区 | 蜜桃av久久久亚洲精品 | 久草综合视频 | 久操97 | 日韩xxxx视频 | 国产中文字幕在线观看 | 久久精品国产免费 | 欧美色噜噜噜 | 怡红院av| 国产成人精品一二三区 | 亚洲成色 | 99福利片| 日本精品二区 | 成人黄大片 | caobi视频| 东方av在| 波多野结衣视频在线 | 在线视频 91 | 精品一区二区av | 国产日韩欧美在线免费观看 | 久久影院亚洲 | 黄色国产区| 激情欧美一区二区免费视频 | 最近日本mv字幕免费观看 | 久久综合给合久久狠狠色 | 国产视频黄 | 亚洲国产精品成人女人久久 | 日韩亚洲精品电影 | 视频一区二区在线 | 色av婷婷 | 久草在线一免费新视频 | 久久久久久国产精品免费 | 中文字幕精品一区二区精品 | 热久久免费视频 | 国产精品美女久久久久久久久 | 草久视频在线观看 | 国产成人精品亚洲精品 | av一级片在线观看 | 国内精品一区二区 | 日韩成人欧美 | 婷婷色视频 | 久久久免费毛片 | 日本精品视频免费 | 国产综合婷婷 | 国产一区二区在线观看视频 | 色视频网站在线 | 精品亚洲午夜久久久久91 | 少妇bbbb搡bbbb搡bbbb | 久久伊人国产精品 | 中文国产在线观看 | 国产高清视频 | 最新av在线网站 | 日本精品视频网站 | 国产一级视频 | 97色狠狠 | 日韩影视精品 | 最新日韩在线观看视频 | 天天插日日操 | 国产精品福利在线 | 最近日韩中文字幕中文 | 成人一级在线观看 | 成人免费xxx在线观看 | 亚洲精品乱码久久久久久蜜桃欧美 | 中文字幕久久网 | 国产成人一二片 |