• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    一種基于Stirling圖枚舉算法的分球入盒問題求解

    2017-11-01 17:14:41彭哲也謝民主
    關(guān)鍵詞:枚舉分球個(gè)球

    彭哲也 謝民主

    (湖南師范大學(xué)物理與信息科學(xué)學(xué)院 湖南 長(zhǎng)沙 410081)

    一種基于Stirling圖枚舉算法的分球入盒問題求解

    彭哲也 謝民主*

    (湖南師范大學(xué)物理與信息科學(xué)學(xué)院 湖南 長(zhǎng)沙 410081)

    已有的分球入盒問題解法通常只關(guān)注分球的總方案數(shù),目前尚沒有公開的計(jì)算機(jī)算法來(lái)枚舉出所有具體的分球方案,而方案的枚舉是生物信息學(xué)中一些分區(qū)優(yōu)化算法的基礎(chǔ)。受第二類Stirling數(shù)的遞推公式的啟發(fā),提出一個(gè)新的數(shù)據(jù)結(jié)構(gòu)——Stirling圖。在此基礎(chǔ)上設(shè)計(jì)一個(gè)算法來(lái)枚舉p個(gè)不同球分配到q個(gè)相同盒子里的所有不同的方案。當(dāng)p和q較大,枚舉出所有的方案不可行時(shí),設(shè)計(jì)另一個(gè)算法在整個(gè)方案空間實(shí)現(xiàn)均勻采樣,輸出指定個(gè)數(shù)的方案。測(cè)試結(jié)果表明,這些算法在內(nèi)存為8 GB的普通PC上可在合理的時(shí)間內(nèi)枚舉出上百萬(wàn)組不同的方案。

    分球入盒問題 第二類Stirling數(shù) 枚舉算法 Stirling圖 均勻采樣

    0 引 言

    分球入盒是組合數(shù)學(xué)中的經(jīng)典問題,在組合優(yōu)化中有廣泛的應(yīng)用。把p個(gè)球放入q個(gè)盒子里,根據(jù)球和盒子是否相同可以分為四種不同的模型:(1) 球相同盒子相同;(2) 球相同盒子不同;(3) 球不同盒子不同;(4) 球不同盒子相同。每種模型又進(jìn)一步分為盒子是否允許為空這兩種情形。本文主要討論球不同盒子相同的這類模型。

    將p個(gè)有區(qū)別的球放入q個(gè)相同的盒子中,要求所有盒子都非空的方案數(shù)為第二類Stirling數(shù),記作S(p,q)[1]。例如,將編號(hào)為A、B、C、D的4個(gè)小球放入2個(gè)相同的盒子里的所有不同放法可以表示如下(盒子用集合表示):[{A},{B,C,D}],[{B},{A,C,D}],[{C},{A,B,D}],[{D},{A,B,C}],[{A,B},{C,D}],[{A,C},{B,D}],[{A,D},{B,C}]。因此S(4,2)=7。

    第二類Stirling數(shù)在數(shù)論、組合論、代數(shù)拓?fù)?、算法最?yōu)化等領(lǐng)域都有相關(guān)的應(yīng)用[2]。顯然,當(dāng)pq且q>1時(shí),第二類Stirling可由以下遞歸算法求解:考慮第p個(gè)球,該球不是獨(dú)占一個(gè)盒子就是和其他球在同一個(gè)盒子中。第p個(gè)球獨(dú)占一個(gè)盒子:其他p-1個(gè)球分配到其他q-1個(gè)非空盒子中的所有方案數(shù)為S(p-1,q-1);第p個(gè)球不獨(dú)占一個(gè)盒子的情況則是:其他p-1個(gè)球放入q個(gè)非空的盒子中,第p個(gè)球放入q個(gè)盒子的任意一個(gè)中,這樣共有qS(p-1,q)種方法。這樣就得到了第二類Stirling數(shù)的遞推公式:

    S(p,q)=S(p-1,q-1)+qS(p-1,q)

    (1)

    將p個(gè)有區(qū)別的球放入q個(gè)相同的盒子中, 當(dāng)盒子允許為空時(shí),令其不同的方法數(shù)為F(p,q),顯然:

    F(p,q)=S(p,1)+S(p,2)+…+S(p,q)

    (2)

    第二類Stirling數(shù)的計(jì)算已經(jīng)得到了廣泛深入的研究[3],但是目前尚沒有公開的計(jì)算機(jī)算法來(lái)枚舉出所有具體的分球方案。枚舉是計(jì)算復(fù)雜性理論中的一類基本運(yùn)算,是理論計(jì)算機(jī)科學(xué)中一個(gè)重要分支。自從1979 年 Valiant 首次提出枚舉運(yùn)算以來(lái)[6],人們對(duì)枚舉運(yùn)算進(jìn)行了大量的研究[7],其中最典型的枚舉方法有枚舉問題實(shí)例的所有解或統(tǒng)計(jì)問題實(shí)例中解的個(gè)數(shù)[8]。

    近年來(lái),隨著DNA測(cè)序的快速發(fā)展,利用多倍體基因組的DNA測(cè)序片段如何組裝出構(gòu)成多倍體的多個(gè)單體型[11]和利用環(huán)境群體微生物的深度DNA測(cè)序序列重構(gòu)該群體的微生物基因組[12]等問題的研究成為生物信息學(xué)的研究熱點(diǎn)。這些問題的數(shù)學(xué)建模通常是把DNA測(cè)序序列分配到多個(gè)分區(qū)中以滿足特定的優(yōu)化目標(biāo)。因此,分球入盒方案的枚舉成為一些求解這類問題優(yōu)化算法的基本操作[11]。

    本文受第二類Stirling數(shù)的遞推公式的啟發(fā)提出了一個(gè)新的數(shù)據(jù)結(jié)構(gòu)——Stirling圖。在此基礎(chǔ)上設(shè)計(jì)一個(gè)算法來(lái)枚舉p個(gè)不同球分配到q個(gè)相同盒子里的所有不同的方案;當(dāng)p和q較大時(shí),不可能枚舉出所有的方案,該算法則在整個(gè)方案空間實(shí)現(xiàn)均勻采樣,輸出指定個(gè)數(shù)的方案。在內(nèi)存為8 GB、主頻為2.5 GHz的普通PC上的測(cè)試結(jié)果表明,該算法可以在幾秒內(nèi)枚舉出上百萬(wàn)組的不同的方案。

    1 Stirling圖

    本文用一個(gè)劃分函數(shù)(Partition function)P(i):{1,…,p}→{1,…,q}來(lái)表示p個(gè)不同的球分配到q個(gè)盒子中的一個(gè)具體方案。例如P(i)|i=1,2,3=1表示將編號(hào)為1,2,3的球放入第一個(gè)盒子里,P(4)=2表示將編號(hào)為4的球放入第二個(gè)盒子里。為了簡(jiǎn)便,下面用P(1..k)=1表示P(i)|i=1,…,k=1;P(1..k)=1..k表示P(i)|i=1,…,k=i。令p個(gè)不同的球分配到q個(gè)非空的相同盒子中的所有不同的方案的集合為S(p,q),集合S(p,q)中不同方案的個(gè)數(shù)|S(p,q)|=S(p,q)。

    當(dāng)p0時(shí),S(p,q)=1,即S(p,q)中只有一種分配方案,不同的球放到不同的盒子里,給球i所在盒子貼上標(biāo)簽i,則這種分配方案用P(1..p) = 1..p表示;當(dāng)p>q=1,這時(shí)S(p,q)中同樣只存在一種分配方案,即所有的球放在唯一的盒子中,給該盒子貼上標(biāo)簽1,則這種分配方案用P(1..p) = 1來(lái)表示;當(dāng)p>q>1時(shí),所有不同分配方案的集合S(p,q)可以用式(1)的證明過程構(gòu)造出來(lái):

    S(p,q) ={P(p)=q·S(p-1,q-1)}∪

    {S(p-1,q)·(P(p)=1…q)}

    (3)

    其中{(P(p)=q)·S(p-1,q-1)}表示將球p單獨(dú)放到一個(gè)盒子中,且把該盒子貼上標(biāo)簽q(用P(p)=q表示),然后把其他p-1個(gè)球分配到其他q-1個(gè)相同非空的盒子中(用S(p-1,q-1)表示)的所有不同方案的集合;而{(P(p)=1..q)·S(p-1,q-1)}則表示將前p-1個(gè)球分配到q個(gè)相同的非空盒子中(用S(p-1,q)表示,每個(gè)盒子由此過程獲得了不同的標(biāo)簽)。然后將第p個(gè)球放到標(biāo)簽為i(i=1,…,q)的盒子中(用P(p)=1..q表示)的所有不同方案的集合。為了使用式(3)構(gòu)造出分配方案的集合S(p,q),本文引入如圖1所示的Stirling圖。

    圖1 Stirling圖

    Stirling圖是一個(gè)有向無(wú)環(huán)圖,其結(jié)點(diǎn)用S(j,k)來(lái)標(biāo)記(j≥k≥1)。Stirling圖的結(jié)點(diǎn)分為兩類,一類是終端結(jié)點(diǎn)S(k,k)和S(j,1),一類是非終端結(jié)點(diǎn)S(j,k)(j>k>1)。非終端結(jié)點(diǎn)S(j,k)結(jié)點(diǎn)有左右兩個(gè)孩子結(jié)點(diǎn),其左孩子結(jié)點(diǎn)為S(j-1,k-1),右孩子結(jié)點(diǎn)為S(j-1,k)。圖中用一條標(biāo)簽為P(j)=k的有向邊由非終端結(jié)點(diǎn)S(j,k)指向其左孩子結(jié)點(diǎn),用標(biāo)簽分別為P(j)=1,P(j)=2,…,P(j)=k的k條有向邊指向其右孩子結(jié)點(diǎn)。對(duì)于兩個(gè)結(jié)點(diǎn)S和S′,如果存在一條路徑從S到S′,則S為S′的先驅(qū)結(jié)點(diǎn),S′為S的后繼結(jié)點(diǎn)。

    在Stirling圖中從非終端結(jié)點(diǎn)S(j,k)到其后繼終端結(jié)點(diǎn)的一條路徑代表把p個(gè)不同的球分配到q個(gè)非空的相同盒子的一種分配方案。如從S(3,2)到S(2,1)的路徑代表一個(gè)分配方案:P(1..2)=1,P(3)=2;即把前兩個(gè)球放到一個(gè)盒子里,第三個(gè)球放到單獨(dú)的盒子了。由式(2)可知,從非終端結(jié)點(diǎn)S(j,k)到其所有后繼終端結(jié)點(diǎn)的路徑就代表把j個(gè)不同的球分配到k個(gè)非空相同盒子的所有不同分配方案。

    對(duì)Stirling圖,本文定義如下操作:

    1)initiate(p,q):構(gòu)造出一個(gè)Stirling圖,包含結(jié)點(diǎn)S(p,q)及其所有后繼結(jié)點(diǎn),及相應(yīng)的有向邊。

    2)destroy():刪除Stirling圖,釋放內(nèi)存。

    3)count(j,k): 計(jì)算出S(j,k),即|S(p,q)|。

    4)pathTransverse(j,k):對(duì)Stirling圖遍歷從結(jié)點(diǎn)S(j,k)開始以任意終端節(jié)點(diǎn)結(jié)束的所有路徑,即枚舉出集合S(j,k)中的所有方案。

    5)randomSample(j,k,n): 從集合S(p,q)中按照均勻分布隨機(jī)選擇n個(gè)不同的方案輸出。

    2 數(shù)據(jù)結(jié)構(gòu)及算法實(shí)現(xiàn)

    本文用一個(gè)(p+1)×(q+1)的二維數(shù)組M(假設(shè)下標(biāo)是從0開始)來(lái)存儲(chǔ)Stirling圖的結(jié)點(diǎn)S(j,k):0

    為了支持pathTransverse和count操作,元素Mj,k包含兩個(gè)整型變量index和count,其中count用來(lái)記錄S(j,k)的值,初始化為-1;index的取值范圍為0到k,0表示進(jìn)行pathTransverse操作時(shí),下次訪問的邊是指向左孩子的邊,index=i(1≤i≤k)表示下次訪問的邊是第i條指向其右孩子的邊。

    相關(guān)操作實(shí)現(xiàn)的偽代碼如下所示:

    structnode{

    int index = 0;

    int count = -1;

    }

    node[][] M;

    void initiate(p, q){M = new node[p+1][q+1];}

    void destroy() { delete[] M;}

    intcount(int j, int k) {

    if (M[j][k].count > 0)

    //已經(jīng)計(jì)算過,返回計(jì)算值

    returnM[j][k].count;

    if (j = k || k = 1) {

    //終端結(jié)點(diǎn)

    M[j][k].count = 1; return 1; }

    M[j][k].count=count(j-1,k-1)+k*count(j-1,k);

    returnM[j][k].count; }

    voidpathTransverse(int j,int k){

    // P記錄球分配方案,第i個(gè)球分配到標(biāo)簽為P[i]的盒

    //子中

    int[] P = new int[j+1];

    output_pathes (j, k, P, j); }

    voidoutput_pathes (int j, int k, int[] P, int balls){

    if (j = k) {

    //終端結(jié)點(diǎn),輸出P(1..k) = 1..k

    for(i= 1 .. k) { P[i] = i; }

    輸出分配方案P[1], …, P[balls];

    return; }

    if (k = 1) {

    //終端結(jié)點(diǎn),輸出P(1..k) = 1

    for(i= 1 .. j) { P[i] = 1; }

    輸出分配方案P[1], …, P[balls];

    return; }

    if (M[j][k].index = 0) {

    P[j] = k;

    //訪問指向左孩子的邊

    M[j][k].index++;

    output_pathes(j-1, k-1, P, balls); }

    //訪問左孩子

    while(M[j][k].index ≤ k) {

    //訪問指向右孩子的邊

    P[j] = M[j][k].index; M[j][k].index++;

    output_pathes(j-1, k, P, balls); }

    //訪問右孩子

    M[j][k].index=0;

    //恢復(fù)初始值

    }

    void randomSample(int j, int k, int n){

    //n必須小于count(j, k)

    int[] P = new int[j+1];

    sample_pathes(j, k, P, j, n);

    }

    voidsample_pathes(intj, intk, int[] P, int b, int n){

    if(n< 1){ return; }

    if (j = k) {

    //終端結(jié)點(diǎn),輸出P(1..k) = 1..k

    for(i = 1 .. k) { P[i] = i; }

    輸出分配方案P[1], …, P[b]; return; }

    if (k = 1) {

    //終端結(jié)點(diǎn),輸出P(1..k) = 1

    for(i= 1 .. j) { P[i] = 1; }

    輸出分配方案P[1], …, P[b]; return; }

    intn_left = n;

    int[] n_a = new int[k+1];

    //通過均勻分布隨機(jī)采樣,獲得通過其指向左右孩子結(jié)

    //點(diǎn)的

    // k+1條有向邊的采樣個(gè)數(shù)

    get_sample_number(j, k, n, n_a)

    if (M[j][k].index = 0) { P[j] = k;

    //訪問指向左孩子的邊

    M[j][k].index++;

    sample_pathes(j-1, k-1, P, b, n_a[0]); }

    //訪問左孩子

    while(M[j][k].index < k) {

    //訪問指向右孩子的邊

    P[j] = M[j][k].index;

    M[j][k].index++;

    sample_pathes(j-1, k, P, b, n_a[P[j]]); }

    //訪問右孩子

    M[j][k].index=0;

    //恢復(fù)初始值

    }

    voidget_sample_number(intj, intk, intn,int[] n_a){

    float r = n/count(j, k);

    //采樣率

    令count(j-1, k-1)*r的整數(shù)部分為a1,小數(shù)部分為e1,count(j, k-1)*r的整數(shù)部分為a2,小數(shù)部分為e2;

    n_a[0] = a1;

    floatpl = e1/(e1+k*e2), pr = e2/(e1+k*e2);

    for(i= 1.. k) {

    n_a[i] = a2; }

    intn_left = n - a1 - k*a2;

    while(n_left> 0){

    //得到一個(gè)在[0, 1]之間均勻分布的隨機(jī)數(shù)rand

    floatrand = random();

    if(rand ≤ pl) {n_a[0]++; }

    else { intind = (rand - pl)/pr +1; n_a[ind]++; }

    n_left--;}

    }

    對(duì)于盒子允許為空的情況,利用Stirling圖枚舉出所有的分配方案和按照均勻分布隨機(jī)采樣n個(gè)不同的方案的算法也容易實(shí)現(xiàn)。在盒子允許為空的情況下,令p個(gè)不同的球分配到q個(gè)相同的盒子中的所有不同方案的集合為F(p,q),顯然有:

    F(p,q)=S(p,q)∪S(p,2)∪…∪S(p,q)

    (4)

    枚舉出F(p,q)的算法如下所示:

    F(int p, int q){

    for(i = 1…q) { pathTransverse(p, i); }

    }

    當(dāng)需要按照均勻分布隨機(jī)采樣輸出n個(gè)不同的方案時(shí),可按照get_sample_number的思路把n按照S(p,1),S(p,2),…,S(p,q)集合的大小劃分成n1,n2,…,nq,然后再依次調(diào)用randomSample(p,1,n1),randomSample(p,2,n2),…,randomSample(p,q,nq)即可。

    3 實(shí)驗(yàn)結(jié)果

    本文算法用Java編程實(shí)現(xiàn),測(cè)試計(jì)算機(jī)內(nèi)存為8 GB,CPU為Intel酷睿i5 3210M 2.5 GHz的PC機(jī)。圖2-圖4分別表示的是F(3,2)、F(4,2)、F(4,3)在控制臺(tái)輸出的枚舉結(jié)果。

    圖2 F(3,2)枚舉結(jié)果

    圖3 F(4,2)枚舉結(jié)果

    圖4 F(4,3)枚舉結(jié)果

    由于目前沒有類似的算法可以橫向比較,我們通過幾個(gè)用例來(lái)測(cè)試本文算法的性能。如表1所示,當(dāng)p和q增大時(shí),F(xiàn)(p,q)枚舉出所有分配方案所需的時(shí)間基本上與方案總數(shù)成線性增長(zhǎng);當(dāng)p=12,q=5時(shí),F(xiàn)(p,q)可以在約1秒鐘的時(shí)間內(nèi)枚舉出2百余萬(wàn)種不同的方案。

    表1 F(p,q)枚舉時(shí)間

    當(dāng)方案總數(shù)過大,枚舉所有的方案不可行時(shí),我們?cè)谡麄€(gè)空間進(jìn)行按照均勻分布隨機(jī)采樣100萬(wàn)個(gè)不同的方案進(jìn)行輸出,表2顯示對(duì)于不同的(p,q),隨機(jī)采樣算法randomSample(p,q,1 000 000)的運(yùn)行時(shí)間。

    表2 randomSample(p,q,1 000 000)運(yùn)行時(shí)間

    4 結(jié) 語(yǔ)

    本文提出一種新的數(shù)據(jù)結(jié)構(gòu)Stirling圖,實(shí)現(xiàn)了該數(shù)據(jù)結(jié)構(gòu)的操作算法,能枚舉出不同的球分配到相同盒子中的所有不同方案。當(dāng)方案總數(shù)過多無(wú)法枚舉時(shí),本文實(shí)現(xiàn)了隨機(jī)采樣操作算法,該算法按照均勻分布隨機(jī)采樣輸出指定個(gè)數(shù)的不同的分配方案。測(cè)試結(jié)果表明,這些算法能有效解決球不同盒子相同的“分球入盒”方案枚舉問題和均勻采樣問題。實(shí)際應(yīng)用中有很多優(yōu)化問題如多倍體單體型組裝[11]、宏基因組組裝[12]等均可轉(zhuǎn)化為“分球入盒”問題求解。本文提出的算法可以為這些問題的求解提供有效的分球入盒方案枚舉基本操作的支持。

    [1] Wikipedia. Stirling numbers of the second kind[EB/OL].(2016-12-7)[2016-12-14].https://en.wikipedia.org/wiki/Stirling_numbers_of_the_second_kind.

    [2] 王娟. 第二類Stirling數(shù)及其推廣[D]. 大連理工大學(xué),2009.

    [3] Abramowitz M, Stegun I A. Handbook of mathematical functions with formulas, graphs, and mathematical tables (9th printing)[M]. New York: Dover, 1972:824-825.

    [4] Boyadzhiev K N. Close Encounters with the Stirling Numbers of the Second Kind[J]. Mathematics Magazine,2012, 85(4):252-266.

    [5] Bleick W W. Asymptotics of Stirling Numbers of the Second Kind[J]. Proceedings of the American Mathematical Society,1974,42(2):575-580.

    [6] Valiant L.The complexity of enumeration and reliability problems[J].SIAM Journal on Computing, 1979,8(3):410-421.

    [7] Flum J, Grohe M .The parameterized complexity of counting problems[J]. SIAM Journal on Computing, 2004,33(4):892-922.

    [8] 劉運(yùn)龍, 王建新. 3-維匹配問題的一種固定參數(shù)枚舉算法[J].計(jì)算機(jī)科學(xué),2010,37(5):210-213.

    [9] 謝民主, 羅鋒, 唐烽. 單體型組裝最大片段割參數(shù)化精確算法[J].小型微型計(jì)算機(jī)系統(tǒng),2014,35(2):353-357.

    [10] 張景云. 改進(jìn)的堆的枚舉算法的研究[J]. 計(jì)算機(jī)應(yīng)用與軟件,2012,29(7):264-265,273.

    [11] Xie Minzhu,Wu Qiong,Wang Jianxin, et al. H-PoP and H-PoPG: heuristic partitioning algorithms for single indivi-dual haplotyping of polyploids[J/OL]. Bioinformatics. [2016-08-16] Doi: 10.1093/bioinformatics/btw537.

    [12] Rose R, Constantinides B, Tapinos A, et al. Challenges in the analysis of viral metagenomes[J/OL].Virus Evolution.[2016-08-03].DOI:http://dx.doi.org/10.1093/ve/vew022.

    RESEARCHONDISTRIBUTINGBALLSINTOBOXESENUMERATIONALGORITHM

    Peng Zheye Xie Minzhu*

    (CollegeofPhysicsandInformationScience,HunanNormalUniversity,Changsha410081,Hunan,China)

    The existing researches on the problem of distributing balls into boxes usually focus on the total number of different ways to distribute balls into boxes, but there are no public computer algorithms to enumerate them. However, enumerating them is the foundation to design some partition optimal algorithms in Bioinformatics. Inspired by the recursive formula of the Stirling numbers of the second kind, the paper proposes a new data structure—Stirling diagram, and based on the data structure, designs an algorithm to enumerate all different ways to distribute p different balls into q same boxes. When p and q are larger and none of the schemes is feasible, we design another algorithm to achieve uniform sampling of a given number of different ways. Test results show that these algorithms can enumerate millions of different distributing ways in a reasonable period of time on a PC with 8 GB memory.

    Distributing balls into boxes problem Stirling numbers of the second kind Enumerating algorithm Stirling diagram Uniform sampling

    TP306.1

    A

    10.3969/j.issn.1000-386x.2017.10.044

    2016-12-15。國(guó)家自然科學(xué)基金項(xiàng)目(61370172)。彭哲也,碩士生,主研領(lǐng)域:生物信息學(xué)。謝民主,教授。

    猜你喜歡
    枚舉分球個(gè)球
    基于理解性教學(xué)的信息技術(shù)教學(xué)案例研究
    速讀·上旬(2022年2期)2022-04-10 16:42:14
    一種高效的概率圖上Top-K極大團(tuán)枚舉算法
    NBA比賽3分球技術(shù)運(yùn)用發(fā)展趨勢(shì)
    該你了
    踢你個(gè)球
    踢你個(gè)球
    踢你個(gè)球
    踢你個(gè)球(1)
    基于太陽(yáng)影子定位枚舉法模型的研究
    USB開發(fā)中易混淆的概念剖析
    99热网站在线观看| 91在线精品国自产拍蜜月| 精品一区在线观看国产| 丰满少妇做爰视频| 91精品国产九色| 久久ye,这里只有精品| 91在线精品国自产拍蜜月| 校园人妻丝袜中文字幕| 黄色视频在线播放观看不卡| 国产亚洲av片在线观看秒播厂| 老女人水多毛片| 在线观看一区二区三区激情| 国产欧美日韩精品一区二区| 免费大片黄手机在线观看| 高清午夜精品一区二区三区| 欧美亚洲 丝袜 人妻 在线| 欧美日韩精品成人综合77777| 一级毛片 在线播放| 七月丁香在线播放| 亚洲精品亚洲一区二区| 内射极品少妇av片p| 一级av片app| 久久久久网色| 免费观看a级毛片全部| 黑丝袜美女国产一区| 久久久久久久久久久丰满| 高清在线视频一区二区三区| 人妻 亚洲 视频| 久久精品国产亚洲网站| 国产黄色视频一区二区在线观看| 亚洲欧美一区二区三区国产| 又爽又黄a免费视频| 亚洲丝袜综合中文字幕| 人人妻人人澡人人爽人人夜夜| 制服丝袜香蕉在线| 看非洲黑人一级黄片| 免费不卡的大黄色大毛片视频在线观看| 多毛熟女@视频| 我要看日韩黄色一级片| 中文字幕精品免费在线观看视频 | 人体艺术视频欧美日本| 亚洲一区二区三区欧美精品| 中国三级夫妇交换| 啦啦啦中文免费视频观看日本| av国产久精品久网站免费入址| 视频中文字幕在线观看| 日韩 亚洲 欧美在线| 欧美 日韩 精品 国产| 最近最新中文字幕大全电影3| 亚洲精华国产精华液的使用体验| 久热这里只有精品99| 成人免费观看视频高清| 一级片'在线观看视频| 欧美xxxx黑人xx丫x性爽| 国产中年淑女户外野战色| 国产欧美日韩一区二区三区在线 | 一级毛片 在线播放| 最近最新中文字幕免费大全7| 91久久精品国产一区二区三区| 国产成人a∨麻豆精品| 中文乱码字字幕精品一区二区三区| 亚洲四区av| 亚洲精品中文字幕在线视频 | 亚洲熟女精品中文字幕| 欧美xxⅹ黑人| 精品99又大又爽又粗少妇毛片| 亚洲婷婷狠狠爱综合网| 日韩中文字幕视频在线看片 | 色哟哟·www| 精品一区在线观看国产| 在线免费十八禁| 晚上一个人看的免费电影| 国产淫片久久久久久久久| 精品久久久久久久末码| 国产国拍精品亚洲av在线观看| 婷婷色麻豆天堂久久| 亚洲最大成人中文| 黄色欧美视频在线观看| 亚洲成人手机| 亚洲av免费高清在线观看| 国模一区二区三区四区视频| 免费观看无遮挡的男女| 亚洲综合精品二区| 18禁动态无遮挡网站| 99热这里只有是精品50| 人人妻人人澡人人爽人人夜夜| 亚洲激情五月婷婷啪啪| xxx大片免费视频| 激情五月婷婷亚洲| 午夜老司机福利剧场| 在线亚洲精品国产二区图片欧美 | 九九在线视频观看精品| 亚洲天堂av无毛| 亚洲精品色激情综合| 观看美女的网站| 国语对白做爰xxxⅹ性视频网站| 久久精品久久久久久久性| 国产成人精品一,二区| 亚洲最大成人中文| 两个人的视频大全免费| 久久精品国产自在天天线| 少妇的逼好多水| 美女福利国产在线 | 欧美丝袜亚洲另类| 国产 一区 欧美 日韩| 韩国av在线不卡| 国模一区二区三区四区视频| a级毛片免费高清观看在线播放| 我的女老师完整版在线观看| 亚洲欧洲国产日韩| 亚洲真实伦在线观看| 99精国产麻豆久久婷婷| av免费观看日本| 在线观看三级黄色| 亚洲国产最新在线播放| 毛片女人毛片| 午夜激情福利司机影院| 成人美女网站在线观看视频| 国产av码专区亚洲av| 91aial.com中文字幕在线观看| 成人国产麻豆网| 欧美亚洲 丝袜 人妻 在线| 午夜福利高清视频| 人体艺术视频欧美日本| 一个人看的www免费观看视频| 国产免费又黄又爽又色| 国产片特级美女逼逼视频| 日本猛色少妇xxxxx猛交久久| 日本与韩国留学比较| 国产乱人视频| 国产免费一级a男人的天堂| 欧美精品人与动牲交sv欧美| 在线观看一区二区三区| 免费大片18禁| 国产黄色视频一区二区在线观看| 青春草视频在线免费观看| 亚洲色图av天堂| 一区二区三区精品91| 久久人人爽人人爽人人片va| 久久人人爽av亚洲精品天堂 | 日韩视频在线欧美| 欧美三级亚洲精品| 亚洲精品乱久久久久久| 99久久综合免费| 国产中年淑女户外野战色| 久久久a久久爽久久v久久| 欧美最新免费一区二区三区| 少妇被粗大猛烈的视频| 国国产精品蜜臀av免费| 一本一本综合久久| 涩涩av久久男人的天堂| 免费黄色在线免费观看| 久久久久久久亚洲中文字幕| 欧美日韩综合久久久久久| 一区二区三区精品91| 日本猛色少妇xxxxx猛交久久| 国产免费一区二区三区四区乱码| 国产精品久久久久久久久免| 人妻制服诱惑在线中文字幕| 最黄视频免费看| 国国产精品蜜臀av免费| 老司机影院成人| av线在线观看网站| 麻豆成人av视频| 一个人看的www免费观看视频| 少妇人妻一区二区三区视频| 嫩草影院新地址| a级毛色黄片| 成年女人在线观看亚洲视频| 中文天堂在线官网| 国产国拍精品亚洲av在线观看| 亚洲婷婷狠狠爱综合网| 欧美成人午夜免费资源| 女性生殖器流出的白浆| 国产精品偷伦视频观看了| 男女边摸边吃奶| 国产精品久久久久久久电影| 人妻 亚洲 视频| 一个人免费看片子| 久热久热在线精品观看| 国产精品蜜桃在线观看| 成年免费大片在线观看| 免费看av在线观看网站| 免费观看的影片在线观看| 中文乱码字字幕精品一区二区三区| 亚洲欧美一区二区三区国产| 久久婷婷青草| 寂寞人妻少妇视频99o| 中文字幕久久专区| 日韩一区二区视频免费看| 熟妇人妻不卡中文字幕| 日日啪夜夜爽| 中文字幕制服av| 国产伦理片在线播放av一区| 黑丝袜美女国产一区| 少妇人妻 视频| 色视频在线一区二区三区| 国产视频内射| 亚洲人成网站在线观看播放| 纵有疾风起免费观看全集完整版| 草草在线视频免费看| 啦啦啦视频在线资源免费观看| 欧美日韩亚洲高清精品| 久久久久久久久久久免费av| 久久99热这里只有精品18| 六月丁香七月| 国产大屁股一区二区在线视频| 精品国产露脸久久av麻豆| 高清日韩中文字幕在线| 我要看黄色一级片免费的| 国产 一区 欧美 日韩| 伦理电影大哥的女人| 午夜福利高清视频| 中国国产av一级| 三级国产精品片| 亚洲欧美成人精品一区二区| 国产精品偷伦视频观看了| 成人漫画全彩无遮挡| 秋霞在线观看毛片| 国产av一区二区精品久久 | 国产精品偷伦视频观看了| 激情 狠狠 欧美| 国产黄色视频一区二区在线观看| 久久亚洲国产成人精品v| 久久国产精品大桥未久av | 免费观看性生交大片5| 最近2019中文字幕mv第一页| 80岁老熟妇乱子伦牲交| av在线观看视频网站免费| 黑人高潮一二区| 大香蕉97超碰在线| 亚洲人成网站高清观看| 久久99热这里只有精品18| 色视频在线一区二区三区| 干丝袜人妻中文字幕| 日本免费在线观看一区| 又粗又硬又长又爽又黄的视频| 在线 av 中文字幕| 大又大粗又爽又黄少妇毛片口| 国产精品一区二区性色av| 综合色丁香网| 欧美xxxx性猛交bbbb| 久热这里只有精品99| 国产黄色免费在线视频| 日韩强制内射视频| 美女福利国产在线 | 亚洲第一av免费看| 国产精品久久久久久精品电影小说 | 免费观看在线日韩| 男人添女人高潮全过程视频| 国产成人一区二区在线| 两个人的视频大全免费| 国产精品一区二区在线不卡| 最后的刺客免费高清国语| 一级二级三级毛片免费看| 老司机影院成人| 简卡轻食公司| 女的被弄到高潮叫床怎么办| 亚洲精品国产av蜜桃| 妹子高潮喷水视频| 久久久久国产精品人妻一区二区| 亚洲精品456在线播放app| 黄色欧美视频在线观看| 日日啪夜夜撸| 亚洲国产精品一区三区| av不卡在线播放| 中文在线观看免费www的网站| 99视频精品全部免费 在线| 国产精品偷伦视频观看了| 久久精品久久精品一区二区三区| 午夜老司机福利剧场| 少妇人妻精品综合一区二区| 又黄又爽又刺激的免费视频.| 又爽又黄a免费视频| 欧美日韩一区二区视频在线观看视频在线| 夫妻午夜视频| 如何舔出高潮| 欧美三级亚洲精品| 大码成人一级视频| 久久久久国产精品人妻一区二区| 97精品久久久久久久久久精品| 小蜜桃在线观看免费完整版高清| 五月开心婷婷网| 久久久成人免费电影| 2021少妇久久久久久久久久久| 18禁裸乳无遮挡动漫免费视频| 亚洲中文av在线| 久久久久久久国产电影| 国产av一区二区精品久久 | 国产精品99久久99久久久不卡 | 国产成人freesex在线| 久久人妻熟女aⅴ| 欧美最新免费一区二区三区| 欧美精品一区二区免费开放| 国产欧美另类精品又又久久亚洲欧美| 欧美一区二区亚洲| 久久6这里有精品| 99热网站在线观看| 边亲边吃奶的免费视频| 一区二区av电影网| 国产男女内射视频| 在线观看国产h片| 国产免费视频播放在线视频| 2021少妇久久久久久久久久久| 婷婷色av中文字幕| 毛片女人毛片| 啦啦啦啦在线视频资源| 国产精品熟女久久久久浪| 亚洲色图综合在线观看| 丝袜脚勾引网站| 热re99久久精品国产66热6| 在线观看一区二区三区激情| 大香蕉久久网| 制服丝袜香蕉在线| 久久久久久伊人网av| 国产高清有码在线观看视频| tube8黄色片| 精品人妻偷拍中文字幕| 91在线精品国自产拍蜜月| 成人影院久久| 香蕉精品网在线| av福利片在线观看| 少妇人妻 视频| 啦啦啦中文免费视频观看日本| 亚洲av综合色区一区| 成人午夜精彩视频在线观看| 在线观看免费日韩欧美大片 | 王馨瑶露胸无遮挡在线观看| 欧美亚洲 丝袜 人妻 在线| 一边亲一边摸免费视频| 亚洲欧美日韩卡通动漫| a级毛片免费高清观看在线播放| 亚洲国产欧美人成| a 毛片基地| 日韩欧美一区视频在线观看 | 99久久中文字幕三级久久日本| 在线观看美女被高潮喷水网站| 黄片wwwwww| 精品国产乱码久久久久久小说| 久久99热这里只频精品6学生| 麻豆成人午夜福利视频| 偷拍熟女少妇极品色| 精品久久久久久久末码| 国产黄色视频一区二区在线观看| 最近最新中文字幕大全电影3| 欧美激情极品国产一区二区三区 | 日韩欧美精品免费久久| 亚洲国产精品999| av免费在线看不卡| 欧美另类一区| 性色av一级| 国产男女超爽视频在线观看| 久久久国产一区二区| 99久久精品国产国产毛片| av在线app专区| 国产成人一区二区在线| 免费av不卡在线播放| 一个人看的www免费观看视频| 国产 一区精品| 嫩草影院入口| 五月天丁香电影| 国产视频内射| 久久婷婷青草| 国产精品人妻久久久久久| 欧美亚洲 丝袜 人妻 在线| 色哟哟·www| 日韩亚洲欧美综合| 久久久久久久久久久丰满| 国产成人精品久久久久久| av一本久久久久| 国产无遮挡羞羞视频在线观看| 国产精品成人在线| 在现免费观看毛片| 人妻夜夜爽99麻豆av| 欧美zozozo另类| 人妻制服诱惑在线中文字幕| 亚洲精品国产成人久久av| 免费在线观看成人毛片| 亚洲美女黄色视频免费看| 午夜激情福利司机影院| 久久久久久九九精品二区国产| a 毛片基地| 国产黄片美女视频| 噜噜噜噜噜久久久久久91| av视频免费观看在线观看| www.色视频.com| 亚洲精品一区蜜桃| 深夜a级毛片| 日韩免费高清中文字幕av| 精品熟女少妇av免费看| 啦啦啦中文免费视频观看日本| 九草在线视频观看| 青春草亚洲视频在线观看| 亚洲精华国产精华液的使用体验| 国产精品人妻久久久影院| 久久久久久久大尺度免费视频| 哪个播放器可以免费观看大片| 亚洲欧美中文字幕日韩二区| 国产欧美亚洲国产| 在线观看免费高清a一片| 久久毛片免费看一区二区三区| 日本av免费视频播放| 国产精品嫩草影院av在线观看| 欧美人与善性xxx| 国产精品爽爽va在线观看网站| 国产色爽女视频免费观看| 2022亚洲国产成人精品| 一级av片app| 22中文网久久字幕| a级毛色黄片| 日韩欧美 国产精品| 少妇裸体淫交视频免费看高清| 日韩,欧美,国产一区二区三区| 777米奇影视久久| 国产伦精品一区二区三区四那| 亚洲精品乱久久久久久| 联通29元200g的流量卡| 国产精品国产三级国产av玫瑰| 99热这里只有是精品50| 九九久久精品国产亚洲av麻豆| 国产69精品久久久久777片| 少妇人妻一区二区三区视频| 久久久久国产精品人妻一区二区| 日本与韩国留学比较| 国产片特级美女逼逼视频| 日韩av在线免费看完整版不卡| 春色校园在线视频观看| 久久久久久久久久成人| 国产视频首页在线观看| 久久久久国产网址| 亚洲成人手机| 亚洲色图av天堂| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美日韩东京热| 亚洲高清免费不卡视频| 国产在线视频一区二区| 精品久久久精品久久久| 国产成人精品福利久久| 我要看黄色一级片免费的| 亚洲色图av天堂| 日韩免费高清中文字幕av| 久久av网站| 久久青草综合色| kizo精华| 国产成人免费无遮挡视频| 18禁动态无遮挡网站| 少妇人妻久久综合中文| 晚上一个人看的免费电影| 亚洲国产精品专区欧美| a级毛色黄片| 少妇丰满av| a 毛片基地| 亚洲第一区二区三区不卡| 欧美xxxx性猛交bbbb| 超碰97精品在线观看| 欧美xxⅹ黑人| 国产免费视频播放在线视频| 高清午夜精品一区二区三区| 人人妻人人添人人爽欧美一区卜 | 亚洲国产日韩一区二区| 国产精品av视频在线免费观看| 老女人水多毛片| 少妇丰满av| 国产精品一二三区在线看| 欧美日韩视频精品一区| 久久精品夜色国产| 大陆偷拍与自拍| 一级二级三级毛片免费看| 岛国毛片在线播放| 搡老乐熟女国产| 亚洲三级黄色毛片| 亚洲国产av新网站| 特大巨黑吊av在线直播| 十分钟在线观看高清视频www | 国产乱来视频区| 美女脱内裤让男人舔精品视频| 在线观看免费高清a一片| 亚洲欧美一区二区三区黑人 | 中文字幕人妻熟人妻熟丝袜美| 国产精品成人在线| 春色校园在线视频观看| 国产亚洲5aaaaa淫片| a 毛片基地| 成人黄色视频免费在线看| 久久久久国产精品人妻一区二区| 国产永久视频网站| 干丝袜人妻中文字幕| 婷婷色麻豆天堂久久| 最近2019中文字幕mv第一页| 少妇人妻一区二区三区视频| 国产精品国产三级国产专区5o| 日韩欧美 国产精品| 国产69精品久久久久777片| 男女啪啪激烈高潮av片| 国产黄色视频一区二区在线观看| 在线观看美女被高潮喷水网站| 91狼人影院| 欧美日韩国产mv在线观看视频 | 伦理电影免费视频| 亚洲最大成人中文| 综合色丁香网| 日韩成人伦理影院| 最后的刺客免费高清国语| 久久精品熟女亚洲av麻豆精品| 中文欧美无线码| 久久99热这里只频精品6学生| 99久久人妻综合| 黑丝袜美女国产一区| 亚洲一级一片aⅴ在线观看| 91久久精品国产一区二区成人| 成年免费大片在线观看| 亚洲国产精品999| 91午夜精品亚洲一区二区三区| 中国国产av一级| 成年免费大片在线观看| 十八禁网站网址无遮挡 | 99re6热这里在线精品视频| 国产伦精品一区二区三区视频9| 老司机影院成人| 欧美老熟妇乱子伦牲交| 久久精品久久久久久噜噜老黄| 韩国高清视频一区二区三区| 欧美精品人与动牲交sv欧美| 久热久热在线精品观看| 欧美精品一区二区免费开放| 国产人妻一区二区三区在| 亚洲天堂av无毛| 男人狂女人下面高潮的视频| 国产精品成人在线| 久久精品国产亚洲网站| 亚洲av中文av极速乱| 97热精品久久久久久| 大陆偷拍与自拍| av国产久精品久网站免费入址| 亚洲精品亚洲一区二区| 狂野欧美激情性bbbbbb| 2022亚洲国产成人精品| 久久久久久久久久久丰满| 99久久精品热视频| 纯流量卡能插随身wifi吗| 亚洲国产欧美在线一区| 大又大粗又爽又黄少妇毛片口| 一本—道久久a久久精品蜜桃钙片| 免费看不卡的av| 国国产精品蜜臀av免费| 久久久精品免费免费高清| 亚洲色图综合在线观看| 女人久久www免费人成看片| 一二三四中文在线观看免费高清| 久久亚洲国产成人精品v| 精品少妇久久久久久888优播| 熟女av电影| 免费观看的影片在线观看| 大话2 男鬼变身卡| 狠狠精品人妻久久久久久综合| 一区二区av电影网| 亚洲图色成人| 亚洲精品国产色婷婷电影| 国产色爽女视频免费观看| 我的女老师完整版在线观看| 国产精品99久久久久久久久| 深爱激情五月婷婷| 免费观看av网站的网址| 日韩欧美精品免费久久| 亚洲欧美精品自产自拍| 女的被弄到高潮叫床怎么办| 少妇丰满av| 另类亚洲欧美激情| 美女xxoo啪啪120秒动态图| 中国美白少妇内射xxxbb| 亚洲精品国产av成人精品| 精品熟女少妇av免费看| 蜜桃亚洲精品一区二区三区| 美女主播在线视频| 成年av动漫网址| 国产精品欧美亚洲77777| 人人妻人人看人人澡| 最近手机中文字幕大全| 黑丝袜美女国产一区| 好男人视频免费观看在线| 2022亚洲国产成人精品| 美女视频免费永久观看网站| 国产成人免费观看mmmm| 中文字幕久久专区| 国产一区亚洲一区在线观看| 日韩,欧美,国产一区二区三区| av卡一久久| 边亲边吃奶的免费视频| 久久久久国产精品人妻一区二区| 国产成人91sexporn| 少妇精品久久久久久久| 国产男人的电影天堂91| 欧美激情极品国产一区二区三区 | 在线观看免费日韩欧美大片 | 免费大片18禁| 日韩成人伦理影院| 日韩成人av中文字幕在线观看| 国产免费一级a男人的天堂| 少妇人妻精品综合一区二区| 一本色道久久久久久精品综合| 成人一区二区视频在线观看| a级毛片免费高清观看在线播放| 人人妻人人澡人人爽人人夜夜| 一区二区三区乱码不卡18| 国产白丝娇喘喷水9色精品| 99热6这里只有精品| 哪个播放器可以免费观看大片| 久久久精品免费免费高清| 成人18禁高潮啪啪吃奶动态图 | 菩萨蛮人人尽说江南好唐韦庄| 欧美日本视频| 欧美性感艳星| av网站免费在线观看视频| 秋霞在线观看毛片| 国产成人精品婷婷|