国产无遮挡裸体免费直播视频,久久精品国产蜜臀av,动漫在线视频一区二区,欧亚日韩一区二区三区,久艹在线 免费视频,国产精品美女网站免费,正在播放 97超级视频在线观看,斗破苍穹年番在线观看免费,51最新乱码中文字幕

zip 的壓縮原理與實(shí)現(xiàn)

 更新時(shí)間:2007年02月15日 00:00:00   作者:  
無損數(shù)據(jù)壓縮是一件奇妙的事情,想一想,一串任意的數(shù)據(jù)能夠根據(jù)一定的規(guī)則轉(zhuǎn)換成只有原來 1/2 - 1/5 長(zhǎng)度的數(shù)據(jù),并且能夠按照相應(yīng)的規(guī)則還原到原來的樣子,聽起來真是很酷。
半年前,苦熬過初學(xué) vc 時(shí)那段艱難的學(xué)習(xí)曲線的我,對(duì) MFC、SDK 開始失望和不滿,這些雖然不算易學(xué),但和 DHTML 沒有實(shí)質(zhì)上的區(qū)別,都是調(diào)用微軟提供的各種各樣的函數(shù),不需要你自己去創(chuàng)建一個(gè)窗口,多線程編程時(shí),也不需要你自己去分配 CPU 時(shí)間。我也做過驅(qū)動(dòng),同樣,有DDK(微軟驅(qū)動(dòng)開發(fā)包),當(dāng)然,也有 DDK 的“參考手冊(cè)”,連一個(gè)最簡(jiǎn)單的數(shù)據(jù)結(jié)構(gòu)都不需要你自己做,一切都是函數(shù)、函數(shù)…… 
微軟的高級(jí)程序員編寫了函數(shù)讓我們這些搞應(yīng)用的去調(diào)用,我不想在這里貶低搞應(yīng)用的人,正是這些應(yīng)用工程師連接起了科學(xué)和社會(huì)之間的橋梁,將來可以做銷售,做管理,用自己逐漸積累起來的智慧和經(jīng)驗(yàn)在社會(huì)上打拼。
但是,在技術(shù)上來說,誠(chéng)實(shí)地說,這并不高深,不是嗎?第一流的公司如微軟、Sybase、Oracle 等總是面向社會(huì)大眾的,這樣才能有巨大的市場(chǎng)。但是他們往往也是站在社會(huì)的最頂層的:操作系統(tǒng)、編譯器、數(shù)據(jù)庫(kù)都值得一代代的專家去不斷研究。這些帝國(guó)般的企業(yè)之所以偉大,恐怕不是“有經(jīng)驗(yàn)”、“能吃苦”這些中國(guó)特色的概念所能涵蓋的,艱深的技術(shù)體系、現(xiàn)代的管理哲學(xué)、強(qiáng)大的市場(chǎng)能力都是缺一不可的吧。我們既然有志于技術(shù),并且正在起步階段,何必急不可耐地要轉(zhuǎn)去做“管理”,做“青年才俊”,那些所謂的“成功人士”的根底能有幾何,這樣子浮躁,胸中的規(guī)模和格局能有多大?

在我發(fā)現(xiàn)vc只是一個(gè)用途廣泛的編程工具,并不能代表“知識(shí)”、“技術(shù)”的時(shí)候,我有些失落,無所不能的不是我,而是 MFC、SDK、DDK,是微軟的工程師,他們做的,正是我想做的,或者說,我也想成為那種層次的人,現(xiàn)在我知道了,他們是專家,但這不會(huì)是一個(gè)夢(mèng),有一天我會(huì)做到的,為什么不能說出我的想法呢。
那時(shí)公司做的系統(tǒng)里有一個(gè)壓縮模塊,領(lǐng)導(dǎo)找了一個(gè) zlib 庫(kù),不讓我自己做壓縮算法,站在公司的立場(chǎng)上,我很理解,真的很理解,自己做算法要多久啊。但那時(shí)自己心中隱藏的一份倔強(qiáng)驅(qū)使我去尋找壓縮原理的資料,我完全沒有意識(shí)到,我即將打開一扇大門,進(jìn)入一個(gè)神奇的“數(shù)據(jù)結(jié)構(gòu)”的世界?!坝?jì)算機(jī)藝術(shù)”的第一線陽光,居然也照到了我這樣一個(gè)平凡的人的身上。

上面說到“計(jì)算機(jī)藝術(shù)”,或者進(jìn)一步細(xì)化說“計(jì)算機(jī)編程藝術(shù)”,聽起來很深?yuàn)W,很高雅,但是在將要進(jìn)入專業(yè)的壓縮算法的研究時(shí),我要請(qǐng)大家做的第一件事情是:忘掉自己的年齡、學(xué)歷,忘掉自己的社會(huì)身份,忘掉編程語言,忘掉“面向?qū)ο蟆?、“三層架?gòu)”等一切術(shù)語。把自己當(dāng)作一個(gè)小孩,有一雙求知的眼睛,對(duì)世界充滿不倦的、單純的好奇,唯一的前提是一個(gè)正常的具有人類理性思維能力的大腦。
下面就讓我們開始一段神奇的壓縮算法之旅吧:


1. 原理部分:
  有兩種形式的重復(fù)存在于計(jì)算機(jī)數(shù)據(jù)中,zip 就是對(duì)這兩種重復(fù)進(jìn)行了壓縮。
  一種是短語形式的重復(fù),即三個(gè)字節(jié)以上的重復(fù),對(duì)于這種重復(fù),zip用兩個(gè)數(shù)字:1.重復(fù)位置距當(dāng)前壓縮位置的距離;2.重復(fù)的長(zhǎng)度,來表示這個(gè)重復(fù),假設(shè)這兩個(gè)數(shù)字各占一個(gè)字節(jié),于是數(shù)據(jù)便得到了壓縮,這很容易理解。
  一個(gè)字節(jié)有 0 - 255 共 256 種可能的取值,三個(gè)字節(jié)有 256 * 256 * 256 共一千六百多萬種可能的情況,更長(zhǎng)的短語取值的可能情況以指數(shù)方式增長(zhǎng),出現(xiàn)重復(fù)的概率似乎極低,實(shí)則不然,各種類型的數(shù)據(jù)都有出現(xiàn)重復(fù)的傾向,一篇論文中,為數(shù)不多的術(shù)語傾向于重復(fù)出現(xiàn);一篇小說,人名和地名會(huì)重復(fù)出現(xiàn);一張上下漸變的背景圖片,水平方向上的像素會(huì)重復(fù)出現(xiàn);程序的源文件中,語法關(guān)鍵字會(huì)重復(fù)出現(xiàn)(我們寫程序時(shí),多少次前后copy、paste?),以幾十 K 為單位的非壓縮格式的數(shù)據(jù)中,傾向于大量出現(xiàn)短語式的重復(fù)。經(jīng)過上面提到的方式進(jìn)行壓縮后,短語式重復(fù)的傾向被完全破壞,所以在壓縮的結(jié)果上進(jìn)行第二次短語式壓縮一般是沒有效果的。
  第二種重復(fù)為單字節(jié)的重復(fù),一個(gè)字節(jié)只有256種可能的取值,所以這種重復(fù)是必然的。其中,某些字節(jié)出現(xiàn)次數(shù)可能較多,另一些則較少,在統(tǒng)計(jì)上有分布不均勻的傾向,這是容易理解的,比如一個(gè) ASCII 文本文件中,某些符號(hào)可能很少用到,而字母和數(shù)字則使用較多,各字母的使用頻率也是不一樣的,據(jù)說字母 e 的使用概率最高;許多圖片呈現(xiàn)深色調(diào)或淺色調(diào),深色(或淺色)的像素使用較多(這里順便提一下:png 圖片格式是一種無損壓縮,其核心算法就是 zip 算法,它和 zip 格式的文件的主要區(qū)別在于:作為一種圖片格式,它在文件頭處存放了圖片的大小、使用的顏色數(shù)等信息);上面提到的短語式壓縮的結(jié)果也有這種傾向:重復(fù)傾向于出現(xiàn)在離當(dāng)前壓縮位置較近的地方,重復(fù)長(zhǎng)度傾向于比較短(20字節(jié)以內(nèi))。這樣,就有了壓縮的可能:給 256 種字節(jié)取值重新編碼,使出現(xiàn)較多的字節(jié)使用較短的編碼,出現(xiàn)較少的字節(jié)使用較長(zhǎng)的編碼,這樣一來,變短的字節(jié)相對(duì)于變長(zhǎng)的字節(jié)更多,文件的總長(zhǎng)度就會(huì)減少,并且,字節(jié)使用比例越不均勻,壓縮比例就越大。
  在進(jìn)一步討論編碼的要求以及辦法前,先提一下:編碼式壓縮必須在短語式壓縮之后進(jìn)行,因?yàn)榫幋a式壓縮后,原先八位二進(jìn)制值的字節(jié)就被破壞了,這樣文件中短語式重復(fù)的傾向也會(huì)被破壞(除非先進(jìn)行解碼)。另外,短語式壓縮后的結(jié)果:那些剩下的未被匹配的單、雙字節(jié)和得到匹配的距離、長(zhǎng)度值仍然具有取值分布不均勻性,因此,兩種壓縮方式的順序不能變。
  在編碼式壓縮后,以連續(xù)的八位作為一個(gè)字節(jié),原先未壓縮文件中所具有的字節(jié)取值不均勻的傾向被徹底破壞,成為隨機(jī)性取值,根據(jù)統(tǒng)計(jì)學(xué)知識(shí),隨機(jī)性取值具有均勻性的傾向(比如拋硬幣試驗(yàn),拋一千次,正反面朝上的次數(shù)都接近于 500 次)。因此,編碼式壓縮后的結(jié)果無法再進(jìn)行編碼式壓縮。
  短語式壓縮和編碼式壓縮是目前計(jì)算機(jī)科學(xué)界研究出的僅有的兩種無損壓縮方法,它們都無法重復(fù)進(jìn)行,所以,壓縮文件無法再次壓縮(實(shí)際上,能反復(fù)進(jìn)行的壓縮算法是不可想象的,因?yàn)樽罱K會(huì)壓縮到 0 字節(jié))。
=====================================

(補(bǔ)充)

壓縮文件無法再次壓縮是因?yàn)椋?
1. 短語式壓縮去掉了三個(gè)字節(jié)以上的重復(fù),壓縮后的結(jié)果中包含的是未匹配的單雙字節(jié),和匹配距離、長(zhǎng)度的組合。這個(gè)結(jié)果當(dāng)然仍然可能包含三個(gè)字節(jié)以上的重復(fù),但是概率極低。因?yàn)槿齻€(gè)字節(jié)有 256 * 256 * 256 共一千六百多萬種可能的情況,一千六百萬分之一的概率導(dǎo)致匹配的距離很長(zhǎng),需要二進(jìn)制數(shù)24位來表示這個(gè)匹配距離,再加上匹配長(zhǎng)度就超過了三個(gè)字節(jié),得不償失。所以只能壓縮掉原始文件中“自然存在的,并非隨機(jī)的短語式重復(fù)傾向”。
2.編碼式壓縮利用各個(gè)單字節(jié)使用頻率不一樣的傾向,使定長(zhǎng)編碼變?yōu)椴欢ㄩL(zhǎng)編碼,給使用頻率高的字節(jié)更短的編碼,使用頻率低的字節(jié)更長(zhǎng)的編碼,起到壓縮的效果。如果把編碼式壓縮的“結(jié)果”按照8位作為1字節(jié),重新統(tǒng)計(jì)各字節(jié)的使用頻率,應(yīng)該是大致相等的。因?yàn)樾碌淖止?jié)使用頻率是隨機(jī)的。相等的頻率再去變換字節(jié)長(zhǎng)短是沒有意義的,因?yàn)樽兌痰淖止?jié)沒有比變長(zhǎng)的字節(jié)更多。

=======================================

  短語式重復(fù)的傾向和字節(jié)取值分布不均勻的傾向是可以壓縮的基礎(chǔ),兩種壓縮的順序不能互換的原因也說了,下面我們來看編碼式壓縮的要求及方法:

首先,為了使用不定長(zhǎng)的編碼表示單個(gè)字符,編碼必須符合“前綴編碼”的要求,即較短的編碼決不能是較長(zhǎng)編碼的前綴,反過來說就是,任何一個(gè)字符的編碼,都不是由另一個(gè)字符的編碼加上若干位 0 或 1 組成,否則解壓縮程序?qū)o法解碼。
看一下前綴編碼的一個(gè)最簡(jiǎn)單的例子:


符號(hào) 編碼
A 0
B 10
C 110
D 1110
E 11110

有了上面的碼表,你一定可以輕松地從下面這串二進(jìn)制流中分辨出真正的信息內(nèi)容了:

1110010101110110111100010 - DABBDCEAAB

要構(gòu)造符合這一要求的二進(jìn)制編碼體系,二叉樹是最理想的選擇??疾煜旅孢@棵二叉樹:

        根(root)
       0  |   1
       +-------+--------+
    0  | 1   0  |  1
    +-----+------+  +----+----+
    |     |  |     |
    a      |  d     e
     0  |  1
     +-----+-----+
     |     |
     b     c

要編碼的字符總是出現(xiàn)在樹葉上,假定從根向樹葉行走的過程中,左轉(zhuǎn)為0,右轉(zhuǎn)為1,則一個(gè)字符的編碼就是從根走到該字符所在樹葉的路徑。正因?yàn)樽址荒艹霈F(xiàn)在樹葉上,任何一個(gè)字符的路徑都不會(huì)是另一字符路徑的前綴路徑,符合要求的前綴編碼也就構(gòu)造成功了:

a - 00 b - 010 c - 011 d - 10 e - 11


接下來來看編碼式壓縮的過程:
為了簡(jiǎn)化問題,假定一個(gè)文件中只出現(xiàn)了 a,b,c,d ,e五種字符,它們的出現(xiàn)次數(shù)分別是
a : 6次
b : 15次
c : 2次
d : 9次
e : 1次
如果用定長(zhǎng)的編碼方式為這四種字符編碼: a : 000 b : 001 c : 010 d : 011 e : 100
那么整個(gè)文件的長(zhǎng)度是 3*6 + 3*15 + 3*2 + 3*9 + 3*1 = 99

用二叉樹表示這四種編碼(其中葉子節(jié)點(diǎn)上的數(shù)字是其使用次數(shù),非葉子節(jié)點(diǎn)上的數(shù)字是其左右孩子使用次數(shù)之和):

          根
           |
      +---------33---------+
      |        |
   +----32---+      +----1---+
   |    |      |    |
+-21-+    +-11-+    +--1--+   
|   |    |   |    |   |
6   15  2  9    1   

(如果某個(gè)節(jié)點(diǎn)只有一個(gè)子節(jié)點(diǎn),可以去掉這個(gè)子節(jié)點(diǎn)。)

         根
         |
        +------33------+
       |     |
    +-----32----+     1
    |      |
  +--21--+  +--11--+
  |   |  |   |
  6   15 2    9

現(xiàn)在的編碼是: a : 000 b : 001 c : 010 d : 011 e : 1 仍然符合“前綴編碼”的要求。

第一步:如果發(fā)現(xiàn)下層節(jié)點(diǎn)的數(shù)字大于上層節(jié)點(diǎn)的數(shù)字,就交換它們的位置,并重新計(jì)算非葉子節(jié)點(diǎn)的值。
先交換11和1,由于11個(gè)字節(jié)縮短了一位,1個(gè)字節(jié)增長(zhǎng)了一位,總文件縮短了10位。

           根
            |
       +----------33---------+
       |        |
   +-----22----+     +----11----+
   |      |     |     |
+--21--+    1      2     9
|     |
6   15

再交換15和1、6和2,最終得到這樣的樹:

           根
            |
       +----------33---------+
       |        |
     +-----18----+    +----15----+
    |      |    |     |
  +--3--+    15   6     9
  |   |
  2   1

這時(shí)所有上層節(jié)點(diǎn)的數(shù)值都大于下層節(jié)點(diǎn)的數(shù)值,似乎無法再進(jìn)一步壓縮了。但是我們把每一層的最小的兩個(gè)節(jié)點(diǎn)結(jié)合起來,常會(huì)發(fā)現(xiàn)仍有壓縮余地。

第二步:把每一層的最小的兩個(gè)節(jié)點(diǎn)結(jié)合起來,重新計(jì)算相關(guān)節(jié)點(diǎn)的值。

在上面的樹中,第一、二、四三層都只有一或二個(gè)節(jié)點(diǎn),無法重新組合,但第三層上有四個(gè)節(jié)點(diǎn),我們把最小的3和6結(jié)合起來,并重新計(jì)算相關(guān)節(jié)點(diǎn)的值,成為下面這棵樹。

           根
            |
       +----------33---------+
       |         |
    +------9-----+    +----24----+
    |      |    |     |
   +--3--+    6   15    9
   |   |
  2  1

然后,再重復(fù)做第一步。
這時(shí)第二層的9小于第三層的15,于是可以互換,有9個(gè)字節(jié)增長(zhǎng)了一位,15個(gè)字節(jié)縮短了一位,文件總長(zhǎng)度又縮短了6位。然后重新計(jì)算相關(guān)節(jié)點(diǎn)的值。

           根
            |
       +----------33---------+
       |        |
       15     +----18----+ 
            |    |
         +------9-----+   9
         |      |
         +--3--+   6
         |   |
         2  1

這時(shí)發(fā)現(xiàn)所有的上層節(jié)點(diǎn)都大于下層節(jié)點(diǎn),每一層上最小的兩個(gè)節(jié)點(diǎn)被并在了一起,也不可能再產(chǎn)生比同層其他節(jié)點(diǎn)更小的父節(jié)點(diǎn)了。

這時(shí)整個(gè)文件的長(zhǎng)度是 3*6 + 1*15 + 4*2 + 2*9 + 4*1 = 63

這時(shí)可以看出編碼式壓縮的一個(gè)基本前提:各節(jié)點(diǎn)之間的值要相差比較懸殊,以使某兩個(gè)節(jié)點(diǎn)的和小于同層或下層的另一個(gè)節(jié)點(diǎn),這樣,交換節(jié)點(diǎn)才有利益。
所以歸根結(jié)底,原始文件中的字節(jié)使用頻率必須相差較大,否則將沒有兩個(gè)節(jié)點(diǎn)的頻率之和小于同層或下層其他節(jié)點(diǎn)的頻率,也就無法壓縮。反之,相差得越懸殊,兩個(gè)節(jié)點(diǎn)的頻率之和比同層或下層節(jié)點(diǎn)的頻率小得越多,交換節(jié)點(diǎn)之后的利益也越大。

在這個(gè)例子中,經(jīng)過上面兩步不斷重復(fù),得到了最優(yōu)的二叉樹,但不能保證在所有情況下,都能通過這兩步的重復(fù)得到最優(yōu)二叉樹,下面來看另一個(gè)例子:

                         根
                        ?。?
             ?。保梗?
              |                  ?。?
      +------12------+           ?。?
     ?。             。?
 ?。担     。罚?
 ?。      。     。      。?
+-2-+  ?。常 。常  。矗?
|  ?。  。  。 。  。  。  。?
1  ?。薄  。薄  。病 。薄  。病  。病  。?

這個(gè)例子中,所有上層節(jié)點(diǎn)都大于等于下層節(jié)點(diǎn),每一層最小的兩個(gè)節(jié)點(diǎn)結(jié)合在了一起,但仍然可以進(jìn)一步優(yōu)化:


                         根
                        ?。?
             ?。保梗?
             ?。                  。?
     ?。保玻           。?
     ?。             。?
 ?。矗     。福?
 ?。      。     。      。?
+-2-+  ?。玻 。矗  。矗?
|  ?。  。  。 。  。  。  。?
1  ?。薄  。薄  。薄 。病  。病  。病  。?

通過最低一層的第4第5個(gè)節(jié)點(diǎn)對(duì)換,第3層的8大于第2層的7。
到這里,我們得出這樣一個(gè)結(jié)論:一棵最優(yōu)二叉編碼樹(所有上層節(jié)點(diǎn)都無法和下層節(jié)點(diǎn)交換),必須符合這樣兩個(gè)條件:
1.所有上層節(jié)點(diǎn)都大于等于下層節(jié)點(diǎn)。
2.某節(jié)點(diǎn),設(shè)其較大的子節(jié)點(diǎn)為m,較小的子節(jié)點(diǎn)為n,m下的任一層的所有節(jié)點(diǎn)都應(yīng)大于等于n下的該層的所有節(jié)點(diǎn)。

當(dāng)符合這兩個(gè)條件時(shí),任一層都無法產(chǎn)生更小的節(jié)點(diǎn)去和下層節(jié)點(diǎn)交換,也無法產(chǎn)生更大的節(jié)點(diǎn)去和上層節(jié)點(diǎn)交換。

上面的兩個(gè)例子是比較簡(jiǎn)單的,實(shí)際的文件中,一個(gè)字節(jié)有256種可能的取值,所以二叉樹的葉子節(jié)點(diǎn)多達(dá)256個(gè),需要不斷的調(diào)整樹形,最終的樹形可能非常復(fù)雜,有一種非常精巧的算法可以快速地建起一棵最優(yōu)二叉樹,這種算法由D.Huffman(戴·霍夫曼)提出,下面我們先來介紹霍夫曼算法的步驟,然后再來證明通過這么簡(jiǎn)單的步驟得出的樹形確實(shí)是一棵最優(yōu)二叉樹。

霍夫曼算法的步驟是這樣的:

·從各個(gè)節(jié)點(diǎn)中找出最小的兩個(gè)節(jié)點(diǎn),給它們建一個(gè)父節(jié)點(diǎn),值為這兩個(gè)節(jié)點(diǎn)之和。
·然后從節(jié)點(diǎn)序列中去除這兩個(gè)節(jié)點(diǎn),加入它們的父節(jié)點(diǎn)到序列中。

重復(fù)上面兩個(gè)步驟,直到節(jié)點(diǎn)序列中只剩下唯一一個(gè)節(jié)點(diǎn)。這時(shí)一棵最優(yōu)二叉樹就已經(jīng)建成了,它的根就是剩下的這個(gè)節(jié)點(diǎn)。

仍以上面的例子來看霍夫曼樹的建立過程。
最初的節(jié)點(diǎn)序列是這樣的:
a(6)  b(15)  c(2)  d(9)  e(1)

把最小的c和e結(jié)合起來
                   | (3)
a(6)   b(15)   d(9)   +------+------+
              |      |
              c     e

不斷重復(fù),最終得到的樹是這樣的:

       根
        |
   +-----33-----+
   |     |
   15   +----18----+    
       |       |
       9  +------9-----+
          |      |
         6     +--3--+
              |   |
              2  1

這時(shí)各個(gè)字符的編碼長(zhǎng)度和前面我們說過的方法得到的編碼長(zhǎng)度是相同的,因而文件的總長(zhǎng)度也是相同的: 3*6 + 1*15 + 4*2 + 2*9 + 4*1 = 63

考察霍夫曼樹的建立過程中的每一步的節(jié)點(diǎn)序列的變化:

6  15 2 9 1
6  15 9 3
15 9  9
15 18
33

下面我們用逆推法來證明對(duì)于各種不同的節(jié)點(diǎn)序列,用霍夫曼算法建立起來的樹總是一棵最優(yōu)二叉樹:

對(duì)霍夫曼樹的建立過程運(yùn)用逆推法:
當(dāng)這個(gè)過程中的節(jié)點(diǎn)序列只有兩個(gè)節(jié)點(diǎn)時(shí)(比如前例中的15和18),肯定是一棵最優(yōu)二叉樹,一個(gè)編碼為0,另一個(gè)編碼為1,無法再進(jìn)一步優(yōu)化。
然后往前步進(jìn),節(jié)點(diǎn)序列中不斷地減少一個(gè)節(jié)點(diǎn),增加兩個(gè)節(jié)點(diǎn),在步進(jìn)過程中將始終保持是一棵最優(yōu)二叉樹,這是因?yàn)椋?
1.按照霍夫曼樹的建立過程,新增的兩個(gè)節(jié)點(diǎn)是當(dāng)前節(jié)點(diǎn)序列中最小的兩個(gè),其他的任何兩個(gè)節(jié)點(diǎn)的父節(jié)點(diǎn)都大于(或等于)這兩個(gè)節(jié)點(diǎn)的父節(jié)點(diǎn),只要前一步是最優(yōu)二叉樹,其他的任何兩個(gè)節(jié)點(diǎn)的父節(jié)點(diǎn)就一定都處在它們的父節(jié)點(diǎn)的上層或同層,所以這兩個(gè)節(jié)點(diǎn)一定處在當(dāng)前二叉樹的最低一層。
2.這兩個(gè)新增的節(jié)點(diǎn)是最小的,所以無法和其他上層節(jié)點(diǎn)對(duì)換。符合我們前面說的最優(yōu)二叉樹的第一個(gè)條件。
3.只要前一步是最優(yōu)二叉樹,由于這兩個(gè)新增的節(jié)點(diǎn)是最小的,即使同層有其他節(jié)點(diǎn),也無法和同層其他節(jié)點(diǎn)重新結(jié)合,產(chǎn)生比它們的父節(jié)點(diǎn)更小的上層節(jié)點(diǎn)來和同層的其他節(jié)點(diǎn)對(duì)換。它們的父節(jié)點(diǎn)小于其他節(jié)點(diǎn)的父節(jié)點(diǎn),它們又小于其他所有節(jié)點(diǎn),只要前一步符合最優(yōu)二叉樹的第二個(gè)條件,到這一步仍將符合。

這樣一步步逆推下去,在這個(gè)過程中霍夫曼樹每一步都始終保持著是一棵最優(yōu)二叉樹。

由于每一步都從節(jié)點(diǎn)序列中刪除兩個(gè)節(jié)點(diǎn),新增一個(gè)節(jié)點(diǎn),霍夫曼樹的建立過程共需 (原始節(jié)點(diǎn)數(shù) - 1) 步,所以霍夫曼算法不失為一種精巧的編碼式壓縮算法。


附:對(duì)于 huffman 樹,《計(jì)算機(jī)程序設(shè)計(jì)藝術(shù)》中有完全不同的證明,大意是這樣的:
1.二叉編碼樹的內(nèi)部節(jié)點(diǎn)(非葉子節(jié)點(diǎn))數(shù)等于外部節(jié)點(diǎn)(葉子節(jié)點(diǎn))數(shù)減1。
2.二叉編碼樹的外部節(jié)點(diǎn)的加權(quán)路徑長(zhǎng)度(值乘以路徑長(zhǎng)度)之和,等于所有內(nèi)部節(jié)點(diǎn)值之和。(這兩條都可以通過對(duì)節(jié)點(diǎn)數(shù)運(yùn)用數(shù)學(xué)歸納法來證明,留給大家做練習(xí)。)
3.對(duì) huffman 樹的建立過程運(yùn)用逆推,當(dāng)只有一個(gè)內(nèi)部節(jié)點(diǎn)時(shí),肯定是一棵最優(yōu)二叉樹。
4.往前步進(jìn),新增兩個(gè)最小的外部節(jié)點(diǎn),它們結(jié)合在一起產(chǎn)生一個(gè)新的內(nèi)部節(jié)點(diǎn),當(dāng)且僅當(dāng)原先的內(nèi)部節(jié)點(diǎn)集合是極小化的,加入這個(gè)新的內(nèi)部節(jié)點(diǎn)后仍是極小化的。(因?yàn)樽钚〉膬蓚€(gè)節(jié)點(diǎn)結(jié)合在一起,并處于最低層,相對(duì)于它們分別和其他同層或上層節(jié)點(diǎn)結(jié)合在一起,至少不會(huì)增加加權(quán)路徑長(zhǎng)度。)
5.隨著內(nèi)部節(jié)點(diǎn)數(shù)逐個(gè)增加,內(nèi)部節(jié)點(diǎn)集合總維持極小化。
2.實(shí)現(xiàn)部分
  如果世界上從沒有一個(gè)壓縮程序,我們看了前面的壓縮原理,將有信心一定能作出一個(gè)可以壓縮大多數(shù)格式、內(nèi)容的數(shù)據(jù)的程序,當(dāng)我們著手要做這樣一個(gè)程序的時(shí)候,會(huì)發(fā)現(xiàn)有很多的難題需要我們?nèi)ヒ粋€(gè)個(gè)解決,下面將逐個(gè)描述這些難題,并詳細(xì)分析 zip 算法是如何解決這些難題的,其中很多問題帶有普遍意義,比如查找匹配,比如數(shù)組排序等等,這些都是說不盡的話題,讓我們深入其中,做一番思考。

我們前面說過,對(duì)于短語式重復(fù),我們用“重復(fù)距當(dāng)前位置的距離”和“重復(fù)的長(zhǎng)度”這兩個(gè)數(shù)字來表示這一段重復(fù),以實(shí)現(xiàn)壓縮,現(xiàn)在問題來了,一個(gè)字節(jié)能表示的數(shù)字大小為 0 -255,然而重復(fù)出現(xiàn)的位置和重復(fù)的長(zhǎng)度都可能超過 255,事實(shí)上,二進(jìn)制數(shù)的位數(shù)確定下來后,所能表示的數(shù)字大小的范圍是有限的,n位的二進(jìn)制數(shù)能表示的最大值是2的n次方減1,如果位數(shù)取得太大,對(duì)于大量的短匹配,可能不但起不到壓縮作用,反而增大了最終的結(jié)果。針對(duì)這種情況,有兩種不同的算法來解決這個(gè)問題,它們是兩種不同的思路。一種稱為 lz77 算法,這是一種很自然的思路:限制這兩個(gè)數(shù)字的大小,以取得折衷的壓縮效果。例如距離取 15 位,長(zhǎng)度取 8 位,這樣,距離的最大取值為 32 k - 1,長(zhǎng)度的最大取值為 255,這兩個(gè)數(shù)字占 23 位,比三個(gè)字節(jié)少一位,是符合壓縮的要求的。讓我們?cè)陬^腦中想象一下 lz77 算法壓縮進(jìn)行時(shí)的情況,會(huì)出現(xiàn)有意思的模型:

   最遠(yuǎn)匹配位置->          當(dāng)前處理位置->
───┸─────────────────╂─────────────>壓縮進(jìn)行方向
   已壓縮部分             ┃    未壓縮部分

  在最遠(yuǎn)匹配位置和當(dāng)前處理位置之間是可以用來查找匹配的“字典”區(qū)域,隨著壓縮的進(jìn)行,“字典”區(qū)域從待壓縮文件的頭部不斷地向后滑動(dòng),直到達(dá)到文件的尾部,短語式壓縮也就結(jié)束了。
  解壓縮也非常簡(jiǎn)單:

         ┎────────拷貝────────┒
 匹配位置    ┃          當(dāng)前處理位置  ┃
   ┃<──匹配長(zhǎng)度──>┃       ┠─────∨────┨
───┸──────────┸───────╂──────────┸─>解壓進(jìn)行方向
   已解壓部分              ┃    未解壓部分

  不斷地從壓縮文件中讀出匹配位置值和匹配長(zhǎng)度值,把已解壓部分的匹配內(nèi)容拷貝到解壓文件尾部,遇到壓縮文件中那些壓縮時(shí)未能得到匹配,而是直接保存的單、雙字節(jié),解壓時(shí)只要直接拷貝到文件尾部即可,直到整個(gè)壓縮文件處理完畢。
  lz77算法模型也被稱為“滑動(dòng)字典”模型或“滑動(dòng)窗口”模型。
  另有一種lzw算法對(duì)待壓縮文件中存在大量簡(jiǎn)單匹配的情況進(jìn)行了完全不同的算法設(shè)計(jì),它只用一個(gè)數(shù)字來表示一段短語,下面來描述一下lzw的壓縮解壓過程,然后來綜合比較兩者的適用情況。
  lzw的壓縮過程:
1) 初始化一個(gè)指定大小的字典,把 256 種字節(jié)取值加入字典。
2) 在待壓縮文件的當(dāng)前處理位置尋找在字典中出現(xiàn)的最長(zhǎng)匹配,輸出該匹配在字典中的序號(hào)。
3) 如果字典沒有達(dá)到最大容量,把該匹配加上它在待壓縮文件中的下一個(gè)字節(jié)加入字典。
4) 把當(dāng)前處理位置移到該匹配后。
5) 重復(fù) 2、3、4 直到文件輸出完畢。

  lzw 的解壓過程:
1) 初始化一個(gè)指定大小的字典,把 256 種字節(jié)取值加入字典。
2) 從壓縮文件中順序讀出一個(gè)字典序號(hào),根據(jù)該序號(hào),把字典中相應(yīng)的數(shù)據(jù)拷貝到解壓文件尾部。
3) 如果字典沒有達(dá)到最大容量,把前一個(gè)匹配內(nèi)容加上當(dāng)前匹配的第一個(gè)字節(jié)加入字典。
4) 重復(fù) 2、3 兩步直到壓縮文件處理完畢。

  從 lzw 的壓縮過程,我們可以歸納出它不同于 lz77 算法的一些主要特點(diǎn):
1) 對(duì)于一段短語,它只輸出一個(gè)數(shù)字,即字典中的序號(hào)。(這個(gè)數(shù)字的位數(shù)決定了字典的最大容量,當(dāng)它的位數(shù)取得太大時(shí),比如 24 位以上,對(duì)于短匹配占多數(shù)的情況,壓縮率可能很低。取得太小時(shí),比如 8 位,字典的容量受到限制。所以同樣需要取舍。)
2) 對(duì)于一個(gè)短語,比如 abcd ,當(dāng)它在待壓縮文件中第一次出現(xiàn)時(shí),ab 被加入字典,第二次出現(xiàn)時(shí),abc 被加入字典,第三次出現(xiàn)時(shí),abcd 才會(huì)被加入字典,對(duì)于一些長(zhǎng)匹配,它必須高頻率地出現(xiàn),并且字典有較大的容量,才會(huì)被最終完整地加入字典。相應(yīng)地,lz77 只要匹配在“字典區(qū)域”中存在,馬上就可以直接使用。
3) 設(shè) lzw 的“字典序號(hào)”取 n 位,它的最大長(zhǎng)度可以達(dá)到 2 的 n 次方;設(shè) lz77 的“匹配長(zhǎng)度”取 n 位,“匹配距離”取 d 位,它的最大長(zhǎng)度也是 2 的 n 次方,但還要多輸出 d 位(d 至少不小于 n),從理論上說 lzw 每輸出一個(gè)匹配只要 n 位,不管是長(zhǎng)匹配還是短匹配,壓縮率要比 lz77 高至少一倍,但實(shí)際上,lzw 的字典中的匹配長(zhǎng)度的增長(zhǎng)由于各匹配互相打斷,很難達(dá)到最大值。而且雖然 lz77 每一個(gè)匹配都要多輸出 d 位,但 lzw 每一個(gè)匹配都要從單字節(jié)開始增長(zhǎng)起,對(duì)于種類繁多的匹配,lzw 居于劣勢(shì)。
  可以看出,在多數(shù)情況下,lz77 擁有更高的壓縮率,而在待壓縮文件中占絕大多數(shù)的是些簡(jiǎn)單的匹配時(shí),lzw 更具優(yōu)勢(shì),GIF 就是采用了 lzw 算法來壓縮背景單一、圖形簡(jiǎn)單的圖片。zip 是用來壓縮通用文件的,這就是它采用對(duì)大多數(shù)文件有更高壓縮率的 lz77 算法的原因。

  接下來 zip 算法將要解決在“字典區(qū)域”中如何高速查找最長(zhǎng)匹配的問題。

(注:以下關(guān)于技術(shù)細(xì)節(jié)的描述是以 gzip 的公開源代碼為基礎(chǔ)的,如果需要完整的代碼,可以在 gzip 的官方網(wǎng)站 www.gzip.org 下載。下面提到的每一個(gè)問題,都首先介紹最直觀簡(jiǎn)單的解決方法,然后指出這種方法的弊端所在,最后介紹 gzip 采用的做法,這樣也許能使讀者對(duì) gzip 看似復(fù)雜、不直觀的做法的意義有更好的理解。)
最直觀的搜索方式是順序搜索:以待壓縮部分的第一個(gè)字節(jié)與窗口中的每一個(gè)字節(jié)依次比較,當(dāng)找到一個(gè)相等的字節(jié)時(shí),再比較后續(xù)的字節(jié)…… 遍歷了窗口后得出最長(zhǎng)匹配。gzip 用的是被稱作“哈希表”的方法來實(shí)現(xiàn)較高效的搜索?!肮#╤ash)”是分散的意思,把待搜索的數(shù)據(jù)按照字節(jié)值分散到一個(gè)個(gè)“桶”中,搜索時(shí)再根據(jù)字節(jié)值到相應(yīng)的“桶”中去尋找。短語式壓縮的最短匹配為 3 個(gè)字節(jié),gzip 以 3 個(gè)字節(jié)的值作為哈希表的索引,但 3 個(gè)字節(jié)共有 2 的 24 次方種取值,需要 16M 個(gè)桶,桶里存放的是窗口中的位置值,窗口的大小為 32K,所以每個(gè)桶至少要有大于兩個(gè)字節(jié)的空間,哈希表將大于 32M,作為 90 年代開發(fā)的程序,這個(gè)要求是太大了,而且隨著窗口的移動(dòng),哈希表里的數(shù)據(jù)會(huì)不斷過時(shí),維護(hù)這么大的表,會(huì)降低程序的效率,gzip 定義哈希表為 2 的 15 次方(32K)個(gè)桶,并設(shè)計(jì)了一個(gè)哈希函數(shù)把 16M 種取值對(duì)應(yīng)到 32K 個(gè)桶中,不同的值被對(duì)應(yīng)到相同的桶中是不可避免的,哈希函數(shù)的任務(wù)是 1.使各種取值盡可能均勻地分布到各個(gè)桶中,避免許多不同的值集中到某些桶中,而另一些是空桶,使搜索的效率降低。2.函數(shù)的計(jì)算盡可能地簡(jiǎn)單,因?yàn)槊看巍安迦搿焙汀八褜ぁ惫1矶家獔?zhí)行哈希函數(shù),哈希函數(shù)的復(fù)雜度直接影響程序的執(zhí)行效率,容易想到的哈希函數(shù)是取 3 個(gè)字節(jié)的左邊(或右邊)15 位二進(jìn)制值,但這樣只要左邊(或右邊)2 個(gè)字節(jié)相同,就會(huì)被放到同一個(gè)桶中,而 2 個(gè)字節(jié)相同的概率是比較高的,不符合“平均分布”的要求。gzip 采用的算法是:A(4,5) + A(6,7,8) ^ B(1,2,3) + B(4,5) + B(6,7,8) ^ C(1,2,3) + C(4,5,6,7,8) (說明:A 指 3 個(gè)字節(jié)中的第 1 個(gè)字節(jié),B 指第 2 個(gè)字節(jié),C 指第 3 個(gè)字節(jié),A(4,5) 指第一個(gè)字節(jié)的第 4,5 位二進(jìn)制碼,“^”是二進(jìn)制位的異或操作,“+”是“連接”而不是“加”,“^”優(yōu)先于“+”)這樣使 3 個(gè)字節(jié)都盡量“參與”到最后的結(jié)果中來,而且每個(gè)結(jié)果值 h 都等于 ((前1個(gè)h << 5) ^ c)取右 15 位,計(jì)算也還簡(jiǎn)單。 
哈希表的具體實(shí)現(xiàn)也值得探討,因?yàn)闊o法預(yù)先知道每一個(gè)“桶”會(huì)存放多少個(gè)元素,所以最簡(jiǎn)單的,會(huì)想到用鏈表來實(shí)現(xiàn):哈希表里存放著每個(gè)桶的第一個(gè)元素,每個(gè)元素除了存放著自身的值,還存放著一個(gè)指針,指向同一個(gè)桶中的下一個(gè)元素,可以順著指針鏈來遍歷該桶中的每一個(gè)元素,插入元素時(shí),先用哈希函數(shù)算出該放到第幾個(gè)桶中,再把它掛到相應(yīng)鏈表的最后。這個(gè)方案的缺點(diǎn)是頻繁地申請(qǐng)和釋放內(nèi)存會(huì)降低運(yùn)行速度;內(nèi)存指針的存放占據(jù)了額外的內(nèi)存開銷。有更少內(nèi)存開銷和更快速的方法來實(shí)現(xiàn)哈希表,并且不需要頻繁的內(nèi)存申請(qǐng)和釋放:gzip 在內(nèi)存中申請(qǐng)了兩個(gè)數(shù)組,一個(gè)叫 head[],一個(gè)叫 pre[],大小都為 32K,根據(jù)當(dāng)前位置 strstart 開始的 3 個(gè)字節(jié),用哈希函數(shù)計(jì)算出在 head[] 中的位置 ins_h,然后把 head[ins_h] 中的值記入 pre[strstart],再把當(dāng)前位置 strstart 記入 head[ins_h]。隨著壓縮的進(jìn)行,head[]里記載著最近的可能的匹配的位置(如果有匹配的話,head[ins_h]不為 0),pre[]中的所有位置與原始數(shù)據(jù)的位置相對(duì)應(yīng),但每一個(gè)位置保存的值是前一個(gè)最近的可能的匹配的位置。(“可能的匹配”是指哈希函數(shù)計(jì)算出的 ins_h 相同。)順著 pre[] 中的指示找下去,直到遇到 0,可以得到所有匹配在原始數(shù)據(jù)中的位置,0 表示不再有更遠(yuǎn)的匹配。
  接下來很自然地要觀察 gzip 具體是如何判斷哈希表中數(shù)據(jù)的過時(shí),如何清理哈希表的,因?yàn)?nbsp;pre[] 里只能存放 32K 個(gè)元素,所以這項(xiàng)工作是必須要做的。
  gzip 從原始文件中讀出兩個(gè)窗口大小的內(nèi)容(共 64K 字節(jié))到一塊內(nèi)存中,這塊內(nèi)存也是一個(gè)數(shù)組,稱作 Window[];申請(qǐng) head[]、pre[] 并清零;strstart 置為 0。然后 gzip 邊搜索邊插入,搜索時(shí)通過計(jì)算 ins_h,檢查 head[] 中是否有匹配,如果有匹配,判斷 strstart 減 head[] 中的位置是否大于 1 個(gè)窗口的大小,如果大于 1 個(gè)窗口的大小,就不到 pre[] 中去搜索了,因?yàn)?nbsp;pre[] 中保存的位置更遠(yuǎn)了,如果不大于,就順著 pre[] 的指示到 Window[] 中逐個(gè)匹配位置開始,逐個(gè)字節(jié)與當(dāng)前位置的數(shù)據(jù)比較,以找出最長(zhǎng)匹配,pre[] 中的位置也要判斷是否超出一個(gè)窗口,如遇到超出一個(gè)窗口的位置或者 0 就不再找下去,找不到匹配就輸出當(dāng)前位置的單個(gè)字節(jié)到另外的內(nèi)存(輸出方法在后文中會(huì)介紹),并把 strstart 插入哈希表,strstart 遞增,如果找到了匹配,就輸出匹配位置和匹配長(zhǎng)度這兩個(gè)數(shù)字到另外的內(nèi)存中,并把 strstart 開始的,直到 strstart + 匹配長(zhǎng)度 為止的所有位置都插入哈希表,strstart += 匹配長(zhǎng)度。插入哈希表的方法為:
pre[strstart % 32K] = head[ins_h];
head[ins_h] = strstart;
可以看出,pre[] 是循環(huán)利用的,所有的位置都在一個(gè)窗口以內(nèi),但每一個(gè)位置保存的值不一定是一個(gè)窗口以內(nèi)的。在搜索時(shí),head[] 和 pre[] 中的位置值對(duì)應(yīng)到 pre[] 時(shí)也要 % 32K。當(dāng) Window[] 中的原始數(shù)據(jù)將要處理完畢時(shí),要把 Window[] 中后一窗的數(shù)據(jù)復(fù)制到前一窗,再讀取 32K 字節(jié)的數(shù)據(jù)到后一窗,strstart -= 32K,遍歷 head[],值小于等于 32K 的,置為 0,大于 32K 的,-= 32K;pre[] 同 head[] 一樣處理。然后同前面一樣處理新一窗的數(shù)據(jù)。
  分析:現(xiàn)在可以看到,雖然 3 個(gè)字節(jié)有 16M 種取值,但實(shí)際上一個(gè)窗口只有 32K 個(gè)取值需要插入哈希表,由于短語式重復(fù)的存在,實(shí)際只有 < 32K 種取值插入哈希表的 32K 個(gè)“桶”中,而且哈希函數(shù)又符合“平均分布”的要求,所以哈希表中實(shí)際存在的“沖突”一般不會(huì)多,對(duì)搜索效率的影響不大。可以預(yù)計(jì),在“一般情況”下,每個(gè)“桶”中存放的數(shù)據(jù),正是我們要找的。哈希表在各種搜索算法中,實(shí)現(xiàn)相對(duì)的比較簡(jiǎn)單,容易理解,“平均搜索速度”最快,哈希函數(shù)的設(shè)計(jì)是搜索速度的關(guān)鍵,只要符合“平均分布”和“計(jì)算簡(jiǎn)單”,就常常能成為諸種搜索算法中的首選,所以哈希表是最流行的一種搜索算法。但在某些特殊情況下,它也有缺點(diǎn),比如:1.當(dāng)鍵碼 k 不存在時(shí),要求找出小于 k 的最大鍵碼或大于 k 的最小鍵碼,哈希表無法有效率地滿足這種要求。2.哈希表的“平均搜索速度”是建立在概率論的基礎(chǔ)上的,因?yàn)槭孪炔荒茴A(yù)知待搜索的數(shù)據(jù)集合,我們只能“信賴”搜索速度的“平均值”,而不能“保證”搜索速度的“上限”。在同人類性命攸關(guān)的應(yīng)用中(如醫(yī)療或宇航領(lǐng)域),將是不合適的。這些情況及其他一些特殊情況下,我們必須求助其他“平均速度”較低,但能滿足相應(yīng)的特殊要求的算法。(見《計(jì)算機(jī)程序設(shè)計(jì)藝術(shù)》第3卷 排序與查找)。幸而“在窗口中搜索匹配字節(jié)串”不屬于特殊情況。

時(shí)間與壓縮率的平衡:
gzip 定義了幾種可供選擇的 level,越低的 level 壓縮時(shí)間越快但壓縮率越低,越高的 level 壓縮時(shí)間越慢但壓縮率越高。
不同的 level 對(duì)下面四個(gè)變量有不同的取值:

nice_length
max_chain
max_lazy
good_length

nice_length:前面說過,搜索匹配時(shí),順著 pre[] 的指示到 Window[] 中逐個(gè)匹配位置開始,找出最長(zhǎng)匹配,但在這過程中,如果遇到一個(gè)匹配的長(zhǎng)度達(dá)到或超過 nice_length,就不再試圖尋找更長(zhǎng)的匹配。最低的 level 定義 nice_length 為 8,最高的 level 定義 nice_length 為 258(即一個(gè)字節(jié)能表示的最大短語匹配長(zhǎng)度 3 + 255)。

max_chain:這個(gè)值規(guī)定了順著 pre[] 的指示往前回溯的最大次數(shù)。最低的 level 定義 max_chain 為 4,最高的 level 定義 max_chain 為 4096。當(dāng) max_chain 和 nice_length 有沖突時(shí),以先達(dá)到的為準(zhǔn)。

max_lazy:這里有一個(gè)懶惰匹配(lazy match)的概念,在輸出當(dāng)前位置(strstart)的匹配之前,gzip 會(huì)去找下一個(gè)位置(strstart + 1)的匹配,如果下一個(gè)匹配的長(zhǎng)度比當(dāng)前匹配的長(zhǎng)度更長(zhǎng),gzip 就放棄當(dāng)前匹配,只輸出當(dāng)前位置處的首個(gè)字節(jié),然后再查找 strstart + 2 處的匹配,這樣的方式一直往后找,如果后一個(gè)匹配比前一個(gè)匹配更長(zhǎng),就只輸出前一個(gè)匹配的首字節(jié),直到遇到前一個(gè)匹配長(zhǎng)于后一個(gè)匹配,才輸出前一個(gè)匹配。
gzip 作者的思路是,如果后一個(gè)匹配比前一個(gè)匹配更長(zhǎng),就犧牲前一個(gè)匹配的首字節(jié)來換取后面的大于等于1的額外的匹配長(zhǎng)度。
max_lazy 規(guī)定了,如果匹配的長(zhǎng)度達(dá)到或超過了這個(gè)值,就直接輸出,不再管后一個(gè)匹配是否更長(zhǎng)。最低的4級(jí) level 不做懶惰匹配,第5級(jí) level 定義 max_lazy 為 4,最高的 level 定義 max_lazy 為 258。

good_length:這個(gè)值也和懶惰匹配有關(guān),如果前一個(gè)匹配長(zhǎng)度達(dá)到或超過 good_length,那在尋找當(dāng)前的懶惰匹配時(shí),回溯的最大次數(shù)減小到 max_chain 的 1/4,以減少當(dāng)前的懶惰匹配花費(fèi)的時(shí)間。第5級(jí) level 定義 good_length 為 4(這一級(jí)等于忽略了 good_length),最高的 level 定義 good_length 為 32。

分析:懶惰匹配有必要嗎?可以改進(jìn)嗎?
gzip 的作者是無損壓縮方面的專家,但是世界上沒有絕對(duì)的權(quán)威,吾愛吾師,更愛真理。我覺得 gzip 的作者對(duì)懶惰匹配的考慮確實(shí)不夠周詳。只要是進(jìn)行了認(rèn)真客觀的分析,誰都有權(quán)利提出自己的觀點(diǎn)。
采用懶惰匹配,需要對(duì)原始文件的更多的位置查找匹配,時(shí)間肯定增加了許多倍,但壓縮率的提高在總體上十分有限。在幾種情況下,它反而增長(zhǎng)了短語壓縮的結(jié)果,所以如果一定要用懶惰匹配,也應(yīng)該改進(jìn)一下算法,下面是具體的分析。
1. 連續(xù)3次以上找到了更長(zhǎng)的匹配,就不應(yīng)該單個(gè)輸出前面的那些字節(jié),而應(yīng)該作為匹配輸出。
2. 于是,如果連續(xù)找到更長(zhǎng)的匹配的次數(shù)大于第一個(gè)匹配的長(zhǎng)度,對(duì)于第一個(gè)匹配,相當(dāng)于沒有做懶惰匹配。
3. 如果小于第一個(gè)匹配的長(zhǎng)度但大于2,就沒有必要作懶惰匹配,因?yàn)檩敵龅目偸莾蓚€(gè)匹配。
4. 所以找到一個(gè)匹配后,最多只需要向后做 2 次懶惰匹配,就可以決定是輸出第一個(gè)匹配,還是輸出1(或 2)個(gè)首字節(jié)加后面的匹配了。
5. 于是,對(duì)于一段原始字節(jié)串,如果不做懶惰匹配時(shí)輸出兩個(gè)匹配(對(duì)于每個(gè)匹配,距離占15位二進(jìn)制數(shù),長(zhǎng)度占8位二進(jìn)制數(shù),加起來約占3字節(jié),輸出兩個(gè)匹配約需要6字節(jié)),做了懶惰匹配如果有改進(jìn)的話,將是輸出1或2個(gè)單字節(jié)加上1個(gè)匹配(也就是約4或5字節(jié))。這樣,懶惰匹配可以使某些短語壓縮的結(jié)果再縮短1/3到1/6。
6. 再觀察這樣一個(gè)例子:
1232345145678[當(dāng)前位置]12345678
不用懶惰匹配,約輸出6字節(jié),用懶惰匹配,約輸出7字節(jié),由于使用了懶惰匹配,把更后面的一個(gè)匹配拆成了兩個(gè)匹配。(如果 678 正好能歸入再后面的一個(gè)匹配,那懶惰匹配可能是有益的。)
7. 綜合考慮各種因素(匹配數(shù)和未匹配的單雙字節(jié)在原始文件中所占的比例,后一個(gè)匹配長(zhǎng)度大于前一個(gè)匹配長(zhǎng)度的概率,等等),經(jīng)過改進(jìn)的懶惰匹配算法,對(duì)總的壓縮率即使有貢獻(xiàn),也仍是很小的,而且也仍然很有可能會(huì)降低壓縮率。再考慮到時(shí)間的確定的明顯的增加與壓縮率的不確定的微弱的增益,也許最好的改進(jìn)是果斷地放棄懶惰匹配。
gzip 在完成短語式壓縮后,將轉(zhuǎn)入編碼式壓縮的階段。這個(gè)階段的實(shí)現(xiàn)是很復(fù)雜的,對(duì)最終的壓縮率至關(guān)重要,我會(huì)詳細(xì)解說 gzip 的做法。gzip 是開放源代碼的無損壓縮程序中最著名的,其中的種種技巧很有啟發(fā)意義,但是他是比較早期的程序,現(xiàn)在有很多的程序已經(jīng)在壓縮率上超過了它,所以我會(huì)根據(jù)自己對(duì)無損壓縮的基本規(guī)律的理解提出對(duì)它的改進(jìn)。

編碼式壓縮的幾點(diǎn)考慮:
1. huffman 算法壓縮率的關(guān)鍵是各節(jié)點(diǎn)值的差異要大,這樣就要求分段編碼輸出。因?yàn)槟承┒温渲心承┕?jié)點(diǎn)的出現(xiàn)頻率較高,另一些段落中這些節(jié)點(diǎn)出現(xiàn)頻率較低,如果不分段輸出,頻率的差異會(huì)被彼此抵消,而不同段落中,節(jié)點(diǎn)的出現(xiàn)頻率不同是常有的。
  要決定分段的大小,必須解決一對(duì)矛盾:上面的分析似乎要求段落越小越好,但由于要保存碼表以對(duì) huffman 壓縮結(jié)果進(jìn)行解壓,每個(gè)段落都要保存一份不同的碼表,所以段落取得太小,保存了碼表后得不償失,這樣,似乎又要求段落要盡量大,使碼表的保存份數(shù)盡量少。
  gzip 采取了這樣的策略來確定段落的大小:lz77 壓縮每產(chǎn)生 4k(?。┑臄?shù)據(jù),就判斷現(xiàn)在對(duì)未編碼部分進(jìn)行編碼輸出是否合適,最多積壓到 32k(大)的時(shí)候,必定進(jìn)行強(qiáng)制輸出,因?yàn)槠接沟臄?shù)據(jù)積壓得太多,后面即使有好的數(shù)據(jù),頻率統(tǒng)計(jì)在一起,也會(huì)被平庸化。
  判斷當(dāng)前輸出合適與否的條件是:1)用預(yù)先設(shè)定好的各節(jié)點(diǎn)長(zhǎng)度和各節(jié)點(diǎn)實(shí)際的出現(xiàn)次數(shù),計(jì)算壓縮結(jié)果的大概值,看這個(gè)值是否小于未壓縮時(shí)的 1/2。2)看目前為止的匹配數(shù)是否小于未匹配的字節(jié)數(shù),因?yàn)?nbsp;lz77 壓縮產(chǎn)生的數(shù)據(jù)包括“匹配”和“未匹配的原始字節(jié)”,段落間的節(jié)點(diǎn)頻率差異主要體現(xiàn)在“未匹配的原始字節(jié)”中。
  上面的判斷只是一種“猜測(cè)”,真正的精確的計(jì)算需要花費(fèi)更多的時(shí)間。
  我覺得 gzip 的策略可以改進(jìn),我的策略是:1)輸出的時(shí)機(jī)是壓縮率的關(guān)鍵之一,現(xiàn)在計(jì)算機(jī)的速度和九十年代時(shí)已經(jīng)今非昔比,現(xiàn)在完全有條件采用真正的建 huffman 樹的方法得到各節(jié)點(diǎn)的碼長(zhǎng),作精確的判斷。2)不應(yīng)該與未壓縮的原始數(shù)據(jù)比較,而應(yīng)該與 lz77 輸出的數(shù)據(jù)比較,否則計(jì)算出的壓縮比很大一部分是短語式壓縮的功勞。3)由于采用了真正的建 huffman 樹的方法,不用再去做匹配數(shù)與未匹配的字節(jié)數(shù)的比較,因?yàn)槟侵皇且环N猜測(cè)。4)每 4k 的數(shù)據(jù)都單獨(dú)統(tǒng)計(jì)頻率,如果是合適的,就先輸出之前的積壓(如果有的話),再輸出當(dāng)前的 4k,這樣可以避免當(dāng)前的數(shù)據(jù)被積壓的數(shù)據(jù)平庸化。如果不合適,就把當(dāng)前的頻率歸入到積壓的數(shù)據(jù)(如果有)的頻率中,再判斷是否合適,如仍不合適就暫緩輸出,否則一起輸出,這和 gzip 的作法是一樣的。說明:幾段差的數(shù)據(jù)積壓到一起仍有可能成為好的數(shù)據(jù),比如 01、 02、……積壓在一起,0 的頻率逐漸高出了其他字節(jié)。5)如果愿意付出更多的時(shí)間,在把當(dāng)前的頻率歸入之前的頻率時(shí),可以先和之前 4k 的頻率合并,如果不合適,和之前 8k 的頻率合并,這樣逐漸往前合并 4k,避免前面不好的數(shù)據(jù)拖累合并后的好的數(shù)據(jù)。6)有了前面的機(jī)制,32k 的強(qiáng)制輸出點(diǎn)可以取消。7)進(jìn)一步的改進(jìn):當(dāng)要輸出時(shí),只輸出積壓的不好的部分,好的數(shù)據(jù)先留著,等后面的 4k,如果新的加入后,仍是好的數(shù)據(jù),就再等,如果會(huì)降低壓縮率,才輸出好的部分。這樣,讓好的數(shù)據(jù)大段的輸出,可以減少碼表的保存份數(shù)。8)再進(jìn)一步的改進(jìn):壞的數(shù)據(jù)放在一起可能會(huì)提高壓縮率,好的數(shù)據(jù)放在一起也可能更好,當(dāng)然,兩種情況也都有可能降低壓縮率,所以前面判斷“好”還是“不好”,“合適”還是“不合適”的標(biāo)準(zhǔn)應(yīng)該從某一個(gè)固定的壓縮率標(biāo)準(zhǔn)改變?yōu)椋禾岣吡藟嚎s率還是降低了壓縮率。(提高的幅度應(yīng)該至少抵消多保存一份碼表的損失;降低的幅度也應(yīng)該至少抵消少保存一份碼表的得益)9)綜合前面的分析,確定分段大小的策略最終調(diào)整為:當(dāng)新的數(shù)據(jù)和前面的未切分?jǐn)?shù)據(jù)放在一起時(shí),兩者中任何一方受到損失,都應(yīng)該設(shè)置切分點(diǎn),積累了兩個(gè)分段后,通過計(jì)算,當(dāng)切分帶來的收益大于少保存一份碼表時(shí),才輸出前一段,否則取消它們之間的切分點(diǎn)。這個(gè)策略實(shí)際上可以涵蓋前面提到的所有改進(jìn),因?yàn)槊總€(gè)實(shí)際的分段之中的數(shù)據(jù)或者相互促進(jìn),或者彼此稍有妨害,但好過多保存一份碼表;而每?jī)蓚€(gè)相鄰的分段之間的數(shù)據(jù)彼此妨害,抵消了少保存一份碼表的收益。這個(gè)策略簡(jiǎn)單直觀地體現(xiàn)了我們?cè)O(shè)置分段的初衷:就是分段輸出必須能提高壓縮率。

2. 如果不考慮碼表,huffman 算法能得到最短的編碼式壓縮結(jié)果,但是這種算法必須保存碼表以便解壓縮,所以不能保證結(jié)果是最佳的。gzip 預(yù)先擬定了一套通用的靜態(tài)的編碼,當(dāng)要輸出一個(gè)段落時(shí),比較 huffman 壓縮結(jié)果加碼表的長(zhǎng)度和靜態(tài)編碼的壓縮結(jié)果長(zhǎng)度,再?zèng)Q定用哪種方法輸出這個(gè)段落。靜態(tài)編碼不需要建樹,計(jì)算壓縮結(jié)果長(zhǎng)度時(shí)耗時(shí)很少。如果各節(jié)點(diǎn)的頻率的差異很小,huffman 壓縮結(jié)果加碼表反而增大了結(jié)果,靜態(tài)編碼也不合適,同樣增大了結(jié)果,gzip 就直接保存 lz77 的原始輸出。由于輸出一個(gè)段落時(shí),增加了靜態(tài)編碼的方案,使輸出的實(shí)際長(zhǎng)度和之前確定分段點(diǎn)時(shí)計(jì)算的值可能不同,那么前面計(jì)算出的這個(gè)分段點(diǎn)是否仍是正確的?前面的分段策略是否需要調(diào)整?
  分析:1)靜態(tài)編碼的各節(jié)點(diǎn)編碼是不變的,對(duì)于段落的合并是無所謂的,兩個(gè)連續(xù)段落即使都采用靜態(tài)編碼,也不用合并,因?yàn)楹喜⒑蠼Y(jié)果長(zhǎng)度是不會(huì)變的。2)所以只對(duì)一種情況可能有影響:一個(gè)段落中拆分出一些部分用 huffman 編碼,另一些部分用靜態(tài)編碼,壓縮結(jié)果更好。當(dāng)這種情況發(fā)生時(shí),則必有一些部分的優(yōu)勢(shì)節(jié)點(diǎn)(頻率高的節(jié)點(diǎn))與靜態(tài)編碼預(yù)先擬定的優(yōu)勢(shì)節(jié)點(diǎn)相近,采用靜態(tài)編碼后有稍許改善,其他部分則與靜態(tài)編碼預(yù)先擬定的優(yōu)勢(shì)節(jié)點(diǎn)有一定分歧,采用靜態(tài)編碼后會(huì)有稍許不利。之所以說“稍許”,是因?yàn)槲覀円阎粋€(gè)段落里的各部分?jǐn)?shù)據(jù)或者互相促進(jìn),或者僅有稍許妨害,說明它們的優(yōu)勢(shì)節(jié)點(diǎn)是大致趨同的??紤]到拆分后可能要多保存幾份碼表,拆分帶來收益的可能性和程度是很小的,而且計(jì)算的復(fù)雜度較大,所以前面的拆分策略可以不作調(diào)整。
  至于直接保存 lz77 的原始輸出,可以看作靜態(tài)編碼的一種特殊形式,只不過它假定各節(jié)點(diǎn)的頻率相近,沒有優(yōu)勢(shì)節(jié)點(diǎn)。它可以套用靜態(tài)編碼的分析,來證明不影響前面已經(jīng)制定的分段策略。

3.采用 huffman 編碼,必須深入研究碼表的保存方式。
  只要計(jì)算一下采用簡(jiǎn)單的方式來保存碼表,需要多大的空間,就知道這是一個(gè)挑戰(zhàn)。
  簡(jiǎn)單地保存碼表的方法是順序地保存每一個(gè)值的碼長(zhǎng)和編碼。之所以要保存碼長(zhǎng),是因?yàn)榫幋a是不定長(zhǎng)的,沒有碼長(zhǎng),解壓時(shí)無法正確讀取編碼。碼長(zhǎng)必須是定長(zhǎng)的,也就是說必須限制 huffman 樹的最大層數(shù),使碼長(zhǎng)的位數(shù)能恰好表示這個(gè)層數(shù)。限制 huffman 樹的最大層數(shù)的方法是:如果規(guī)定的最大層數(shù)為 n,則在 n - 1 層找到一個(gè)葉子節(jié)點(diǎn) a(如果 n - 1 層沒有葉子節(jié)點(diǎn),就逐層地往上尋找,直到找到一個(gè)葉子節(jié)點(diǎn)),在節(jié)點(diǎn) a 的位置放一個(gè)非葉子節(jié)點(diǎn) A,使 a 成為 A 的子節(jié)點(diǎn),把某個(gè)超過 n 層的葉子節(jié)點(diǎn) b 提上來作為 A 的另一個(gè)子節(jié)點(diǎn),此時(shí) b 的父節(jié)點(diǎn) B 只剩下一個(gè)子節(jié)點(diǎn) c,取消 B,把 c 放在 B 的位置,重復(fù)這樣的過程,直到所有 n 層以下的節(jié)點(diǎn)都被提上來。之所以要從 n - 1 層開始逐層往上找,是因?yàn)橄聦拥墓?jié)點(diǎn)頻率小,碼長(zhǎng)變化后的影響小。假設(shè)每一層節(jié)點(diǎn)的頻率相近,那么上層父節(jié)點(diǎn)的頻率是其下層子節(jié)點(diǎn)的兩倍,第 11 層節(jié)點(diǎn)的頻率只有第一層節(jié)點(diǎn)頻率的 1 / 1024,所以應(yīng)該從下往上找。
  現(xiàn)在就開始計(jì)算碼表大?。?
  對(duì)于 256 個(gè)原始字節(jié)值,限制它的 huffman 樹的層數(shù)為 0 - 15,碼長(zhǎng)就需要 4 位,256 個(gè)碼長(zhǎng)需要 4 bit * 256 = 128 字節(jié);而 256 個(gè)新編碼需要至少 256 字節(jié)。(當(dāng)二叉樹的所有葉子節(jié)點(diǎn)都放在第 8 層 —— 不算根節(jié)點(diǎn)一層,正好能放下 2 的 8 次方 = 256 個(gè)葉子節(jié)點(diǎn),其中任何一個(gè)葉子節(jié)點(diǎn)往上升,至少造成兩個(gè)葉子節(jié)點(diǎn)往下降。換一個(gè)角度說,如果在第 8 層以上存在一個(gè)葉子節(jié)點(diǎn) a,在節(jié)點(diǎn) a 的位置放一個(gè)非葉子節(jié)點(diǎn) A,使 a 成為 A 的子節(jié)點(diǎn),把某個(gè)超過 8 層的葉子節(jié)點(diǎn) b 提上來作為 A 的另一個(gè)子節(jié)點(diǎn),此時(shí) b 的父節(jié)點(diǎn) B 只剩下一個(gè)子節(jié)點(diǎn) c,取消 B,把 c 放在 B 的位置,此時(shí) a 增長(zhǎng)了一位,c 縮短了一位,b 縮短了至少一位,編碼的平均位長(zhǎng)縮短。所以,當(dāng)?shù)?nbsp;8 層以上不存在葉子節(jié)點(diǎn),所有葉子節(jié)點(diǎn)都放在第 8 層時(shí),編碼的平均位長(zhǎng)達(dá)到最短 —— 8位。)這套碼表共需至少 128 + 256 = 384 字節(jié)。
  256 個(gè)“匹配長(zhǎng)度”的情況與原始字節(jié)值相同,兩套碼表共需至少 384 * 2 = 768 字節(jié)。
  對(duì)于 32k 個(gè)“匹配距離”,如果限制該 huffman 樹的層數(shù)為 0 - 31,保存每個(gè)值的碼長(zhǎng)需要 5 位,新編碼的平均長(zhǎng)度超過 15 位。(因?yàn)樗腥~子節(jié)點(diǎn)都放在第 15 層 —— 不算根節(jié)點(diǎn)一層,正好能放下 2 的 15 次方 = 32k 個(gè)葉子節(jié)點(diǎn)。)這套碼表要超過80k 字節(jié)( (5 + 15) * 32k / 8 = 80k )。
  前面討論分段策略時(shí)已經(jīng)說過,為了避免個(gè)段落間節(jié)點(diǎn)頻率差異被互相抵消,要求段落劃分盡量細(xì)致、準(zhǔn)確,最小的段落可以僅為 4k,而采用上面這種簡(jiǎn)單的方式,碼表要超過 80k,顯然是無法接受的。
  對(duì)碼表的保存方式的深入研究,確實(shí)是個(gè)無法繞開的挑戰(zhàn),如果不攻克這個(gè)難關(guān),編碼式壓縮無法進(jìn)行下去!挑戰(zhàn)會(huì)帶來樂趣,困難會(huì)激發(fā)豪情。我們所要做的是:觀察 gzip 如何一步步地通過繁復(fù)但又巧妙的做法解決這個(gè)難題,對(duì)其中的做法的道理務(wù)求知其然、知其所以然,通過觀察、思考,把握無損壓縮內(nèi)在的、深層的、本質(zhì)的規(guī)律!事實(shí)上,對(duì) gzip 的這些做法進(jìn)行閱讀(源代碼)、分析、挖掘其中的智慧,本身就是一個(gè)對(duì)智慧、耐力、乃至決心的長(zhǎng)期的挑戰(zhàn),我接受了這個(gè)挑戰(zhàn),并把它描述、解釋出來,讀者面對(duì)的挑戰(zhàn)是花費(fèi)較長(zhǎng)期的時(shí)間去閱讀、理解,希望讀者完全有耐力、豪情、興趣來接受這個(gè)挑戰(zhàn),深化自己的技術(shù)層次、思維層次。

3.1 只保存碼長(zhǎng),并增加一些特殊的值。

3.1.1 把 huffman 樹的每一層上的葉子節(jié)點(diǎn)都換到該層的左邊,按照其原始值從小到大依次排列,非葉子節(jié)點(diǎn)則集中在該層右邊,這時(shí)仍是一棵二叉樹,得到的編碼仍符合前綴編碼的要求。每個(gè)葉子節(jié)點(diǎn)的編碼長(zhǎng)度不變,所以壓縮率也不變。僅需要按照原始值從小到大依次保存每個(gè)值的碼長(zhǎng),解壓時(shí)就可以還原這套編碼表,還原方法是:碼長(zhǎng)為 n 的第一個(gè)值的編碼是碼長(zhǎng)為 n - 1 的最后一個(gè)值的編碼加 1,并左移一位(也就是說,在編碼最后加個(gè) 0),而碼長(zhǎng)為 n 的其他值的編碼是前一個(gè)碼長(zhǎng)為 n 的值的編碼加 1。從上面所說的樹的角度來解釋,每一層的第一個(gè)葉子節(jié)點(diǎn)是其上層最后一個(gè)葉子節(jié)點(diǎn)的右邊一個(gè)節(jié)點(diǎn)的左子節(jié)點(diǎn),所以它的編碼是上層最后一個(gè)葉子節(jié)點(diǎn)的編碼加 1 并左移一位,而每一層上的葉子節(jié)點(diǎn)都緊密排列,所以除了第一個(gè)葉子節(jié)點(diǎn)外,其他葉子節(jié)點(diǎn)的編碼都是前一個(gè)葉子節(jié)點(diǎn)的編碼加 1。編程上的實(shí)現(xiàn)方法是:遍歷碼表,得到每個(gè)碼長(zhǎng)(n)上有多少個(gè)值,計(jì)算出每個(gè)碼長(zhǎng)上第一個(gè)值的編碼,放在數(shù)組 bit_len[]中,再次遍歷碼表,依次根據(jù)每個(gè)值的碼長(zhǎng)(n),賦予它的編碼為該碼長(zhǎng)上的前一個(gè)值的編碼 (bit_len[n]) 加 1,bit_len[n] ++。
  由于只需要保存碼長(zhǎng),現(xiàn)在碼表由超過 80k 字節(jié)減小到約 20k 字節(jié)。

3.1.2 如何只保存在段落中出現(xiàn)過的節(jié)點(diǎn)(有效節(jié)點(diǎn))的編碼?
  一個(gè) ASCⅡ文本,128 以后的值是不會(huì)在文件中出現(xiàn)的,按照 3.1.1 的方法,碼表中后半部分(都是 0)在解壓縮時(shí)是用不到的。為了避免這類浪費(fèi),只保存有效節(jié)點(diǎn)(碼長(zhǎng)不為 0 的節(jié)點(diǎn)),一種方法是保存有效節(jié)點(diǎn)的原始值和新編碼的碼長(zhǎng),當(dāng)有效節(jié)點(diǎn)超過所有節(jié)點(diǎn)的1/4,這種方法保存的碼表的大小會(huì)超過 3.1.1 的方法。
  gzip 采用的方法是:在 3.1.1 的基礎(chǔ)上,于若干種碼長(zhǎng)之外,增加一些特殊的值,他們表示當(dāng)前為之前一個(gè)碼長(zhǎng)或 0 碼長(zhǎng)(無效節(jié)點(diǎn))的重復(fù),遇到這種值,那后面的一個(gè)數(shù)字表示重復(fù)的次數(shù)。第一種值代表當(dāng)前為之前一個(gè)碼長(zhǎng)的重復(fù) 3 - 6 次,后面跟著 2 bit 為具體的重復(fù)次數(shù);第二種值代表當(dāng)前為 0 碼長(zhǎng)的重復(fù) 3 - 10 次,后面跟著 3 bit 為具體的重復(fù)次數(shù);第三種值代表當(dāng)前為 0 碼長(zhǎng)的重復(fù) 11 - 138 次,后面跟著 7 bit 為具體的重復(fù)次數(shù)。限制最小重復(fù)次數(shù)為 3,可以確保這種方法得到的碼表不會(huì)大過 3.1.1。第一種值限制最大重復(fù)次數(shù)為 6,是因?yàn)檫B續(xù) 6 個(gè)值以上的碼長(zhǎng)相等(說明頻率十分接近)的情況不常見,做這個(gè)限制可以節(jié)省附加 bit;第二第三種值區(qū)分重復(fù)次數(shù)的范圍,也是為了節(jié)省附加 bit。在只有少數(shù)有效節(jié)點(diǎn)的情況下,這種方法只需要保存較少的數(shù)據(jù),同時(shí)也具有簡(jiǎn)單的去重復(fù)的作用。
  如果最大碼長(zhǎng)是 15,0 - 15 共 16 種值,一個(gè)碼長(zhǎng)需要 4 位,加上上面 3 種值,共 19 種值,需要 5 位,在重復(fù)不多時(shí),加了這 3 種值,是不是會(huì)增大碼表?其實(shí)不用擔(dān)心,gzip 會(huì)對(duì)碼表再進(jìn)行一次 huffman 壓縮,根據(jù)這 19 種值的頻率分配給它們可變碼長(zhǎng)的編碼,不會(huì)造成浪費(fèi),由于涉及到一些其他情況,對(duì)碼表的再編碼壓縮在后面還會(huì)詳細(xì)介紹!

3.2 把原始字節(jié)值和匹配長(zhǎng)度值建在一棵樹上。
  現(xiàn)在先考慮另一個(gè)問題:如何使解壓時(shí)能區(qū)分當(dāng)前是一個(gè)未匹配的字節(jié),還是一個(gè)匹配?未匹配字節(jié)值和匹配長(zhǎng)度、匹配距離是三棵不同的 huffman 樹,它們的編碼互相不符合前綴編碼的要求,部分節(jié)點(diǎn)甚至可能編碼相同,解壓時(shí)如何區(qū)分?
  第一種方法是用標(biāo)志位。輸出壓縮結(jié)果時(shí),除了輸出每一段的碼表、重新編碼后的數(shù)據(jù)流,還要保存對(duì)應(yīng)于這一段數(shù)據(jù)的標(biāo)志位流,流中的每一位 0 或 1 表示當(dāng)前是一個(gè)未匹配的字節(jié),還是一個(gè)匹配。
  第二種方法是給原始字節(jié)值和匹配長(zhǎng)度值不同的編碼,并符合前綴編碼的要求。最好的做法是把它們建在一棵樹上,以確保它們符合前綴編碼的要求,并由它們的頻率來確定各自的碼長(zhǎng)。
  第一種方法相當(dāng)于原始字節(jié)值和匹配長(zhǎng)度值的編碼都增長(zhǎng)一位。
  第二種方法中這兩套節(jié)點(diǎn)的碼長(zhǎng)變化要根據(jù)具體節(jié)點(diǎn)各自的頻率而定。
  經(jīng)過分析,第二種方法更好,因?yàn)榈谝环N方法可以看作是第二種方法的變種,相當(dāng)于簡(jiǎn)單地在兩棵 huffman 樹的根節(jié)點(diǎn)上再加一個(gè)父節(jié)點(diǎn),這樣顯然是不能保證最佳的結(jié)果的。

3.3 把匹配長(zhǎng)度、匹配距離變?yōu)殚L(zhǎng)度范圍、距離范圍,減少節(jié)點(diǎn)。
  經(jīng)過上面對(duì)保存碼表的方法的改進(jìn)后,現(xiàn)在碼表還有多大?
  由于有了上面介紹的去重復(fù)機(jī)制,碼表的實(shí)際大小和節(jié)點(diǎn)的重復(fù)情況有關(guān),如果有很多連續(xù) 3 個(gè)以上節(jié)點(diǎn)的碼長(zhǎng)相等的情況出現(xiàn),或有很多連續(xù) 3 個(gè)以上的無效節(jié)點(diǎn)的情況出現(xiàn),碼表可能是很小的,但作為通用的無損壓縮算法,必須考慮重復(fù)不多的情況。“匹配距離”是碼表中最主要的部分,我們來分析一下它的重復(fù)情況,“匹配距離”共有 32k 個(gè)取值,如果一個(gè)段落不到 32k,“匹配距離”的有效節(jié)點(diǎn)數(shù)當(dāng)然是不可能到 32k 的,思考一下,可以知道,它的有效節(jié)點(diǎn)數(shù)和這樣幾個(gè)因素有關(guān):一段有多長(zhǎng),段落中匹配數(shù)和未匹配數(shù)的比例,決定了它有多少個(gè)值,再加上這些值的重復(fù)性,決定了它有多少個(gè)有效節(jié)點(diǎn)。再分析一下這些值的重復(fù)性:不同于原始字節(jié)和“匹配長(zhǎng)度”都只有 256 個(gè)取值,它有 32k 個(gè)取值,相同的匹配有相同的匹配長(zhǎng)度但不一定有相同的匹配距離,所以它的去值范圍廣,重復(fù)率低,有效節(jié)點(diǎn)多。雖然實(shí)際的情況無法預(yù)測(cè),但我們可以做一些“大致合理”的假設(shè),以便對(duì)碼表的大小有一個(gè)基本的概念,假如短語式壓縮的輸出段落的大小為 90k 字節(jié),其中未匹配字節(jié)數(shù)和匹配數(shù)的比例為 3 : 1,每個(gè)未匹配字節(jié)占 8 位;每個(gè)匹配中,長(zhǎng)度占 8 位,距離占 15 位,共 23 位,約為未匹配字節(jié)的 3 倍,所以匹配占了 90k 字節(jié)中的約 45k 字節(jié),匹配數(shù)約 15k 個(gè),也就是說有 15k 個(gè)距離值,假如距離值的平均節(jié)點(diǎn)頻率為 3,那么去掉重復(fù)后有 5k 個(gè)有效距離值節(jié)點(diǎn),保存到碼表時(shí)每個(gè)碼長(zhǎng)需要 5 位,保存 5k 個(gè)碼長(zhǎng)需要 5k * 5 / 8 約 3k 字節(jié),算上無效節(jié)點(diǎn)、碼長(zhǎng)的重復(fù)的因素,原始字節(jié)值、匹配長(zhǎng)度的保存,最終碼表約 5k 字節(jié),為 90k 的 18 分之一。當(dāng)段落減小時(shí),有效節(jié)點(diǎn)趨于稀疏,無效節(jié)點(diǎn)容易連成片,去重復(fù)機(jī)制能發(fā)揮更大的作用;當(dāng)段落增大時(shí),無效節(jié)點(diǎn)密度減小,可能無法大片連接,去重復(fù)機(jī)制的效用降低,碼表的比例可能會(huì)增大。一旦“匹配距離”需要保存的碼長(zhǎng)數(shù)達(dá)到了 32k個(gè),碼表達(dá)到最大,之后段落再增大也不會(huì)增大碼表,于是碼表的比例又會(huì)逐漸下降。當(dāng)然段落通常不會(huì)達(dá)到這么大,使得“匹配距離”需要保存的碼長(zhǎng)數(shù)能有機(jī)會(huì)達(dá)到 32k。
  gzip 以犧牲壓縮率的代價(jià)來換取碼表的進(jìn)一步的大幅度減小。我們先描述一下它的具體做法,再來分析其利弊。
  gzip 把匹配長(zhǎng)度劃成 29 個(gè)范圍,把匹配距離劃成 30 個(gè)范圍,根據(jù)每個(gè)范圍中節(jié)點(diǎn)的總頻率,為 29 個(gè)長(zhǎng)度范圍加 258 個(gè)字節(jié)值建 huffman 樹:l_tree,為 30 個(gè)距離范圍建 huffman 樹:d_tree。輸出一個(gè)值時(shí)先輸出該值所在范圍的編碼,再輸出附加碼,即它是該范圍中的第幾個(gè)值。這樣碼表中只需保存范圍的碼長(zhǎng)。范圍的大小都是 2 的乘方,所以范圍大小和附加碼的位長(zhǎng)是互相決定的。
29 個(gè)長(zhǎng)度范圍的附加碼位長(zhǎng)是:
{0,0,0,0,0,0,0,0,1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5,0};
30 個(gè)距離范圍的附加碼位長(zhǎng)是:
{0,0,0,0,1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13};
可以看出:范圍的劃分是從小到大的。為什么不平均劃分呢?
  如果仍以單個(gè)節(jié)點(diǎn)的角度來看,被分到同一范圍的節(jié)點(diǎn)相當(dāng)于被賦予了相同的碼長(zhǎng):范圍編碼的碼長(zhǎng)加附加碼的碼長(zhǎng)。若頻率差別很大的節(jié)點(diǎn)因劃分入同一個(gè)范圍而擁有相同的碼長(zhǎng),就不符合 huffman 編碼的初衷,會(huì)對(duì)壓縮率產(chǎn)生不良影響。因此要求劃分后,范圍里的節(jié)點(diǎn)頻率相近,以盡量降低同一個(gè)范圍里不同節(jié)點(diǎn)間的相互影響。
  “匹配長(zhǎng)度”從短到長(zhǎng),頻率會(huì)逐漸衰減,而且衰減的幅度有從大到小的特點(diǎn),這個(gè)特點(diǎn)是在大多數(shù)原始文件中“自然存在”的。比如在 google 上搜索,2 個(gè)字的短語和 22 個(gè)字的短語,搜到的結(jié)果數(shù)差別巨大,200 個(gè)字和 220 個(gè)字,搜到的結(jié)果數(shù)差別就沒有那么大。頻率大致上單向地逐步變化,所以劃分范圍后,范圍內(nèi)節(jié)點(diǎn)的頻率較接近;變化速度由大到小,所以范圍的劃分應(yīng)該從小到大。
  “匹配距離”也有類似的特點(diǎn),對(duì)大多數(shù)文件來說,匹配發(fā)生在 1k 以內(nèi)比發(fā)生在 5k 左右的可能性要大得多,但發(fā)生在 28k 處附近的可能性和發(fā)生在 32k 處附近的可能性的差別就沒那么明顯。所以范圍劃分也應(yīng)該是從小到大。
  “未匹配的原始字節(jié)”不具有頻率衰減或遞增的單向變化的規(guī)律,它們的頻率分布往往是參差不齊、難以預(yù)測(cè)的,不可能用預(yù)先設(shè)定的范圍表對(duì)它們進(jìn)行大致合理的劃分,就像“匹配長(zhǎng)度”和“匹配距離”那樣。雖然也可以通過計(jì)算分析,對(duì)它們進(jìn)行不設(shè)定范圍數(shù)量和大小的劃分,以求每個(gè)范圍中的各節(jié)點(diǎn)頻率大致相近,但 1) “匹配距離”的劃分已經(jīng)大幅度地縮小了碼表的大??;2) 由于不具有頻率單向變化的趨向,要強(qiáng)行劃出節(jié)點(diǎn)頻率相近并且節(jié)點(diǎn)數(shù)是 2 的乘方的范圍太勉強(qiáng),難度也大;3) 未匹配的字節(jié)數(shù)一般要大于“匹配數(shù)”(注意:不是“匹配字節(jié)數(shù)”),強(qiáng)行劃分造成的不良反應(yīng)較大。所以 gzip 保留了這套節(jié)點(diǎn),沒去拆分。
  長(zhǎng)度范圍的最后一個(gè)附加碼位長(zhǎng)是 0,是因?yàn)殚L(zhǎng)度大于 258 的匹配都被截?cái)嗟?nbsp;258,所以 258 的頻率可能會(huì)高出前面的節(jié)點(diǎn),單獨(dú)劃為一個(gè)范圍。
  如果一個(gè)范圍里的節(jié)點(diǎn)頻率相同,節(jié)點(diǎn)數(shù)是 2 的乘方,且沒有無效節(jié)點(diǎn),那么這個(gè)范圍可以看作 huffman 樹中的一棵子樹,范圍的編碼可以看作這棵子樹的根的編碼,這樣的劃分是不會(huì)影響壓縮率的。
  對(duì)壓縮率的損害來自頻率不一致,以及無效節(jié)點(diǎn)的存在。范圍里的有效節(jié)點(diǎn)如果沒有過半,“附加碼”的位數(shù)就至少有一位浪費(fèi)了,也就是說,范圍里所有有效節(jié)點(diǎn)的碼長(zhǎng)無端增長(zhǎng)了一位,如果有效節(jié)點(diǎn)沒有過 1/4,至少就有 2 位附加碼浪費(fèi)。
  劃分范圍的收益是使碼表減小到不足 0.2k,加上后面會(huì)介紹的對(duì)碼表的第二次壓縮,碼表的最終大小是微不足道的。
  現(xiàn)在我們來近似地估計(jì)一下劃分范圍在“一般情況”下對(duì)壓縮率的損害的情況,以便有一個(gè)大致的概念,仍舉前面的例子:段落大小為 90k,設(shè)其中未匹配字節(jié)數(shù)和匹配數(shù)的比例為 3:1,未匹配字節(jié)有 45k 個(gè),匹配距離值和匹配長(zhǎng)度值各 15k 個(gè),有效距離值節(jié)點(diǎn)為 5k個(gè)(節(jié)點(diǎn)平均頻率為 3),無效距離值節(jié)點(diǎn)為 32k - 5k = 27k 個(gè),有效距離值節(jié)點(diǎn)的平均密度為 5/32,不到 1/6。范圍的劃分是前小后大,有效節(jié)點(diǎn)頻率是前大后小,無效節(jié)點(diǎn)是前少后多。距離值有 15k 個(gè),設(shè)前面有效節(jié)點(diǎn)頻率高、密度較大的部分占一半,約 7k 個(gè)值,這個(gè)部分中無效節(jié)點(diǎn)帶來的損害較小,而且范圍劃分細(xì),節(jié)點(diǎn)間頻率不一致帶來的損害也小,姑且不去計(jì)算。后面的范圍劃分大、有效節(jié)點(diǎn)密度小的部分損害較大,這個(gè)部分占了約 7k 個(gè)值,由于前面的部分有效節(jié)點(diǎn)密度大,所以假設(shè)這個(gè)部分有效節(jié)點(diǎn)密度為 1/8(也就是說,約一半的匹配發(fā)生在 1k 距離以內(nèi),且 1k 以內(nèi)無效節(jié)點(diǎn)很少,那么 4k / 31k 約等于 1/8),附加碼浪費(fèi)了 3 位,7k 個(gè)值浪費(fèi) 3 位,共浪費(fèi)了 21k bit 約等于 3k 字節(jié)。
  再看頻率不一致帶來的損害:huffman 編碼如果要達(dá)到 50% 的壓縮率,需要節(jié)點(diǎn)間頻率的差異達(dá)到幾百倍。讀者可以虛擬一些節(jié)點(diǎn)頻率,試著建一下 huffman 樹,會(huì)發(fā)現(xiàn)當(dāng)節(jié)點(diǎn)頻率差異在幾十倍甚至只有幾倍的時(shí)候,壓縮率其實(shí)微乎其微。經(jīng)過上面這樣合理地劃分范圍,范圍內(nèi)的節(jié)點(diǎn)頻率差異一般不會(huì)那么大,所以我們假設(shè)頻率不一致造成的損害為 1k - 2k。
  匹配長(zhǎng)度值的取值范圍只有 258 個(gè),而且匹配長(zhǎng)度可能很少會(huì)超過 20 字節(jié),而前 20 字節(jié)的范圍劃分是很細(xì)的,所以無效節(jié)點(diǎn)的損害和頻率不一致的損害都較小。
  這樣,在這個(gè)例子中,劃分范圍帶來的損害約在 5k - 6k,和不劃分范圍時(shí)碼表的大小非常相似,至少也是在一個(gè)數(shù)量級(jí)上。 
  再來看看損害比例變化的趨勢(shì):當(dāng)段落很小時(shí),范圍中的有效值稀疏,損害比例會(huì)加大。而不劃分范圍時(shí),碼表的去重復(fù)機(jī)制會(huì)有更大作用,無效節(jié)點(diǎn)連成片,損害比例減小。反之,段落增大,范圍里有效節(jié)點(diǎn)密度大,損害比例降低,而不劃分范圍時(shí),無效節(jié)點(diǎn)可能無法大片連接,去重復(fù)機(jī)制的效用降低,損害比例增大。
  由于劃分范圍能使 huffman 樹的節(jié)點(diǎn)從最多 32k 減到不足 320 個(gè),從而使壓縮速度顯著改善。綜上所述,段落?。ū热绮坏?nbsp;10k),不宜劃分范圍,否則劃分范圍是有益的。
3.4 對(duì)碼表進(jìn)行第二次壓縮。
  目前為止,碼表中只需要保存各個(gè)節(jié)點(diǎn)經(jīng)過 huffman 編碼后的新編碼的碼長(zhǎng)。共兩棵樹,l_tree: 256 個(gè)原始字節(jié)值加 29 個(gè)長(zhǎng)度范圍值加 1 個(gè)段落中止符,共 286 個(gè)節(jié)點(diǎn),段落中止符用來在解壓時(shí)標(biāo)示一個(gè)段落的終結(jié)。d_tree: 30 個(gè)距離范圍值。也就是說,共需要保存 286 + 30 = 316 個(gè)編碼的碼長(zhǎng)。gzip 限制 huffman 樹的最大層數(shù)為 15,這樣,碼長(zhǎng)就有 0 - 15 共 16 種值,再加上前面介紹過的去重復(fù)機(jī)制使用的 3 種特殊值,共 19 種值,如果就這樣保存碼表的話,每個(gè)碼長(zhǎng)都需要 5 位,才能表示 19 種值。我們觀察一下,316 個(gè)碼長(zhǎng),一共只有 19 種值,碼長(zhǎng)值的重復(fù)是必然的,而且由于 huffman 樹上每層的節(jié)點(diǎn)數(shù)不同,所以各個(gè)碼長(zhǎng)值的頻率也不一樣。所以還可以為這 19 種值再建 huffman 樹,進(jìn)行第二次編碼。這棵樹只有 19 個(gè)節(jié)點(diǎn),限制它的層數(shù)為 0 - 7,可以用 3 個(gè) bit 表示這 19 個(gè)節(jié)點(diǎn)的“長(zhǎng)度”。這樣,用新的“碼長(zhǎng)的編碼”來保存 316 個(gè)碼長(zhǎng),另需額外保存 3 * 19 = 57 bit,就可以解壓出這 19 個(gè)“碼長(zhǎng)的編碼”。(至于這 57 bit,就沒有必要再作第 3 次編碼了)


4. 解決了碼表的問題,現(xiàn)在再回過頭來看靜態(tài)編碼。
  靜態(tài)編碼是 gzip 預(yù)先設(shè)定的編碼方案,它的碼表是固定的。
  該如何合理設(shè)計(jì)這套編碼?作為 huffman 編碼的補(bǔ)助,它的耗時(shí)應(yīng)盡量少,前面說過,lz77 輸出一個(gè)分段之前,要比較 huffman 編碼和靜態(tài)編碼的壓縮結(jié)果,為了直接利用 lz77 輸出時(shí)做的匹配長(zhǎng)度范圍、匹配距離范圍的頻率的統(tǒng)計(jì),靜態(tài)編碼采用了同樣的范圍-附加碼的方案,這樣可以快速得到靜態(tài)編碼的壓縮結(jié)果大小。
  靜態(tài)編碼的碼長(zhǎng)的分配是這樣的:29 個(gè)長(zhǎng)度范圍中前 24 個(gè)范圍的碼長(zhǎng)為 7,后 5 個(gè)范圍的碼長(zhǎng)為8。原始字節(jié)值中 0 - 143 的碼長(zhǎng)為 8,144 - 255 的碼長(zhǎng)為 9。而 30 個(gè)距離范圍的碼長(zhǎng)為 5。根據(jù)這些預(yù)先設(shè)定的碼長(zhǎng)建立靜態(tài)的 l_tree 和 d_tree,編碼也就產(chǎn)生了。結(jié)合前面提到的附加碼位數(shù)的定義:
29 個(gè)長(zhǎng)度范圍的附加碼位長(zhǎng):
{0,0,0,0,0,0,0,0,1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5,0};
30 個(gè)距離范圍的附加碼位長(zhǎng):
{0,0,0,0,1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13};
讀者可以知道每一個(gè)值的實(shí)際碼長(zhǎng)。長(zhǎng)度范圍值和原始字節(jié)值建在一棵樹上,節(jié)點(diǎn)多所以碼長(zhǎng)較長(zhǎng),30 個(gè)距離范圍值只需要 5 位二進(jìn)制數(shù)表示。短匹配的長(zhǎng)度范圍值位長(zhǎng)較短,字節(jié)值 0 - 143 的位長(zhǎng)中等,其他字節(jié)值和長(zhǎng)匹配的長(zhǎng)度范圍值較長(zhǎng)。這樣的分配反映了 gzip 作者對(duì)“大多數(shù)”文件中各種值的頻率的粗略估計(jì)。作為一個(gè)通用的壓縮算法,無法預(yù)先知道一個(gè)文件的實(shí)際情況,不可能做精確的估計(jì)。
  進(jìn)一步的思考:靜態(tài)編碼有必要嗎?靜態(tài)編碼采用了和 huffman 編碼相同的范圍-附加碼的方案,在碼長(zhǎng)的分配上不可能超過 huffman 編碼,如果能“獲勝”,那就是勝在不需要保存碼表上,而前面分析過,碼表是很小的,對(duì)壓縮率沒有多大影響,所以 gzip 設(shè)計(jì)的這個(gè)靜態(tài)編碼方案應(yīng)該是可有可無的。

5. 關(guān)于堆排序算法。
  似乎已經(jīng)解決了所有的難題,但是對(duì)于沒有學(xué)過數(shù)據(jù)結(jié)構(gòu)的讀者,仍然有一個(gè)會(huì)對(duì)程序效率產(chǎn)生影響的問題需要關(guān)注,那就是“排序”。
  已經(jīng)講過,huffman 算法就是從一個(gè)節(jié)點(diǎn)序列中,不斷找出兩個(gè)最小的節(jié)點(diǎn),為它們建一個(gè)父節(jié)點(diǎn),值為這兩個(gè)節(jié)點(diǎn)之和,然后從節(jié)點(diǎn)序列中去除這兩個(gè)節(jié)點(diǎn),加入它們的父節(jié)點(diǎn)到序列中,不斷重復(fù)這樣的步驟,直到節(jié)點(diǎn)序列中只剩下一個(gè)節(jié)點(diǎn)。如何快速地找出最小的元素呢?
  在普通的線性羅列的數(shù)據(jù)結(jié)構(gòu)中,從 N 個(gè)元素中找出最小的元素的時(shí)間和 N 成正比,如果數(shù)據(jù)以我們所要介紹的“堆”的結(jié)構(gòu)存儲(chǔ),時(shí)間和 lg N 成正比(注:lg 以 2 為底數(shù),如 lg 256 = 8,lg 1024 = 10 ...)。 集合中的元素越多,堆排序算法的優(yōu)勢(shì)越突出,而且堆排序非常適合于在數(shù)據(jù)序列中不斷地取走最小的元素并加入新的元素。

5.1 什么是堆?
  堆首先是一棵“完全二叉樹”,即所有的葉子節(jié)點(diǎn)都在樹的最低二層,最低一層的節(jié)點(diǎn)依次靠左排列的二叉樹。如圖:

                       完全二叉樹
                        ?。?
              +----------○----------+
              |                    ?。?
     ?。穑         。穑?
     ?。             。         。      。?
  +---○---+     ?。穑   。穑  。穑?
 ?。      。     。      。   。  。  。  。?
+-○-+  ?。穑 。穑  。穑  觥   觥   觥   ?
|  ?。  。  。 。  。  。  。?
■   ■   ■   ■  ■   ■   ■   ■


  堆分大根堆和小根堆,大根堆的所有子節(jié)點(diǎn)都小于它的父節(jié)點(diǎn),小根堆的所有子節(jié)點(diǎn)都大于它的父節(jié)點(diǎn)。下面就是一個(gè)小根堆:

                         小根堆
                         ?。?
             ?。玻?
             ?。                     。?
     ?。常         。福?
     ?。             。         。       。?
 ?。叮     。矗   。保担  。保福?
 ?。      。     。      。   。   。  。   。?
+-8-+  ?。梗 。担  。担 。保丁 。玻啊  。保埂  。玻?
|  ?。  。  。 。  。  。  。?
9  ?。埂 。保薄 。保场 。丁  。浮  。丁  。?

5.2 堆如何在內(nèi)存中存儲(chǔ)?
  堆存放在一個(gè)數(shù)組中,存放的順序是:從根開始,依次存放每一層從左至右的節(jié)點(diǎn)。
5.3 如何尋找任意節(jié)點(diǎn)的子節(jié)點(diǎn)和父節(jié)點(diǎn)?
  數(shù)組中第 k 個(gè)元素,它的左子節(jié)點(diǎn)是第 2k 個(gè)元素,右子節(jié)點(diǎn)是第 2k + 1 個(gè)元素。它的父節(jié)點(diǎn)是┖ k/2 ┚(注:┖ X ┚表示小于等于 X 的最大整數(shù))。
5.4 如何建立堆?
  先把 n 個(gè)元素依次放入數(shù)組中,令變量 k = ┖ n/2 ┚,這時(shí)第 k 個(gè)元素是最后一個(gè)元素的父節(jié)點(diǎn),從第 k 個(gè)元素的兩個(gè)子節(jié)點(diǎn)中找出較小的一個(gè)與 k 元素比較,如果小于 k 元素,就和 k 元素交換一下位置,換位后的原先的 k 元素再和新的子節(jié)點(diǎn)比較(如果有子節(jié)點(diǎn)的話),直到它不再小于新的子節(jié)點(diǎn)或沒有子節(jié)點(diǎn)。令 k = k - 1。再重復(fù)上面的做法直到 k < 1,一個(gè)堆就建成了。
5.5 如何從堆中找出第二個(gè)最小的元素?
  把堆中第一個(gè)元素(最小的元素)存放到其他地方,把第 n 個(gè)元素(最后一個(gè))放到第一個(gè)的位置,再用前面的方法和下層節(jié)點(diǎn)交換直到它放到合適的位置,這時(shí)數(shù)組仍然是一個(gè)堆,第一個(gè)元素是最小的節(jié)點(diǎn),數(shù)組的最后一個(gè)有效節(jié)點(diǎn)是第 n - 1 個(gè)元素。
  花費(fèi)的時(shí)間和交換的次數(shù)成正比,最大的可能的交換次數(shù)是: 堆的層數(shù) - 1 =┏ lg (元素?cái)?shù) + 1) ┒- 1(注:┏ X ┒表示大于等于 X 的最小整數(shù))。
  現(xiàn)在可以看到,堆之所以采用完全二叉樹的形式,是為了樹的層數(shù)盡可能少。
  而抽出最后一個(gè)元素放到樹根,而不是抽出第二層的元素,是為了維持完全二叉樹的結(jié)構(gòu)!
5.6 如何加入新的元素到堆中?
  把第一個(gè)元素存放到其他地方,把新的元素放到第一個(gè)的位置,再用前面的方法和下層節(jié)點(diǎn)交換,直到它被放到合適的位置,此時(shí)數(shù)組中仍然是一個(gè)堆。

6. 建 huffman 樹和編碼的算法:
  如果現(xiàn)在有 n 個(gè)待編碼的節(jié)點(diǎn),按照原始數(shù)值從小到大存放在數(shù)組 tree[n] 中,那么,將要建立的 huffman 樹總共會(huì)有 2n -1 節(jié)點(diǎn),包括葉子節(jié)點(diǎn)和非葉子節(jié)點(diǎn)。申請(qǐng)一塊內(nèi)存,大小是能放下 huffman 樹的所有節(jié)點(diǎn),先把 n 個(gè)待編碼節(jié)點(diǎn)放入這塊內(nèi)存的左端,然后用“堆排序”算法先把它們建成一個(gè)堆。
  然后不斷用“堆排序”算法取出頻率最小的節(jié)點(diǎn),把它們從右到左、從小到大排放在內(nèi)存塊的右端,每當(dāng)取出兩個(gè)節(jié)點(diǎn),給它們生成一個(gè)父節(jié)點(diǎn),頻率等于它們之和,加入堆中。這樣直到堆中只剩下一個(gè)根節(jié)點(diǎn),這時(shí),內(nèi)存中從左到右存儲(chǔ)的是頻率從大到小的所有節(jié)點(diǎn),一棵 huffman 樹其實(shí)也就建成了,層數(shù)小的節(jié)點(diǎn)在前,層數(shù)大的節(jié)點(diǎn)在后,每一層的節(jié)點(diǎn)又是按頻率從大到小依次排列。
  申請(qǐng)兩個(gè)數(shù)組:bl_count[],bl_base[]。置根節(jié)點(diǎn)的碼長(zhǎng)為 0,從左至右,所有節(jié)點(diǎn)的碼長(zhǎng)(len)為它的父節(jié)點(diǎn)的碼長(zhǎng) + 1,如果是葉子節(jié)點(diǎn),bl_count[len]++,得到了每一層上的葉子節(jié)點(diǎn)數(shù)目。令變量 code = 0,然后根據(jù) bl_count[] 生成 bl_base[]:碼長(zhǎng) len 從 1 開始遞增,bl_base[len] = code = (code + bl_count[len - 1]) << 1,得到了每一層上第一個(gè)葉子節(jié)點(diǎn)的編碼。
  現(xiàn)在所有待編碼節(jié)點(diǎn)都被賦予了碼長(zhǎng),遍歷待編碼節(jié)點(diǎn),根據(jù)它們的碼長(zhǎng)得到它們的編碼:序號(hào) n 遞增,tree[n].code = bl_base[ tree[n].len ] ++。
  注意:我們前面討論碼表的時(shí)候說過,gzip 對(duì) huffman 編碼進(jìn)行了改進(jìn),只需要得到每一個(gè)葉子節(jié)點(diǎn)(待編碼節(jié)點(diǎn))的碼長(zhǎng),就可以進(jìn)行編碼,而不需要關(guān)心它的父節(jié)點(diǎn)的編碼是什么。而保存碼表時(shí),只需要保存碼長(zhǎng)。

動(dòng)態(tài) huffman 壓縮和解壓的整個(gè)流程:
壓縮:
  lz77 的壓縮過程中輸出未匹配的單雙字節(jié),和匹配,并統(tǒng)計(jì)各字節(jié)值和匹配長(zhǎng)度范圍、匹配距離范圍的頻率,根據(jù)這些頻率建立兩棵 huffman 樹:ltree、dtree,得到這兩棵樹上所有節(jié)點(diǎn)的長(zhǎng)度和編碼。
  統(tǒng)計(jì)這兩棵樹節(jié)點(diǎn)長(zhǎng)度的使用頻率,對(duì)各節(jié)點(diǎn)長(zhǎng)度建立 huffman 樹:bl_tree,得到 bl_tree 的長(zhǎng)度和編碼。
  存儲(chǔ) bl_tree 的節(jié)點(diǎn)長(zhǎng)度數(shù)組。
  再用 bl_tree 的編碼存儲(chǔ) ltree、dtree 的節(jié)點(diǎn)長(zhǎng)度數(shù)組。
  再用 ltree 的編碼存儲(chǔ)各字節(jié)值和匹配長(zhǎng)度范圍(及附加碼)的流;用 dtree 的編碼存儲(chǔ)匹配距離范圍(及附加碼)的流。
解壓:
  先根據(jù) bl_tree 的節(jié)點(diǎn)長(zhǎng)度數(shù)組得到 bl_tree 的編碼。
  再用這些編碼得到 ltree、dtree 的節(jié)點(diǎn)長(zhǎng)度數(shù)組,進(jìn)而得到 ltree、dtree 的編碼。
  再根據(jù) ltree、dtree 的編碼及附加碼的定義,得到 lz77 的輸出的原始結(jié)果:各字節(jié)值和匹配長(zhǎng)度的流,匹配距離的流。
后記:
  寫作本文花費(fèi)了超過一年的業(yè)余時(shí)間,其實(shí)看懂 gzip 源碼只用了一個(gè)半月的業(yè)余時(shí)間,等真正開始寫這篇文章的時(shí)候,發(fā)現(xiàn)深入分析無損壓縮算法要投入的心力會(huì)遠(yuǎn)超過我原來的想象,不是光靠“毅然決然的態(tài)度”和“拼搏精神”就可以完成的。只有耐心地去付出。
  經(jīng)過了一年多的時(shí)間,終于有了現(xiàn)在這樣質(zhì)量的這篇文章。這期間,我的工作已經(jīng)從應(yīng)用工程師轉(zhuǎn)變到了研究員,應(yīng)該說,寫作這篇文章對(duì)促成我把今后的工作轉(zhuǎn)變?yōu)楦阊芯渴怯杏绊懙?。所以這篇文章對(duì)我自己的人生道路當(dāng)然是有重要的意義,我也希望它會(huì)促成讀者投身研究的決心。
  巴甫洛夫說:“科學(xué)研究需要的是偉大的熱情和艱苦的勞作”,從看到這句話起,我就一直很喜歡它,常常會(huì)想起這句話。希望這篇文章能使讀者聯(lián)想到“長(zhǎng)久的熱情和耐心的勞作”,并在生活和工作中貫徹這種精神。
  一篇文章發(fā)布以后,它的全部?jī)r(jià)值就在于讀者的閱讀,感謝讀者諸君。

相關(guān)文章

最新評(píng)論

亚洲国产美女一区二区三区软件| 蜜桃久久久久久久人妻| sw137 中文字幕 在线| 精品视频国产在线观看| 性感美女福利视频网站| 粉嫩av蜜乳av蜜臀| 人人爱人人妻人人澡39| 日韩北条麻妃一区在线| rct470中文字幕在线| 国产1区,2区,3区| 亚洲一区二区三区五区| 黄色视频在线观看高清无码| 欧美 亚洲 另类综合| 日韩精品中文字幕福利| 在线观看免费视频网| 2020av天堂网在线观看| 男人天堂最新地址av| 国产成人午夜精品福利| 亚洲成人午夜电影在线观看| av手机在线观播放网站| 大学生A级毛片免费视频| 亚洲午夜高清在线观看| 99久久超碰人妻国产| 97人人妻人人澡人人爽人人精品 | 免费男阳茎伸入女阳道视频| chinese国产盗摄一区二区| 大香蕉伊人国产在线| 99精品国产免费久久| 亚洲成人国产综合一区| 男人操女人的逼免费视频| 18禁网站一区二区三区四区| 欧美视频中文一区二区三区| 真实国产乱子伦一区二区| 欧美亚洲国产成人免费在线 | 亚洲av无乱一区二区三区性色| 日韩欧美国产一区ab| 国产伊人免费在线播放| 影音先锋女人av噜噜色| 成人国产影院在线观看| 3D动漫精品啪啪一区二区下载| 亚洲1069综合男同| 青青草在观免费国产精品| 免费岛国喷水视频在线观看| 国产97视频在线精品| 青草亚洲视频在线观看| 亚洲av日韩av网站| av老司机亚洲一区二区| 精品久久婷婷免费视频| 日韩人妻xxxxx| 欧美性感尤物人妻在线免费看| 99精品久久久久久久91蜜桃| 又色又爽又黄又刺激av网站| 香港一级特黄大片在线播放| 日本熟妇喷水xxx| 亚洲久久午夜av一区二区| 一区国内二区日韩三区欧美| 亚洲av极品精品在线观看| 无忧传媒在线观看视频| 一区国内二区日韩三区欧美| 五十路在线观看完整版| 老司机在线精品福利视频| 家庭女教师中文字幕在线播放| 精品av国产一区二区三区四区| 青青热久免费精品视频在线观看| 嫩草aⅴ一区二区三区| 久久久精品精品视频视频| 国产女孩喷水在线观看| 超碰在线中文字幕一区二区| 欧美黑人性暴力猛交喷水| 午夜激情久久不卡一区二区| 精品黑人巨大在线一区| 中文字幕一区二区自拍| 操日韩美女视频在线免费看 | 91九色porny国产在线| 欧美精品国产综合久久| 丰满的子国产在线观看| 一区二区熟女人妻视频| 人妻熟女中文字幕aⅴ在线| 污污小视频91在线观看| 3344免费偷拍视频| 婷婷久久久综合中文字幕| 在线视频精品你懂的| 日韩欧美制服诱惑一区在线| chinese国产盗摄一区二区| 大香蕉伊人国产在线| 欧美亚洲国产成人免费在线| 欧美在线偷拍视频免费看| 2022精品久久久久久中文字幕| 天天夜天天日天天日| 激情伦理欧美日韩中文字幕| 成年人黄视频在线观看| 噜噜色噜噜噜久色超碰| 黄色成人在线中文字幕| 老司机99精品视频在线观看| 视频 一区二区在线观看| 青青色国产视频在线| 久久久制服丝袜中文字幕| 经典亚洲伊人第一页| 亚洲最大黄了色网站| 亚洲 中文 自拍 另类 欧美| 久久久91蜜桃精品ad| 国产va精品免费观看| 在线观看免费av网址大全| 午夜久久久久久久精品熟女| 99精品视频在线观看免费播放| 国产性色生活片毛片春晓精品| 色秀欧美视频第一页| 亚洲1卡2卡三卡4卡在线观看| weyvv5国产成人精品的视频| 干逼又爽又黄又免费的视频| av老司机精品在线观看| 日本精品视频不卡一二三| 人妻丝袜精品中文字幕| 中文人妻AV久久人妻水| 精品一区二区三区三区色爱| 中文字幕第1页av一天堂网| 日本免费视频午夜福利视频| 亚洲欧美国产综合777| 国产又粗又猛又爽又黄的视频美国| 97超碰人人搞人人| 天堂va蜜桃一区入口| 国产精品人妻熟女毛片av久| 神马午夜在线观看视频| 区一区二区三国产中文字幕| 99精品国产自在现线观看| 国产精品3p和黑人大战| 神马午夜在线观看视频| 绝顶痉挛大潮喷高潮无码| 久久久久久国产精品| av天堂中文免费在线| 视频在线免费观看你懂得| 人妻丝袜榨强中文字幕| 久久久久久cao我的性感人妻| 亚洲天堂精品久久久| 97精品成人一区二区三区| 午夜久久久久久久99| 搡老妇人老女人老熟女| 91亚洲手机在线视频播放| 99精品视频之69精品视频| 精品黑人一区二区三区久久国产 | av网址在线播放大全| 午夜蜜桃一区二区三区| 一级a看免费观看网站| 亚洲精品在线资源站| 欧美一级视频一区二区| 9色精品视频在线观看| v888av在线观看视频| 18禁精品网站久久| 天天躁日日躁狠狠躁av麻豆| 日韩av熟妇在线观看| 日韩a级精品一区二区| 午夜激情精品福利视频| 日本人妻欲求不满中文字幕| 任我爽精品视频在线播放| jiujiure精品视频在线| 青草亚洲视频在线观看| 激情五月婷婷免费视频| 激情人妻校园春色亚洲欧美| 亚洲综合色在线免费观看| 亚洲日本一区二区久久久精品| 青青青激情在线观看视频| 欧美区一区二区三视频| 人妻少妇精品久久久久久| 91精品国产观看免费| 91福利视频免费在线观看| 欧美专区第八页一区在线播放| 国产九色91在线观看精品| 国产午夜亚洲精品不卡在线观看| 久久久久久cao我的性感人妻| 少妇被强干到高潮视频在线观看 | 五月天色婷婷在线观看视频免费| 狍和女人的王色毛片| 久草视频在线一区二区三区资源站| 自拍偷拍亚洲精品第2页| 天天色天天操天天透| 免费一级特黄特色大片在线观看| 91精品一区二区三区站长推荐| 热99re69精品8在线播放| 亚洲 色图 偷拍 欧美| 啊啊好慢点插舔我逼啊啊啊视频| 美女av色播在线播放| av俺也去在线播放| 免费看国产av网站| 蜜桃视频在线欧美一区| 日本韩国在线观看一区二区| 自拍偷拍亚洲精品第2页| 1024久久国产精品| 一级黄片久久久久久久久| 日韩欧美中文国产在线| 欧美一区二区三区久久久aaa| 国产又粗又黄又硬又爽| 丝袜亚洲另类欧美变态| 久草视频首页在线观看| 91片黄在线观看喷潮| 大鸡巴插入美女黑黑的阴毛| 国产亚洲欧美45p| 欧美日韩高清午夜蜜桃大香蕉| 国产视频网站国产视频| 亚洲 清纯 国产com| 55夜色66夜色国产精品站| 一个人免费在线观看ww视频| 99国产精品窥熟女精品| 一区二区在线视频中文字幕| 蜜桃色婷婷久久久福利在线| 日本男女操逼视频免费看| 亚洲 中文 自拍 无码| 国产精品大陆在线2019不卡| 国产 在线 免费 精品| 免费在线福利小视频| 蜜桃视频在线欧美一区| 欧美精品黑人性xxxx| 国产精品久久久久网| 在线免费观看国产精品黄色| 少妇人妻久久久久视频黄片| 中文字幕人妻三级在线观看| 日本熟妇一区二区x x| 久久久久久久精品老熟妇| 老熟妇xxxhd老熟女| 亚洲特黄aaaa片| 天天操天天干天天日狠狠插 | 黑人解禁人妻叶爱071| 亚洲va国产va欧美va在线| 在线观看国产网站资源| 欧洲亚洲欧美日韩综合| 日本熟女精品一区二区三区| 女同互舔一区二区三区| 精品91自产拍在线观看一区| 五月婷婷在线观看视频免费| 欧美日韩熟女一区二区三区| 中文字幕人妻av在线观看| 中文字幕人妻一区二区视频 | 91精品国产91久久自产久强| 亚洲一区二区三区av网站| 亚洲欧美一区二区三区爱爱动图| 白嫩白嫩美女极品国产在线观看| 深田咏美亚洲一区二区| 国产自拍黄片在线观看| 亚洲精品国品乱码久久久久| 直接观看免费黄网站| 在线免费91激情四射| 在线观看的a站 最新| 日本性感美女写真视频| 成人免费毛片aaaa| 久久久久久9999久久久久| 66久久久久久久久久久| 天天日夜夜干天天操| 亚洲美女高潮喷浆视频| 日本最新一二三区不卡在线| 视频久久久久久久人妻| 日本脱亚入欧是指什么| 老鸭窝日韩精品视频观看| 大香蕉伊人国产在线| 小泽玛利亚视频在线观看| 欧美国产亚洲中英文字幕| 在线观看视频 你懂的| 懂色av之国产精品| 99精品免费观看视频| 亚洲综合自拍视频一区| 大鸡巴操b视频在线| 国产精品视频一区在线播放| 亚洲av色图18p| av中文字幕在线导航| 丝袜肉丝一区二区三区四区在线| 国产精品视频欧美一区二区 | 欧美色呦呦最新网址| 青青在线视频性感少妇和隔壁黑丝| 日韩精品电影亚洲一区| 男人和女人激情视频| 日美女屁股黄邑视频| 97超碰免费在线视频| 日韩精品激情在线观看| 早川濑里奈av黑人番号| 国产精品自拍偷拍a| 亚洲一区二区激情在线| 亚洲中文精品字幕在线观看 | 好吊视频—区二区三区| 亚洲精品亚洲人成在线导航| 沈阳熟妇28厘米大战黑人| 天天操天天干天天插| 日韩加勒比东京热二区| 天天日天天舔天天射进去| 亚洲区欧美区另类最新章节| 91九色porny国产蝌蚪视频| 97成人免费在线观看网站| 欧美少妇性一区二区三区| 久久丁香婷婷六月天| 男女第一次视频在线观看| 影音先锋女人av噜噜色| 亚洲一区制服丝袜美腿| 亚洲av无码成人精品区辽| 色婷婷精品大在线观看| 熟女人妻在线观看视频| 人人妻人人爱人人草| 日韩伦理短片在线观看| 日韩成人性色生活片| 三级av中文字幕在线观看| 精品一区二区三区欧美| 少妇高潮一区二区三区| 干逼又爽又黄又免费的视频| 青青青青青青青青青国产精品视频| 成年人免费看在线视频| 2021最新热播中文字幕| 国产精彩福利精品视频| 女人精品内射国产99| 日本高清在线不卡一区二区| 宅男噜噜噜666国产| 四川乱子伦视频国产vip| av男人天堂狠狠干| 干逼又爽又黄又免费的视频| www天堂在线久久| 日本一区精品视频在线观看| 人妻爱爱 中文字幕| 好男人视频在线免费观看网站| 大尺度激情四射网站| 午夜精品福利91av| 大陆胖女人与丈夫操b国语高清| 在线观看国产免费麻豆| 国产亚洲欧美另类在线观看| 在线国产精品一区二区三区| 黑人乱偷人妻中文字幕| 丝袜美腿欧美另类 中文字幕| 丝袜长腿第一页在线| 强行扒开双腿猛烈进入免费版| 97人妻人人澡爽人人精品| 蜜臀av久久久久久久| 91国内精品久久久久精品一| av一区二区三区人妻| 一个色综合男人天堂| 无码日韩人妻精品久久| 中文字幕在线观看极品视频| 久草电影免费在线观看| 亚洲自拍偷拍综合色| 秋霞午夜av福利经典影视| 99亚洲美女一区二区三区| 91精品国产黑色丝袜| 97超碰最新免费在线观看| 在线免费观看亚洲精品电影| 久久久久久九九99精品| 成年人中文字幕在线观看| 中文 成人 在线 视频| 在线 中文字幕 一区| 国产中文精品在线观看| 久草免费人妻视频在线| 日韩北条麻妃一区在线| 中文字幕av一区在线观看| 黄色中文字幕在线播放| 亚洲中文精品人人免费| 狠狠鲁狠狠操天天晚上干干| 天天摸天天亲天天舔天天操天天爽| 亚洲综合在线视频可播放| 9国产精品久久久久老师| 天堂av狠狠操蜜桃| 亚洲男人让女人爽的视频| 亚洲国产成人最新资源| 少妇系列一区二区三区视频| 亚洲在线免费h观看网站| 国产va在线观看精品| 国产亚州色婷婷久久99精品| 午夜频道成人在线91| 日韩写真福利视频在线观看| 九色精品视频在线播放| 人人妻人人澡欧美91精品| 国产精品人久久久久久| 国产综合高清在线观看| 丝袜亚洲另类欧美变态| 成人亚洲精品国产精品| 亚洲av天堂在线播放| 色哟哟在线网站入口| 97国产在线av精品| 亚洲福利天堂久久久久久| 精品亚洲在线免费观看| 亚洲视频乱码在线观看| 黄色黄色黄片78在线| 亚洲伊人色一综合网| 美女 午夜 在线视频| 无码精品一区二区三区人| 老司机99精品视频在线观看| 久久久精品精品视频视频| 国产一区av澳门在线观看| 制服丝袜在线人妻中文字幕| 91免费观看国产免费| 午夜精品福利91av| 天天通天天透天天插| 亚洲最大免费在线观看| 国产精品自拍偷拍a| 五月色婷婷综合开心网4438| 91成人精品亚洲国产| 亚洲熟妇无码一区二区三区| 青青操免费日综合视频观看| 少妇被强干到高潮视频在线观看| 最新91精品视频在线| 亚洲粉嫩av一区二区三区| av一区二区三区人妻| 日本人妻少妇18—xx| 在线免费观看视频一二区| 男人的天堂一区二区在线观看| 亚洲欧美综合另类13p| 亚洲免费成人a v| 精品人妻伦一二三区久| 色花堂在线av中文字幕九九| 插小穴高清无码中文字幕| 人妻少妇亚洲精品中文字幕| 夜色17s精品人妻熟女| 中文字幕人妻三级在线观看| 91在线视频在线精品3| 久久久久久久久久久免费女人| 91久久综合男人天堂| 香蕉aⅴ一区二区三区| 2019av在线视频| 99精品视频在线观看婷婷| 欧美成人猛片aaaaaaa| 午夜蜜桃一区二区三区| 中文字幕在线视频一区二区三区| 国产精品一区二区三区蜜臀av| 偷拍美女一区二区三区| 91av精品视频在线| 亚洲欧美成人综合在线观看| 天天操天天爽天天干| 一区二区三区的久久的蜜桃的视频 | brazzers欧熟精品系列| 成人福利视频免费在线| 欧美地区一二三专区| 亚洲一区二区激情在线| 国产精品视频欧美一区二区| 久久精品美女免费视频| 午夜场射精嗯嗯啊啊视频| 99精品国产免费久久| 成年人黄色片免费网站| 亚洲午夜高清在线观看| 91啪国自产中文字幕在线| okirakuhuhu在线观看| 国产精品三级三级三级| 色偷偷伊人大杳蕉综合网| 欧美黄片精彩在线免费观看| 午夜毛片不卡免费观看视频| 国产精品熟女久久久久浪潮| 91精品高清一区二区三区| 边摸边做超爽毛片18禁色戒| 婷婷综合蜜桃av在线| 特级欧美插插插插插bbbbb| 亚洲一区二区三区久久午夜| 中文乱理伦片在线观看| 99精品视频在线观看免费播放| 日视频免费在线观看| 国产精品自拍视频大全| 专门看国产熟妇的网站| 欧美另类一区二区视频| 日本又色又爽又黄又粗| 精品一区二区三区在线观看| sw137 中文字幕 在线| 一区二区三区综合视频| 综合激情网激情五月天| 成人蜜桃美臀九一一区二区三区| 高潮视频在线快速观看国家快速| 香港一级特黄大片在线播放| 久久精品久久精品亚洲人| 东京热男人的av天堂| 亚洲成a人片777777| 国产三级片久久久久久久| 国产九色91在线视频| 超碰97人人澡人人| 亚洲综合一区成人在线| 日韩av免费观看一区| 免费费一级特黄真人片| free性日本少妇| 99精品免费久久久久久久久a| 在线观看av亚洲情色| 青青草在观免费国产精品| 免费手机黄页网址大全| 黑人解禁人妻叶爱071| 黑人3p华裔熟女普通话| 中文字幕在线观看国产片| 女同性ⅹxx女同h偷拍| 自拍偷拍vs一区二区三区| 精品国产污污免费网站入口自| 韩国男女黄色在线观看| 91国语爽死我了不卡| 亚洲中文字幕人妻一区| 午夜精品久久久久麻豆影视| 国产V亚洲V天堂无码欠欠| 真实国产乱子伦一区二区| 午夜久久香蕉电影网| 九色porny九色9l自拍视频| 日韩特级黄片高清在线看| 国产97视频在线精品| 日本裸体熟妇区二区欧美| 人妻少妇一区二区三区蜜桃| chinese国产盗摄一区二区| 国产真实乱子伦a视频| 天堂av在线播放免费| 欧美特级特黄a大片免费| 美女张开两腿让男人桶av| 真实国产乱子伦一区二区| 大香蕉伊人中文字幕| 久久久久久久精品成人热| 亚洲精品国产久久久久久| 欧美性受xx黑人性猛交| 黄色资源视频网站日韩| 亚洲精品 日韩电影| 亚洲一区二区激情在线| 国产伊人免费在线播放| 日韩在线视频观看有码在线| 唐人色亚洲av嫩草| 真实国产乱子伦一区二区| 99热久久极品热亚洲| 开心 色 六月 婷婷| 不卡日韩av在线观看| 美女视频福利免费看| 综合精品久久久久97| 久久99久久99精品影院| 天美传媒mv视频在线观看| 骚货自慰被发现爆操| 韩国一级特黄大片做受| 久久三久久三久久三久久| 99热99re在线播放| 黄色成年网站午夜在线观看 | 久久久久久九九99精品| 97超碰国语国产97超碰| 午夜dv内射一区区| 久草视频在线看免费| 亚洲av无乱一区二区三区性色| 久久精品国产23696| 动漫美女的小穴视频| 日本丰满熟妇BBXBBXHD| 超级福利视频在线观看| 男人操女人逼逼视频网站| 男人操女人逼逼视频网站| 无码日韩人妻精品久久| 青青擦在线视频国产在线| 岳太深了紧紧的中文字幕| 天天日天天操天天摸天天舔| 亚洲色偷偷综合亚洲AV伊人| jiujiure精品视频在线| 日韩欧美高清免费在线| asmr福利视频在线观看| 国产精品久久综合久久| 国产午夜福利av导航| 馒头大胆亚洲一区二区| 激情国产小视频在线| 亚洲国产成人在线一区| 国产午夜男女爽爽爽爽爽视频| 人妻丝袜诱惑我操她视频| 亚洲 中文 自拍 另类 欧美| 成人动漫大肉棒插进去视频| 亚洲国产香蕉视频在线播放 | 任我爽精品视频在线播放| 人妻丝袜av在线播放网址| 青娱乐蜜桃臀av色| 亚洲欧美一区二区三区爱爱动图| 亚洲欧美另类自拍偷拍色图| 成年人啪啪视频在线观看| 人人妻人人爽人人添夜| 成人资源在线观看免费官网| 欧美精品激情在线最新观看视频| 日本人竟这样玩学生妹| 日韩欧美亚洲熟女人妻| 日比视频老公慢点好舒服啊| 青青擦在线视频国产在线| 漂亮 人妻被中出中文| 成人国产小视频在线观看| 亚洲天堂有码中文字幕视频| 国产欧美日韩第三页| 一区二区视频在线观看视频在线| 人妻激情图片视频小说| 久草视频在线看免费| 狠狠操操操操操操操操操| 久久丁香婷婷六月天| 91精品国产综合久久久蜜 | 亚洲2021av天堂| 欧美亚洲牲夜夜综合久久| 97欧洲一区二区精品免费| 大鸡巴操娇小玲珑的女孩逼| 日本丰满熟妇大屁股久久| 一区二区三区欧美日韩高清播放| 白白操白白色在线免费视频 | 小穴多水久久精品免费看| 国产av欧美精品高潮网站| 99热这里只有国产精品6| 日韩人妻丝袜中文字幕| 成人av在线资源网站| 亚洲精品国产综合久久久久久久久| 亚洲国产最大av综合| 亚洲精品中文字幕下载| 久久久精品999精品日本| 黄色视频成年人免费观看| 91亚洲手机在线视频播放| 亚洲天天干 夜夜操| 亚洲国际青青操综合网站| 抽查舔水白紧大视频| 真实国产乱子伦一区二区| 一级黄片大鸡巴插入美女| 在线播放 日韩 av| 日本少妇人妻xxxxxhd| 欧美男同性恋69视频| 欧美一区二区三区在线资源| 国产91久久精品一区二区字幕| 婷婷综合蜜桃av在线| 制服丝袜在线人妻中文字幕| 天天色天天爱天天爽| 国产精品成人xxxx| 亚洲精品亚洲人成在线导航| 91精品啪在线免费| 中文字幕之无码色多多| 日本熟妇色熟妇在线观看| 天堂av狠狠操蜜桃| 99国产精品窥熟女精品| 丰满的子国产在线观看| 女生自摸在线观看一区二区三区 | 啊用力插好舒服视频| 中文字幕—97超碰网| 日本精品美女在线观看| 免费黄页网站4188| 美女日逼视频免费观看| 人妻最新视频在线免费观看| 福利视频网久久91| 亚洲一级av大片免费观看| 91国偷自产一区二区三区精品| 亚洲超碰97人人做人人爱| 久久精品国产亚洲精品166m| 老司机你懂得福利视频| 东游记中文字幕版哪里可以看到| 国产三级片久久久久久久 | 亚洲av黄色在线网站| 一区二区三区综合视频| 免费观看国产综合视频| 天天摸天天干天天操科普| 男人在床上插女人视频| 欧美一区二区三区乱码在线播放 | 社区自拍揄拍尻屁你懂的| 91综合久久亚洲综合| 国产精彩福利精品视频| 岛国一区二区三区视频在线| 社区自拍揄拍尻屁你懂的| 成人sm视频在线观看| 欧美日韩熟女一区二区三区| 国产亚洲欧美另类在线观看| 99精品视频在线观看婷婷| 欧美另类重口味极品在线观看| 午夜毛片不卡在线看| 免费av岛国天堂网站| 黄片大全在线观看观看| 国产麻豆国语对白露脸剧情| 亚洲男人的天堂a在线| 午夜蜜桃一区二区三区| 国产妇女自拍区在线观看| 黄色黄色黄片78在线| 天天日天天操天天摸天天舔 | sejizz在线视频| 一色桃子久久精品亚洲| 色婷婷综合激情五月免费观看| 成人av电影免费版| 在线制服丝袜中文字幕| 人妻少妇亚洲精品中文字幕| 日本韩国亚洲综合日韩欧美国产| 在线免费91激情四射| 粗大的内捧猛烈进出爽大牛汉子 | 欧美老鸡巴日小嫩逼| 农村胖女人操逼视频| 哥哥姐姐综合激情小说| 2022中文字幕在线| 青青青青青青青在线播放视频| 在线观看911精品国产| 亚洲第17页国产精品| 免费手机黄页网址大全| 18禁污污污app下载| 欧美久久久久久三级网| 久久久久五月天丁香社区| 国产精品国产精品一区二区| 人妻3p真实偷拍一二区| 2020中文字幕在线播放| 男女之间激情网午夜在线| 97人妻色免费视频| 91极品新人『兔兔』精品新作| 亚洲国产在人线放午夜| 久久精品国产999| 大香蕉伊人中文字幕| 国产超码片内射在线| 99精品视频之69精品视频| 国产妇女自拍区在线观看| 国产黄色片蝌蚪九色91| 欧洲日韩亚洲一区二区三区| 91 亚洲视频在线观看| 大鸡巴操娇小玲珑的女孩逼| 午夜av一区二区三区| 日韩欧美高清免费在线| 亚洲图片欧美校园春色| 一区二区三区毛片国产一区| 一本一本久久a久久精品综合不卡| 国产在线拍揄自揄视频网站| 十八禁在线观看地址免费| 狠狠躁夜夜躁人人爽天天久天啪| 亚洲成人黄色一区二区三区 | 人人人妻人人澡人人| 日本丰满熟妇大屁股久久| 人妻少妇性色欲欧美日韩| 日韩成人性色生活片| 国产又大又黄免费观看| 婷婷综合蜜桃av在线| 精品久久久久久久久久中文蒉| 亚洲熟妇无码一区二区三区| 蜜桃色婷婷久久久福利在线 | 欧美黑人与人妻精品| 视频 一区二区在线观看| 中国把吊插入阴蒂的视频| 久久免看30视频口爆视频| 免费人成黄页网站在线观看国产| 啊用力插好舒服视频| 天天操夜夜操天天操天天操| 午夜青青草原网在线观看| 亚洲av日韩精品久久久| 婷婷午夜国产精品久久久| 免费69视频在线看| 337p日本大胆欧美人| yellow在线播放av啊啊啊| 国产视频在线视频播放| 亚洲高清自偷揄拍自拍| 99av国产精品欲麻豆| 一区国内二区日韩三区欧美| 青青青青爽手机在线| 亚洲高清国产拍青青草原| 欧美一级片免费在线成人观看| 美女日逼视频免费观看| 日韩二区视频一线天婷婷五| 久久久极品久久蜜桃| 免费观看理论片完整版| 精品人妻每日一部精品| 国产精品视频一区在线播放| 中文亚洲欧美日韩无线码 | 日韩少妇人妻精品无码专区| 日日日日日日日日夜夜夜夜夜夜| 国产高清精品极品美女| 久久尻中国美女视频| 中文字幕在线观看极品视频| 欧美特级特黄a大片免费| 丰满少妇翘臀后进式| 又粗又长 明星操逼小视频| 亚洲av无码成人精品区辽| 中文字幕日本人妻中出| 黄色中文字幕在线播放| 国产91精品拍在线观看| 国产大学生援交正在播放| 99亚洲美女一区二区三区| 亚洲自拍偷拍综合色| 夏目彩春在线中文字幕| 欧美偷拍亚洲一区二区| 日本熟妇喷水xxx| aⅴ五十路av熟女中出| 婷婷午夜国产精品久久久| 欧美香蕉人妻精品一区二区| 阿v天堂2014 一区亚洲| 大胆亚洲av日韩av| 一级黄色片夫妻性生活| 亚洲另类综合一区小说| 午夜久久久久久久精品熟女| 中文字幕高清在线免费播放 | 青青青青青操视频在线观看| 99国内精品永久免费视频| 五月婷婷在线观看视频免费| 青青草亚洲国产精品视频| 中文字幕奴隷色的舞台50| 欧美一区二区三区激情啪啪啪| 91人妻精品一区二区久久| 欧亚乱色一区二区三区| 国产日本欧美亚洲精品视| 又粗又长 明星操逼小视频| 狠狠躁狠狠爱网站视频| 偷拍自拍亚洲美腿丝袜| 国产麻豆乱子伦午夜视频观看| 国产污污污污网站在线| 视频一区二区三区高清在线| 久久久久久久久久久免费女人| 无码国产精品一区二区高潮久久4 日韩欧美一级精品在线观看 | 欧美激情精品在线观看| 国产视频网站国产视频| 青青擦在线视频国产在线| 中国视频一区二区三区| 国产视频在线视频播放| 亚洲欧美成人综合在线观看| 欧美特色aaa大片| 日韩a级精品一区二区| av黄色成人在线观看| 亚洲蜜臀av一区二区三区九色| 成人乱码一区二区三区av| 五十路在线观看完整版| 国产卡一卡二卡三乱码手机| 蝴蝶伊人久久中文娱乐网| 美女在线观看日本亚洲一区| 11久久久久久久久久久| 大陆精品一区二区三区久久| 国产乱弄免费视频观看| AV无码一区二区三区不卡| 亚洲视频在线观看高清| 欧美美女人体视频一区| 激情人妻校园春色亚洲欧美 | 超碰在线观看免费在线观看| 久久久久五月天丁香社区| 老司机欧美视频在线看| 天天日天天干天天干天天日| 大陆av手机在线观看| av手机在线免费观看日韩av| 2021久久免费视频| 18禁网站一区二区三区四区| 青青草原色片网站在线观看| av乱码一区二区三区| 国产精品人妻熟女毛片av久| 3344免费偷拍视频| 久久久久91精品推荐99| 青青青艹视频在线观看| 亚洲成人情色电影在线观看| 2019av在线视频| 亚洲天堂有码中文字幕视频| 91精品视频在线观看免费| 美女被肏内射视频网站| 亚洲熟色妇av日韩熟色妇在线| 人妻爱爱 中文字幕| 老鸭窝日韩精品视频观看| 免费十精品十国产网站| 亚洲精品午夜aaa久久| 在线免费观看黄页视频| 成人国产小视频在线观看| 天天操天天爽天天干| 偷偷玩弄新婚人妻h视频| 97人妻人人澡爽人人精品| 亚洲一区自拍高清免费视频| 国产卡一卡二卡三乱码手机| 亚洲精品国偷自产在线观看蜜桃| 美女福利视频导航网站| 77久久久久国产精产品| 国产精品久久9999| 三级黄色亚洲成人av| 亚洲高清国产自产av| 快点插进来操我逼啊视频| av线天堂在线观看| 最新国产亚洲精品中文在线| 一二三区在线观看视频| 国产一区自拍黄视频免费观看| 亚洲国产精品中文字幕网站| 久久免费看少妇高潮完整版| 中文字幕1卡1区2区3区| 很黄很污很色的午夜网站在线观看| 97国产福利小视频合集| 日韩人妻xxxxx| 偷拍自拍福利视频在线观看| 亚洲欧美成人综合视频| 东京热男人的av天堂| 天堂女人av一区二区| 成人影片高清在线观看| 中文字幕高清免费在线人妻 | 美女福利视频网址导航| 国产成人精品午夜福利训2021| 欧美va不卡视频在线观看| av高潮迭起在线观看| 青青草原色片网站在线观看| 国产 在线 免费 精品| 精品av国产一区二区三区四区 | 亚洲人一区二区中文字幕| 亚洲一级av大片免费观看| 亚洲精品中文字幕下载| 18禁污污污app下载| 国产片免费观看在线观看| 2020国产在线不卡视频| 日韩激情文学在线视频| 国产日韩精品电影7777| 在线观看亚洲人成免费网址| 午夜精品久久久久久99热| 久草视频在线免播放| 大鸡八强奸视频在线观看| 国产精品午夜国产小视频| 中文字幕一区二区自拍| 国产又粗又黄又硬又爽| 午夜毛片不卡免费观看视频| 亚洲无线观看国产高清在线| 欧美国产亚洲中英文字幕| 熟女国产一区亚洲中文字幕| 亚洲成人线上免费视频观看| 一本久久精品一区二区| www天堂在线久久| 99久久激情婷婷综合五月天| 日本最新一二三区不卡在线| 国产精品女邻居小骚货| 天天插天天狠天天操| 在线视频免费观看网| 中文字幕日韩精品就在这里| 在线观看成人国产电影| 日韩写真福利视频在线观看| 不卡精品视频在线观看| 天天干天天日天天谢综合156 | 可以在线观看的av中文字幕| 免费岛国喷水视频在线观看| 自拍偷拍,中文字幕| 57pao国产一区二区| 欧美成人猛片aaaaaaa| 精品首页在线观看视频| 天天日天天做天天日天天做| 日本熟女50视频免费| 天天操天天干天天插| 夜色撩人久久7777| www久久久久久久久久久| 自拍偷拍亚洲另类色图| 久草免费人妻视频在线| 2018在线福利视频| 99久久成人日韩欧美精品| gogo国模私拍视频| 国产揄拍高清国内精品对白| 日视频免费在线观看| 沙月文乃人妻侵犯中文字幕在线 | 亚洲免费福利一区二区三区| 熟女在线视频一区二区三区| 欧美3p在线观看一区二区三区| 人妻无码中文字幕专区| 亚洲自拍偷拍精品网| 国产福利小视频免费观看| 在线国产中文字幕视频| 国产97在线视频观看| 中国黄片视频一区91| 午夜在线观看岛国av,com| 青青青aaaa免费| 伊人日日日草夜夜草| 黑人进入丰满少妇视频| 边摸边做超爽毛片18禁色戒| 国产精品伦理片一区二区| 色综合久久久久久久久中文| 精品亚洲中文字幕av| 精品欧美一区二区vr在线观看 | 国产视频在线视频播放| 18禁精品网站久久| 午夜精彩视频免费一区| 香蕉aⅴ一区二区三区| 国产精品女邻居小骚货| 人妻激情图片视频小说| 99婷婷在线观看视频| 2020av天堂网在线观看| 久久香蕉国产免费天天| 亚洲伊人久久精品影院一美女洗澡| 精品一区二区三区三区色爱| 天天日天天干天天插舔舔| 亚洲精品三级av在线免费观看| 久久永久免费精品人妻专区| 99精品亚洲av无码国产另类| 亚洲av天堂在线播放| 国产日韩av一区二区在线| 好太好爽好想要免费| 天天日夜夜干天天操| av俺也去在线播放| 亚洲 中文字幕在线 日韩| 国产三级片久久久久久久| 色秀欧美视频第一页| 亚洲视频在线视频看视频在线| 一区国内二区日韩三区欧美| 亚洲第17页国产精品| 免费观看污视频网站| 超黄超污网站在线观看| 日韩写真福利视频在线观看| 国产精品国产三级麻豆| 成年人午夜黄片视频资源| 一二三中文乱码亚洲乱码one | 在线国产中文字幕视频| 亚洲高清国产一区二区三区| 日韩美在线观看视频黄| 黄色无码鸡吧操逼视频| 欧美视频综合第一页| 伊人成人在线综合网| 视频 一区二区在线观看| 日本女人一级免费片| aⅴ五十路av熟女中出| 精品一区二区亚洲欧美| 在线不卡日韩视频播放| 91天堂精品一区二区| 亚洲一区二区三区偷拍女厕91| 国产成人无码精品久久久电影| 日韩伦理短片在线观看| 97国产精品97久久| 国产精品日韩欧美一区二区| 青青青青在线视频免费观看| 丰满少妇人妻xxxxx| 精品一区二区三区在线观看| 99re久久这里都是精品视频| 2o22av在线视频| 成人av免费不卡在线观看| 一色桃子人妻一区二区三区| 免费在线看的黄片视频| 欧美麻豆av在线播放| 中文字幕国产专区欧美激情| aiss午夜免费视频| 顶级尤物粉嫩小尤物网站| 亚洲国产成人最新资源| 国产va在线观看精品| 亚洲2021av天堂| 中文字幕日本人妻中出| 成人福利视频免费在线| 天天日天天操天天摸天天舔| 久草视频在线一区二区三区资源站| japanese日本熟妇另类| ka0ri在线视频| 91老师蜜桃臀大屁股| 97人妻色免费视频| 国产精品久久久久国产三级试频| 2022中文字幕在线| 亚洲av日韩精品久久久久久hd| 蜜桃专区一区二区在线观看| 人妻av无码专区久久绿巨人| 亚洲2021av天堂| 亚洲另类在线免费观看| 黑人乱偷人妻中文字幕| 男人的天堂在线黄色| 国产亚洲四十路五十路| 无码精品一区二区三区人| 日视频免费在线观看| av网址在线播放大全| 任你操任你干精品在线视频 | 欧美黑人性猛交xxxxⅹooo| 青草久久视频在线观看| 国产福利在线视频一区| 日本xx片在线观看| 欧美精品一区二区三区xxxx| 九一传媒制片厂视频在线免费观看| 国产精品国色综合久久 | 国产三级影院在线观看| 黄色三级网站免费下载| 成年美女黄网站18禁久久| 国内精品在线播放第一页| 一本久久精品一区二区| 成人18禁网站在线播放| 91久久人澡人人添人人爽乱| 中文人妻AV久久人妻水| 真实国模和老外性视频| 婷婷色国产黑丝少妇勾搭AV| 青青伊人一精品视频| 欧美精品久久久久久影院| 欧美久久久久久三级网| 一区二区麻豆传媒黄片 | 欧洲精品第一页欧洲精品亚洲| 亚洲av一妻不如妾| 成人影片高清在线观看 | 懂色av蜜桃a v| 大香蕉大香蕉在线看| 日本av在线一区二区三区| 中文字幕人妻熟女在线电影| 在线观看av观看av| 色哟哟在线网站入口| 日本三极片视频网站观看| 97黄网站在线观看| 97精品综合久久在线| 日本性感美女写真视频| 搡老熟女一区二区在线观看| 午夜福利人人妻人人澡人人爽| 超碰97人人做人人爱| 亚洲精品 日韩电影| aiss午夜免费视频| 日本人妻欲求不满中文字幕| 日韩不卡中文在线视频网站| 青草久久视频在线观看| 青青青青操在线观看免费| 337p日本大胆欧美人| 日韩少妇人妻精品无码专区 | av黄色成人在线观看| 欧美黄色录像免费看的| 淫秽激情视频免费观看| 天天做天天干天天舔| 亚洲va国产va欧美va在线| 久草视频在线一区二区三区资源站 | 99re国产在线精品| 免费费一级特黄真人片 | 18禁免费av网站| 日本少妇在线视频大香蕉在线观看| 久草视频福利在线首页| 日本一道二三区视频久久 | 91国产在线视频免费观看| 美洲精品一二三产区区别| 99热99这里精品6国产| 午夜免费观看精品视频| 黄色av网站免费在线| 熟女妇女老妇一二三区| 97精品人妻一区二区三区精品| 一区二区熟女人妻视频| 亚洲国产香蕉视频在线播放| 日韩欧美在线观看不卡一区二区| 男人的天堂av日韩亚洲| 精品91自产拍在线观看一区| 水蜜桃一区二区三区在线观看视频| 一区二区在线视频中文字幕 | 狠狠躁夜夜躁人人爽天天天天97| 护士特殊服务久久久久久久| 亚洲乱码中文字幕在线| 青青草国内在线视频精选| 老熟妇凹凸淫老妇女av在线观看| 11久久久久久久久久久| 国产视频一区二区午夜| 亚洲图片欧美校园春色| 精品人人人妻人人玩日产欧| av中文字幕福利网| 亚洲国产免费av一区二区三区| 99热久久极品热亚洲| 日韩人妻在线视频免费| av黄色成人在线观看| 国产亚洲四十路五十路| 视频一区 视频二区 视频| 在线观看黄色成年人网站| 欧美成人猛片aaaaaaa| 红杏久久av人妻一区| 欧美一区二区三区乱码在线播放 | 国产精品午夜国产小视频| 中文字幕视频一区二区在线观看 | 亚洲最大免费在线观看| 人人超碰国字幕观看97| 玖玖一区二区在线观看| 综合精品久久久久97| av在线免费中文字幕| 天天干天天操天天玩天天射| av高潮迭起在线观看| 男女之间激情网午夜在线| 班长撕开乳罩揉我胸好爽| 欧美日韩激情啪啪啪| 一区二区久久成人网| 一区二区在线观看少妇| 亚洲高清国产一区二区三区| 亚洲 清纯 国产com| av乱码一区二区三区| 中文字幕人妻三级在线观看| 婷婷久久久久深爱网| 曰本无码人妻丰满熟妇啪啪| 91亚洲精品干熟女蜜桃频道| 大香蕉日本伊人中文在线| 亚洲高清国产自产av| 日韩精品中文字幕播放| 中国黄片视频一区91| 超碰公开大香蕉97| 午夜毛片不卡免费观看视频| 精内国产乱码久久久久久| 馒头大胆亚洲一区二区| 天堂av在线官网中文| 亚洲精品ww久久久久久| 久久久久久久99精品| 日韩午夜福利精品试看| 日韩美av高清在线| 久久精品国产999| 精品av久久久久久久| 日日摸夜夜添夜夜添毛片性色av| 国产亚洲视频在线二区| sw137 中文字幕 在线| 99一区二区在线观看| 一级黄片大鸡巴插入美女| 18禁美女无遮挡免费| caoporm超碰国产| 亚洲丝袜老师诱惑在线观看| av大全在线播放免费| 插逼视频双插洞国产操逼插洞| 日韩人妻在线视频免费| 这里有精品成人国产99| 亚洲 清纯 国产com| 免费成人va在线观看| 国产免费高清视频视频| 熟女妇女老妇一二三区| 最近中文字幕国产在线| 亚洲 国产 成人 在线| 亚洲精品国产综合久久久久久久久| 色av色婷婷人妻久久久精品高清| 青青草人人妻人人妻| 伊人综合aⅴ在线网| 中文字幕无码日韩专区免费| 93视频一区二区三区| 亚洲午夜电影之麻豆| 日韩欧美一级黄片亚洲| 国产黄色大片在线免费播放 | 午夜在线观看一区视频| 日日夜夜大香蕉伊人| 日韩国产乱码中文字幕| 成人av天堂丝袜在线观看 | 肏插流水妹子在线乐播下载| av天堂资源最新版在线看| 欧美精品免费aaaaaa| 欧美交性又色又爽又黄麻豆| 国产一区二区久久久裸臀| 欧美在线精品一区二区三区视频 | 日本熟妇色熟妇在线观看| 成人av久久精品一区二区| 福利国产视频在线观看| 午夜毛片不卡在线看| 成人蜜桃美臀九一一区二区三区| 男人和女人激情视频| 蜜臀av久久久久蜜臀av麻豆| 在线免费观看欧美小视频| 日韩美女搞黄视频免费| 亚洲精品国偷自产在线观看蜜桃| 蜜桃色婷婷久久久福利在线| 五月天色婷婷在线观看视频免费| 日本少妇人妻xxxxxhd| 天天色天天爱天天爽| 亚洲一区二区三区偷拍女厕91| 人妻丰满熟妇综合网| 视频二区在线视频观看| 国产成人午夜精品福利| 91片黄在线观看喷潮| 亚洲成人黄色一区二区三区| 啪啪啪啪啪啪啪啪av| 热久久只有这里有精品| 欧美一区二区三区久久久aaa| 最新国产精品网址在线观看| 四川乱子伦视频国产vip| 日本高清成人一区二区三区 | 亚洲综合在线视频可播放| 久久香蕉国产免费天天| 欧美美女人体视频一区| 青青草在观免费国产精品| 人妻久久久精品69系列| 日韩av大胆在线观看| 动漫av网站18禁| 久久久久久性虐视频| 欧美性感尤物人妻在线免费看| 肏插流水妹子在线乐播下载| 中文字幕乱码av资源| 91中文字幕免费在线观看| 久草视频中文字幕在线观看| 在线观看的黄色免费网站| 2021最新热播中文字幕| 亚洲美女高潮喷浆视频| 99久久激情婷婷综合五月天| 亚洲国产精品美女在线观看| 欧美伊人久久大香线蕉综合| 40道精品招牌菜特色| 亚洲一区二区久久久人妻| 日韩欧美国产精品91| 夜夜操,天天操,狠狠操| 中国熟女@视频91| 97超碰人人搞人人| 久久久久久久久久性潮| av老司机亚洲一区二区| 天天做天天爽夜夜做少妇| 亚洲高清免费在线观看视频| 不卡日韩av在线观看| 91九色porny国产在线| 男人天堂av天天操| 亚洲高清视频在线不卡| 最新中文字幕乱码在线| 超级福利视频在线观看| 日韩av有码一区二区三区4| 亚洲成人国产综合一区| 精品91高清在线观看| 国产美女精品福利在线| 97超碰人人搞人人| 国产精品视频一区在线播放| 亚洲成人激情av在线| 亚洲1069综合男同| 日韩美女搞黄视频免费| av手机在线观播放网站| 88成人免费av网站| 国产九色91在线视频| 在线观看欧美黄片一区二区三区| 把腿张开让我插进去视频| 大尺度激情四射网站| 免费观看成年人视频在线观看| 老鸭窝在线观看一区| 偷拍自拍 中文字幕| 欧美偷拍自拍色图片| mm131美女午夜爽爽爽| 非洲黑人一级特黄片| 亚洲卡1卡2卡三卡四老狼| 日韩欧美亚洲熟女人妻| 欧美精品 日韩国产| 国产高清97在线观看视频| 在线观看视频污一区| 日本黄在免费看视频| 日本高清在线不卡一区二区| 91九色国产porny蝌蚪| 中国老熟女偷拍第一页| 动漫精品视频在线观看| 超碰在线观看免费在线观看| 在线免费91激情四射| 亚洲精品亚洲人成在线导航| 视频二区在线视频观看| 福利在线视频网址导航| 免费在线黄色观看网站| 99久久激情婷婷综合五月天| 涩涩的视频在线观看视频| 人妻熟女在线一区二区| 亚洲精品成人网久久久久久小说| 91亚洲国产成人精品性色| 啊啊好慢点插舔我逼啊啊啊视频| 中文字幕之无码色多多| 国产精品人妻66p| 色秀欧美视频第一页| 一个色综合男人天堂| 大胆亚洲av日韩av| 97年大学生大白天操逼| 91麻豆精品传媒国产黄色片| 精品区一区二区三区四区人妻| 91九色porny蝌蚪国产成人| 最新97国产在线视频| 深田咏美亚洲一区二区| 天天日天天爽天天干| 大尺度激情四射网站| 国产欧美日韩第三页| 最新黄色av网站在线观看| 天堂av在线播放免费| 国产亚洲成人免费在线观看| 99精品视频在线观看免费播放| 亚洲蜜臀av一区二区三区九色| av无限看熟女人妻另类av| 亚洲推理片免费看网站| 天天做天天爽夜夜做少妇| 久青青草视频手机在线免费观看| 一二三中文乱码亚洲乱码one | 无码日韩人妻精品久久| 亚洲精品在线资源站| 亚洲免费国产在线日韩| 一区二区麻豆传媒黄片| 久久精品36亚洲精品束缚| 亚洲欧美一卡二卡三卡| 91精品免费久久久久久| 91啪国自产中文字幕在线| 农村胖女人操逼视频| 大陆精品一区二区三区久久| 天天操夜夜操天天操天天操| 欧美国产亚洲中英文字幕| 老司机在线精品福利视频| 337p日本大胆欧美人| 中文字幕乱码av资源| 国产成人自拍视频播放 | 大香蕉大香蕉在线看| 国产精品熟女久久久久浪潮| 亚洲国产欧美一区二区丝袜黑人| 午夜影院在线观看视频羞羞羞| 女同性ⅹxx女同h偷拍| 亚洲第一伊人天堂网| 中文字幕最新久久久| 日韩精品电影亚洲一区| 欧美黑人与人妻精品| 亚欧在线视频你懂的| 欧美国产亚洲中英文字幕| 亚洲 人妻 激情 中文| 边摸边做超爽毛片18禁色戒| 亚洲少妇人妻无码精品| 韩国三级aaaaa高清视频 | 国产美女精品福利在线| 成人国产激情自拍三区| 一区二区久久成人网| 精品老妇女久久9g国产| 自拍偷拍亚洲另类色图| 精品人妻一二三区久久| 黄色视频在线观看高清无码| 爱爱免费在线观看视频| 91国内精品自线在拍白富美| 日日操综合成人av| 免费福利av在线一区二区三区| 亚洲成人情色电影在线观看| 国产乱子伦精品视频潮优女| 亚洲人妻av毛片在线| 最新的中文字幕 亚洲| 2020av天堂网在线观看| 日本阿v视频在线免费观看| 自拍偷区二区三区麻豆| 偷拍美女一区二区三区| 97超碰国语国产97超碰| 男人靠女人的逼视频| 一区二区三区激情在线| 国产精品一区二区三区蜜臀av | 国产av欧美精品高潮网站| 天天干天天操天天插天天日| 国产va精品免费观看 | 精品一线二线三线日本| 黄片三级三级三级在线观看| 淫秽激情视频免费观看| 东游记中文字幕版哪里可以看到| 国产日本欧美亚洲精品视| 欧美偷拍自拍色图片| 久久久噜噜噜久久熟女av| 国产精品福利小视频a| 亚洲av无乱一区二区三区性色| 成人久久精品一区二区三区| 51精品视频免费在线观看| 99国内精品永久免费视频| av在线免费资源站| 精品人人人妻人人玩日产欧| 在线观看亚洲人成免费网址| 日日夜夜大香蕉伊人| 日本少妇在线视频大香蕉在线观看| 最新国产亚洲精品中文在线| 国产欧美日韩在线观看不卡| 亚洲免费va在线播放| 中文字幕亚洲久久久| 精品久久久久久久久久久久人妻| 日本乱人一区二区三区| 亚洲少妇人妻无码精品| 日本高清在线不卡一区二区| AV天堂一区二区免费试看| 色偷偷伊人大杳蕉综合网| 99热99这里精品6国产| 亚洲精品久久视频婷婷| 性色蜜臀av一区二区三区| 天天日夜夜操天天摸| 天天操天天爽天天干| 男女啪啪视频免费在线观看| 欧美视频中文一区二区三区| 特级无码毛片免费视频播放| 人人超碰国字幕观看97| 国产黄色片蝌蚪九色91| sw137 中文字幕 在线| 18禁污污污app下载| 欧洲黄页网免费观看| 又粗又硬又猛又爽又黄的| 成人色综合中文字幕| 欧美日韩在线精品一区二区三| 欧洲亚洲欧美日韩综合| japanese五十路熟女熟妇| 欧美日韩v中文在线| 男生舔女生逼逼视频| 国产视频网站国产视频| 亚洲美女美妇久久字幕组| 三上悠亚和黑人665番号| 人妻无码中文字幕专区| 香港三日本三韩国三欧美三级| 青青青aaaa免费| 国产欧美精品免费观看视频| 国产精品sm调教视频| 99热这里只有精品中文| 日韩欧美一级精品在线观看| 天堂av在线官网中文| 2022国产综合在线干| 国产av福利网址大全| 中文字幕之无码色多多| 日韩精品中文字幕福利| 欧美亚洲少妇福利视频| av日韩在线免费播放| 2021久久免费视频| 一区二区久久成人网| 啪啪啪啪啪啪啪啪av| 久草视频首页在线观看| 人人超碰国字幕观看97| 天干天天天色天天日天天射| 天天艹天天干天天操| 天天操天天干天天日狠狠插| 啪啪啪啪啪啪啪免费视频| 99久久久无码国产精品性出奶水 | av中文字幕电影在线看| 国产精品系列在线观看一区二区| 一区二区三区精品日本| 天天日天天干天天插舔舔| 99精品视频之69精品视频| 日韩av免费观看一区| 亚洲天堂精品福利成人av| 自拍偷拍 国产资源| 天天摸天天亲天天舔天天操天天爽 | 青青青青在线视频免费观看| 亚洲高清国产拍青青草原| 大鸡八强奸视频在线观看| 日韩黄色片在线观看网站| 黄色片一级美女黄色片| 亚洲天堂精品久久久| 亚洲欧美激情中文字幕| 亚洲福利精品视频在线免费观看| 又粗又硬又猛又黄免费30| 亚洲一区二区三区精品乱码| 亚洲欧美清纯唯美另类| 日本裸体熟妇区二区欧美| jul—619中文字幕在线| av中文字幕电影在线看| 国产老熟女伦老熟妇ⅹ| 粉嫩小穴流水视频在线观看| 亚洲午夜在线视频福利| 经典av尤物一区二区| 1024久久国产精品| 中文亚洲欧美日韩无线码| 亚洲高清视频在线不卡| 久久久久久久久久性潮| 9色精品视频在线观看| 亚洲偷自拍高清视频| 国产综合视频在线看片| 中文字幕高清在线免费播放| 成人sm视频在线观看| 日日爽天天干夜夜操| 亚洲一级 片内射视正片| 日韩a级黄色小视频| 亚洲精品色在线观看视频| 92福利视频午夜1000看 | 99热99这里精品6国产| 亚洲综合自拍视频一区| 91福利视频免费在线观看| 农村胖女人操逼视频| 中文亚洲欧美日韩无线码| 欧美日韩一区二区电影在线观看| 93人妻人人揉人人澡人人| 国产熟妇一区二区三区av| 亚洲 欧美 精品 激情 偷拍| 亚洲欧美综合另类13p| 韩国爱爱视频中文字幕| 1区2区3区4区视频在线观看| 午夜精品久久久久久99热| 国产视频精品资源网站| 四虎永久在线精品免费区二区| 精品国产成人亚洲午夜| 亚洲av午夜免费观看| 1000小视频在线| 亚洲嫩模一区二区三区| 含骚鸡巴玩逼逼视频| 色综合色综合色综合色| 国产美女精品福利在线| 欧美一级片免费在线成人观看| 国产精品久久久久久久精品视频| 自拍偷拍一区二区三区图片| caoporn蜜桃视频| 激情综合治理六月婷婷| 一区二区视频视频视频| 同居了嫂子在线播高清中文| 做爰视频毛片下载蜜桃视频1| 内射久久久久综合网| 国产男女视频在线播放| 在线观看亚洲人成免费网址| 91九色porny国产蝌蚪视频| 999九九久久久精品| av黄色成人在线观看| 亚国产成人精品久久久| 亚洲图片偷拍自拍区| 四川五十路熟女av| 香蕉av影视在线观看| 久久久久久久久久性潮| 欧亚日韩一区二区三区观看视频| 中文字幕第一页国产在线| 毛片一级完整版免费| 91chinese在线视频| 成人亚洲精品国产精品| 亚洲最大免费在线观看| 宅男噜噜噜666免费观看| 3344免费偷拍视频| 少妇人妻真实精品视频| 精品久久久久久高潮| 色哟哟在线网站入口| 久久精品国产999| 中文字幕在线免费第一页| 国产久久久精品毛片| 美味人妻2在线播放| 久久精品36亚洲精品束缚| 日本啪啪啪啪啪啪啪| 91chinese在线视频| 久久久久91精品推荐99| 国产91久久精品一区二区字幕| 小泽玛利亚视频在线观看| 影音先锋女人av噜噜色| 视频在线免费观看你懂得| 国产使劲操在线播放| 成人30分钟免费视频| 免费在线黄色观看网站| 国产日韩精品电影7777| 又色又爽又黄的美女裸体| 四川乱子伦视频国产vip| 国产男女视频在线播放| 农村胖女人操逼视频| 国产精品污污污久久| 女同性ⅹxx女同h偷拍| 91桃色成人网络在线观看| 欧美viboss性丰满| 中文字幕在线观看国产片| 黄色av网站免费在线| 91中文字幕免费在线观看| 91试看福利一分钟| 蜜桃精品久久久一区二区| 又色又爽又黄又刺激av网站| 亚洲老熟妇日本老妇| 成年人黄色片免费网站| 性感美女福利视频网站| 亚洲在线观看中文字幕av| 欧洲亚洲欧美日韩综合| 巨乳人妻日下部加奈被邻居中出| 国产黑丝高跟鞋视频在线播放 | 91社福利《在线观看| av破解版在线观看| nagger可以指黑人吗| 免费国产性生活视频| 国产精品国产三级国产午| 人妻少妇中文有码精品| 人妻久久无码中文成人| 一级黄片大鸡巴插入美女| 啊啊啊想要被插进去视频| 欲满人妻中文字幕在线| 2020久久躁狠狠躁夜夜躁 | 久久香蕉国产免费天天| 日本黄色三级高清视频| 成年午夜影片国产片| 国产午夜男女爽爽爽爽爽视频| 玖玖一区二区在线观看| 2020av天堂网在线观看| 国产激情av网站在线观看| 日韩欧美一级黄片亚洲| 日韩精品二区一区久久| 在线观看黄色成年人网站| 久久午夜夜伦痒痒想咳嗽P| 人妻丝袜诱惑我操她视频| 亚洲少妇高潮免费观看| 2o22av在线视频| 9色精品视频在线观看| 欧洲黄页网免费观看| 日本特级片中文字幕| 玖玖一区二区在线观看| 一区二区三区视频,福利一区二区| 在线观看的黄色免费网站| 福利午夜视频在线观看| 男人在床上插女人视频| weyvv5国产成人精品的视频| 蜜桃色婷婷久久久福利在线| 精品高潮呻吟久久av| 国产精品久久9999| 精品一区二区三区欧美| 久久午夜夜伦痒痒想咳嗽P| 一色桃子久久精品亚洲| 日韩欧美制服诱惑一区在线| 中文字幕欧美日韩射射一| 欧美男人大鸡吧插女人视频| 午夜dv内射一区区| 激情五月婷婷综合色啪| 天天艹天天干天天操| 欧洲日韩亚洲一区二区三区| 日韩欧美一级精品在线观看| 熟妇一区二区三区高清版| 国产不卡av在线免费| 国产亚州色婷婷久久99精品| 熟女在线视频一区二区三区| 91社福利《在线观看| 888欧美视频在线| 99re国产在线精品| 国产日韩欧美视频在线导航| 91福利视频免费在线观看| 欧美一区二区中文字幕电影| 成人精品在线观看视频| h国产小视频福利在线观看| 九一传媒制片厂视频在线免费观看 | 亚洲av可乐操首页| 国产中文精品在线观看| 中国无遮挡白丝袜二区精品 | 一区二区在线观看少妇| 午夜在线精品偷拍一区二| 538精品在线观看视频| 欧美日韩不卡一区不区二区| 青青草人人妻人人妻| 91精品国产综合久久久蜜| 日本人竟这样玩学生妹| 姐姐的朋友2在线观看中文字幕| 1024久久国产精品| 888亚洲欧美国产va在线播放| 绯色av蜜臀vs少妇| 亚洲一级 片内射视正片| 久久这里只有精品热视频| 狠狠鲁狠狠操天天晚上干干| 午夜成午夜成年片在线观看| 日本成人不卡一区二区| 国产高清女主播在线| 日本又色又爽又黄又粗| 欧美视频中文一区二区三区| 夜夜嗨av蜜臀av| 粉嫩欧美美人妻小视频| 国产黄色高清资源在线免费观看| 黄页网视频在线免费观看| 欧美一区二区三区乱码在线播放 | 加勒比视频在线免费观看| 亚洲高清国产自产av| 成年人啪啪视频在线观看| 欧美日本在线视频一区| 国产超码片内射在线| 午夜蜜桃一区二区三区| 国产午夜激情福利小视频在线| 成年女人免费播放视频| 丰满的继坶3中文在线观看| 亚洲男人让女人爽的视频| 欧美成人黄片一区二区三区| 六月婷婷激情一区二区三区| 日韩熟女系列一区二区三区| 日本少妇精品免费视频| av俺也去在线播放| 骚逼被大屌狂草视频免费看| av亚洲中文天堂字幕网| 日本少妇人妻xxxxx18| 免费观看成年人视频在线观看| 亚洲av成人网在线观看| 少妇人妻100系列| 天天综合天天综合天天网| 自拍偷区二区三区麻豆| 日韩av免费观看一区| 亚洲激情偷拍一区二区| 亚洲一区二区人妻av| 黄色黄色黄片78在线| 漂亮 人妻被中出中文| av在线免费观看亚洲天堂| 久久久久久久一区二区三| 亚洲一级 片内射视正片| av天堂中文字幕最新| 天天干夜夜操啊啊啊| 97资源人妻免费在线视频| 亚洲成人激情视频免费观看了| 91久久国产成人免费网站| 成人性爱在线看四区| 免费在线观看视频啪啪| 不卡日韩av在线观看| 99精品国产aⅴ在线观看| 久草视频在线看免费| 在线观看的a站 最新| 绝色少妇高潮3在线观看| 国产在线一区二区三区麻酥酥| 一二三区在线观看视频| 免费在线观看视频啪啪| av在线资源中文字幕| 操操网操操伊剧情片中文字幕网| 亚洲国产成人最新资源| 男人的天堂一区二区在线观看| av日韩在线观看大全| 中文字幕在线欧美精品| 亚洲精品国产在线电影| 亚洲在线一区二区欧美| 国产精品成人xxxx| 天天操,天天干,天天射| 亚洲欧美日韩视频免费观看| 免费69视频在线看| 成人在线欧美日韩国产| 福利午夜视频在线合集| 欧美日韩熟女一区二区三区| 青青草亚洲国产精品视频| 久久久久久久99精品| 91精品资源免费观看| 黑人借宿ntr人妻的沦陷2| 亚洲国产精品黑丝美女| av老司机精品在线观看| 夜色福利视频在线观看| 中文字幕乱码av资源| 熟女人妻一区二区精品视频| 在线成人日韩av电影| 五十路在线观看完整版| 精品老妇女久久9g国产| 精品亚洲在线免费观看| 亚洲人妻av毛片在线| 美女日逼视频免费观看| 蜜桃色婷婷久久久福利在线| 九色porny九色9l自拍视频| 黄色成人在线中文字幕| 99久久成人日韩欧美精品| 无码中文字幕波多野不卡| 亚洲精品午夜久久久久| 欧美视频综合第一页| 国产又粗又硬又大视频| 香蕉片在线观看av| 中文字幕第1页av一天堂网| 欧美精产国品一二三区| 亚洲一区自拍高清免费视频| 青青青国产片免费观看视频| 在线免费观看靠比视频的网站| 97香蕉碰碰人妻国产樱花| 亚洲熟妇x久久av久久| 午夜青青草原网在线观看| 久久久久久久久久一区二区三区| 国产品国产三级国产普通话三级| 91国内视频在线观看| 视频二区在线视频观看| 久久久久久国产精品| 福利国产视频在线观看| 一区二区三区视频,福利一区二区| 精品老妇女久久9g国产| 最新欧美一二三视频| 国产janese在线播放| 欧美精品黑人性xxxx| 中文字幕高清免费在线人妻| 亚洲成av人无码不卡影片一| 2019av在线视频| 97精品视频在线观看| 亚洲欧美激情国产综合久久久| 天天色天天操天天透| 亚洲激情偷拍一区二区| 啪啪啪啪啪啪啪免费视频| 亚洲高清一区二区三区视频在线| 欧美视频一区免费在线| 成人18禁网站在线播放| 欧美视频一区免费在线| 在线免费91激情四射 | 都市激情校园春色狠狠| av天堂中文字幕最新| 首之国产AV医生和护士小芳| 天天爽夜夜爽人人爽QC| 国产极品美女久久久久久| 国产精品黄色的av| 日本女人一级免费片| 日韩特级黄片高清在线看| 国产视频一区二区午夜| 一区二区三区av高清免费| 人妻少妇亚洲精品中文字幕| aiss午夜免费视频| 国产自拍在线观看成人| 日本韩国免费一区二区三区视频| 男生舔女生逼逼的视频| 亚欧在线视频你懂的| 91精品国产综合久久久蜜| 小穴多水久久精品免费看| 亚洲国产精品久久久久久6| 亚洲成人av一区在线| 扒开让我视频在线观看| 国产卡一卡二卡三乱码手机| 亚洲国产在线精品国偷产拍| 天天摸天天亲天天舔天天操天天爽 | 1区2区3区4区视频在线观看| 丰满的子国产在线观看| sw137 中文字幕 在线| 国产精品久久综合久久| 中文字幕在线乱码一区二区| 2021国产一区二区| 在线观看免费视频网| 亚洲免费福利一区二区三区| 手机看片福利盒子日韩在线播放| 99国产精品窥熟女精品| 青青草人人妻人人妻| 黄色大片免费观看网站| 制服丝袜在线人妻中文字幕| 老司机99精品视频在线观看| 中文字幕日韩精品日本| 岳太深了紧紧的中文字幕| 综合精品久久久久97| 一色桃子人妻一区二区三区| 日韩欧美在线观看不卡一区二区 | 超碰公开大香蕉97| 国产一区二区久久久裸臀| 一区二区三区蜜臀在线| 亚洲一区二区三区偷拍女厕91| 大陆av手机在线观看| av大全在线播放免费| 中文字幕人妻三级在线观看| 亚洲日本一区二区久久久精品| 韩国爱爱视频中文字幕| 国产剧情演绎系列丝袜高跟| 熟女妇女老妇一二三区| 午夜福利人人妻人人澡人人爽| 晚上一个人看操B片| 1024久久国产精品| 亚洲国产成人最新资源| 婷婷色国产黑丝少妇勾搭AV| 国内自拍第一页在线观看| 老司机免费视频网站在线看| 一区二区视频在线观看视频在线| 日韩二区视频一线天婷婷五| 日美女屁股黄邑视频| 亚洲精品乱码久久久久久密桃明| 97国产精品97久久| 91国语爽死我了不卡| rct470中文字幕在线| 亚洲图片偷拍自拍区| 日本精品一区二区三区在线视频。| 美女福利视频导航网站| 日本少妇在线视频大香蕉在线观看| av线天堂在线观看| 首之国产AV医生和护士小芳| 中国产一级黄片免费视频播放| 在线制服丝袜中文字幕| 三上悠亚和黑人665番号| 国产精品精品精品999| 免费在线福利小视频| 日韩精品中文字幕播放| 激情人妻校园春色亚洲欧美| 亚洲丝袜老师诱惑在线观看| 亚洲中文字幕乱码区| 欧美一级视频一区二区| 2020韩国午夜女主播在线| av成人在线观看一区| 国产福利小视频二区| 在线观看免费视频色97| 一区二区三区欧美日韩高清播放| 人妻丝袜榨强中文字幕| 色吉吉影音天天干天天操 | 天天日天天鲁天天操| 午夜毛片不卡免费观看视频| 18禁美女黄网站色大片下载| av中文字幕福利网| 日韩av大胆在线观看| 国产丰满熟女成人视频| 国产美女一区在线观看| 久久精品国产亚洲精品166m| 熟女91pooyn熟女| 亚洲图库另类图片区| 国产又粗又黄又硬又爽| 亚洲av色图18p| 人妻久久久精品69系列| 国产成人精品一区在线观看| 国产日韩一区二区在线看| 国产又色又刺激在线视频| 亚洲人人妻一区二区三区| 亚洲欧美自拍另类图片| 国产成人一区二区三区电影网站| 在线观看免费视频网| 国产精品欧美日韩区二区| 日韩中文字幕精品淫| 免费黄色成人午夜在线网站| 蜜桃色婷婷久久久福利在线| 国产乱弄免费视频观看| 欧美在线精品一区二区三区视频| 中文字幕在线乱码一区二区| 黑人变态深video特大巨大| 国产亚洲国产av网站在线| 欧美va亚洲va天堂va| 99久久中文字幕一本人| 国产精品久久久久久久女人18| 在线观看免费视频色97| 亚洲一级特黄特黄黄色录像片| 91麻豆精品秘密入口在线观看 | 男女第一次视频在线观看| 偷偷玩弄新婚人妻h视频| 亚洲国产第一页在线观看| 国语对白xxxx乱大交| 美女 午夜 在线视频| 91精品激情五月婷婷在线| 久久久91蜜桃精品ad| tube69日本少妇| 99热99re在线播放| 搡老妇人老女人老熟女| 国产亚洲精品视频合集| 97超碰免费在线视频| 亚洲欧美成人综合视频| 色爱av一区二区三区| 偷拍自拍视频图片免费| 欧美亚洲自偷自拍 在线| 韩国爱爱视频中文字幕| 天天日天天干天天要| 中文字幕亚洲中文字幕| 青青草精品在线视频观看| 天天躁日日躁狠狠躁躁欧美av | 黄片大全在线观看观看| 欧洲欧美日韩国产在线| 亚洲最大免费在线观看| 操操网操操伊剧情片中文字幕网| 91精品啪在线免费| 亚洲图库另类图片区| 欧美在线精品一区二区三区视频| 成人精品在线观看视频| 日韩av熟妇在线观看| 久久麻豆亚洲精品av| 91免费放福利在线观看| 国产亚洲欧美45p| 最新欧美一二三视频| 老司机99精品视频在线观看| 国产高清精品一区二区三区| 中文字幕中文字幕人妻| 亚洲老熟妇日本老妇| 操人妻嗷嗷叫视频一区二区| 亚洲欧美福利在线观看| 好吊操视频这里只有精品| 91成人在线观看免费视频| 久久机热/这里只有| 人人爱人人妻人人澡39| 欧美日韩v中文在线| 91极品大一女神正在播放| 日本韩国亚洲综合日韩欧美国产| 国产精品人妻熟女毛片av久| 亚洲的电影一区二区三区| 国产janese在线播放| 11久久久久久久久久久| 少妇人妻久久久久视频黄片| 精品高潮呻吟久久av| 中国视频一区二区三区| 中文乱理伦片在线观看| 黄色男人的天堂视频| 在线国产精品一区二区三区| 国产精品福利小视频a| 人人妻人人爽人人澡人人精品| 亚洲高清自偷揄拍自拍| 欧美专区日韩专区国产专区| 99热久久极品热亚洲| 亚洲中文精品人人免费| 五十路熟女人妻一区二区9933| 11久久久久久久久久久| 黄片大全在线观看观看| 国产黄色a级三级三级三级| 亚洲少妇高潮免费观看| 大鸡吧插逼逼视频免费看| 97人妻人人澡爽人人精品| 日本免费视频午夜福利视频| 夫妻在线观看视频91| 2022天天干天天操| 大香蕉伊人国产在线| 免费啪啪啪在线观看视频| 久久久精品999精品日本| 亚洲成人情色电影在线观看| 中文字幕日韩精品就在这里| 国产日韩欧美美利坚蜜臀懂色| 热久久只有这里有精品| www日韩毛片av| 福利在线视频网址导航| 91she九色精品国产| 亚洲综合一区成人在线| www天堂在线久久| 2o22av在线视频| 国产精品手机在线看片| 午夜精品福利一区二区三区p| 成熟丰满熟妇高潮xx×xx | 蜜臀av久久久久蜜臀av麻豆| 青青热久免费精品视频在线观看| 天天干天天操天天摸天天射| 9国产精品久久久久老师| 阴茎插到阴道里面的视频| 91破解版永久免费| 中文字幕欧美日韩射射一| 天天综合天天综合天天网| 2o22av在线视频| 粗大的内捧猛烈进出爽大牛汉子| 一区二区三区四区视频在线播放| 好太好爽好想要免费| 男人插女人视频网站| 午夜极品美女福利视频| 欧美天堂av无线av欧美| 最新国产精品网址在线观看| 国产+亚洲+欧美+另类| 亚洲欧美综合另类13p| 国产精品手机在线看片| 久草视频 久草视频2| 男人的天堂在线黄色| 欧美日韩人妻久久精品高清国产| 特黄老太婆aa毛毛片| 亚洲av日韩精品久久久| 日本美女性生活一级片| 18禁网站一区二区三区四区| 一区二区三区四区中文| 99亚洲美女一区二区三区| 国产精品自拍视频大全| 天天躁日日躁狠狠躁av麻豆| 国产女人露脸高潮对白视频| 人妻自拍视频中国大陆| 欧美成人猛片aaaaaaa| 亚洲区欧美区另类最新章节| 91色九色porny| 国产精品熟女久久久久浪潮| 黄色片年轻人在线观看| 人人妻人人爽人人添夜| 伊人日日日草夜夜草| 国产使劲操在线播放| 偷拍自拍亚洲美腿丝袜| av网址在线播放大全| 超碰中文字幕免费观看| 老司机99精品视频在线观看 | 日本少妇在线视频大香蕉在线观看 | 欧洲亚洲欧美日韩综合| 亚洲 自拍 色综合图| 99热99re在线播放| 亚洲最大免费在线观看| xxx日本hd高清| 精品久久久久久高潮| 大鸡巴后入爆操大屁股美女| 欧洲日韩亚洲一区二区三区| 自拍偷拍 国产资源| 国产在线观看黄色视频| 久久h视频在线观看| 五月激情婷婷久久综合网| 亚洲欧美一区二区三区爱爱动图| 国产亚洲欧美45p| 一区二区三区综合视频| 欧美精品中文字幕久久二区| 日日爽天天干夜夜操| 护士特殊服务久久久久久久| 久久久久久久99精品| 欧美日韩不卡一区不区二区| 亚洲 色图 偷拍 欧美| 青青青爽视频在线播放| 青青青青在线视频免费观看| 极品性荡少妇一区二区色欲| 女生自摸在线观看一区二区三区 | 国产精品视频资源在线播放| 亚洲成av人无码不卡影片一| 国产黄色a级三级三级三级 | 天天日天天天天天天天天天天| 日本高清撒尿pissing| 日视频免费在线观看| 国产露脸对白在线观看| 亚洲一区久久免费视频| 熟女视频一区,二区,三区| 天天操天天射天天操天天天| 欧美另类一区二区视频| 亚洲熟妇无码一区二区三区| 又粗又硬又猛又爽又黄的| 91福利在线视频免费观看| 一级黄色片夫妻性生活| 天天日天天鲁天天操| 福利视频广场一区二区| 久草电影免费在线观看| 国产成人精品久久二区91| 久草视频在线看免费| 女同性ⅹxx女同h偷拍| 999热精品视频在线| 国产片免费观看在线观看| 日本黄色三级高清视频| 欧美日韩一区二区电影在线观看| 激情图片日韩欧美人妻| 欧洲国产成人精品91铁牛tv| 在线网站你懂得老司机| 国产午夜亚洲精品不卡在线观看| 日韩少妇人妻精品无码专区| 中文字幕第一页国产在线| 综合页自拍视频在线播放| 韩国AV无码不卡在线播放| 欧洲亚洲欧美日韩综合| av中文字幕福利网| 在线播放国产黄色av| 欧美女同性恋免费a| 神马午夜在线观看视频| 亚洲的电影一区二区三区| 一区二区三区久久久91| 亚洲va国产va欧美精品88| 欧美国品一二三产区区别| 91在线免费观看成人| 国产视频在线视频播放| 国产三级影院在线观看| 日韩av有码一区二区三区4| 男人的天堂一区二区在线观看| 亚洲精品一区二区三区老狼| 精品国产高潮中文字幕| 丝袜肉丝一区二区三区四区在线看| 神马午夜在线观看视频| 日本免费视频午夜福利视频| 中文字幕+中文字幕| 午夜精品九一唐人麻豆嫩草成人 | 加勒比视频在线免费观看| 美女大bxxxx内射| 青青青青操在线观看免费| 亚洲av日韩av第一区二区三区| 18禁无翼鸟成人在线| 97超碰免费在线视频| 91精品国产91久久自产久强| 2012中文字幕在线高清| 精品久久久久久高潮| 夜色撩人久久7777| 天天日天天干天天要| 亚洲男人的天堂a在线| 久久热这里这里只有精品| 11久久久久久久久久久| 岳太深了紧紧的中文字幕| 国产精品黄色的av| 亚洲免费视频欧洲免费视频| 亚洲中文字字幕乱码| 久久这里只有精品热视频 | 色爱av一区二区三区| 国产成人精品av网站| 欧美另类z0z变态| 日韩成人综艺在线播放| 亚洲欧美一卡二卡三卡| 国产精品探花熟女在线观看| 国产视频网站国产视频| 久草视频中文字幕在线观看| 九一传媒制片厂视频在线免费观看| 久久久久国产成人精品亚洲午夜| 在线视频自拍第三页| av视屏免费在线播放| 中文字幕成人日韩欧美| 一区二区熟女人妻视频| 午夜青青草原网在线观看| 成人av天堂丝袜在线观看| 亚洲精品乱码久久久本| 免费在线福利小视频| 丰满少妇人妻xxxxx| 亚洲中文精品字幕在线观看| 欧亚日韩一区二区三区观看视频| 无码精品一区二区三区人| 色婷婷综合激情五月免费观看| 硬鸡巴动态操女人逼视频| 姐姐的朋友2在线观看中文字幕 | 久久香蕉国产免费天天| 国产精选一区在线播放| 大香蕉玖玖一区2区| 欧美黄片精彩在线免费观看| 天天日天天干天天爱| 麻豆性色视频在线观看| 亚洲中文精品字幕在线观看| 欧亚日韩一区二区三区观看视频 | 自拍偷区二区三区麻豆| 国产97在线视频观看| 亚洲在线一区二区欧美| 国产亚洲精品品视频在线| 国产刺激激情美女网站| 91p0rny九色露脸熟女| 亚洲一级特黄特黄黄色录像片| 国产高清在线观看1区2区| 一级a看免费观看网站| 激情五月婷婷综合色啪| 婷婷久久久久深爱网| 日本少妇人妻xxxxxhd| 国产麻豆国语对白露脸剧情| 动漫av网站18禁| 国产精品视频资源在线播放| 欧美亚洲一二三区蜜臀| 日本韩国亚洲综合日韩欧美国产 | 亚洲第一伊人天堂网| 中文字幕在线一区精品| 可以免费看的www视频你懂的| 亚洲av人人澡人人爽人人爱| 蜜桃臀av蜜桃臀av| aiss午夜免费视频| 特一级特级黄色网片| 精品人人人妻人人玩日产欧| 91免费放福利在线观看| 成人24小时免费视频| 午夜频道成人在线91| 激情五月婷婷免费视频| 青娱乐极品视频青青草| 成人av电影免费版| 人人爱人人妻人人澡39| av老司机精品在线观看| 国产精品三级三级三级| 亚洲va国产va欧美精品88| 国产黄色大片在线免费播放| 日本av高清免费网站| 影音先锋女人av噜噜色| 精品人妻伦一二三区久| 精品高潮呻吟久久av| 青青草视频手机免费在线观看| 精品欧美一区二区vr在线观看| 亚洲黄色av网站免费播放| 97精品人妻一区二区三区精品| 国产刺激激情美女网站| 国产精品一二三不卡带免费视频| 亚洲av极品精品在线观看| aiss午夜免费视频| 欧美日韩人妻久久精品高清国产| 国产女人被做到高潮免费视频| 久久久久久久久久久久久97| 国产精品免费不卡av| 亚洲综合乱码一区二区| 91综合久久亚洲综合| 色伦色伦777国产精品| 99精品一区二区三区的区| 视频在线亚洲一区二区| 美女大bxxxx内射| 中文字幕乱码av资源| 9久在线视频只有精品| 老司机在线精品福利视频| 韩国三级aaaaa高清视频| 日韩av熟妇在线观看| 蜜桃视频入口久久久| 国产高潮无码喷水AV片在线观看| 一区二区视频视频视频| 无码精品一区二区三区人| 精彩视频99免费在线| 天天夜天天日天天日| 五十路人妻熟女av一区二区| 青青青青视频在线播放| 青青青青青青青在线播放视频| 亚洲国产美女一区二区三区软件 | 啪啪啪18禁一区二区三区| 啪啪啪18禁一区二区三区| 欧美专区第八页一区在线播放| 91综合久久亚洲综合| 人人爽亚洲av人人爽av| 国产又粗又猛又爽又黄的视频美国| 天堂va蜜桃一区入口| 爱爱免费在线观看视频| 国产激情av网站在线观看| 亚洲推理片免费看网站| 精品高跟鞋丝袜一区二区| 少妇人妻真实精品视频| 欧美xxx成人在线| 91精品国产高清自在线看香蕉网| 中文字幕国产专区欧美激情| 亚洲va国产va欧美va在线| 国产亚洲视频在线二区| av老司机亚洲一区二区| 亚洲综合乱码一区二区| 亚洲熟女久久久36d| 精品一区二区亚洲欧美| 中文字幕 亚洲av| 国产一线二线三线的区别在哪| 老有所依在线观看完整版| 99久久成人日韩欧美精品| 亚洲中文字字幕乱码| 97超碰最新免费在线观看| 日本www中文字幕| 2o22av在线视频| 97a片免费在线观看| 玩弄人妻熟妇性色av少妇| huangse网站在线观看| 丁香花免费在线观看中文字幕 | 国产av福利网址大全| 4个黑人操素人视频网站精品91 | 久久久久久9999久久久久| 最新的中文字幕 亚洲| 日视频免费在线观看| 9国产精品久久久久老师| 2022中文字幕在线| 国产久久久精品毛片| 成年人啪啪视频在线观看| 成年午夜影片国产片| 丝袜国产专区在线观看| 精品久久婷婷免费视频| 青青青青青免费视频| 成人资源在线观看免费官网| 天天操天天射天天操天天天| 1区2区3区不卡视频| 夜色17s精品人妻熟女| 欧美精品亚洲精品日韩在线| 中文字幕欧美日韩射射一| 中文字幕av第1页中文字幕| 精彩视频99免费在线| 中文字幕 人妻精品| 99精品一区二区三区的区| 中文字幕日韩91人妻在线| 久草极品美女视频在线观看| 91片黄在线观看喷潮| 女同性ⅹxx女同h偷拍| 亚洲精品国品乱码久久久久| 激情伦理欧美日韩中文字幕| 欧美男同性恋69视频| 直接观看免费黄网站| 91人妻精品久久久久久久网站| 成年人午夜黄片视频资源| 中文字幕一区二区亚洲一区| 青青青青操在线观看免费| 亚洲国产欧美一区二区三区久久 | 女生被男生插的视频网站| 日韩特级黄片高清在线看| 人妻av无码专区久久绿巨人| 91精品激情五月婷婷在线| 亚洲欧美一区二区三区电影| 99亚洲美女一区二区三区| 日本女人一级免费片| 青青青青青操视频在线观看| 欧美一区二区三区久久久aaa| 亚洲精品 日韩电影| 中文字幕国产专区欧美激情| 日韩a级精品一区二区| 国产高潮无码喷水AV片在线观看| 免费一级黄色av网站| 国产女人露脸高潮对白视频| 国产中文精品在线观看| 久草视频首页在线观看| 天天日天天干天天舔天天射| 国产亚洲四十路五十路| 欧美一区二区三区激情啪啪啪 | 中英文字幕av一区| 青青青青青手机视频| 亚洲av无乱一区二区三区性色| 在线观看的黄色免费网站| 亚洲国产欧美一区二区丝袜黑人 | 久久久人妻一区二区| 国产视频网站国产视频| 色狠狠av线不卡香蕉一区二区| 日韩写真福利视频在线观看| 五十路熟女av天堂| 亚洲av色图18p| 婷婷色中文亚洲网68| 视频 一区二区在线观看| 五月色婷婷综合开心网4438| 亚洲美女自偷自拍11页| 国产又色又刺激在线视频| 五月天中文字幕内射| 国产亚洲精品品视频在线| 色哟哟在线网站入口| 四川五十路熟女av| av天堂加勒比在线| 精品久久久久久久久久久久人妻 | 国产午夜男女爽爽爽爽爽视频| 亚洲午夜电影之麻豆| 中文字幕在线乱码一区二区| 青青青青爽手机在线| 亚洲va欧美va人人爽3p| 亚洲成人av一区在线| 日本av熟女在线视频| 日本三极片视频网站观看| 神马午夜在线观看视频| 午夜av一区二区三区| 999久久久久999| 亚洲少妇人妻无码精品| av资源中文字幕在线观看| 亚洲一区二区三区久久受 | 91成人精品亚洲国产| aiss午夜免费视频| 专门看国产熟妇的网站| 91福利在线视频免费观看| 清纯美女在线观看国产| 福利片区一区二体验区| 夜夜嗨av一区二区三区中文字幕| 久久久精品999精品日本 | 成人资源在线观看免费官网| 韩国男女黄色在线观看| 一个色综合男人天堂| 色呦呦视频在线观看视频| 一区二区三区精品日本| 精品久久久久久久久久中文蒉| 日本高清在线不卡一区二区| 91大屁股国产一区二区| 在线观看国产网站资源| 亚洲高清视频在线不卡| jiuse91九色视频| 精品国产午夜视频一区二区| 午夜美女福利小视频| 超碰97免费人妻麻豆| 免费看国产av网站| 岛国免费大片在线观看| 伊人精品福利综合导航| 国产高清在线观看1区2区| 日本黄色三级高清视频| 欧美视频不卡一区四区| 骚逼被大屌狂草视频免费看| 在线观看视频 你懂的| 激情综合治理六月婷婷| 亚洲激情唯美亚洲激情图片| 人人妻人人人操人人人爽| 超级碰碰在线视频免费观看| 欧洲黄页网免费观看| 日本一二三中文字幕| 521精品视频在线观看| 亚洲蜜臀av一区二区三区九色| 午夜大尺度无码福利视频| 伊人精品福利综合导航| 精品黑人巨大在线一区| 亚洲av日韩精品久久久久久hd| 大香蕉伊人国产在线| 日韩精品中文字幕播放| 日本人妻少妇18—xx| 久久久久久久99精品| 日韩中文字幕在线播放第二页 | 非洲黑人一级特黄片| 人妻av无码专区久久绿巨人| 熟女俱乐部一二三区| 一区二区三区的久久的蜜桃的视频| 操操网操操伊剧情片中文字幕网| 涩爱综合久久五月蜜臀| 五十路在线观看完整版| 超碰公开大香蕉97| 中文字幕av熟女人妻| 适合午夜一个人看的视频| 小穴多水久久精品免费看| 日韩熟女系列一区二区三区| 91精品激情五月婷婷在线| 18禁美女黄网站色大片下载| 成人免费公开视频无毒 | 在线观看av观看av| 一区二区三区日本伦理| 亚洲精品精品国产综合| 91色老99久久九九爱精品| 巨乳人妻日下部加奈被邻居中出 | 欧美日韩高清午夜蜜桃大香蕉| 午夜大尺度无码福利视频| av视网站在线观看| 91she九色精品国产| 亚洲精品乱码久久久本|