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

    多點(diǎn)種子預(yù)劃分的二階段社區(qū)發(fā)現(xiàn)算法

    2021-10-07 06:24:32佟帥陳德運(yùn)楊海陸

    佟帥 陳德運(yùn) 楊海陸

    摘 要:社區(qū)發(fā)現(xiàn)是在線社交網(wǎng)絡(luò)研究領(lǐng)域中的重要內(nèi)容,基于種子擴(kuò)張的社區(qū)發(fā)現(xiàn)算法具有時(shí)間復(fù)雜度低、識(shí)別精度高以及不受社區(qū)形態(tài)限制等特點(diǎn),近年來(lái)在網(wǎng)絡(luò)局部社區(qū)發(fā)現(xiàn)任務(wù)中得到了廣泛的應(yīng)用。然而,該方法在種子選取時(shí)沒(méi)有考慮種子之間的關(guān)聯(lián)性,因此識(shí)別出的社區(qū)結(jié)構(gòu)個(gè)數(shù)較多、結(jié)構(gòu)松散。針對(duì)這一問(wèn)題,提出一種基于多點(diǎn)種子預(yù)劃分的二階段社區(qū)發(fā)現(xiàn)算法。首先識(shí)別網(wǎng)絡(luò)中的高影響力節(jié)點(diǎn),利用K-means算法將高影響力節(jié)點(diǎn)加以聚合,得到高影響力社區(qū)簇。然后提出一種吸引力度量函數(shù),選擇性的將網(wǎng)絡(luò)中的剩余節(jié)點(diǎn)合并到社區(qū)簇以完成社區(qū)識(shí)別任務(wù)。實(shí)驗(yàn)結(jié)果表明,二階段社區(qū)發(fā)現(xiàn)方法能夠發(fā)現(xiàn)尺寸較大,個(gè)數(shù)較少的社區(qū)結(jié)構(gòu),進(jìn)而在中觀層面捕捉群組之間的關(guān)聯(lián)性。

    關(guān)鍵詞:復(fù)雜網(wǎng)絡(luò);局部社區(qū)發(fā)現(xiàn);種子擴(kuò)張;節(jié)點(diǎn)影響力;K-means算法

    DOI:10.15938/j.jhust.2021.04.011

    中圖分類(lèi)號(hào):TP391.4

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1007-2683(2021)04-0078-09

    Abstract:Community detection is an important content in the field of online social networks research. The community detection algorithm based on seed expansion has the characteristics of low time complexity, high recognition accuracy, and is not restricted by the shape of the community. In recent years, it has been widely used in local community discovery of complex networks. However, this method does not consider the correlation between seeds when selecting seeds, so the number of identified community structures is large and the structure is loose. Aiming at this problem, a two-stage community discovery algorithm for multi-point seed prepartition was proposed. First, high-impact nodes in the network are identified, and high-impact nodes are aggregated using the k-means algorithm to obtain high-impact community clusters. Next, an attractiveness measurement function is proposed to selectively merge the remaining nodes in the network into the community cluster to complete the community identification task. The experimental results show that the two-stage community discovery method can find community structures with larger sizes and fewer numbers, and then capture the association between groups at the meso level.

    Keywords:complex network; local community detection; seed expansion; node influence; K-means algorithm

    0 引 言

    現(xiàn)實(shí)世界的復(fù)雜網(wǎng)絡(luò)通??梢猿橄鬄閳D模型,圖中節(jié)點(diǎn)代表現(xiàn)實(shí)世界的實(shí)體,兩個(gè)節(jié)點(diǎn)之間的鏈接代表實(shí)體之間的關(guān)系。社區(qū)[1]是復(fù)雜網(wǎng)絡(luò)中的稠密子圖,保證了社區(qū)內(nèi)部節(jié)點(diǎn)之間的鏈接較為緊密,社區(qū)之間節(jié)點(diǎn)的鏈接較為稀疏[2-5]。探索社區(qū)結(jié)構(gòu)有助于人們理解復(fù)雜網(wǎng)絡(luò)的自組織以及群聚特性,是復(fù)雜網(wǎng)絡(luò)中觀層次最重要的屬性之一。

    從社區(qū)的層次化角度來(lái)看,已有的社區(qū)識(shí)別算法可分為全局優(yōu)化算法和局部?jī)?yōu)化算法兩種。前者在宏觀角度尋找社區(qū)結(jié)構(gòu),描述了宏觀層面網(wǎng)絡(luò)節(jié)點(diǎn)的自組織特性,適合以數(shù)理統(tǒng)計(jì)為目的的網(wǎng)絡(luò)特征分析[1-2]。后者從微觀層面刻畫(huà)節(jié)點(diǎn)的局部?jī)A向性,更適合發(fā)現(xiàn)社區(qū)的演化規(guī)律及形成過(guò)程。基于種子節(jié)點(diǎn)擴(kuò)張的社區(qū)發(fā)現(xiàn)算法[6],就是局部?jī)?yōu)化算法的典型代表。

    2009年,Lancichinetti等[7]提出重疊社區(qū)發(fā)現(xiàn)算法LFM,根據(jù)定義的自適應(yīng)函數(shù)Fitness進(jìn)行基于種子擴(kuò)張的社區(qū)識(shí)別。Coscia等[8]提出的DEMON算法以整個(gè)網(wǎng)絡(luò)節(jié)點(diǎn)為起始種子,利用標(biāo)簽傳播算法識(shí)別網(wǎng)絡(luò)中的局部社區(qū)結(jié)構(gòu)。Baumes等[9]提出一種鏈接聚合方法。先按照遞增或遞減原則對(duì)節(jié)點(diǎn)進(jìn)行排序,作為初始社區(qū),如果節(jié)點(diǎn)的添加不能提高任何社區(qū)的密度,則該節(jié)點(diǎn)將作為新的種子節(jié)點(diǎn),并生成新社區(qū)。李婕等[10]采用基于派系過(guò)濾的算法選擇種子節(jié)點(diǎn)進(jìn)行社區(qū)識(shí)別。Su等[11]基于隨機(jī)游走算法,使用緊密連接子圖作為社區(qū)識(shí)別的初始種子社區(qū)。

    提高基于種子擴(kuò)張的社區(qū)識(shí)別性能的關(guān)鍵在于種子節(jié)點(diǎn)的選取,近年來(lái)有研究者提出借助影響力分析方法,增加初始種子節(jié)點(diǎn)位于社區(qū)內(nèi)核的概率。Clauset等[12]依次探索網(wǎng)絡(luò)節(jié)點(diǎn),推斷給定節(jié)點(diǎn)所在社區(qū)。Luo等[13]將度的概念從單節(jié)點(diǎn)擴(kuò)展到子圖,在此基礎(chǔ)上給出了網(wǎng)絡(luò)模塊化的定義,以此進(jìn)行社區(qū)識(shí)別。Chen等[14]提出一種新的局部社區(qū)測(cè)度,先提取所有可能的候選社區(qū),然后對(duì)社區(qū)層次進(jìn)行優(yōu)化。吳英俊等[15]提出LS算法,通過(guò)分析社區(qū)與節(jié)點(diǎn)之間的鏈接相似性來(lái)尋找局部社區(qū)結(jié)構(gòu)。Fanrong等[16]提出一種基于最大團(tuán)擴(kuò)展的局部社區(qū)檢測(cè)算法LCDMC,首先找到包含源節(jié)點(diǎn)的最大派系集合,然后利用貪婪優(yōu)化方法擴(kuò)展社區(qū)。Yao等[17]分析了高影響節(jié)點(diǎn)在社區(qū)檢測(cè)中的作用。齊金山等[18]提出了一種結(jié)合Jaccard系數(shù)的節(jié)點(diǎn)影響力計(jì)算公式,提高了算法對(duì)星形社區(qū)的匹配性。

    上述社區(qū)識(shí)別方法雖然在特定領(lǐng)域具有一定的性能優(yōu)勢(shì),但普遍存在以下兩方面問(wèn)題。首先,上述方法在種子選取時(shí)沒(méi)有考慮種子之間的關(guān)聯(lián)性,導(dǎo)致識(shí)別出社區(qū)結(jié)構(gòu)個(gè)數(shù)較多、穩(wěn)定性較差,不利于在中觀層面上挖掘群組之間的關(guān)聯(lián)性。其次,由于社會(huì)網(wǎng)絡(luò)結(jié)構(gòu)的多樣性,上述方法在影響力識(shí)別時(shí)忽略了網(wǎng)絡(luò)的結(jié)構(gòu)特性,導(dǎo)致網(wǎng)絡(luò)中的高影響力節(jié)點(diǎn)呈現(xiàn)出較低的影響力評(píng)分的假象,這使得識(shí)別出的社區(qū)結(jié)構(gòu)較為松散,內(nèi)聚性較差。

    為了解決上述問(wèn)題,提出一種多點(diǎn)種子預(yù)劃分的二階段社區(qū)發(fā)現(xiàn)算法(two-stage community detection algorithm,TSCDA),其創(chuàng)新之處主要體現(xiàn)在以下三方面。首先,在節(jié)點(diǎn)的影響力計(jì)算時(shí)融入邊介數(shù)屬性,增強(qiáng)節(jié)點(diǎn)影響力的結(jié)構(gòu)因素;其次,通過(guò)預(yù)劃分高影響力節(jié)點(diǎn)生成骨干網(wǎng)絡(luò),增強(qiáng)種子的結(jié)構(gòu)密度;最后,以稠密子圖代替節(jié)點(diǎn)作為種子進(jìn)行擴(kuò)張,增加了社區(qū)結(jié)構(gòu)的穩(wěn)定性。仿真結(jié)果表明,TSCDA能夠挖掘出數(shù)量較少、規(guī)模較大的社區(qū)結(jié)構(gòu),并在模塊度以及F-Score等指標(biāo)上表現(xiàn)出一定的性能優(yōu)勢(shì)。

    1 基于節(jié)點(diǎn)影響力識(shí)別的種子節(jié)點(diǎn)選取方法

    用無(wú)向圖G=(V,E)表示社交網(wǎng)絡(luò),其中V代表網(wǎng)絡(luò)中個(gè)數(shù)為n的節(jié)點(diǎn)集合,E代表網(wǎng)絡(luò)中個(gè)數(shù)為m的鏈接關(guān)系集合。圖1給出了一種具有強(qiáng)聚合特性的社交網(wǎng)絡(luò)結(jié)構(gòu)(社區(qū)網(wǎng)絡(luò))。

    社會(huì)網(wǎng)絡(luò)中,節(jié)點(diǎn)u的鄰居集合定義為:

    這里u和v代表網(wǎng)絡(luò)節(jié)點(diǎn),V代表網(wǎng)絡(luò)中的節(jié)點(diǎn)集,如果u和v之間存在邊(u,v),則(u,v)屬于邊集E。節(jié)點(diǎn)u的節(jié)點(diǎn)度D(u)定義為與u直接相連的邊的個(gè)數(shù)(或鄰居元素的個(gè)數(shù)),滿足D(u)=|N(u)|。

    為了在影響力計(jì)算中引入結(jié)構(gòu)的相似屬性,提出從節(jié)點(diǎn)相似性以及鏈接相似性雙重角度構(gòu)造影響力評(píng)價(jià)函數(shù)。

    Jaccard系數(shù)Juv用于比較社會(huì)網(wǎng)絡(luò)中節(jié)點(diǎn)之間的結(jié)構(gòu)相似性。其定義為:

    N(u)和N(v)分別代表節(jié)點(diǎn)u和節(jié)點(diǎn)v的鄰居節(jié)點(diǎn)。JuvJuv取值范圍為[0,1],其值越大,節(jié)點(diǎn)u和節(jié)點(diǎn)v的共有鄰居節(jié)點(diǎn)就越多,鄰域結(jié)構(gòu)也就相對(duì)越稠密。

    邊介數(shù)用來(lái)度量社會(huì)網(wǎng)絡(luò)中所有節(jié)點(diǎn)間的最短路徑中經(jīng)過(guò)該邊的路徑的數(shù)目占最短路徑總數(shù)的比例。邊介數(shù)值越大,該邊越有可能成為連接社區(qū)內(nèi)部節(jié)點(diǎn)的重要途徑。其定義為:

    其中g(shù)uvst表示網(wǎng)絡(luò)中任意節(jié)點(diǎn)s與任意節(jié)點(diǎn)t之間的最短路徑同時(shí)經(jīng)過(guò)節(jié)點(diǎn)u和v(即邊euv)的條數(shù)。根據(jù)式(3),可以計(jì)算出圖1各邊的邊介數(shù)指數(shù),例如節(jié)點(diǎn)8和節(jié)點(diǎn)10,即(e8,10)之間邊強(qiáng)度為EBET(e8,9)=11;節(jié)點(diǎn)9和節(jié)點(diǎn)10,即(e9,10)之間邊強(qiáng)度為EBET(e9,10)=11。邊介數(shù)反映了社區(qū)內(nèi)部鏈接的稠密程度,是增加社區(qū)穩(wěn)定性的關(guān)鍵參數(shù)。

    受牛頓萬(wàn)有引力定律的啟發(fā),提出一種基于Jaccard系數(shù)以及邊介數(shù)的影響力計(jì)算方法。社會(huì)網(wǎng)絡(luò)節(jié)點(diǎn)u對(duì)節(jié)點(diǎn)v的影響力評(píng)分定義為:

    這里D(u)和D(v)分別表示節(jié)點(diǎn)u和節(jié)點(diǎn)v的度。由于節(jié)點(diǎn)度在局部環(huán)境反映了節(jié)點(diǎn)鏈接能力,因此在本文中被用來(lái)衡量社會(huì)網(wǎng)絡(luò)節(jié)點(diǎn)的質(zhì)量。節(jié)點(diǎn)之間的距離d(u,v)用節(jié)點(diǎn)間的相異程度加以衡量,滿足d(u,v)=1-Juv。其中Juv代表節(jié)點(diǎn)u和節(jié)點(diǎn)v的Jaccard相似度。G是影響力常量,在實(shí)驗(yàn)中通常取值為1。

    節(jié)點(diǎn)在局部社區(qū)的影響力評(píng)分Iscore(u)定義為該節(jié)點(diǎn)對(duì)其所有鄰居的節(jié)點(diǎn)影響力評(píng)分之和。進(jìn)而有:

    式(5)中,v∈N(u)表明節(jié)點(diǎn)v是節(jié)點(diǎn)u的鄰居節(jié)點(diǎn)。Iscore(u)的值越大,說(shuō)明節(jié)點(diǎn)u在局部環(huán)境的影響強(qiáng)度越大,表明越有可能成為種子節(jié)點(diǎn)。根據(jù)式(5),可以計(jì)算出圖1各節(jié)點(diǎn)的影響力分?jǐn)?shù)。例如節(jié)點(diǎn)10的影響力分?jǐn)?shù):

    因此節(jié)點(diǎn)10在網(wǎng)絡(luò)中影響力分?jǐn)?shù)Iscore(u10)=275。

    下面給出基于節(jié)點(diǎn)影響力識(shí)別的種子節(jié)點(diǎn)選取過(guò)程。算法通過(guò)式(5)計(jì)算社會(huì)網(wǎng)絡(luò)每個(gè)節(jié)點(diǎn)的影響力評(píng)分。其偽代碼如下:

    算法1 具有高影響力的初始種子選取算法

    輸入:社會(huì)網(wǎng)絡(luò)G=(V,E)

    輸出:高質(zhì)量節(jié)點(diǎn)集合S

    1)初始化集合S為空集;

    2)遍歷社會(huì)網(wǎng)絡(luò)的所有節(jié)點(diǎn),根據(jù)公式(5)計(jì)算社會(huì)網(wǎng)絡(luò)所有節(jié)點(diǎn)的節(jié)點(diǎn)影響力評(píng)分;

    3)對(duì)于任意節(jié)點(diǎn),如果它的影響力評(píng)分大于其任一鄰居的影響力評(píng)分,則為高影響力種子節(jié)點(diǎn),將其存儲(chǔ)到集合S;

    4)輸出集合S。

    根據(jù)算法1,計(jì)算圖1中節(jié)點(diǎn)1到節(jié)點(diǎn)10的影響力評(píng)分,具體如表1所示。在本文中,如果一個(gè)節(jié)點(diǎn)的影響力分?jǐn)?shù)大于其任一鄰居的影響力評(píng)分,則為高影響種子節(jié)點(diǎn)。例如,節(jié)點(diǎn)8的鄰居節(jié)點(diǎn)為6,7,9,10。通過(guò)表1可知,節(jié)點(diǎn)8的影響力評(píng)分大于節(jié)點(diǎn)6,7,9,10的影響力評(píng)分,因此節(jié)點(diǎn)8為高影響力節(jié)點(diǎn)。

    2 多點(diǎn)種子預(yù)劃分的二階段社區(qū)發(fā)現(xiàn)算法

    本節(jié)給出多點(diǎn)種子預(yù)劃分的二階段社區(qū)發(fā)現(xiàn)算法。首先,根據(jù)算法1選取網(wǎng)絡(luò)中的高影響力種子節(jié)點(diǎn),對(duì)這些節(jié)點(diǎn)利用K-means聚類(lèi)算法進(jìn)行初始劃分,識(shí)別出由高質(zhì)量節(jié)點(diǎn)組成的種子社區(qū)即骨干社區(qū);然后提出一種基于吸引力度量的社區(qū)識(shí)別算法將社會(huì)網(wǎng)絡(luò)中其余節(jié)點(diǎn)按特定規(guī)則有選擇性的加入骨干社區(qū)中,從而完成社區(qū)劃分。該方法的優(yōu)勢(shì)在于充分考慮了種子節(jié)點(diǎn)之間稠密性以及關(guān)聯(lián)性,有助于提高社區(qū)的穩(wěn)定程度。

    2.1 基于K-means聚類(lèi)的骨干社區(qū)識(shí)別

    本文采用K-means算法作為骨干社區(qū)的初始化算法,其距離度量公式在式(7)給出。

    式中:minDis(A,B)為在社會(huì)網(wǎng)絡(luò)中節(jié)點(diǎn)A,B間最短路徑長(zhǎng)度;comNeighbor(A,B)為節(jié)點(diǎn)A,B間共同鄰居數(shù)量;C為迭代中心,初始值為0;Dis(A,B)為A、B兩點(diǎn)的距離,該值與A、B兩點(diǎn)的最短路徑長(zhǎng)度呈正相關(guān),與A、B兩點(diǎn)的共同鄰居數(shù)量呈負(fù)相關(guān)。

    K-means算法在迭代時(shí)依靠簇內(nèi)所有節(jié)點(diǎn)的平均坐標(biāo)設(shè)定新中心。然而社區(qū)不存在坐標(biāo),因此需要重新定義中心選取方式。本文選取到簇內(nèi)其他節(jié)點(diǎn)的距離之和最小的節(jié)點(diǎn)作為新的中心節(jié)點(diǎn),具體如式(8)所示。

    式中:N(i)表示第i個(gè)社區(qū)內(nèi)的節(jié)點(diǎn)集合。K-means算法的最終輸出1棵骨干社區(qū)層次化樹(shù),采用模塊度函數(shù)在層次化樹(shù)上進(jìn)行分割,獲取質(zhì)量最高的社區(qū)劃分。模塊度函數(shù)Q是Newman[19]在2004年提出的概念,用以評(píng)價(jià)社區(qū)劃分的質(zhì)量,其定義為:

    式中:eii代表社區(qū)i中邊的期望;ai代表鏈接到社區(qū)i中的邊數(shù)的期望。模塊度的形式化定義是社區(qū)內(nèi)的邊數(shù)減去隨機(jī)產(chǎn)生的邊的期望,因此模塊度越高,社區(qū)內(nèi)外邊比例就越大,社區(qū)劃分結(jié)果就越好。

    我們發(fā)現(xiàn),算法在層次聚類(lèi)時(shí)模塊度函數(shù)Q的取值并不是單調(diào)遞增或單調(diào)遞減的,因此每當(dāng)有社區(qū)進(jìn)行合并操作,就需要重新計(jì)算Q值,能對(duì)產(chǎn)生最大Q值的兩個(gè)社區(qū)合并。為簡(jiǎn)化計(jì)算,每次僅計(jì)

    算Q的變化部分,即向Q值增大最多或者減少最小的方向進(jìn)行社區(qū)合并。Q函數(shù)的變化量計(jì)算公式如下所示。

    對(duì)能產(chǎn)生最大Q值的兩社區(qū)進(jìn)行合并,直到所有節(jié)點(diǎn)合并至同一社區(qū)。算法的偽代碼如下:

    算法2 利用K-means聚類(lèi)算法實(shí)現(xiàn)種子節(jié)點(diǎn)初始劃分

    輸入:3個(gè)種子節(jié)點(diǎn)集合,社會(huì)網(wǎng)絡(luò)G

    輸出:骨干社區(qū)集合BC={BC1,BC2,…BCn}

    1)選取這3個(gè)種子節(jié)點(diǎn)進(jìn)行初始聚類(lèi);

    2)根據(jù)距離度量式(7)計(jì)算任意兩節(jié)點(diǎn)之間的距離;

    3)比較任意兩點(diǎn)間的距離,將節(jié)點(diǎn)劃分至距離最近的節(jié)點(diǎn)簇內(nèi);

    4)對(duì)于每個(gè)簇內(nèi),分別計(jì)算簇內(nèi)每個(gè)節(jié)點(diǎn)與簇內(nèi)的其他節(jié)點(diǎn)距離之和,取和最小的節(jié)點(diǎn)作為新的中心節(jié)點(diǎn);

    5)迭代步驟4)、5),直到中心節(jié)點(diǎn)不發(fā)生改變;

    6)對(duì)K-means算法產(chǎn)生的社區(qū)按照式(10)計(jì)算兩兩之間合并產(chǎn)生的ΔQ,將使ΔQ增加或減少的兩個(gè)社區(qū)進(jìn)行合并,并按照式(9)計(jì)算當(dāng)前Q值,保存Q值及當(dāng)前社區(qū)劃分;

    7)對(duì)步驟7)進(jìn)行迭代,直到所有節(jié)點(diǎn)都包含在一個(gè)社區(qū)或不能進(jìn)行社區(qū)合并;

    8)比較每次合并的Q值大小,取Q值最大的社區(qū)作為最終的結(jié)果。

    2.2 基于吸引力度量的社區(qū)識(shí)別過(guò)程

    算法2結(jié)束后,會(huì)對(duì)社會(huì)網(wǎng)絡(luò)中影響力較高的前3個(gè)節(jié)點(diǎn)進(jìn)行預(yù)劃分,識(shí)別出由高影響力節(jié)點(diǎn)組成的種子社區(qū)即為骨干社區(qū)。而網(wǎng)絡(luò)中影響力較低的n-3個(gè)零散節(jié)點(diǎn)可能處于孤立狀態(tài),因此需要將其合并到骨干社區(qū)之中。提出一種基于影響力分析的吸引力度量函數(shù),將網(wǎng)絡(luò)中零散節(jié)點(diǎn)添加到對(duì)其吸引力較高的骨干社區(qū)。對(duì)于任意骨干社區(qū)BC和節(jié)點(diǎn)u,BC對(duì)u的吸引力度量函數(shù)attract(BC,u)的定義為:

    式中:∑v∈BC∧v∈N(u)NG(u,v)表示節(jié)點(diǎn)u對(duì)當(dāng)前所在骨干社區(qū)的所有鄰居節(jié)點(diǎn)的影響力評(píng)分?jǐn)?shù)之和;Iscore(u)表示節(jié)點(diǎn)u的影響力評(píng)分。attract(BC,u)的值越大,節(jié)點(diǎn)u受骨干社區(qū)BC吸引而加入社區(qū)BC的可能性就越大,反之節(jié)點(diǎn)u成為骨干社區(qū)BC內(nèi)部節(jié)點(diǎn)的概率較低。

    設(shè)定了吸引力閾值β,對(duì)于節(jié)點(diǎn)u以及任意社區(qū)BC,若attract(BC,u)>β,則節(jié)點(diǎn)u屬于骨干社區(qū)BC。節(jié)點(diǎn)u可能同時(shí)加入多個(gè)社區(qū)結(jié)構(gòu),因此社區(qū)具有重疊性。在圖1中,假設(shè)骨干社區(qū)BC1包含的節(jié)點(diǎn)集合為{1,2,3,4,5},骨干社區(qū)BC1對(duì)節(jié)點(diǎn)6的吸引力函數(shù)為:

    假設(shè)BC2的骨干社區(qū)包含節(jié)點(diǎn)集合為{7,8,9,10},骨干社區(qū)BC2對(duì)節(jié)點(diǎn)6的吸引力為:

    若β=0.4,則節(jié)點(diǎn)6不屬于骨干社區(qū)BC1,屬于骨干社區(qū)BC2。

    下面給出基于種子擴(kuò)張的社區(qū)發(fā)現(xiàn)算法,算法的偽代碼如下:

    算法3 基于種子擴(kuò)張的社區(qū)發(fā)現(xiàn)算法

    輸入:高影響力節(jié)點(diǎn)構(gòu)成的社區(qū)集合BC,社區(qū)網(wǎng)絡(luò)G

    輸出:社區(qū)劃分結(jié)構(gòu)

    1)將網(wǎng)絡(luò)中高影響力節(jié)點(diǎn)標(biāo)記為true,其余n-3個(gè)節(jié)點(diǎn)標(biāo)記為false;

    2)將網(wǎng)絡(luò)中標(biāo)記為false的節(jié)點(diǎn)按照式(11)依次對(duì)骨干社區(qū)集合BC={BC1,BC2…BCn}進(jìn)行計(jì)算;

    3)對(duì)每個(gè)標(biāo)記為false的節(jié)點(diǎn)按照不同的骨干社區(qū)的計(jì)算,選取attract(BC,u)最大的值并將該節(jié)點(diǎn)加入到對(duì)應(yīng)的骨干網(wǎng)絡(luò);

    4)將加入骨干社區(qū)網(wǎng)絡(luò)的節(jié)點(diǎn)標(biāo)記為true;

    5)重復(fù)執(zhí)行步驟2),步驟3)和步驟4),直到社區(qū)網(wǎng)絡(luò)中所有節(jié)點(diǎn)都標(biāo)記為true;

    6)輸出社區(qū)劃分結(jié)構(gòu)。

    為了挖掘網(wǎng)絡(luò)的社區(qū)結(jié)構(gòu),首先將網(wǎng)絡(luò)除3個(gè)高影響力節(jié)點(diǎn)外的節(jié)點(diǎn)都標(biāo)記為false,若節(jié)點(diǎn)劃分到一個(gè)社區(qū)中,則該節(jié)點(diǎn)標(biāo)記為ture。

    在基于吸引力度量的種子擴(kuò)張過(guò)程的偽代碼中,第一步為遍歷n個(gè)高影響力評(píng)分節(jié)點(diǎn)社區(qū)集合。第二步以及第三步為計(jì)算集合BCi中對(duì)n-3吸引度量函數(shù)值attract(BCi,v)不小于給定的閾值β的節(jié)點(diǎn)。若節(jié)點(diǎn)attract(BCi,v)值大于或等于閾值,將該節(jié)點(diǎn)合并到當(dāng)前的社區(qū)中,若一個(gè)節(jié)點(diǎn)屬于多個(gè)社區(qū),則稱該節(jié)點(diǎn)為重疊節(jié)點(diǎn)。將該節(jié)點(diǎn)加入到attract(BCi,v)值最大的骨干社區(qū)中。最后將加入骨干社區(qū)的節(jié)點(diǎn)標(biāo)記為true,得到最終的社區(qū)結(jié)構(gòu)劃分結(jié)果。社區(qū)合并本質(zhì)上是一種貪心策略,然而由于網(wǎng)絡(luò)本身存在稀疏性,在實(shí)際的種子擴(kuò)張中零散節(jié)點(diǎn)鄰接社區(qū)要遠(yuǎn)小于其鄰接的社區(qū)數(shù),這有效的控制了算法的時(shí)間開(kāi)銷(xiāo)。

    2.3 算法的時(shí)間復(fù)雜度分析

    假設(shè)社會(huì)網(wǎng)絡(luò)G包含有n個(gè)節(jié)點(diǎn)和m條邊,算法1在求解節(jié)點(diǎn)的影響力強(qiáng)度時(shí),時(shí)間復(fù)雜度為O(n)。在算法2中基于K-means聚類(lèi)算法將O(l)個(gè)節(jié)點(diǎn)進(jìn)行合并得出骨干網(wǎng)絡(luò),其時(shí)間復(fù)雜度為O(kl),最壞的情況下將O(n)個(gè)節(jié)點(diǎn)進(jìn)行合并,故算法2的時(shí)間復(fù)雜度為O(kn)。在算法3中,在社區(qū)擴(kuò)展階段將O(n)個(gè)節(jié)點(diǎn)加入到O(k)個(gè)骨干社區(qū)中,該過(guò)程重復(fù)執(zhí)行O(n)次,所以找到所有的最終社區(qū)的時(shí)間復(fù)雜度為O(kn2),若k<

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

    本節(jié)給出算法在人工合成網(wǎng)絡(luò)以及真實(shí)數(shù)據(jù)集上的運(yùn)行結(jié)果。實(shí)驗(yàn)的運(yùn)行環(huán)境為Intel Core i5-7300HQ 2.5GHz處理器,8GB內(nèi)存,Windows 10操作系統(tǒng),算法采用Python與Matlab R2016a混合編程。

    3.1 NMI指數(shù)

    Danon等[20]提出歸一化互信息(normalized mutual information,NMI),度量社區(qū)集合之間的相似性。NMI基于混淆矩陣N,其中行表示真實(shí)的社區(qū)結(jié)構(gòu),列表示算法生成的社區(qū)。Ni表示矩陣N第i行所有元素的總和,Nj表示矩陣N第j列所有元素的總和。Nij是矩陣N的元素,表示同時(shí)屬于真實(shí)社區(qū)i和算法生成社區(qū)j的節(jié)點(diǎn)數(shù)量。NMI公式定義為:

    式中:CA代表網(wǎng)絡(luò)中真實(shí)社區(qū)的數(shù)量;CB代表由算法生成社區(qū)的數(shù)量。如果探測(cè)出的社區(qū)結(jié)構(gòu)和真實(shí)網(wǎng)絡(luò)的社區(qū)結(jié)構(gòu)相同則NMI=1;反之,如果所探測(cè)的社區(qū)結(jié)構(gòu)和真實(shí)網(wǎng)絡(luò)完全不同,則NMI=0。傳統(tǒng)意義上的歸一化互信息被用來(lái)量化兩個(gè)分布之間的差異性,為了使之適用于社區(qū)差異性量化,實(shí)驗(yàn)中采用的了NMI的簡(jiǎn)化版本[21]。

    3.2 F-SCORE

    F-Score也是常用的社區(qū)發(fā)現(xiàn)算法度量指標(biāo)。首先定義算法生成社區(qū)在全部社區(qū)中的比例。

    式中,CR代表高質(zhì)量節(jié)點(diǎn)所在的真實(shí)社區(qū),CF代表高質(zhì)量節(jié)點(diǎn)所在的探測(cè)社區(qū)。

    Precision是正確劃分節(jié)點(diǎn)的數(shù)量除以節(jié)點(diǎn)在CF中的數(shù)量,其定義為:

    結(jié)合式(15)以及式(16),F(xiàn)-Score定義為:

    3.3 人工合成網(wǎng)絡(luò)社區(qū)識(shí)別結(jié)果

    利用LFR-Benchmark提供的MATLAB數(shù)據(jù)生成器生成網(wǎng)絡(luò),網(wǎng)絡(luò)中的部分重要參數(shù)設(shè)置如下。首先,生成網(wǎng)絡(luò)規(guī)模為N=1 000的復(fù)雜網(wǎng)絡(luò),其中最小社區(qū)尺寸為10、最大社區(qū)尺寸為100。接下來(lái)設(shè)置網(wǎng)絡(luò)的平均度為5、最大節(jié)點(diǎn)度為25、混合參數(shù)u的取值范圍為0.1~0.8。在LFR Benchmark中,混合參數(shù)是一個(gè)獨(dú)立參數(shù),其值越大社區(qū)越難以發(fā)現(xiàn),因此通常用來(lái)刻畫(huà)算法的魯棒性。

    首先驗(yàn)證了不同混合參數(shù)下K-means聚類(lèi)算法中的k取值不同時(shí)NMI指標(biāo)變化情況,實(shí)驗(yàn)結(jié)果如圖2所示??梢钥闯觯跏紩r(shí)社區(qū)質(zhì)量會(huì)隨種子社區(qū)的不斷增多而不斷上升,但達(dá)到某一極值時(shí),社區(qū)質(zhì)量開(kāi)始下降。這是因?yàn)槌跏挤N子過(guò)多會(huì)導(dǎo)致社區(qū)個(gè)數(shù)增加,進(jìn)而降低社區(qū)內(nèi)部的緊密程度。另一個(gè)發(fā)現(xiàn)是,社區(qū)結(jié)構(gòu)較為模糊的網(wǎng)絡(luò)中(u=0.8),需要更多的種子社區(qū)才能達(dá)到較好的發(fā)現(xiàn)效果。

    接下來(lái),驗(yàn)證TSCDA算法的性能。與Clauset[12]、LWP[13]、Chen[14]、LS[15]、LCDMC[16]以及VI[17]6種社區(qū)識(shí)別算法在人工合成網(wǎng)絡(luò)社區(qū)進(jìn)行比較,為了保證算法結(jié)果的客觀性,將相關(guān)參數(shù)按原有文獻(xiàn)加以設(shè)定,算法的NMI指標(biāo)以及F-Score在圖3以及圖4給出。

    如圖2及圖3所示,當(dāng)混合參數(shù)u等于0.1時(shí),TSCDA在NMI以及F-Score的性能指標(biāo)明顯優(yōu)于其它算法,具有一定的性能優(yōu)勢(shì)。隨著混合參數(shù)u的不斷增加,TSCDA算法和其他6種算法的NMI指標(biāo)以及F-Score均呈現(xiàn)明顯的下降趨勢(shì),當(dāng)混合參數(shù)u等于0.5以及0.6時(shí),TSCDA算法略優(yōu)于VI算法和LCDMC算法,優(yōu)于LWP以及LS算法;當(dāng)混合參數(shù)u等于0.7以及0.8時(shí),各算法性能較為接近,但總體來(lái)看,TSCDA算法仍然要優(yōu)于其他6種社區(qū)識(shí)別方法,這表明TSCDA算法在結(jié)構(gòu)指標(biāo)上更貼近于內(nèi)嵌社區(qū)。

    3.4 真實(shí)網(wǎng)絡(luò)社區(qū)識(shí)別結(jié)果

    為了驗(yàn)證TSCDA算法在真實(shí)網(wǎng)絡(luò)社區(qū)識(shí)別的性能,本節(jié)選取了在5種真實(shí)網(wǎng)絡(luò)上,對(duì)比TSCDA算法與前文6種社區(qū)識(shí)別算法的性能差異。表2列出了所選擇的真實(shí)網(wǎng)絡(luò)的網(wǎng)絡(luò)特征,其中Node代表節(jié)點(diǎn)個(gè)數(shù)、Link代表鏈接條數(shù)、d—表示網(wǎng)絡(luò)的平均度、|C—|表示真實(shí)社區(qū)的平均大小、Community表示社區(qū)個(gè)數(shù)。我們從Amazon社區(qū)中移除了前5 000個(gè)子社區(qū),保留了它的不同的社區(qū)。

    還在DBLP計(jì)算機(jī)學(xué)科文獻(xiàn)數(shù)據(jù)網(wǎng)絡(luò)上測(cè)試了TSCDA算法與Clauset算法、LWP算法、Chen算法、LS算法、LCDMC算法以及VI算法的NMI指標(biāo)以及F-Score指標(biāo)。DBLP文獻(xiàn)網(wǎng)絡(luò)的特性如表3所示。|C|代表社區(qū)規(guī)模的范圍,例如,(0,10]代表樣本網(wǎng)絡(luò)包含的社區(qū)大小大于0,小于等于10。我們同樣移除了DBLP網(wǎng)絡(luò)前5 000個(gè)社區(qū)中的子社區(qū),并保留了它們不同的社區(qū)結(jié)構(gòu)。

    表4給出了各算法在真實(shí)網(wǎng)絡(luò)上的NMI指標(biāo)以及F-Score指標(biāo)。相比其它6種方法,TSCDA具有一定的性能優(yōu)勢(shì)。一方面,TSCDA基于節(jié)點(diǎn)影響力擴(kuò)展社區(qū),因此在識(shí)別社區(qū)成員時(shí)準(zhǔn)確度更高。另一方面,本文在節(jié)點(diǎn)影響力計(jì)算中加入了邊介數(shù)屬性,不僅考慮了節(jié)點(diǎn)自身的重要程度,還考慮了鏈接在局部區(qū)域中的重要程度。LWP和LS算法在個(gè)別網(wǎng)絡(luò)中顯示出比TSCDA更好的性能,例如在Football網(wǎng)絡(luò)中LWP算法的NMI指標(biāo)高于TSCDA約0.09、在Amazon網(wǎng)絡(luò)中LS算法的F-Score高于TSCDA約0.02。一個(gè)可能的原因是Football網(wǎng)絡(luò)的平均節(jié)點(diǎn)度較高,使得各算法在社區(qū)識(shí)別時(shí)更容易發(fā)現(xiàn)結(jié)構(gòu)較為緊湊的局部結(jié)構(gòu)。TSCDA算法在Dolpins網(wǎng)的準(zhǔn)確性有所下降,這是由于Dolpins網(wǎng)絡(luò)中的鏈接關(guān)系依靠于節(jié)點(diǎn)之間的接觸歷史,因此影響力識(shí)別時(shí)準(zhǔn)確度一般。但是相比其它6種算法TSCDA算法仍然表現(xiàn)出一定的穩(wěn)定性,進(jìn)一步說(shuō)明了TSCDA算法具有較高的魯棒性。

    圖5和圖6給出了各算法在DBLP引文網(wǎng)絡(luò)上的性能指標(biāo)。實(shí)驗(yàn)結(jié)果顯示,當(dāng)社區(qū)規(guī)模在0到10時(shí)(Label 0),TSCDA在社區(qū)識(shí)別性能上略低于其他6種算法表現(xiàn)的差。然而當(dāng)社區(qū)規(guī)模大于10時(shí)(Label 1~10),TSCDA的表現(xiàn)性能明顯優(yōu)于Clauset、LWP、Chen、LS、LCDMC和VI算法,這表明TSCDA善于挖掘個(gè)數(shù)較多的社區(qū)結(jié)果,社區(qū)粒度更加細(xì)致。

    4 結(jié) 論

    為了解決傳統(tǒng)社區(qū)識(shí)別方法社區(qū)穩(wěn)定性較差這一問(wèn)題,從節(jié)點(diǎn)相似性以及鏈接相似性兩個(gè)角度入手,提出一種多點(diǎn)種子預(yù)劃分的社區(qū)發(fā)現(xiàn)算法TSCDA。首先,在種子擴(kuò)張算法的基礎(chǔ)上,提出利用高質(zhì)量節(jié)點(diǎn)構(gòu)建骨干社區(qū)進(jìn)行社區(qū)識(shí)別,提高社區(qū)結(jié)構(gòu)的致密性;其次,將邊介數(shù)屬性加入到節(jié)點(diǎn)影響力的計(jì)算中,提高影響力計(jì)算中的結(jié)構(gòu)相關(guān)性。人工合成網(wǎng)絡(luò)和真實(shí)網(wǎng)絡(luò)數(shù)據(jù)集上的仿真結(jié)果驗(yàn)證了TSCDA算法的效率和識(shí)別性能。

    參 考 文 獻(xiàn):

    [1] SANTO Fortunato. Community Detection in Graphs[J]. Physics Reports,2009,486(3):26.

    [2] JAVED M A, YOUNIS M S, LATIF S, et al. Community Detection in Networks:A Multidisciplinary Review[J]. Journal of Network and Computer Applications,2018, 108:87.

    [3] AMOR B, VUIK S, CALLAHAN R, et al. Community Detection and Role Identification in Directed Networks: Understanding the Twitter Network of the Caredata Debate[J]. Dynamic Networks and Cyber-security, 2015, 25(2):19.

    [4] BAZZI M, PORTER M A, WILLIAMS S, et al. Community Detection in Temporal Multilayer Networks, with an Application to Correlation Networks[J]. Multiscale Modeling & Simulation, 2016, 14(1):1.

    [5] SHARMA S, SINGH A. An Efficient Method for Link Prediction in Complex Multiplex Networks[C]// 2015 11th International Conference on Signal-Image Technology & Internet-Based Systems(SITIS). IEEE, 2015:452.

    [6] DING X, ZHANG J, YANG J. A Robust Two-stage Algorithm for Local Community Detection[J]. Knowledge Based Systems, 2018, 14(7):62.

    [7] LANCICHINETTI A, FORTUNATO S, KERTSZ, JNOS. Detecting the Overlapping and Hierarchical Community Structure in Complex Networks[J]. New Journal of Physics, 2009, 11(3):15.

    [8] COSCIA M, ROSSETTI G, GIANNOTTI F, et al. DEMON:a Local-First Discovery Method for Overlapping Communities[C]// Proceedings of the 18th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, 2012:19.

    [9] BAUMES J, GOLDBERG M, MAGDON-ISMAIL M. Efficient Identification of Overlapping Communities[C]// International Conference on Intelligence and Security Informatics, 2005:11.

    [10]李婕, 王興偉, 郭靜, 等. 面向移動(dòng)通信網(wǎng)絡(luò)的局部擴(kuò)張群組構(gòu)造方法[J]. 東北大學(xué)學(xué)報(bào)(自然科學(xué)版),2017,38(12):1691.

    LI Jie, WANG Xingwei, GUO Jing, et al. Clique Percolation Based Local Fitness Method for User Clustering in Telecommunication Network[J]. Journal of Northeastern University Natural(Science), 2017,38(12):1691.

    [11]SU Y, WANG B, CHENG F, et al. An Algorithm Based on Positive and Negative Links for Community Detection in Signed Networks[J]. Scientific Reports, 2017, 7(1):10874.

    [12]CLAUSET, ARON. Finding Local Community Structure in Networks[J]. Physical Review E Statistical Nonlinear & Soft Matter Physics, 2005, 72(2):26.

    [13]LUO F , WANG J Z , PROMISLOW E. Exploring Local Community Structures in Large Networks[J]. Web Intelligence & Agent Systems, 2008, 6(4):387.

    [14]JIYANG Chen, OSMAR Zaane, RANDY Goebel. Local Community Identification in Social Networks[C]// International Conference on Advances in Social Network Analysis & Mining. IEEE Computer Society, 2009:79.

    [15]YING JUN, WU, HAN, et al. Local Community Detection Using Link Similarity[J]. Journal of Computer Science & Technology, 2012, 23(3):69.

    [16]FANRONG M , MU Z , YONG Z , et al. Local Community Detection in Complex Networks Based on Maximum Cliques Extension[J]. Mathematical Problems in Engineering, 2014, 36(6):1.

    [17]YAOY , WU W , LEI M , et al. Community Detection Based on Variable Vertex Influence[C]// IEEE International Conference on Data Science in Cyberspace. IEEE, 2016:145.

    [18]齊金山, 梁循, 王怡. 基于種子節(jié)點(diǎn)選擇的重疊社區(qū)發(fā)現(xiàn)算法[J]. 計(jì)算機(jī)應(yīng)用研究,2017(12):20.

    QI Jinshan, LIANG Xun, WANG Yi. Overlapping Community Detection Algorithm Based on Selection of Seed Nodes[J]. Application Research of Computers,2017(12):20.

    [19]NEWMANM E J. Fast Algorithm for Detecting Community Structure in Networks[J]. Physical Review E, 2004, 69(6):133.

    [20]DANON, LEON, DUCH, et al. Comparing Community Structure Identification[J]. Journal of Statistical Mechanics, 2005, 206(9):69.

    [21]BAGROW, JAMES P. Evaluating Local Community Methods in Networks[J]. Journal of Statistical Mechanics:Theory and Experiment, 2008, 8(5):26.

    (編輯:溫澤宇)

    叶爱在线成人免费视频播放| 亚洲av中文字字幕乱码综合| 九色成人免费人妻av| 日本黄色视频三级网站网址| 可以在线观看的亚洲视频| 在线视频色国产色| 精品国产乱子伦一区二区三区| 亚洲无线在线观看| 国产精品国产高清国产av| 在线观看日韩欧美| 成人三级黄色视频| 757午夜福利合集在线观看| 精品一区二区三区四区五区乱码| 亚洲精品美女久久av网站| 91av网一区二区| 亚洲精品在线美女| 亚洲欧美精品综合久久99| 亚洲精品美女久久av网站| 午夜福利欧美成人| 午夜成年电影在线免费观看| 久久亚洲精品不卡| 亚洲成av人片在线播放无| 久久九九热精品免费| 久久精品夜夜夜夜夜久久蜜豆| 成人一区二区视频在线观看| 桃红色精品国产亚洲av| 88av欧美| 国产av在哪里看| 日韩欧美国产在线观看| 男女床上黄色一级片免费看| 女同久久另类99精品国产91| 两人在一起打扑克的视频| 国产精品一区二区三区四区久久| 精品福利观看| x7x7x7水蜜桃| 国产精品日韩av在线免费观看| 久久香蕉精品热| 亚洲avbb在线观看| 日韩欧美免费精品| 亚洲熟女毛片儿| 中文字幕久久专区| 成人鲁丝片一二三区免费| 免费看十八禁软件| 久久久久久久久久黄片| 淫妇啪啪啪对白视频| 在线观看日韩欧美| 国产三级黄色录像| 日本五十路高清| 国产精品 欧美亚洲| 成年女人永久免费观看视频| 日韩中文字幕欧美一区二区| 国产av一区在线观看免费| 亚洲黑人精品在线| 搡老妇女老女人老熟妇| 亚洲国产看品久久| 首页视频小说图片口味搜索| 免费电影在线观看免费观看| 大型黄色视频在线免费观看| 国产精品一区二区三区四区免费观看 | 色精品久久人妻99蜜桃| 国产精品九九99| 身体一侧抽搐| 成人永久免费在线观看视频| 制服人妻中文乱码| cao死你这个sao货| 国产精品一区二区三区四区久久| 男人舔女人下体高潮全视频| 伊人久久大香线蕉亚洲五| 国产精品国产高清国产av| 国产精品亚洲av一区麻豆| 国产高清视频在线观看网站| 91av网一区二区| 国产激情欧美一区二区| 免费在线观看成人毛片| 亚洲18禁久久av| 999久久久国产精品视频| 黑人操中国人逼视频| 在线播放国产精品三级| 一个人观看的视频www高清免费观看 | 国产成人一区二区三区免费视频网站| 午夜a级毛片| 无限看片的www在线观看| 在线看三级毛片| 精品国产乱码久久久久久男人| 色哟哟哟哟哟哟| 欧美成人免费av一区二区三区| 制服丝袜大香蕉在线| 久久久国产欧美日韩av| 成人永久免费在线观看视频| 成人三级做爰电影| 日日干狠狠操夜夜爽| 国产一区在线观看成人免费| 日本熟妇午夜| 美女cb高潮喷水在线观看 | 国产成人影院久久av| 日韩av在线大香蕉| 视频区欧美日本亚洲| 人妻夜夜爽99麻豆av| 熟女少妇亚洲综合色aaa.| 日韩av在线大香蕉| 一夜夜www| 999久久久精品免费观看国产| 蜜桃久久精品国产亚洲av| 欧美激情在线99| 亚洲熟妇中文字幕五十中出| 精品电影一区二区在线| 欧美中文综合在线视频| 日韩高清综合在线| 天天添夜夜摸| 国产黄a三级三级三级人| 亚洲专区国产一区二区| 久久久久精品国产欧美久久久| 久久午夜亚洲精品久久| 亚洲av美国av| 日韩人妻高清精品专区| 精品久久久久久久人妻蜜臀av| 午夜激情福利司机影院| 91av网站免费观看| 午夜a级毛片| 亚洲精品在线美女| 国产在线精品亚洲第一网站| 男插女下体视频免费在线播放| 久久香蕉国产精品| 丁香欧美五月| 亚洲中文字幕一区二区三区有码在线看 | 香蕉国产在线看| 天堂动漫精品| 大型黄色视频在线免费观看| 国产精品九九99| 成年女人看的毛片在线观看| 国产人伦9x9x在线观看| 国产精华一区二区三区| 亚洲中文av在线| 亚洲自拍偷在线| 国产精品乱码一区二三区的特点| 国产av一区在线观看免费| 在线免费观看不下载黄p国产 | 精品无人区乱码1区二区| 免费看十八禁软件| 日本 av在线| 国产午夜福利久久久久久| 亚洲一区二区三区不卡视频| 性欧美人与动物交配| 母亲3免费完整高清在线观看| 亚洲乱码一区二区免费版| 国产精品亚洲一级av第二区| 国产精品九九99| 99精品在免费线老司机午夜| aaaaa片日本免费| 久久午夜亚洲精品久久| 午夜免费激情av| 日本 欧美在线| 亚洲熟妇熟女久久| 久久久国产欧美日韩av| 国产三级中文精品| 一级毛片高清免费大全| 中文在线观看免费www的网站| 少妇的丰满在线观看| 美女 人体艺术 gogo| 久久久久久久精品吃奶| 热99在线观看视频| 18禁观看日本| 日韩国内少妇激情av| 久久亚洲真实| 中文字幕人妻丝袜一区二区| 成人特级黄色片久久久久久久| 看黄色毛片网站| 色尼玛亚洲综合影院| 日本a在线网址| 久久精品国产亚洲av香蕉五月| 国产高清视频在线播放一区| av天堂中文字幕网| 久久精品国产综合久久久| e午夜精品久久久久久久| 午夜精品久久久久久毛片777| 久久久国产欧美日韩av| 欧美日韩精品网址| 少妇的逼水好多| 高清在线国产一区| 女警被强在线播放| 91久久精品国产一区二区成人 | 狠狠狠狠99中文字幕| 中文字幕人妻丝袜一区二区| www日本黄色视频网| 欧美日韩黄片免| 国产精品亚洲美女久久久| 国产精品一区二区免费欧美| 女人被狂操c到高潮| 美女cb高潮喷水在线观看 | 精品一区二区三区四区五区乱码| 亚洲国产精品成人综合色| 久久久久久大精品| 亚洲激情在线av| 老熟妇仑乱视频hdxx| 亚洲中文字幕日韩| av天堂中文字幕网| 久久香蕉精品热| 看黄色毛片网站| 一区二区三区国产精品乱码| 亚洲国产欧美网| 色综合婷婷激情| 亚洲第一电影网av| 成年女人毛片免费观看观看9| 性色avwww在线观看| 美女cb高潮喷水在线观看 | 免费电影在线观看免费观看| 亚洲欧美精品综合一区二区三区| 国产亚洲精品久久久com| 欧美性猛交黑人性爽| 亚洲美女黄片视频| 一级黄色大片毛片| 亚洲色图av天堂| 99久久成人亚洲精品观看| 免费在线观看成人毛片| 又爽又黄无遮挡网站| 精品日产1卡2卡| 人妻丰满熟妇av一区二区三区| 中文在线观看免费www的网站| 国产单亲对白刺激| 美女cb高潮喷水在线观看 | 制服丝袜大香蕉在线| 国产成人福利小说| 亚洲激情在线av| 日本 av在线| 18禁黄网站禁片免费观看直播| 悠悠久久av| 国产乱人伦免费视频| 欧美日韩精品网址| 亚洲电影在线观看av| 在线观看免费视频日本深夜| 亚洲一区二区三区不卡视频| 久久国产精品人妻蜜桃| 国产在线精品亚洲第一网站| 99国产精品99久久久久| 韩国av一区二区三区四区| 久久午夜综合久久蜜桃| 国产亚洲欧美98| 国产麻豆成人av免费视频| 国产欧美日韩一区二区精品| 免费看十八禁软件| 国产精品精品国产色婷婷| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av美国av| 久久天堂一区二区三区四区| 精品久久久久久,| 国产成+人综合+亚洲专区| 色综合婷婷激情| 可以在线观看毛片的网站| 亚洲av成人精品一区久久| 老司机福利观看| 久久久国产成人精品二区| 久久久久性生活片| 非洲黑人性xxxx精品又粗又长| 欧美色视频一区免费| 亚洲av第一区精品v没综合| 一夜夜www| 亚洲av第一区精品v没综合| 亚洲中文av在线| 19禁男女啪啪无遮挡网站| 黄色视频,在线免费观看| 亚洲av第一区精品v没综合| 嫩草影院入口| 天天添夜夜摸| 国产三级黄色录像| 欧美性猛交╳xxx乱大交人| 男人舔女人的私密视频| 亚洲精品一卡2卡三卡4卡5卡| 日日夜夜操网爽| 亚洲欧美日韩高清在线视频| xxx96com| 伦理电影免费视频| 色综合欧美亚洲国产小说| 桃色一区二区三区在线观看| 国内精品久久久久久久电影| 欧美日韩一级在线毛片| 国产淫片久久久久久久久 | 99视频精品全部免费 在线 | 大型黄色视频在线免费观看| 人人妻,人人澡人人爽秒播| 一a级毛片在线观看| 99国产极品粉嫩在线观看| 成人特级黄色片久久久久久久| 亚洲专区国产一区二区| 亚洲欧美精品综合一区二区三区| 国产精品亚洲一级av第二区| 欧美日韩黄片免| 精品久久久久久久人妻蜜臀av| 日韩欧美在线乱码| 亚洲专区国产一区二区| 国产精品久久久久久亚洲av鲁大| 国内精品一区二区在线观看| 欧美一级a爱片免费观看看| 丰满人妻熟妇乱又伦精品不卡| 怎么达到女性高潮| 免费无遮挡裸体视频| 久久精品91蜜桃| 国产乱人视频| 这个男人来自地球电影免费观看| 久久久久久国产a免费观看| 最近视频中文字幕2019在线8| 久久久久久大精品| 欧美性猛交╳xxx乱大交人| 国产伦人伦偷精品视频| 露出奶头的视频| 国产精品av久久久久免费| 观看美女的网站| 岛国在线免费视频观看| 久久久久久久精品吃奶| 国产精品久久久久久人妻精品电影| 亚洲熟妇熟女久久| 色吧在线观看| 久久国产乱子伦精品免费另类| 亚洲乱码一区二区免费版| 九色成人免费人妻av| 久久久久九九精品影院| 久久久久精品国产欧美久久久| 国产亚洲精品一区二区www| 欧美精品啪啪一区二区三区| 免费在线观看影片大全网站| 他把我摸到了高潮在线观看| 90打野战视频偷拍视频| 两个人视频免费观看高清| 亚洲色图 男人天堂 中文字幕| 在线免费观看不下载黄p国产 | 99热这里只有是精品50| 日日摸夜夜添夜夜添小说| 老司机深夜福利视频在线观看| 国产精品99久久99久久久不卡| 一区二区三区激情视频| 日韩欧美 国产精品| 又黄又爽又免费观看的视频| 九九在线视频观看精品| 欧美又色又爽又黄视频| 久久久国产欧美日韩av| 亚洲片人在线观看| 久久久国产精品麻豆| 久久久久九九精品影院| 久久精品国产99精品国产亚洲性色| 国产高清videossex| 亚洲男人的天堂狠狠| 伦理电影免费视频| xxx96com| 亚洲专区字幕在线| 日本黄大片高清| 悠悠久久av| 国产av在哪里看| 男人舔女人下体高潮全视频| 免费电影在线观看免费观看| 久久久久久大精品| 亚洲精品中文字幕一二三四区| 一本久久中文字幕| 日本撒尿小便嘘嘘汇集6| 中文字幕精品亚洲无线码一区| 欧美丝袜亚洲另类 | 成人欧美大片| 国产av一区在线观看免费| 日日夜夜操网爽| 亚洲欧美日韩东京热| 亚洲国产高清在线一区二区三| 欧美日韩瑟瑟在线播放| 欧美三级亚洲精品| 久久久久久人人人人人| 精品久久久久久,| 亚洲精品粉嫩美女一区| 91久久精品国产一区二区成人 | 国产视频内射| 国产单亲对白刺激| 亚洲人成网站高清观看| 天堂av国产一区二区熟女人妻| 精品午夜福利视频在线观看一区| 亚洲国产精品久久男人天堂| 99国产精品99久久久久| 三级男女做爰猛烈吃奶摸视频| www.自偷自拍.com| 草草在线视频免费看| 欧美性猛交黑人性爽| 嫩草影院入口| aaaaa片日本免费| 欧美中文日本在线观看视频| 久久久久久久久免费视频了| 无限看片的www在线观看| 国产一区二区三区视频了| 亚洲成人中文字幕在线播放| 美女高潮喷水抽搐中文字幕| 亚洲av第一区精品v没综合| 1024手机看黄色片| 亚洲avbb在线观看| 色综合婷婷激情| 日韩精品青青久久久久久| 久久国产乱子伦精品免费另类| 少妇人妻一区二区三区视频| 亚洲成人精品中文字幕电影| 哪里可以看免费的av片| 美女cb高潮喷水在线观看 | 精品午夜福利视频在线观看一区| 亚洲第一电影网av| 国产精品久久视频播放| 欧美日韩福利视频一区二区| 亚洲精品456在线播放app | 国产免费男女视频| 999久久久精品免费观看国产| 欧美av亚洲av综合av国产av| 国产av不卡久久| 亚洲精华国产精华精| 成年女人毛片免费观看观看9| 99国产精品99久久久久| 国语自产精品视频在线第100页| 大型黄色视频在线免费观看| 99久久精品热视频| 亚洲av美国av| 99国产精品一区二区蜜桃av| 国产精华一区二区三区| а√天堂www在线а√下载| 欧美成狂野欧美在线观看| 狂野欧美激情性xxxx| 女警被强在线播放| 国产在线精品亚洲第一网站| 中文字幕高清在线视频| 欧美日韩综合久久久久久 | 99热这里只有精品一区 | 女警被强在线播放| 欧美大码av| 一二三四在线观看免费中文在| 亚洲片人在线观看| 我要搜黄色片| 久久久久九九精品影院| 欧美另类亚洲清纯唯美| 久久精品国产综合久久久| 国产午夜精品论理片| 日韩国内少妇激情av| 波多野结衣高清无吗| 日韩免费av在线播放| 国产精品一及| 国产成人精品久久二区二区91| 桃色一区二区三区在线观看| 国产午夜精品论理片| 精品电影一区二区在线| 18禁裸乳无遮挡免费网站照片| 午夜福利免费观看在线| 夜夜夜夜夜久久久久| 精品不卡国产一区二区三区| 亚洲乱码一区二区免费版| 18美女黄网站色大片免费观看| 国产精品 欧美亚洲| 精品国产亚洲在线| 久久99热这里只有精品18| 午夜免费激情av| av黄色大香蕉| 亚洲av五月六月丁香网| 午夜福利视频1000在线观看| 日本 欧美在线| 国产高清videossex| 999久久久精品免费观看国产| 偷拍熟女少妇极品色| 国产v大片淫在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | 天堂av国产一区二区熟女人妻| 午夜福利高清视频| 亚洲狠狠婷婷综合久久图片| 99热6这里只有精品| 看片在线看免费视频| 日韩欧美精品v在线| 欧美日韩乱码在线| 久久久久精品国产欧美久久久| 网址你懂的国产日韩在线| www日本黄色视频网| 国产成人影院久久av| 最新美女视频免费是黄的| 国产爱豆传媒在线观看| 麻豆国产97在线/欧美| 精品欧美国产一区二区三| 亚洲欧美日韩高清在线视频| 97碰自拍视频| 两个人视频免费观看高清| 男女那种视频在线观看| 亚洲欧美精品综合一区二区三区| 波多野结衣高清作品| 最近在线观看免费完整版| 香蕉久久夜色| 国产成人系列免费观看| 日韩欧美免费精品| 99久久成人亚洲精品观看| 亚洲人成网站在线播放欧美日韩| 国产欧美日韩精品亚洲av| 久久热在线av| 国产精品影院久久| 国产真实乱freesex| 女同久久另类99精品国产91| 国产美女午夜福利| 老汉色∧v一级毛片| 免费观看的影片在线观看| 亚洲av成人不卡在线观看播放网| av女优亚洲男人天堂 | 亚洲五月天丁香| 欧美不卡视频在线免费观看| 真人一进一出gif抽搐免费| 亚洲精品一区av在线观看| 真人一进一出gif抽搐免费| 国产男靠女视频免费网站| 欧洲精品卡2卡3卡4卡5卡区| 99久久精品热视频| 国产激情偷乱视频一区二区| 国产伦人伦偷精品视频| 欧美日韩瑟瑟在线播放| 午夜精品在线福利| 午夜影院日韩av| 欧美中文日本在线观看视频| 美女黄网站色视频| 黑人操中国人逼视频| 国产一区二区三区在线臀色熟女| 岛国在线观看网站| 毛片女人毛片| 亚洲精品在线美女| 国产69精品久久久久777片 | 十八禁网站免费在线| 国产熟女xx| 欧美日韩中文字幕国产精品一区二区三区| 曰老女人黄片| 99精品久久久久人妻精品| 岛国在线观看网站| 国产淫片久久久久久久久 | 很黄的视频免费| 色在线成人网| 亚洲人成网站高清观看| 亚洲熟妇中文字幕五十中出| 国产亚洲精品av在线| 欧美成人一区二区免费高清观看 | 一边摸一边抽搐一进一小说| 色哟哟哟哟哟哟| 亚洲色图 男人天堂 中文字幕| 无遮挡黄片免费观看| 在线免费观看不下载黄p国产 | or卡值多少钱| 久久久久久久久中文| 天天躁日日操中文字幕| 亚洲专区中文字幕在线| aaaaa片日本免费| 一a级毛片在线观看| 国产精品一区二区免费欧美| 免费高清视频大片| 岛国在线观看网站| 午夜精品在线福利| 伊人久久大香线蕉亚洲五| 欧美中文综合在线视频| 波多野结衣高清作品| 免费在线观看日本一区| 亚洲av第一区精品v没综合| 久久精品亚洲精品国产色婷小说| 色精品久久人妻99蜜桃| 人人妻,人人澡人人爽秒播| 丰满人妻熟妇乱又伦精品不卡| 欧美黄色片欧美黄色片| 日韩国内少妇激情av| 一本一本综合久久| 97碰自拍视频| 级片在线观看| 午夜福利欧美成人| 国产欧美日韩精品一区二区| 欧美激情在线99| 国产乱人伦免费视频| 三级毛片av免费| 他把我摸到了高潮在线观看| 两个人看的免费小视频| 日本a在线网址| av在线蜜桃| 国产成人aa在线观看| 日本 av在线| 色综合婷婷激情| 国产欧美日韩精品亚洲av| 一个人看视频在线观看www免费 | 亚洲成a人片在线一区二区| 欧美在线黄色| 亚洲国产中文字幕在线视频| 19禁男女啪啪无遮挡网站| 白带黄色成豆腐渣| 99riav亚洲国产免费| 97碰自拍视频| 国产人伦9x9x在线观看| 亚洲熟妇熟女久久| 国产精品98久久久久久宅男小说| 亚洲欧美激情综合另类| 18禁黄网站禁片免费观看直播| 中文字幕久久专区| 啪啪无遮挡十八禁网站| 精品一区二区三区视频在线观看免费| 亚洲欧美精品综合久久99| 亚洲中文av在线| 香蕉久久夜色| 国产淫片久久久久久久久 | 久久久国产成人精品二区| 国产一区在线观看成人免费| 看片在线看免费视频| 中国美女看黄片| 一区二区三区高清视频在线| 国内久久婷婷六月综合欲色啪| 最近视频中文字幕2019在线8| 婷婷亚洲欧美| 三级毛片av免费| 1024香蕉在线观看| 99国产精品99久久久久| 男人的好看免费观看在线视频| 啦啦啦免费观看视频1| 制服人妻中文乱码| 欧美成人性av电影在线观看| 成人国产一区最新在线观看| 欧美色欧美亚洲另类二区| 五月玫瑰六月丁香| 亚洲色图av天堂| 亚洲欧美日韩高清专用| 久久婷婷人人爽人人干人人爱| 国产午夜精品论理片| 色精品久久人妻99蜜桃| 亚洲成人精品中文字幕电影| 国产激情欧美一区二区| 香蕉丝袜av|