青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

posts - 100,  comments - 15,  trackbacks - 0
一、引言

計算機的出現使得很多原本十分繁瑣的工作得以大幅度簡化,但是也有一些在人們直觀看來很容易的問題卻需要拿出一套并不簡單的通用解決方案,比如幾何問題。作為計算機科學的一個分支,計算幾何主要研究解決幾何問題的算法。在現代工程和數學領域,計算幾何在圖形學、機器人技術、超大規模集成電路設計和統計等諸多領域有著十分重要的應用。在本文中,我們將對計算幾何常用的基本算法做一個全面的介紹,希望對您了解并應用計算幾何的知識解決問題起到幫助。

二、目錄

  本文整理的計算幾何基本概念和常用算法包括如下內容:

  矢量的概念

  矢量加減法

  矢量叉積

  折線段的拐向判斷

  判斷點是否在線段上

  判斷兩線段是否相交

  判斷線段和直線是否相交

  判斷矩形是否包含點

  判斷線段、折線、多邊形是否在矩形中

  判斷矩形是否在矩形中

  判斷圓是否在矩形中

  判斷點是否在多邊形中

  判斷線段是否在多邊形內

  判斷折線是否在多邊形內

  判斷多邊形是否在多邊形內

  判斷矩形是否在多邊形內

  判斷圓是否在多邊形內

  判斷點是否在圓內

  判斷線段、折線、矩形、多邊形是否在圓內

  判斷圓是否在圓內

  計算點到線段的最近點

  計算點到折線、矩形、多邊形的最近點

  計算點到圓的最近距離及交點坐標

  計算兩條共線的線段的交點

  計算線段或直線與線段的交點

  求線段或直線與折線、矩形、多邊形的交點

  求線段或直線與圓的交點

  凸包的概念

  凸包的求法

三、算法介紹

矢量的概念:

如果一條線段的端點是有次序之分的,我們把這種線段成為有向線段(directed segment)。如果有向線段p1p2的起點p1在坐標原點,我們可以把它稱為矢量(vector)p2。

矢量叉積:

計算矢量叉積是與直線和線段相關算法的核心部分。設矢量P = (x1,y1) ,Q = (x2,y2),則矢量叉積定義為由(0,0)、p1、p2和p1+p2所組成的平行四邊形的帶符號的面積,即:P × Q = x1*y2 - x2*y1,其結果是一個標量。顯然有性質 P × Q = - ( Q × P ) 和 P × ( - Q ) = - ( P × Q )。一般在不加說明的情況下,本文下述算法中所有的點都看作矢量,兩點的加減法就是矢量相加減,而點的乘法則看作矢量叉積。

  叉積的一個非常重要性質是可以通過它的符號判斷兩矢量相互之間的順逆時針關系:

  若 P × Q > 0 , 則P在Q的順時針方向。
  若 P × Q < 0 , 則P在Q的逆時針方向。
  若 P × Q = 0 , 則P與Q共線,但可能同向也可能反向。

折線段的拐向判斷:

  折線段的拐向判斷方法可以直接由矢量叉積的性質推出。對于有公共端點的線段p0p1和p1p2,通過計算(p2 - p0) × (p1 - p0)的符號便可以確定折線段的拐向:

  若(p2 - p0) × (p1 - p0) > 0,則p0p1在p1點拐向右側后得到p1p2。

  若(p2 - p0) × (p1 - p0) < 0,則p0p1在p1點拐向左側后得到p1p2。

  若(p2 - p0) × (p1 - p0) = 0,則p0、p1、p2三點共線。

  具體情況可參照下圖:

判斷點是否在線段上:

  設點為Q,線段為P1P2 ,判斷點Q在該線段上的依據是:( Q - P1 ) × ( P2 - P1 ) = 0 且 Q 在以 P1,P2為對角頂點的矩形內。前者保證Q點在直線P1P2上,后者是保證Q點不在線段P1P2的延長線或反向延長線上,對于這一步驟的判斷可以用以下過程實現:

  ON-SEGMENT(pi,pj,pk)

  if min(xi,xj)<=xk<=max(xi,xj) and min(yi,yj)<=yk<=max(yi,yj)

  then return true;

  else return false;

  特別要注意的是,由于需要考慮水平線段和垂直線段兩種特殊情況,min(xi,xj)<=xk<=max(xi,xj)和min(yi,yj)<=yk<=max(yi,yj)兩個條件必須同時滿足才能返回真值。

判斷兩線段是否相交:

我們分兩步確定兩條線段是否相交:

  (1)快速排斥試驗

  設以線段 P1P2 為對角線的矩形為R, 設以線段 Q1Q2 為對角線的矩形為T,如果R和T不相交,顯然兩線段不會相交。

  (2)跨立試驗

  如果兩線段相交,則兩線段必然相互跨立對方。若P1P2跨立Q1Q2 ,則矢量 ( P1 - Q1 ) 和( P2 - Q1 )位于矢量( Q2 - Q1 ) 的兩側,即( P1 - Q1 ) × ( Q2 - Q1 ) * ( P2 - Q1 ) × ( Q2 - Q1 ) < 0。上式可改寫成( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) > 0。當 ( P1 - Q1 ) × ( Q2 - Q1 ) = 0 時,說明 ( P1 - Q1 ) 和 ( Q2 - Q1 )共線,但是因為已經通過快速排斥試驗,所以 P1 一定在線段 Q1Q2上;同理,( Q2 - Q1 ) ×(P2 - Q1 ) = 0 說明 P2 一定在線段 Q1Q2上。所以判斷P1P2跨立Q1Q2的依據是:( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) >= 0。同理判斷Q1Q2跨立P1P2的依據是:( Q1 - P1 ) × ( P2 - P1 ) * ( P2 - P1 ) × ( Q2 - P1 ) >= 0。具體情況如下圖所示:

在相同的原理下,對此算法的具體的實現細節可能會與此有所不同,除了這種過程外,大家也可以參考《算法導論》上的實現。

判斷線段和直線是否相交:

有了上面的基礎,這個算法就很容易了。如果線段P1P2和直線Q1Q2相交,則P1P2跨立Q1Q2,即:( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) >= 0。

判斷矩形是否包含點:

只要判斷該點的橫坐標和縱坐標是否夾在矩形的左右邊和上下邊之間。

判斷線段、折線、多邊形是否在矩形中:

因為矩形是個凸集,所以只要判斷所有端點是否都在矩形中就可以了。

判斷矩形是否在矩形中:

只要比較左右邊界和上下邊界就可以了。

判斷圓是否在矩形中:

很容易證明,圓在矩形中的充要條件是:圓心在矩形中且圓的半徑小于等于圓心到矩形四邊的距離的最小值。

判斷點是否在多邊形中:

  判斷點P是否在多邊形中是計算幾何中一個非常基本但是十分重要的算法。以點P為端點,向左方作射線L,由于多邊形是有界的,所以射線L的左端一定在多邊形外,考慮沿著L從無窮遠處開始自左向右移動,遇到和多邊形的第一個交點的時候,進入到了多邊形的內部,遇到第二個交點的時候,離開了多邊形,……所以很容易看出當L和多邊形的交點數目C是奇數的時候,P在多邊形內,是偶數的話P在多邊形外。

  但是有些特殊情況要加以考慮。如圖下圖(a)(b)(c)(d)所示。在圖(a)中,L和多邊形的頂點相交,這時候交點只能計算一個;在圖(b)中,L和多邊形頂點的交點不應被計算;在圖(c)和(d) 中,L和多邊形的一條邊重合,這條邊應該被忽略不計。如果L和多邊形的一條邊重合,這條邊應該被忽略不計。

為了統一起見,我們在計算射線L和多邊形的交點的時候,1。對于多邊形的水平邊不作考慮;2。對于多邊形的頂點和L相交的情況,如果該頂點是其所屬的邊上縱坐標較大的頂點,則計數,否則忽略;3。對于P在多邊形邊上的情形,直接可判斷P屬于多邊行。由此得出算法的偽代碼如下:

  count ← 0;
  以P為端點,作從右向左的射線L;
  for 多邊形的每條邊s
  do if P在邊s上
  then return true;
  if s不是水平的
  then if s的一個端點在L上
  if 該端點是s兩端點中縱坐標較大的端點
  then count ← count+1
  else if s和L相交
  then count ← count+1;
  if count mod 2 = 1
  then return true;
  else return false;

  其中做射線L的方法是:設P'的縱坐標和P相同,橫坐標為正無窮大(很大的一個正數),則P和P'就確定了射線L。

  判斷點是否在多邊形中的這個算法的時間復雜度為O(n)。

  另外還有一種算法是用帶符號的三角形面積之和與多邊形面積進行比較,這種算法由于使用浮點數運算所以會帶來一定誤差,不推薦大家使用。

判斷線段是否在多邊形內:

  線段在多邊形內的一個必要條件是線段的兩個端點都在多邊形內,但由于多邊形可能為凹,所以這不能成為判斷的充分條件。如果線段和多邊形的某條邊內交(兩線段內交是指兩線段相交且交點不在兩線段的端點),因為多邊形的邊的左右兩側分屬多邊形內外不同部分,所以線段一定會有一部分在多邊形外(見圖a)。于是我們得到線段在多邊形內的第二個必要條件:線段和多邊形的所有邊都不內交。

  線段和多邊形交于線段的兩端點并不會影響線段是否在多邊形內;但是如果多邊形的某個頂點和線段相交,還必須判斷兩相鄰交點之間的線段是否包含于多邊形內部(反例見圖b)。

因此我們可以先求出所有和線段相交的多邊形的頂點,然后按照X-Y坐標排序(X坐標小的排在前面,對于X坐標相同的點,Y坐標小的排在前面,這種排序準則也是為了保證水平和垂直情況的判斷正確),這樣相鄰的兩個點就是在線段上相鄰的兩交點,如果任意相鄰兩點的中點也在多邊形內,則該線段一定在多邊形內。

  證明如下:

  命題1:

  如果線段和多邊形的兩相鄰交點P1 ,P2的中點P' 也在多邊形內,則P1, P2之間的所有點都在多邊形內。

  證明:

  假設P1,P2之間含有不在多邊形內的點,不妨設該點為Q,在P1, P'之間,因為多邊形是閉合曲線,所以其內外部之間有界,而P1屬于多邊行內部,Q屬于多邊性外部,P'屬于多邊性內部,P1-Q-P'完全連續,所以P1Q和QP'一定跨越多邊形的邊界,因此在P1,P'之間至少還有兩個該線段和多邊形的交點,這和P1P2是相鄰兩交點矛盾,故命題成立。證畢。

  由命題1直接可得出推論:

  推論2:

  設多邊形和線段PQ的交點依次為P1,P2,……Pn,其中Pi和Pi+1是相鄰兩交點,線段PQ在多邊形內的充要條件是:P,Q在多邊形內且對于i =1, 2,……, n-1,Pi ,Pi+1的中點也在多邊形內。

  在實際編程中,沒有必要計算所有的交點,首先應判斷線段和多邊形的邊是否內交,倘若線段和多邊形的某條邊內交則線段一定在多邊形外;如果線段和多邊形的每一條邊都不內交,則線段和多邊形的交點一定是線段的端點或者多邊形的頂點,只要判斷點是否在線段上就可以了。

  至此我們得出算法如下:

  if 線端PQ的端點不都在多邊形內
  then return false;
  點集pointSet初始化為空;
  for 多邊形的每條邊s
  do if 線段的某個端點在s上
  then 將該端點加入pointSet;
  else if s的某個端點在線段PQ上
  then 將該端點加入pointSet;
  else if s和線段PQ相交 // 這時候已經可以肯定是內交了
  then return false;
  將pointSet中的點按照X-Y坐標排序;
  for pointSet中每兩個相鄰點 pointSet[i] , pointSet[ i+1]
  do if pointSet[i] , pointSet[ i+1] 的中點不在多邊形中
  then return false;
  return true;

  這個過程中的排序因為交點數目肯定遠小于多邊形的頂點數目n,所以最多是常數級的復雜度,幾乎可以忽略不計。因此算法的時間復雜度也是O(n)。

判斷折線是否在多邊形內:

  只要判斷折線的每條線段是否都在多邊形內即可。設折線有m條線段,多邊形有n個頂點,則該算法的時間復雜度為O(m*n)。

 判斷多邊形是否在多邊形內:

  只要判斷多邊形的每條邊是否都在多邊形內即可。判斷一個有m個頂點的多邊形是否在一個有n個頂點的多邊形內復雜度為O(m*n)。

  判斷矩形是否在多邊形內:

將矩形轉化為多邊形,然后再判斷是否在多邊形內。

判斷圓是否在多邊形內:

  只要計算圓心到多邊形的每條邊的最短距離,如果該距離大于等于圓半徑則該圓在多邊形內。計算圓心到多邊形每條邊最短距離的算法在后文闡述。

判斷點是否在圓內:

計算圓心到該點的距離,如果小于等于半徑則該點在圓內。

判斷線段、折線、矩形、多邊形是否在圓內:

因為圓是凸集,所以只要判斷是否每個頂點都在圓內即可。

判斷圓是否在圓內:

 設兩圓為O1,O2,半徑分別為r1, r2,要判斷O2是否在O1內。先比較r1,r2的大小,如果r1<r2則O2不可能在O1內;否則如果兩圓心的距離大于r1 - r2 ,則O2不在O1內;否則O2在O1內。

計算點到線段的最近點:

  如果該線段平行于X軸(Y軸),則過點point作該線段所在直線的垂線,垂足很容易求得,然后計算出垂足,如果垂足在線段上則返回垂足,否則返回離垂足近的端點;如果該線段不平行于X軸也不平行于Y軸,則斜率存在且不為0。設線段的兩端點為pt1和pt2,斜率為:k = ( pt2.y - pt1. y ) / (pt2.x - pt1.x );該直線方程為:y = k* ( x - pt1.x) + pt1.y。其垂線的斜率為 - 1 / k,垂線方程為:y = (-1/k) * (x - point.x) + point.y 。

  聯立兩直線方程解得:x = ( k^2 * pt1.x + k * (point.y - pt1.y ) + point.x ) / ( k^2 + 1) ,y = k * ( x - pt1.x) + pt1.y;然后再判斷垂足是否在線段上,如果在線段上則返回垂足;如果不在則計算兩端點到垂足的距離,選擇距離垂足較近的端點返回。

計算點到折線、矩形、多邊形的最近點:

只要分別計算點到每條線段的最近點,記錄最近距離,取其中最近距離最小的點即可。

計算點到圓的最近距離及交點坐標:

  如果該點在圓心,因為圓心到圓周任一點的距離相等,返回UNDEFINED。

  連接點P和圓心O,如果PO平行于X軸,則根據P在O的左邊還是右邊計算出最近點的橫坐標為centerPoint.x - radius 或 centerPoint.x + radius。如果PO平行于Y軸,則根據P在O的上邊還是下邊計算出最近點的縱坐標為 centerPoint.y -+radius或 centerPoint.y - radius。如果PO不平行于X軸和Y軸,則PO的斜率存在且不為0,這時直線PO斜率為k = ( P.y - O.y )/ ( P.x - O.x )。直線PO的方程為:y = k * ( x - P.x) + P.y。設圓方程為:(x - O.x ) ^2 + ( y - O.y ) ^2 = r ^2,聯立兩方程組可以解出直線PO和圓的交點,取其中離P點較近的交點即可。

計算兩條共線的線段的交點:

  對于兩條共線的線段,它們之間的位置關系有下圖所示的幾種情況。圖(a)中兩條線段沒有交點;圖 (b) 和 (d) 中兩條線段有無窮焦點;圖 (c) 中兩條線段有一個交點。設line1是兩條線段中較長的一條,line2是較短的一條,如果line1包含了line2的兩個端點,則是圖(d)的情況,兩線段有無窮交點;如果line1只包含line2的一個端點,那么如果line1的某個端點等于被line1包含的line2的那個端點,則是圖(c)的情況,這時兩線段只有一個交點,否則就是圖(b)的情況,兩線段也是有無窮的交點;如果line1不包含line2的任何端點,則是圖(a)的情況,這時兩線段沒有交點。

 

計算線段或直線與線段的交點:

  設一條線段為L0 = P1P2,另一條線段或直線為L1 = Q1Q2 ,要計算的就是L0和L1的交點。

  1. 首先判斷L0和L1是否相交(方法已在前文討論過),如果不相交則沒有交點,否則說明L0和L1一定有交點,下面就將L0和L1都看作直線來考慮。

  2. 如果P1和P2橫坐標相同,即L0平行于Y軸

  a) 若L1也平行于Y軸,

  i. 若P1的縱坐標和Q1的縱坐標相同,說明L0和L1共線,假如L1是直線的話他們有無窮的交點,假如L1是線段的話可用"計算兩條共線線段的交點"的算法求他們的交點(該方法在前文已討論過);

  ii. 否則說明L0和L1平行,他們沒有交點;

  b) 若L1不平行于Y軸,則交點橫坐標為P1的橫坐標,代入到L1的直線方程中可以計算出交點縱坐標;

  3. 如果P1和P2橫坐標不同,但是Q1和Q2橫坐標相同,即L1平行于Y軸,則交點橫坐標為Q1的橫坐標,代入到L0的直線方程中可以計算出交點縱坐標;

  4. 如果P1和P2縱坐標相同,即L0平行于X軸

  a) 若L1也平行于X軸,

  i. 若P1的橫坐標和Q1的橫坐標相同,說明L0和L1共線,假如L1是直線的話他們有無窮的交點,假如L1是線段的話可用"計算兩條共線線段的交點"的算法求他們的交點(該方法在前文已討論過);

  ii. 否則說明L0和L1平行,他們沒有交點;

  b) 若L1不平行于X軸,則交點縱坐標為P1的縱坐標,代入到L1的直線方程中可以計算出交點橫坐標;

  5. 如果P1和P2縱坐標不同,但是Q1和Q2縱坐標相同,即L1平行于X軸,則交點縱坐標為Q1的縱坐標,代入到L0的直線方程中可以計算出交點橫坐標;

  6. 剩下的情況就是L1和L0的斜率均存在且不為0的情況

  a) 計算出L0的斜率K0,L1的斜率K1 ;

  b) 如果K1 = K2

  i. 如果Q1在L0上,則說明L0和L1共線,假如L1是直線的話有無窮交點,假如L1是線段的話可用"計算兩條共線線段的交點"的算法求他們的交點(該方法在前文已討論過);

  ii. 如果Q1不在L0上,則說明L0和L1平行,他們沒有交點。

  c) 聯立兩直線的方程組可以解出交點來

  這個算法并不復雜,但是要分情況討論清楚,尤其是當兩條線段共線的情況需要單獨考慮,所以在前文將求兩條共線線段的算法單獨寫出來。另外,一開始就先利用矢量叉乘判斷線段與線段(或直線)是否相交,如果結果是相交,那么在后面就可以將線段全部看作直線來考慮。需要注意的是,我們可以將直線或線段方程改寫為ax+by+c=0的形式,這樣一來上述過程的部分步驟可以合并,縮短了代碼長度,但是由于先要求出參數,這種算法將花費更多的時間。

求線段或直線與折線、矩形、多邊形的交點:

分別求與每條邊的交點即可。

求線段或直線與圓的交點:

  設圓心為O,圓半徑為r,直線(或線段)L上的兩點為P1,P2。

  1. 如果L是線段且P1,P2都包含在圓O內,則沒有交點;否則進行下一步。

  2. 如果L平行于Y軸,

  a) 計算圓心到L的距離dis;
  b) 如果dis > r 則L和圓沒有交點;
  c) 利用勾股定理,可以求出兩交點坐標,但要注意考慮L和圓的相切情況。

  3. 如果L平行于X軸,做法與L平行于Y軸的情況類似;

  4. 如果L既不平行X軸也不平行Y軸,可以求出L的斜率K,然后列出L的點斜式方程,和圓方程聯立即可求解出L和圓的兩個交點;

  5. 如果L是線段,對于2,3,4中求出的交點還要分別判斷是否屬于該線段的范圍內。

凸包的概念:

  點集Q的凸包(convex hull)是指一個最小凸多邊形,滿足Q中的點或者在多邊形邊上或者在其內。下圖中由紅色線段表示的多邊形就是點集Q={p0,p1,...p12}的凸包。

 

凸包的求法:

  現在已經證明了凸包算法的時間復雜度下界是O(n*logn),但是當凸包的頂點數h也被考慮進去的話,Krikpatrick和Seidel的剪枝搜索算法可以達到O(n*logh),在漸進意義下達到最優。最常用的凸包算法是Graham掃描法和Jarvis步進法。本文只簡單介紹一下Graham掃描法,其正確性的證明和Jarvis步進法的過程大家可以參考《算法導論》。

  對于一個有三個或以上點的點集Q,Graham掃描法的過程如下:

  令p0為Q中Y-X坐標排序下最小的點

  設<p1,p2,...pm>為對其余點按以p0為中心的極角逆時針排序所得的點集(如果有多個點有相同的極角,除了距p0最遠的點外全部移除

  壓p0進棧S
  壓p1進棧S
  壓p2進棧S
  for i ← 3 to m
  do while 由S的棧頂元素的下一個元素、S的棧頂元素以及pi構成的折線段不拐向左側
  對S彈棧
  壓pi進棧S
  return S;

  此過程執行后,棧S由底至頂的元素就是Q的凸包頂點按逆時針排列的點序列。需要注意的是,我們對點按極角逆時針排序時,并不需要真正求出極角,只需要求出任意兩點的次序就可以了。而這個步驟可以用前述的矢量叉積性質實現。

四、結語

  盡管人類對幾何學的研究從古代起便沒有中斷過,但是具體到借助計算機來解決幾何問題的研究,還只是停留在一個初級階段,無論從應用領域還是發展前景來看,計算幾何學都值得我們認真學習、加以運用,希望這篇文章能帶你走進這個豐富多彩的世界。

posted on 2009-10-06 14:38 wyiu 閱讀(172) 評論(0)  編輯 收藏 引用 所屬分類: 算法
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            国产视频在线观看一区| 亚洲一区二区精品视频| 日韩午夜在线播放| 亚洲日本在线视频观看| 亚洲国产精品一区二区第一页| 国产视频一区欧美| 国产一区二区三区久久| 激情久久影院| 日韩视频在线播放| 亚洲一区二区在线免费观看视频 | 另类亚洲自拍| 欧美黄在线观看| 一区二区不卡在线视频 午夜欧美不卡'| 亚洲美女在线观看| 午夜欧美大尺度福利影院在线看| 久久精品一本久久99精品| 免费中文日韩| 国产精品影音先锋| 最新日韩欧美| 亚洲欧美三级在线| 美女被久久久| 一本一本久久a久久精品综合妖精 一本一本久久a久久精品综合麻豆 | 欧美福利一区二区| 一本大道久久a久久综合婷婷| 亚洲欧美日韩精品久久久| 久久尤物电影视频在线观看| 欧美剧在线观看| 国产一区二区久久精品| 在线视频日韩| 欧美激情成人在线| 销魂美女一区二区三区视频在线| 欧美激情国产日韩| 国产一区二区三区四区hd| 一本色道久久综合狠狠躁篇怎么玩 | 久色成人在线| 国产精品一区二区女厕厕| 亚洲一级在线观看| 久久综合色88| 国产手机视频一区二区| 在线中文字幕一区| 欧美激情乱人伦| 久久国产精品久久国产精品| 欧美视频久久| 亚洲免费激情| 欧美电影打屁股sp| 久久久夜夜夜| 国内一区二区三区| 欧美一区二区三区视频| 亚洲精品日本| 欧美精品在线极品| 亚洲免费不卡| 亚洲日本电影| 久久综合伊人77777蜜臀| 国产精品一二一区| 欧美一区在线看| 亚洲午夜激情网站| 国产精品乱码一区二区三区| 中文国产成人精品| 亚洲精品久久久久久久久久久久久 | 一区二区三区在线看| 亚洲性感激情| 艳女tv在线观看国产一区| 欧美成人免费播放| 亚洲精品久久久久久久久久久久| 免费高清在线视频一区·| 久久精品最新地址| 亚洲高清视频在线观看| 欧美国产日韩在线观看| 麻豆91精品91久久久的内涵| 亚洲国产日本| 亚洲理伦电影| 国产精品啊啊啊| 欧美一区二区三区视频免费播放| 午夜精品三级视频福利| 国产欧美亚洲视频| 久久婷婷国产综合精品青草| 久久性色av| 久热国产精品| 中文一区二区| 欧美伊人久久| 亚洲国产专区| 亚洲国产日韩欧美综合久久 | 国内久久精品视频| 欧美成人免费观看| 欧美日韩成人一区| 香蕉亚洲视频| 久久久久久噜噜噜久久久精品| 久久不射网站| 亚洲美女少妇无套啪啪呻吟| 一区二区三区欧美激情| 国产亚洲一区二区三区在线播放| 麻豆freexxxx性91精品| 欧美日韩黄视频| 久久精品2019中文字幕| 裸体丰满少妇做受久久99精品 | 久久精品二区三区| 欧美成人黄色小视频| 亚洲网站在线| 久久精品伊人| 亚洲综合日韩在线| 米奇777在线欧美播放| 亚洲自拍偷拍视频| 久久综合亚州| 久久精品在线播放| 欧美精品videossex性护士| 欧美在线视频播放| 欧美精品三级在线观看| 久热国产精品| 国产美女精品一区二区三区| 91久久精品网| 在线观看亚洲视频啊啊啊啊| 一区二区三区四区五区精品| 亚洲国产一区二区三区高清| 午夜日韩在线| 午夜精品一区二区三区在线播放 | 在线一区欧美| 欧美高清视频www夜色资源网| 欧美一区二区三区电影在线观看| 麻豆精品视频在线| 久久精品一区二区国产| 欧美午夜一区二区福利视频| 欧美高清成人| 亚洲福利视频三区| 欧美在线视频导航| 久久99伊人| 国产农村妇女毛片精品久久莱园子| 91久久精品日日躁夜夜躁欧美| 在线播放亚洲一区| 久久精品久久综合| 麻豆精品一区二区综合av| 国产日产欧产精品推荐色 | 亚洲国产精品视频一区| 性做久久久久久久久| 羞羞视频在线观看欧美| 国产精品成人v| 一本色道久久99精品综合 | 亚洲国产日韩一区| 亚洲欧洲综合| 美女999久久久精品视频| 国产亚洲成精品久久| 亚洲欧美激情视频在线观看一区二区三区 | 亚洲国产精品久久精品怡红院| 久久免费视频一区| 欧美电影在线免费观看网站| 在线观看国产成人av片| 久久久国产一区二区| 免费成人网www| 亚洲精品视频免费在线观看| 欧美精品电影在线| 99热免费精品| 午夜在线精品偷拍| 国内在线观看一区二区三区| 久久久噜噜噜久久中文字免| 欧美激情一区二区三区| 在线一区二区三区做爰视频网站| 国产精品扒开腿爽爽爽视频| 亚洲欧美自拍偷拍| 免费91麻豆精品国产自产在线观看| 激情久久综合| 欧美成人福利视频| 一本色道久久综合亚洲精品按摩 | 亚洲激情电影在线| 欧美激情亚洲自拍| 亚洲一区区二区| 麻豆精品视频在线观看| 日韩视频一区二区在线观看 | 浪潮色综合久久天堂| 91久久国产自产拍夜夜嗨| 亚洲欧美电影院| 国产亚洲精久久久久久| 久久久久久尹人网香蕉| 99视频精品全部免费在线| 欧美日韩网址| 欧美一区视频在线| 亚洲第一级黄色片| 性感少妇一区| 亚洲国产日韩欧美在线图片| 欧美精品色综合| 欧美专区日韩专区| 亚洲美女视频在线观看| 狂野欧美激情性xxxx| 中文欧美在线视频| 伊人成人开心激情综合网| 欧美精品一区二区在线播放| 亚洲欧美日韩在线一区| 亚洲国产精品久久久| 欧美一区二区三区视频在线| 亚洲免费黄色| 亚洲国产一区二区精品专区| 国产日韩欧美高清免费| 欧美日一区二区在线观看 | av成人动漫| 亚洲国产99| 国产自产女人91一区在线观看| 欧美日韩一区二区三区在线| 老牛国产精品一区的观看方式| 亚洲一区二区精品视频| 日韩视频―中文字幕| 91久久久久| 亚洲欧洲在线看|