• <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ā)中易混淆的概念剖析
    一个人看的www免费观看视频| 草草在线视频免费看| 亚洲欧美精品专区久久| av播播在线观看一区| 亚洲成人精品中文字幕电影| 久久久午夜欧美精品| 亚洲一区二区三区欧美精品 | 一级毛片我不卡| 插阴视频在线观看视频| 一区二区三区乱码不卡18| 一区二区三区乱码不卡18| 高清欧美精品videossex| 性色avwww在线观看| 国产亚洲av嫩草精品影院| 高清欧美精品videossex| 国产精品熟女久久久久浪| 国产亚洲精品久久久com| 好男人在线观看高清免费视频| 日韩免费高清中文字幕av| 国产精品久久久久久av不卡| 久久99热这里只频精品6学生| 日韩成人av中文字幕在线观看| 亚洲aⅴ乱码一区二区在线播放| 免费观看无遮挡的男女| 69av精品久久久久久| 欧美成人a在线观看| 国模一区二区三区四区视频| 丰满少妇做爰视频| 毛片女人毛片| 简卡轻食公司| 韩国高清视频一区二区三区| 97精品久久久久久久久久精品| 亚洲成人一二三区av| 中文字幕制服av| 亚洲不卡免费看| 老女人水多毛片| 成年女人在线观看亚洲视频 | 国产av不卡久久| 国产成人福利小说| 国产视频首页在线观看| 嫩草影院精品99| 观看美女的网站| 91久久精品国产一区二区成人| 午夜老司机福利剧场| 中文乱码字字幕精品一区二区三区| 国产精品久久久久久精品电影小说 | 精品人妻视频免费看| 婷婷色综合大香蕉| 亚洲国产av新网站| 高清av免费在线| 欧美另类一区| 偷拍熟女少妇极品色| av国产精品久久久久影院| 六月丁香七月| 高清欧美精品videossex| 国产人妻一区二区三区在| 国产精品99久久99久久久不卡 | 欧美xxxx黑人xx丫x性爽| 国产亚洲av嫩草精品影院| 精品人妻熟女av久视频| 天堂俺去俺来也www色官网| 如何舔出高潮| 中文在线观看免费www的网站| 联通29元200g的流量卡| 久久久欧美国产精品| 特级一级黄色大片| 国产高清不卡午夜福利| 亚洲精品aⅴ在线观看| 日本与韩国留学比较| 亚洲欧洲国产日韩| 亚洲精品乱久久久久久| 97在线视频观看| 在线观看一区二区三区| 免费av毛片视频| 2018国产大陆天天弄谢| 久久人人爽av亚洲精品天堂 | 毛片一级片免费看久久久久| 欧美日韩国产mv在线观看视频 | 久久99热6这里只有精品| 如何舔出高潮| 久久久久久久亚洲中文字幕| 国产又色又爽无遮挡免| 伊人久久精品亚洲午夜| 精品国产一区二区三区久久久樱花 | 韩国av在线不卡| 国产欧美日韩一区二区三区在线 | www.av在线官网国产| 九九在线视频观看精品| 久久久午夜欧美精品| 精品久久久久久电影网| 国产男人的电影天堂91| 青青草视频在线视频观看| 六月丁香七月| 51国产日韩欧美| 欧美极品一区二区三区四区| 亚洲在线观看片| 熟女av电影| 99精国产麻豆久久婷婷| 交换朋友夫妻互换小说| 精品人妻视频免费看| 成人国产av品久久久| 免费黄频网站在线观看国产| 一级毛片aaaaaa免费看小| 国产伦在线观看视频一区| 成人亚洲欧美一区二区av| 男女边摸边吃奶| 成人国产麻豆网| 国产国拍精品亚洲av在线观看| 亚洲av一区综合| 亚洲国产最新在线播放| 韩国av在线不卡| 精品亚洲乱码少妇综合久久| 联通29元200g的流量卡| 国产精品一区二区三区四区免费观看| 日本免费在线观看一区| 亚洲av国产av综合av卡| 高清日韩中文字幕在线| 嫩草影院精品99| 婷婷色综合www| 香蕉精品网在线| 国产色爽女视频免费观看| 亚洲精华国产精华液的使用体验| 女人久久www免费人成看片| 国产亚洲5aaaaa淫片| 免费不卡的大黄色大毛片视频在线观看| 青青草视频在线视频观看| 国内精品宾馆在线| 蜜桃久久精品国产亚洲av| 免费看不卡的av| 在线天堂最新版资源| 亚洲成人中文字幕在线播放| 久久人人爽人人爽人人片va| 91精品伊人久久大香线蕉| 国产精品99久久99久久久不卡 | 99热全是精品| 观看免费一级毛片| 在线精品无人区一区二区三 | 亚洲精品亚洲一区二区| 免费av毛片视频| 综合色丁香网| 国产精品99久久99久久久不卡 | 国产精品国产三级国产专区5o| 观看免费一级毛片| 久久鲁丝午夜福利片| 欧美精品一区二区大全| 亚洲欧美成人综合另类久久久| 国产高清国产精品国产三级 | 乱系列少妇在线播放| 97在线人人人人妻| 久久精品国产亚洲av涩爱| av网站免费在线观看视频| 菩萨蛮人人尽说江南好唐韦庄| 欧美性感艳星| 亚洲欧美中文字幕日韩二区| 久久久久久久久久久丰满| 色5月婷婷丁香| 男的添女的下面高潮视频| 亚洲欧美成人综合另类久久久| 夫妻午夜视频| 五月玫瑰六月丁香| av国产久精品久网站免费入址| 日韩一区二区三区影片| 国产探花在线观看一区二区| 欧美日韩综合久久久久久| av天堂中文字幕网| 黄片无遮挡物在线观看| 欧美国产精品一级二级三级 | 亚洲无线观看免费| 久久久久国产网址| 亚洲精品视频女| 精品国产露脸久久av麻豆| 哪个播放器可以免费观看大片| 丰满乱子伦码专区| 国产精品爽爽va在线观看网站| 欧美少妇被猛烈插入视频| 久久久久九九精品影院| 大香蕉97超碰在线| 成人鲁丝片一二三区免费| 极品教师在线视频| 亚洲精品一二三| 国产女主播在线喷水免费视频网站| 91久久精品电影网| 久久久久久久久大av| 观看美女的网站| 少妇猛男粗大的猛烈进出视频 | 天天一区二区日本电影三级| 国产男女超爽视频在线观看| 少妇人妻一区二区三区视频| 国产伦精品一区二区三区视频9| 中文精品一卡2卡3卡4更新| 91久久精品国产一区二区成人| 亚洲精品影视一区二区三区av| 亚洲国产最新在线播放| 热99国产精品久久久久久7| 久久久色成人| 久久久久久久久久久丰满| 最近最新中文字幕大全电影3| 极品少妇高潮喷水抽搐| 中文字幕人妻熟人妻熟丝袜美| 极品少妇高潮喷水抽搐| 免费大片黄手机在线观看| 最近中文字幕2019免费版| 永久免费av网站大全| 亚洲精品,欧美精品| 日本猛色少妇xxxxx猛交久久| 一级二级三级毛片免费看| 女人十人毛片免费观看3o分钟| 女人十人毛片免费观看3o分钟| 一区二区三区免费毛片| 国产av国产精品国产| 毛片一级片免费看久久久久| 亚洲精品影视一区二区三区av| 国产精品久久久久久久久免| 亚洲人成网站在线播| 国产v大片淫在线免费观看| 国产91av在线免费观看| 人人妻人人爽人人添夜夜欢视频 | 看黄色毛片网站| 成人午夜精彩视频在线观看| 日韩欧美 国产精品| 国产精品福利在线免费观看| 99热国产这里只有精品6| 中国三级夫妇交换| 女人十人毛片免费观看3o分钟| 嫩草影院入口| 国产精品熟女久久久久浪| 韩国高清视频一区二区三区| 高清在线视频一区二区三区| 男插女下体视频免费在线播放| 国产黄色视频一区二区在线观看| 极品教师在线视频| 男女那种视频在线观看| 人人妻人人看人人澡| 久久久精品免费免费高清| 另类亚洲欧美激情| 亚洲av二区三区四区| 亚洲av日韩在线播放| 欧美日韩在线观看h| 成人国产av品久久久| 一区二区三区免费毛片| 一本久久精品| 又爽又黄a免费视频| 另类亚洲欧美激情| 色综合色国产| 在线观看一区二区三区| 2018国产大陆天天弄谢| 精品少妇久久久久久888优播| 欧美日韩国产mv在线观看视频 | 午夜福利视频精品| 国产欧美日韩一区二区三区在线 | 天天躁日日操中文字幕| 亚洲四区av| 观看免费一级毛片| 国产成人一区二区在线| 欧美最新免费一区二区三区| 1000部很黄的大片| 国产成人精品婷婷| 亚洲,一卡二卡三卡| 国产在视频线精品| 免费黄色在线免费观看| 一级毛片aaaaaa免费看小| 人妻一区二区av| 看十八女毛片水多多多| 日韩人妻高清精品专区| a级毛片免费高清观看在线播放| 大片免费播放器 马上看| 亚洲人与动物交配视频| 久久精品国产亚洲网站| 在线a可以看的网站| 亚洲精品一二三| 亚洲精品乱码久久久v下载方式| 欧美bdsm另类| 日韩成人伦理影院| 老师上课跳d突然被开到最大视频| 麻豆精品久久久久久蜜桃| 亚洲天堂国产精品一区在线| 成人黄色视频免费在线看| 久久国产乱子免费精品| 又黄又爽又刺激的免费视频.| 2021少妇久久久久久久久久久| 免费观看的影片在线观看| 精品少妇久久久久久888优播| 国产一区二区亚洲精品在线观看| 免费观看无遮挡的男女| 亚洲av免费在线观看| 激情五月婷婷亚洲| 身体一侧抽搐| 综合色av麻豆| 国产精品一及| 国产精品99久久99久久久不卡 | 高清视频免费观看一区二区| 波多野结衣巨乳人妻| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品第二区| 亚洲精品色激情综合| 欧美日韩在线观看h| 亚洲经典国产精华液单| 国产老妇女一区| 亚洲天堂av无毛| 男女边吃奶边做爰视频| 亚洲av不卡在线观看| 亚洲,欧美,日韩| 麻豆久久精品国产亚洲av| 男女国产视频网站| 一区二区三区免费毛片| 中文字幕制服av| 国产午夜精品一二区理论片| 国产黄频视频在线观看| 男女边吃奶边做爰视频| 亚洲精品日本国产第一区| 亚洲欧美成人精品一区二区| 三级国产精品欧美在线观看| 夫妻性生交免费视频一级片| 日韩人妻高清精品专区| 国国产精品蜜臀av免费| 国产免费视频播放在线视频| 超碰97精品在线观看| 久久女婷五月综合色啪小说 | 日韩一区二区视频免费看| h日本视频在线播放| 久久久久久久久久久免费av| 91精品国产九色| 99热这里只有精品一区| 九草在线视频观看| 亚洲精品成人久久久久久| 有码 亚洲区| 国产精品熟女久久久久浪| 热re99久久精品国产66热6| 亚洲人成网站在线观看播放| 女人久久www免费人成看片| 三级国产精品欧美在线观看| 免费av毛片视频| 亚洲最大成人av| 特大巨黑吊av在线直播| 波多野结衣巨乳人妻| 久久精品国产亚洲av涩爱| 欧美丝袜亚洲另类| 欧美丝袜亚洲另类| 欧美xxxx性猛交bbbb| 一级毛片我不卡| 天天躁夜夜躁狠狠久久av| 精品久久国产蜜桃| 国产成年人精品一区二区| 好男人在线观看高清免费视频| av线在线观看网站| 欧美日本视频| 欧美三级亚洲精品| 少妇的逼好多水| 成人漫画全彩无遮挡| 男插女下体视频免费在线播放| 内地一区二区视频在线| av福利片在线观看| 亚洲真实伦在线观看| 大码成人一级视频| 亚洲最大成人av| 婷婷色综合www| 精品久久久精品久久久| 亚洲欧美日韩另类电影网站 | 日韩电影二区| 色5月婷婷丁香| 国产精品久久久久久久电影| 少妇人妻久久综合中文| 国产精品女同一区二区软件| 一级毛片我不卡| 国产精品蜜桃在线观看| av在线app专区| 免费观看av网站的网址| 国产精品无大码| 亚洲国产欧美人成| 欧美激情在线99| 狂野欧美激情性bbbbbb| 午夜精品国产一区二区电影 | 国精品久久久久久国模美| 亚洲精品国产av成人精品| 国产精品99久久99久久久不卡 | 国产69精品久久久久777片| 欧美3d第一页| 日本一本二区三区精品| 中文字幕久久专区| 少妇的逼水好多| h日本视频在线播放| 日本黄大片高清| 国产永久视频网站| 两个人的视频大全免费| 中文欧美无线码| 九色成人免费人妻av| 六月丁香七月| 久久久久久久精品精品| 性色avwww在线观看| 青春草国产在线视频| 男女无遮挡免费网站观看| 亚洲电影在线观看av| av免费观看日本| 亚洲伊人久久精品综合| 亚洲婷婷狠狠爱综合网| 美女高潮的动态| 熟妇人妻不卡中文字幕| 欧美人与善性xxx| 日产精品乱码卡一卡2卡三| 一区二区三区四区激情视频| 香蕉精品网在线| 国产亚洲av片在线观看秒播厂| 内射极品少妇av片p| 国产黄色免费在线视频| 欧美日韩视频高清一区二区三区二| 少妇 在线观看| 精品久久久噜噜| 亚洲av二区三区四区| 午夜福利高清视频| 国产亚洲午夜精品一区二区久久 | 久久久久久久久大av| 国产免费一区二区三区四区乱码| av女优亚洲男人天堂| 在线精品无人区一区二区三 | 亚洲色图综合在线观看| av线在线观看网站| 久久亚洲国产成人精品v| 精品人妻偷拍中文字幕| 一本色道久久久久久精品综合| 深爱激情五月婷婷| 国产大屁股一区二区在线视频| 免费av不卡在线播放| 亚洲成人中文字幕在线播放| 日本猛色少妇xxxxx猛交久久| 秋霞在线观看毛片| 美女被艹到高潮喷水动态| 全区人妻精品视频| 日韩人妻高清精品专区| 亚洲一级一片aⅴ在线观看| 综合色丁香网| 午夜精品国产一区二区电影 | 97在线人人人人妻| 日韩成人伦理影院| 久久热精品热| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 自拍偷自拍亚洲精品老妇| 亚洲欧洲日产国产| 一区二区三区免费毛片| 亚洲,欧美,日韩| 精品久久久久久久久亚洲| 青春草亚洲视频在线观看| 久久精品综合一区二区三区| 亚洲av电影在线观看一区二区三区 | 91精品一卡2卡3卡4卡| 日本免费在线观看一区| 日日撸夜夜添| 18禁裸乳无遮挡免费网站照片| 欧美日韩亚洲高清精品| 国产一区二区三区综合在线观看 | 亚洲欧美一区二区三区国产| 国产日韩欧美在线精品| 99久久九九国产精品国产免费| 老司机影院成人| 水蜜桃什么品种好| 真实男女啪啪啪动态图| 夫妻性生交免费视频一级片| 少妇的逼好多水| 好男人在线观看高清免费视频| 亚洲精品乱码久久久v下载方式| 白带黄色成豆腐渣| 少妇人妻一区二区三区视频| 97热精品久久久久久| 成年免费大片在线观看| 成人午夜精彩视频在线观看| 成人毛片a级毛片在线播放| 国产精品久久久久久精品电影| 好男人在线观看高清免费视频| 亚洲av中文字字幕乱码综合| 亚洲精品视频女| 少妇裸体淫交视频免费看高清| 久久人人爽人人片av| 18+在线观看网站| av又黄又爽大尺度在线免费看| 2021少妇久久久久久久久久久| 久久久久久久亚洲中文字幕| 亚洲成人av在线免费| 亚洲精品国产成人久久av| 国产片特级美女逼逼视频| 青春草视频在线免费观看| 亚洲欧美一区二区三区黑人 | 亚洲精品一二三| 午夜福利在线在线| 少妇猛男粗大的猛烈进出视频 | 国产毛片a区久久久久| 少妇人妻一区二区三区视频| 99九九线精品视频在线观看视频| 新久久久久国产一级毛片| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久精品电影| 一区二区三区四区激情视频| 色视频www国产| 国产高清不卡午夜福利| 久久久久精品久久久久真实原创| 免费观看性生交大片5| 国产精品蜜桃在线观看| av在线亚洲专区| 国产午夜精品久久久久久一区二区三区| 久久影院123| 91精品国产九色| 久久亚洲国产成人精品v| 777米奇影视久久| 久久精品夜色国产| 中国美白少妇内射xxxbb| 自拍偷自拍亚洲精品老妇| 天美传媒精品一区二区| 男人爽女人下面视频在线观看| 日韩欧美精品免费久久| 日本三级黄在线观看| 女人久久www免费人成看片| 国产精品久久久久久精品电影| 91在线精品国自产拍蜜月| 极品少妇高潮喷水抽搐| 亚洲av成人精品一二三区| 在线观看免费高清a一片| 国产欧美日韩一区二区三区在线 | 国产精品一及| 蜜桃久久精品国产亚洲av| 人人妻人人澡人人爽人人夜夜| 视频区图区小说| 插阴视频在线观看视频| 日本三级黄在线观看| 99精国产麻豆久久婷婷| 亚洲自偷自拍三级| 国产片特级美女逼逼视频| 国产日韩欧美在线精品| 久久99热这里只有精品18| 色播亚洲综合网| 国产在线一区二区三区精| 91久久精品电影网| 欧美xxxx黑人xx丫x性爽| 日日啪夜夜爽| 欧美xxⅹ黑人| 777米奇影视久久| 97精品久久久久久久久久精品| 午夜视频国产福利| 亚洲色图av天堂| 精品人妻视频免费看| 日韩强制内射视频| 搡老乐熟女国产| 欧美xxxx性猛交bbbb| 国产综合懂色| 成人高潮视频无遮挡免费网站| 女人久久www免费人成看片| 日本与韩国留学比较| 久久久成人免费电影| 熟妇人妻不卡中文字幕| 国产伦精品一区二区三区四那| 大陆偷拍与自拍| 国产黄频视频在线观看| 一本色道久久久久久精品综合| 嫩草影院新地址| 国产精品国产三级国产专区5o| 97在线人人人人妻| 日韩国内少妇激情av| 久久精品人妻少妇| 97超视频在线观看视频| 毛片女人毛片| kizo精华| 欧美日韩在线观看h| 观看美女的网站| 精品人妻一区二区三区麻豆| 中文字幕av成人在线电影| 一级二级三级毛片免费看| 亚洲国产精品成人久久小说| 免费观看av网站的网址| 久久久久久国产a免费观看| 热re99久久精品国产66热6| 大香蕉久久网| 久久久久久久久大av| 美女内射精品一级片tv| 国产视频内射| 亚洲人成网站高清观看| 熟女电影av网| 男男h啪啪无遮挡| 国产日韩欧美亚洲二区| 国产黄频视频在线观看| 能在线免费看毛片的网站| 精品酒店卫生间| 中国三级夫妇交换| 国产视频内射| 99精国产麻豆久久婷婷| 日韩强制内射视频| 国产欧美日韩一区二区三区在线 | 国产精品蜜桃在线观看| 国产黄频视频在线观看| 性色avwww在线观看| 男女无遮挡免费网站观看| 男人舔奶头视频| 卡戴珊不雅视频在线播放| 日韩国内少妇激情av| 亚洲综合色惰| 黄色视频在线播放观看不卡| 人人妻人人澡人人爽人人夜夜| 91久久精品国产一区二区成人| 亚洲av成人精品一区久久| 2022亚洲国产成人精品| 午夜免费鲁丝| 99re6热这里在线精品视频| 我的女老师完整版在线观看| 秋霞在线观看毛片| 日本熟妇午夜| 国产精品三级大全| 极品教师在线视频| 交换朋友夫妻互换小说| 亚洲精品日本国产第一区| kizo精华| av国产免费在线观看| 成人鲁丝片一二三区免费| 一级毛片黄色毛片免费观看视频| 日韩亚洲欧美综合| 国产白丝娇喘喷水9色精品| 欧美zozozo另类| 精品视频人人做人人爽| 亚洲国产欧美在线一区| 久久精品国产亚洲网站|