日韩av黄I国产麻豆传媒I国产91av视频在线观看I日韩一区二区三区在线看I美女国产在线I麻豆视频国产在线观看I成人黄色短片

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

ZIP压缩算法详细分析及解压实例

發布時間:2023/12/8 编程问答 45 豆豆
生活随笔 收集整理的這篇文章主要介紹了 ZIP压缩算法详细分析及解压实例 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

最近自己實現了一個ZIP壓縮數據的解壓程序,覺得有必要把ZIP壓縮格式進行一下詳細總結,數據壓縮是一門通信原理和計算機科學都會涉及到的學科,在通信原理中,一般稱為信源編碼,在計算機科學里,一般稱為數據壓縮,兩者本質上沒啥區別,在數學家看來,都是映射。一方面在進行通信的時候,有必要將待傳輸的數據進行壓縮,以減少帶寬需求;另一方面,計算機存儲數據的時候,為了減少磁盤容量需求,也會將文件進行壓縮,盡管現在的網絡帶寬越來越高,壓縮已經不像90年代初那個時候那么迫切,但在很多場合下仍然需要,其中一個原因是壓縮后的數據容量減小后,磁盤訪問IO的時間也縮短,盡管壓縮和解壓縮過程會消耗CPU資源,但是CPU計算資源增長得很快,但是磁盤IO資源卻變化得很慢,比如目前主流的SATA硬盤仍然是7200轉,如果把磁盤的IO壓力轉化到CPU上,總體上能夠提升系統運行速度。壓縮作為一種非常典型的技術,會應用到很多很多場合下,比如文件系統、數據庫、消息傳輸、網頁傳輸等等各類場合。盡管壓縮里面會涉及到很多術語和技術,但無需擔心,博主盡量將其描述得通俗易懂。另外,本文涉及的壓縮算法非常主流并且十分精巧,理解了ZIP的壓縮過程,對理解其它相關的壓縮算法應該就比較容易了。

?

1、引子

壓縮可以分為無損壓縮和有損壓縮,有損,指的是壓縮之后就無法完整還原原始信息,但是壓縮率可以很高,主要應用于視頻、話音等數據的壓縮,因為損失了一點信息,人是很難察覺的,或者說,也沒必要那么清晰照樣可以看可以聽;無損壓縮則用于文件等等必須完整還原信息的場合,ZIP自然就是一種無損壓縮,在通信原理中介紹數據壓縮的時候,往往是從信息論的角度出發,引出香農所定義的熵的概念,這方面的介紹實在太多,這里換一種思路,從最原始的思想出發,為了達到壓縮的目的,需要怎么去設計算法。而ZIP為我們提供了相當好的案例。

盡管我們不去探討信息論里面那些復雜的概念,不過我們首先還是要從兩位信息論大牛談起。因為是他們奠基了今天大多數無損數據壓縮的核心,包括ZIP、RAR、GZIP、GIF、PNG等等大部分無損壓縮格式。這兩位大牛的名字分別是Jacob Ziv和Abraham Lempel,是兩位以色列人,在1977年的時候發表了一篇論文《A Universal Algorithm for Sequential Data Compression》,從名字可以看出,這是一種通用壓縮算法,所謂通用壓縮算法,指的是這種壓縮算法沒有對數據的類型有什么限定。不過論文我覺得不用仔細看了,因為博主作為一名通信專業的PHD,看起來也焦頭爛額,不過我們后面可以看到,它的思想還是很簡單的,之所以看起來復雜,主要是因為IEEE的某些雜志就是這個特點,需要從數學上去證明,這種壓縮算法到底有多優,比如針對一個各態歷經的隨機序列(不用追究什么叫各態歷經隨機序列),經過這樣的壓縮算法后,是否可以接近信息論里面的極限(也就是前面說的熵的概念)等等,不過在理解其思想之前,個人認為沒必要深究這些東西,除非你要發論文。這兩位大牛提出的這個算法稱為LZ77,兩位大牛過了一年又提了一個類似的算法,稱為LZ78,思想類似,ZIP這個算法就是基于LZ77的思想演變過來的,但ZIP對LZ77編碼之后的結果又繼續進行壓縮,直到難以壓縮為止。除了LZ77、LZ78,還有很多變種的算法,基本都以LZ開頭,如LZW、LZO、LZMA、LZSS、LZR、LZB、LZH、LZC、LZT、LZMW、LZJ、LZFG等等,非常多,LZW也比較流行,GIF那個動畫格式記得用了LZW。我也寫過解碼程序,以后有時間可以再寫一篇,但感覺跟LZ77這些類似,寫的必要性不大。

ZIP的作者是一個叫Phil Katz的人,這個人算是開源界的一個具有悲劇色彩的傳奇人物。雖然二三十年前,開源這個詞還沒有現在這樣風起云涌,但是總有一些具有黑客精神的牛人,內心里面充滿了自由,無論他處于哪個時代。Phil Katz這個人是個牛逼程序員,成名于DOS時代,我個人也沒有經歷過那個時代,我是從Windows98開始接觸電腦的,只是從書籍中得知,那個時代網速很慢,撥號使用的是只有幾十Kb(比特不是字節)的貓,56Kb實際上是這種貓的最高速度,在ADSL出現之后,這種技術被迅速淘汰。當時記錄文件的也是硬盤,但是在電腦之間拷貝文件的是軟盤,這個東西我大一還用過,最高容量記得是1.44MB,這還是200X年的軟盤,以前的軟盤容量具體多大就不知道了,Phil Katz上網的時候還不到1990年,WWW實際上就沒出現,瀏覽器當然是沒有的,當時上網干嘛呢?基本就是類似于網管敲各種命令,這樣實際上也可以聊天、上論壇不是嗎,傳個文件不壓縮的話肯定死慢死慢的,所以壓縮在那個時代很重要。當時有個商業公司提供了一種稱為ARC的壓縮軟件,可以讓你在那個時代聊天更快,當然是要付費的,Phil Katz就感覺到不爽,于是寫了一個PKARC,免費的,看名字知道是兼容ARC的,于是網友都用PKARC了,ARC那個公司自然就不爽,把哥們告上了法庭,說牽涉了知識產權等等,結果Phil Katz坐牢了。。。牛人就是牛人, 在牢里面冥思苦想,決定整一個超越ARC的牛逼算法出來,牢里面就是適合思考,用了兩周就整出來的,稱為PKZIP,不僅免費,而且這次還開源了,直接公布源代碼,因為算法都不一樣了,也就不涉及到知識產權了,于是ZIP流行開來,不過Phil Katz這個人沒有從里面賺到一分錢,還是窮困潦倒,因為喝酒過多等眾多原因,2000年的時候死在一個汽車旅館里。英雄逝去,精神永存,現在我們用UE打開ZIP文件,我們能看到開頭的兩個字節就是PK兩個字符的ASCII碼。

?

2、一個案例的入門思考

好了,Phil Katz在牢里面到底思考了什么?用什么樣的算法來壓縮數據呢?我們想一個簡單的例子:

生,容易。活,容易。生活,不容易。

上面這句話假如不壓縮,如果使用Unicode編碼,每個字會用2個字節表示。為什么是兩個字節呢?Unicode是一種國際標準,把常見各國的字符,比如英文字符、日文字符、韓文字符、中文字符、拉丁字符等等全部制定了一個標準,顯然,用2個字節可以最多表示2^16=65536個字符,那么65536就夠了嗎?生僻字其實是很多的,比如光康熙字典里面收錄的漢字就好幾萬,所以實際上是不夠的,那么是不是擴到4個字節?也可以,這樣空間倒是變大了,可以收錄更多字符,但一方面擴到4個字節就一定保證夠嗎?另一方面,4個字節是不是太浪費空間了,就為了那些一般情況都不會出現的生僻字?所以,一般情況下,使用2個字節表示,當出現生僻字的時候,再使用4個字節表示。這實際上就體現了信息論中數據壓縮基本思想,出現頻繁的那些字符,表示得短一些;出現稀少的,可以表示得長些(反正一般情況下也不會出現),這樣整體長度就會減小。除了Unicode,ASCII編碼是針對英文字符的編碼方案,用1個字節即可,除了這兩種編碼方案,還有很多地區性的編碼方案,比如GB2312可以對中文簡體字進行編碼,Big5可以對中文繁體字進行編碼。兩個文件如果都使用一種編碼方案,那是沒有問題的,不過考慮到國際化,還是盡量使用Unicode這種國際標準吧。不過這個跟ZIP沒啥關系,純屬背景介紹。

好了,回到我們前面說的例子,一共有17個字符(包括標點符號),如果用普通Unicode表示,一共是17*2=34字節。可不可以壓縮呢?所有人一眼都可以看出里面出現了很多重復的字符,比如里面出現了好多次容易(實際上是容易加句號三個字符)這個詞,第一次出現的時候用普通的Unicode,第二次出現的“容易?!眲t可以用(距離、長度)表示,距離的意思是接下來的字符離前面重復的字符隔了幾個,長度則表示有幾個重復字符,上面的例子的第二個“容易?!本捅硎緸?#xff08;5,3),就是距離為5個字符,長度是3,在解壓縮的時候,解到這個地方的時候,往前跳5個字符,把這個位置的連續3個字符拷貝過來就完成了解壓縮,這實際上不就是指針的概念?沒有錯,跟指針很類似,不過在數據壓縮領域,一般稱為字典編碼,為什么叫字典呢,當我們去查一個字的時候,我們總是先去目錄查找這個字在哪一頁,再翻到那一頁去看,指針不也是這樣,指針不就是內存的地址,要對一個內存進行操作,我們先拿到指針,然后去那塊內存去操作。所謂的指針、字典、索引、目錄等等術語,不同的背景可能稱呼不同,但我們要理解他們的本質。如果使用(5,3)這種表示方法,原來需要用6個字節表示,現在只需要記錄5和3即可。那么,5和3怎么記錄呢?一種方法自然還是可以用Unicode,那么就相當于節省了2個字節,但是有兩個問題,第一個問題是解壓縮的時候怎么知道是正常的5和3這兩個字符,還是這只是一個特殊標記呢?所以前面還得加一個標志來區分一下,到底接下來的Unicode碼是指普通字符,還是指距離和長度,如果是普通Unicode,則直接查Unicode碼表,如果是距離和長度,則往前面移動一段距離,拷貝即可。第二個問題,還是壓縮程度不行,這么一弄,感覺壓縮不了多少,如果重復字符比較長那倒是比較劃算,因為反正“距離+長度”就夠了,但比如這個例子,如果5和3前面加一個特殊字節,豈不是又是3個字節,那還不如不壓縮。咋辦呢?能不能對(5,3)這種整數進行再次壓縮?這里就利用了我們前面說的一個基本原則:出現的少的整數多編一些比特,出現的多的整數少編一些比特。那么,比如3、4、5、6、7、8、9這些距離誰出現得多?誰出現的少呢?誰知道?

壓縮之前當然不知道,不過掃描一遍不就知道了?比如,后面那個重復的字符串“容易。”按照前面的規則可以表示為(7,3),即離前面重復的字符串距離為7,長度為3。(7,3)指著前面跟自己一樣那個字符串。那么,為什么不指著第一個“容易?!币钢诙€“容易。”呢?如果指著第一個,那就不是(7,3)了,就是(12,3)了。當然,表示為(12,3)也可以解壓縮,但是有一個問題,就是12這個值比7大,大了又怎么了?我們在生活中會發現一些普遍規律,重復現象往往具有局部性。比如,你跟一個人說話,你說了一句話以后,往往很快會重復一遍,但是你不會隔了5個小時又重復這句話,這種特點在文件里面也存在著,到處都是這種例子,比如你在編程的時候,你定義了一個變量int nCount,這個nCount一般你很快就會用到,不會離得很遠。我們前面所說的距離代表了你隔了多久再說這句話,這個距離一般不大,既然如此,應該以離當前字符串距離最近的那個作為記錄的依據(也就是指向離自己最近那個重復字符串),這樣的話,所有的標記都是一些短距離,比如都是3、4、5、6、7而不會是3、5、78、965等等,如果大多數都是一些短距離,那么這些短距離就可以用短一些的比特表示,長一些的距離不太常見,則用一些長一些的比特表示。這樣, 總體的表示長度就會減少。好了,我們前面得到了(5,3)、(7、3)這種記錄重復的表示,距離有兩種:5、7;長度只有1種:3。咋編碼?越短越好。

既然表示的比特越短越好,3表示為0、5表示為10、7表示為11,行不行?這樣(5,3),(7,3)就只需要表示為100、110,這樣豈不是很短?貌似可以,貌似很高效。

但解壓縮遇到10這兩個比特的時候,怎么知道10表示5呢?這種表示方法是一個映射表,稱為碼表。我們設計的上面這個例子的碼表如下:

3–>0

5–>10

7–>11

這個碼表也得傳過去或者記錄在壓縮文件里才行啊,否則無法解壓縮,但會不會記錄了碼表以后整體空間又變大了,會不會起不到壓縮的作用?而且一個碼表怎么記錄?碼表記錄下來也是一堆數據,是不是也需要編碼?碼表是否可以繼續壓縮?那豈不是又需要新的碼表?壓縮會不會是一個永無止境的過程?作為一個入門級的同學,大概想到這兒就不容易想下去了。

?

3、ZIP中的LZ編碼思想

上面我們說的重復字符串用指針標記記錄下來,這種方法就是LZ這兩個人提出來的,理解起來比較簡單。后面分析(5,3)這種指針標記應該怎么編碼的時候,就涉及到一種非常廣泛的編碼方式,Huffman編碼,Huffman大致和香農是一個時代的人,這種編碼方式是他在MIT讀書的時候提出來的。接下來,我們來看看ZIP是怎么做的。

以上面的例子,一個很簡單的示意圖如下:

可以看出,ZIP中使用的LZ77算法和前面分析的類似,當然,如果仔細對比的話,ZIP中使用的算法和LZ提出來的LZ77算法其實還是有差異的,不過我建議不用仔細去扣里面的差異,思想基本是相同的,我們后面會簡要分析一下兩者的差異。LZ77算法一般稱為“滑動窗口壓縮”,我們前面說過,該算法的核心是在前面的歷史數據中尋找重復字符串,但如果要壓縮的文件有100MB,是不是從文件頭開始找?不是,這里就涉及前面提過的一個規律,重復現象是具有局部性的,它的基本假設是,如果一個字符串要重復,那么也是在附近重復,遠的地方就不用找了,因此設置了一個滑動窗口,ZIP中設置的滑動窗口是32KB,那么就是往前面32KB的數據中去找,這個32KB隨著編碼不斷進行而往前滑動。當然,理論上講,把滑動窗口設置得很大,那樣就有更大的概率找到重復的字符串,壓縮率不就更高了?初看起來如此,找的范圍越大,重復概率越大,不過仔細想想,可能會有問題,一方面,找的范圍越大,計算量會增大,不顧一切地增大滑動窗口,甚至不設置滑動窗口,那樣的軟件可能不可用,你想想,現在這種方式,我們在壓縮一個大文件的時候,速度都已經很慢了,如果增大滑動窗口,速度就更慢,從工程實現角度來說,設置滑動窗口是必須的;另一方面,找的范圍越大,距離越遠,出現的距離很多,也不利于對距離進行進一步壓縮吧,我們前面說過,距離和長度最好出現的值越少越好,那樣更好壓縮,如果出現的很多,如何記錄距離和長度可能也存在問題。不過,我相信滑動窗口設置得越大,最終的結果應該越好一些,不過應該不會起到特別大的作用,比如壓縮率提高了5%,但計算量增加了10倍,這顯然有些得不償失。

在第一個圖中,“容易?!笔且粋€重復字符串,距離distance=5,字符串長度length=3。當對這三個字符壓縮完畢后,接下來滑動窗口向前移動3個字符,要壓縮的是“我…”這個字符串,但這個串在滑動窗口內沒找到,所以無法使用distance+length的方式記錄。這種結果稱為literal。literal的中文含義是原義的意思,表示沒有使用distance+length的方式記錄的那些普通字符。literal是不是就用原始的編碼方式,比如Unicode方式表示?ZIP里不是這么做的,ZIP把literal認為也是一個數,盡管不能用distance+length表示,但不代表不可以繼續壓縮。另外,如果“我”出現在了滑動窗口內,是不是就可以用distance+length的方式表示?也不是,因為一個字出現重復,不值得用這種方式表示,兩個字呢?distance+length就是兩個整數,看起來也不一定值得,ZIP中確實認為2個字節如果在滑動窗口內找到重復,也不管,只有3個字節以上的重復字符串,才會用distance+length表示,重復字符串的長度越長越好,因為不管多長,都用distance+length表示就行了。

這樣的話,一段字符串最終就可以表示成literal、distance+length這兩種形式了。LZ系列算法的作用到此為止,下面,Phil Katz考慮使用Huffman對上面的這些LZ壓縮后的結果進行二次壓縮。個人認為接下來的過程才是ZIP的核心,所以我們要熟悉一下Huffman編碼。

?

4、ZIP中的Huffman編碼思想

上面LZ壓縮結果有三類(literal、distance、length),我們拿出distance一類來舉例。distance代表重復字符串離前一個一模一樣的字符串之間的距離,是一個大于0的整數。如何對一個整數進行編碼呢?一種方法是直接用固定長度表示,比如采用計算機里面描述一個4字節整數那樣去記錄,這也是可以的,主要問題當然是浪費存儲空間,在ZIP中,distance這個數表示的是重復字符串之間的距離,顯然,一般而言,越小的距離,出現的數量可能越多,而越大的距離,出現的數量可能越少,那么,按照我們前面所說的原則,小的值就用較少比特表示,大的值就用較多比特表示,在我們這個場景里,distance當然也不會無限大,比如不會超過滑動窗口的最大長度,假如對一個文件進行LZ壓縮后,得到的distance值為:

3、6、4、3、4、3、4、3、5

這個例子里,3出現了4次,4出現了3次,5出現了1次,6出現了1次。當然,不同的文件得到的結果不同,這里只是舉一個典型的例子,因為只有4種值,所以我們沒有必要對其它整數編碼。只需要對這4個整數進行編碼即可。

那么,怎么設計一個算法,符合3的編碼長度最短?6的編碼長度最長這種直觀上可行的原則(我們并沒有說這是理論上最優的方式)呢?

看起來似乎很難想出來。我們先來簡化一下,用固定長度表示。這里有4個整數,只要使用2個比特表示即可。于是這樣表示就很簡單:

00–>3; 01–>4; 10–>5;? 11–>6。

00、01這種編碼結果稱為碼字,碼字的平均長度是2。上面這個對應關系即為碼表,在壓縮時,需要將碼表放在最前面,后面的數字就用碼字表示,解碼時,先把碼表記錄在內存里,比如用一個哈希表記錄下來,解壓縮時,遇到00,就解碼為3等等。

因為出現了9個數,所以全部碼字總長度為18個比特。(我們暫時不考慮記錄碼表到底要占多少空間)

想要編碼結果更短,因為3出現的最多,所以考慮把3的碼字縮短點,比如3是不是可以用1個比特表示,這樣才算縮短吧,因為0和1只是二進制的一個標志,所以用0還是1沒有本質區別,那么,我們暫定把3用比特0表示。那么,4、5、6還能用0開頭的碼字表示呢?

這樣會存在問題,因為4、5、6的編碼結果如果以0開頭,那么,在解壓縮的時候,遇到比特0,就不知道是表示3還是表示4、5、6了,就無法解碼,當然,似乎理論上也不是不可以,比如可以往后解解看,比如假定0表示3的條件下往后解,如果無效則說明這個假設不對,但這種方式很容易出現兩個字符串的編碼結果是一樣的,這個誰來保證?所以,4、5、6都得以1開頭才行,那么,按照這個原則,4用1個比特也不行,因為5、6要么以0開頭,要么以1開頭,就無法編碼了,所以我們將4的碼字增加至2個比特,比如10,于是我們得到了部分碼表:

0–>3;10–>4。

按照這個道理,5、6既不能以0開頭,也不能以10開頭了,因為同樣存在無法解碼的問題,所以5應該以11開頭,就定為11行不行呢?也不行,因為6就不知道怎么編碼了,6也不能以0開頭,也不能以10、11開頭,那就無法表示了,所以,迫不得已,我們必須把5擴展一位,比如110,那么,6顯然就可以用111表示了,反正也沒有其他數了。于是我們得到了最終的碼表:

0–>3;10–>4;110–>5;111–>6。

看起來,編碼結果只能是這樣了,我們來算一下,碼字的總長度減少了沒有,原來的9個數是3、6、4、3、4、3、4、3、5,分別對應的碼字是:

0、111、10、0、10、0、10、0、110

算一下,總共16個比特,果然比前面那種方式變短了。我們在前面的設計過程中,是按照這些值出現次數由高到底的順序去找碼字的,比如先確定3,再確定4、5、6等等。按照一個碼字不能是另一個碼字的前綴這一規則,逐步獲得所有的碼字。這種編碼規則有一個專用術語,稱為前綴碼。Huffman編碼基本上就是這么做的,把出現頻率排個序,然后逐個去找,這個逐個去找的過程,就引入了二叉樹。不過Huffman的算法一般是從頻率由低到高排序,從樹的下面依次往上合并,不過本質上沒區別,理解思想即可。上面的結果可以用一顆二叉樹表示為下圖:

這棵樹也稱為碼樹,其實就是碼表的一種形式化描述,每個節點(除了葉子節點)都會分出兩個分支,左分支代表比特0,右邊分支代表1,從根節點到葉子節點的一個比特序列就是碼字。因為只有葉子節點可以是碼字,所以這樣也符合一個碼字不能是另一個碼字的前綴這一原則。說到這里,可以說一下另一個話題,就是一個映射表map在內存中怎么存儲,沒有相關經驗的可以跳過,map實現的是key–>value這樣的一個表,map的存儲一般有哈希表和樹形存儲兩類,樹形存儲就可以采用上面這棵樹,樹的中間節點并沒有什么含義,葉子節點的值表示value,從根節點到葉子節點上分支的值就是key,這樣比較適合存儲那些key由多個不等長字符組成的場合,比如key如果是字符串,那么把二叉樹的分支擴展很多,成為多叉樹,每個分支就是a,b,c,d這種字符,這棵樹也就是Trie樹,是一種很好使的數據結構。利用樹的遍歷算法,就實現了一個有序Map。

好了,我們理解了Huffman編碼的思想,我們來看看distance的實際情況。ZIP中滑動窗口大小固定為32KB,也就是說,distance的值范圍是1-32768。那么,通過上面的方式,統計頻率后,就得到32768個碼字,按照上面這種方式可以構建出來。于是我們會遇到一個最大的問題,那就是這棵樹太大了,怎么記錄呢?

好了,個人認為到了ZIP的核心了,那就是碼樹應該怎么縮小,以及碼樹怎么記錄的問題。

?

5、ZIP中Huffman碼樹的記錄方式

分析上面的例子,看看這個碼表:

0–>3;10–>4;110–>5;111–>6。

我們之前提過,0和1就是二進制的一個標志,互換一下其實根本不影響編碼長度,所以,下面的碼表其實是一樣的:

1–>3;00–>4;010–>5;011–>6。

1–>3;01–>4;000–>5;001–>6。

0–>3;11–>4;100–>5;101–>6。

。。。。。

這些都可以,而且編碼長度完全一樣,只是碼字不同而已。

對比一下第一個和第二個例子,對應的碼樹是這個樣子:

也就是說,我們把碼樹的任意節點的左右分支旋轉(0、1互換),也可以稱為樹的左右互換,其實不影響編碼長度,也就是說,這些碼表其實都是一樣好的,使用哪個都可以。

這個規律暗示了什么信息呢?暗示了碼表可以怎么記錄呢?Phil Katz當年在牢里也深深地思考了這一問題。

為了體會Phil Katz當時的心情,我們有必要盯著這兩棵樹思考幾分鐘:怎么把一顆樹用最少的比特記錄下來?

Phil Katz當時思考的邏輯我猜是這樣的,既然這些樹的壓縮程度都一樣,那干脆使用最特殊的那棵樹,反正壓縮程度都一樣,只要ZIP規定了這棵樹的特殊性,那么我記錄的信息就可以最少,這種特殊化的思想在后面還會看到。不同的樹當然有不同的特點,比如數據結構里面常見的平衡樹也是一類特殊的樹,他選的樹就是左邊那棵,這棵樹有一個特點,越靠左邊越淺,越往右邊越深,是這些樹中最不平衡的樹。ZIP里的壓縮算法稱為Deflate算法,這棵樹也稱為Deflate樹,對應的解壓縮算法稱為Inflate,Deflate的大致意思是把輪胎放氣了,意為壓縮;Inflate是給輪胎打氣的意思,意為解壓。那么,Deflate樹的特殊性又帶來什么了?

揭曉答案吧,Phil Katz認為換來換去只有碼字長度不變,如果規定了一類特殊的樹,那么就只需要記錄碼字長度即可。比如,一個有效的碼表是0–>3;10–>4;110–>5;111–>6。但只需要記錄這個對應關系即可:

3  4  5  6

1  2  3  3

也就是說,把1、2、3、3記錄下來,解壓一邊照著左邊那棵樹的形狀構造一顆樹,然后只需要1、2、3、3這個信息自然就知道是0、10、110、111。這就是Phil Katz想出來的ZIP最核心的一點:這棵碼樹用碼字長度序列記錄下來即可。

當然,只把1、2、3、3這個序列記錄下來還不行,比如不知道111對應5還是對應6?

所以,構造出樹來只是知道了有哪些碼字了,但是這些碼字到底對應哪些整數還是不知道。

Phil Katz于是又出現了一個想法:記錄1、2、3、3還是記錄1、3、2、3,或者3、3、2、1,其實都能構造出這棵樹來,那么,為什么不按照一個特殊的順序記錄呢?這個順序就是整數的大小順序,比如上面的3、4、5、6是整數大小順序排列的,那么,記錄的順序就是1、2、3、3。而不是2、3、3、1。

好了,根據1、2、3、3這個信息構造出了碼字,這些碼字對應的整數一個比一個大,假如我們知道編碼前的整數就是3、4、5、6這四個數,那就能對應起來了,不過到底是哪四個還是不知道啊?這個整數可以表示距離啊,距離不知道怎么去解碼LZ?

Phil Katz又想了,既然distance的范圍是1-32768,那么就按照這個順序記錄。上面的例子1和2沒有,那就記錄長度0。所以記錄下來的碼字長度序列為:

0、0、1、2、3、3、0、0、0、0、0、。。。。。。。。。。。。

這樣就知道構造出來的碼字對應哪個整數了吧,但因為distance可能的值很多(32768個),但實際出現的往往不多,中間會出現很多0(也就是根本就沒出現這個距離),不過這個問題倒是可以對連續的0做個特殊標記,這樣是不是就行了呢?還有什么問題?

我們還是要站在時代的高度來看待這個問題。我們明白,每個distance肯定對應唯一一個碼字,使用Huffman編碼可以得到所有碼字,但是因為distance可能非常多,雖然一般不會有32768這么多,但對一個大些的文件進行LZ編碼,distance上千還是很正常的,所以這棵樹很大,計算量、消耗的內存都容易超越了那個時代的硬件條件,那么怎么辦呢?這里再次體現了Phil Katz對Huffman編碼掌握的深度,他把distance劃分成多個區間,每個區間當做一個整數來看,這個整數稱為Distance Code。當一個distance落到某個區間,則相當于是出現了那個Code,多個distance對應于一個Distance Code,Distance雖然很多,但Distance Code可以劃分得很少,只要我們對Code進行Huffman編碼,得到Code的編碼后,Distance Code再根據一定規則擴展出來。那么,劃分多少個區間?怎么劃分區間呢?我們分析過,越小的距離,出現的越多;越大的距離,出現的越少,所以這種區間劃分不是等間隔的,而是越來越稀疏的,類似于下面的劃分:

1、2、3、4這四個特殊distance不劃分,或者說1個Distance就是1個區間;5,6作為一個區間;7、8作為一個區間等等,基本上,區間的大小都是1、2、4、8、16、32這么遞增的,越往后,涵蓋的距離越多。為什么這么分呢?首先自然是距離越小出現頻率越高,所以距離值小的時候,劃分密一些,這樣相當于一個放大鏡,可以對小的距離進行更精細地編碼,使得其編碼長度與其出現次數盡量匹配;對于距離較大那些,因為出現頻率低,所以可以適當放寬一些。另一個原因是,只要知道這個區間Code的碼字,那么對于這個區間里面的所有distance,后面追加相應的多個比特即可,比如,17-24這個區間的Huffman碼字是110,因為17-24這個區間有8個整數,于是按照下面的規則即可獲得其distance對應的碼字:

17–>110 000

18–>110 001

19–>110 010

20–>110 011

21–>110 100

22–>110 101

23–>110 110

24–>110 111

這樣計算復雜度和內存消耗是不是很小了,因為需要進行Huffman編碼的整數一下字變少了,這棵樹不會多大,計算起來時間和空間復雜度降低,擴展起來也比較簡單。當然,從理論上來說,這樣的編碼方式實際上將編碼過程分為了兩級,并不是理論上最優的,把所有distance當作一個大空間去編碼才可能得到最優結果,不過還是那句話,工程實現的限制,在壓縮軟件實現上,我們不能用壓縮率作為衡量一個算法優劣的唯一指標,其實耗費的時間和空間同樣是指標,所以需要看綜合指標。很多其他軟件也一樣,擴展性、時間空間復雜度、穩定性、移植性、維護的方便性等等是工程上很重要的東西。我沒有看過RAR是如何壓縮的,有可能是在類似的地方進行了改進,如果如此,那也是站在巨人的肩膀上,而且硬件條件不同,進行一些改進也并不奇怪。

具體來說,Phil Katz把distance劃分為30個區間,如下圖:

這個圖是我從David Salomon的《Data Compression The Complete Reference》這本書(第四版)中拷貝出來的,下面的有些圖也是,如果需要對數據壓縮進行全面的了解,這本書幾乎是最全的了,強烈推薦。

當然,你要問為什么是30個區間,我也沒分析過,也許是復雜度和壓縮率經過試驗之后的一種折中吧。

其中,左邊的Code表示區間的編號,是0-29,共30個區間,這只是個編號,沒有特別的含義,但Huffman就是對0-29這30個Code進行編碼的,得到區間的碼字;

bits表示distance的碼字需要在Code的碼字基礎上擴展幾位,比如0就表示不擴展,最大的13表示要擴展13位,因此,最大的區間包含的distance數量為8192個。

Distance一列則表示這個區間涵蓋的distance范圍。

理解了碼樹如何有效記錄,以及如何縮小碼樹的過程,我覺得就理解了ZIP的精髓。

?

6、ZIP中literal和length的壓縮方式

說完了distance,LZ編碼結果還有兩類:literal和length。這兩類也利用了類似于distance的方式進行壓縮。

前面分析過,literal表示未匹配的字符,我們前面之所以拿漢字來舉例,完全是為了便于理解,ZIP之所以是通用壓縮,它實際上是針對字節作為基本字符來編碼的,所以一個literal至多有256種可能。

length表示重復字符串長度,length=1當然不會出現,因為一個字符不值得用distance+length去記錄,重復字符串當然越長越好,Phil Katz(下面還是簡稱PK了,拷貝太麻煩)認為,length=2也不值得用這種方式記錄,還是太短了,所以PK把length最小值認為是3,必須3個以上字符的字符串出現重復才用distance+length記錄。那么,最大的length是多少呢?理論上當然可以很長很長,比如一個文件就是連續的0,這個重復字符串長度其實接近于這個文件的實際長度。但是PK把length的范圍做了限制,限定length的個數跟literal一樣,也只有256個,因為PK認為,一個重復字符串達到了256個已經很長了,概率非常小;另外,其實哪怕超過了256,我還是認為是一段256再加上另外一段,增加一個distance+length就行了嘛,并不影響結果。而且這樣做,我想同樣也考慮了硬件條件吧。

初看有點奇怪的在于,將literal和length二者合二為一,什么意思呢?就是對這兩種整數(literal本質上是一個字節)共用一個Huffman碼表,一會兒解釋為什么。PK對Huffman的理解我覺得達到了爐火純青的地步,前面已經看到,后面還會看到。他認為Huffman編碼的輸入反正說白了就是一個集合的元素就行,無論這個元素是啥,所以多個集合看做一個集合當作Huffman編碼的輸入沒啥問題。literal用整數0-255表示,256是一個結束標志,解碼以后結果是256表示解碼結束;從257開始表示length,所以257這個數表示length=3,258這個數表示length=4等等,但PK也不是一直這么一一對應,和distance一樣,也是把length(總共256個值)劃分為29個區間,其結果如下圖:

其中的含義和distance類似,不再贅述,所以literal/length這個Huffman編碼的輸入元素一共285個,其中256表示解碼結束標志。為什么要把二者合二為一呢?因為當解碼器接收到一個比特流的時候,首先可以按照literal/length這個碼表來解碼,如果解出來是0-255,就表示未匹配字符,如果是256,那自然就結束,如果是257-285之間,則表示length,把后面擴展比特加上形成length后,后面的比特流肯定就表示distance,因此,實際上通過一個Huffman碼表,對各類情況進行了統一,而不是通過加一個什么標志來區分到底是literal還是重復字符串。

好了,理解了上面的過程,就理解了ZIP壓縮的第二步,第一步是LZ編碼,第二步是對LZ編碼后結果(literal、distance、length)進行的再編碼,因為literal/length是一個碼表,我稱其為Huffman碼表1,distance那個碼表稱為Huffman碼表2。前面我們已經分析了,Huffman碼樹用一個碼字長度序列表示,稱為CL(Code Length),記錄兩個碼表的碼字長度序列分別記為CL1、CL2。碼樹記錄下來,對literal/length的編碼比特流稱為LIT比特流;對distance的編碼比特流稱為DIST比特流。

按照上面的方法,LZ的編碼結果就變成四塊:CL1、CL2、LIT比特流、DIST比特流。CL1、CL2是碼字長度的序列,這個序列說白了就是一堆正整數,因此,PK繼續深挖,認為這個序列還應該繼續壓縮,也就是說,對碼表進行壓縮。

?

7、ZIP中對CL進行再次壓縮的方法

這里仍然沿用Huffman的想法,因為CL也是一堆整數,那么當然可以再次應用Huffman編碼。不過在這之前,PK對CL序列進行了一點處理。這個處理也是很精巧的。

CL序列表示一系列整數對應的碼字長度,對于literal/length來說,總共有0-285這么多符號,所以這個序列長度為286,每個符號都有一個碼字長度,當然,這里面可能會出現大段連續的0,因為某些字符或長度不存在,尤其是對英文文本編碼的時候,非ASCII字符就根本不會出現,length較大的值出現概率也很小,所以出現大段的0是很正常的;對于distance也類似,也可能出現大段的0。PK于是先進行了一下游程編碼。在說什么是游程編碼之前,我們談談PK對CL序列的認識。

literal/length的編碼符號總共286個(回憶:256個Literal+1個結束標志+29個length區間),distance的編碼符號總共30個(回憶:30個區間),所以這顆碼樹不會特別深,Huffman編碼后的碼字長度不會特別長,PK認為最長不會超過15,也就是樹的深度不會超過15,這個是否是理論證明我還沒有分析,有興趣的同學可以分析一下。因此,CL1和CL2這兩個序列的任意整數值的范圍是0-15。0表示某個整數沒有出現(比如literal=0x12, length Code=8, distance Code=15等等)。

什么叫游程呢?就是一段完全相同的數的序列。什么叫游程編碼呢?說起來原理更簡單,就是對一段連續相同的數,記錄這個數一次,緊接著記錄出現了多少個即可。David的書中舉了這個例子,比如CL序列如下:

4, 4, 4, 4, 4, 3, 3, 3, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 0, 0, 0, 0, 0, 0, 2, 2, 2, 2
那么,游程編碼的結果為:

4, 16, 01(二進制), 3, 3, 3, 6, 16, 11(二進制), 16, 00(二進制), 17,011(二進制), 2, 16, 00(二進制)
這是什么意思呢?因為CL的范圍是0-15,PK認為重復出現2次太短就不用游程編碼了,所以游程長度從3開始。用16這個特殊的數表示重復出現3、4、5、6個這樣一個游程,分別后面跟著00、01、10、11表示(實際存儲的時候需要低比特優先存儲,需要把比特倒序來存,博文的一些例子有時候會忽略這點,實際寫程序的時候一定要注意,否則會得到錯誤結果)。于是4,4,4,4,4,這段游程記錄為4,16,01,也就是說,4這個數,后面還會連續出現了4次。6,16,11,16,00表示6后面還連續跟著6個6,再跟著3個6;因為連續的0出現的可能很多,所以用17、18這兩個特殊的數專門表示0游程,17后面跟著3個比特分別記錄長度為3-10(總共8種可能)的游程;18后面跟著7個比特表示11-138(總共128種可能)的游程。17,011(二進制)表示連續出現6個0;18,0111110(二進制)表示連續出現62個0??傊涀?#xff0c;0-15是CL可能出現的值,16表示除了0以外的其它游程;17、18表示0游程。因為二進制實際上也是個整數,所以上面的序列用整數表示為:

4, 16, 1, 3, 3, 3, 6, 16, 3, 16, 0, 17, 3, 2, 16, 0

我們又看到了一串整數,這串整數的值的范圍是0-18。這個序列稱為SQ(Sequence的意思)。因為有兩個CL1、CL2,所以對應的有兩個SQ1、SQ2。

針對SQ1、SQ2,PK用了第三個Huffman碼表來對這兩個序列進行編碼。通過統計各個整數(0-18范圍內)的出現次數,按照相同的思路,對SQ1和SQ2進行了Huffman編碼,得到的碼流記為SQ1 bits和SQ2 bits。同時,這里又需要記錄第三個碼表,稱為Huffman碼表3。同理,這個碼表也用相同的方法記錄,也等效為一個碼長序列,稱為CCL,因為至多有0-18個,PK認為樹的深度至多為7,于是CCL的范圍是0-7。

當得到了CCL序列后,PK決定不再折騰,對這個序列用普通的3比特定長編碼記錄下來即可,即000代表0,111代表7。但實際上還有一點小折騰,就是最后這個序列如果全部記錄,那就需要19*3=57個比特,PK認為CL序列里面CL范圍為0-15,特殊的幾個值是16、17、18,如果把CCL序列位置置換一下,把16、17、18這些放前面,那么這個CCL序列就很可能最后面跟著一串0(因為CL=14,15這些很可能沒有),所以最后還引入了一個置換,其示意圖如下,分別表示置換前的CCL序列和置換后的CCL??梢钥闯?#xff0c;16、17、18對應的CCL被放到了前面,這樣如果尾部出現了一些0,就只需要記錄CCL長度即可,后面的0不記錄。可以繼續節省一些比特,不過這個例子尾部置換后只有1個0:

不過粗看起來,這個置換效果并不好,我一開始接觸這個置換的時候,我覺得應該按照16、17、18、0、1、2、3、。。。這樣的順序來存儲,如果按照我理解的,那么置換后的結果如下:

2、4、0、4、5、5、1、5、0、6、0、0、0、0、0、0、0、0、0

這樣后面的一大串0直接截斷,比PK的方法更短。但PK卻按照上面的順序。我總是認為,我覺得牛人可能出錯了的時候,往往是我自己錯了,所以我又仔細想了一下,上面的順序特點比較明顯,直觀上看,PK認為CL為0和中間的值出現得比較多(放在了前面),但CL比較小的和比較大的出現得比較少(1、15、2、14這些放在了后面,你看,后面交叉著放),在文件比較小的時候,這種方法效果不算好,上面就是一個典型的例子,但文件比較大了以后,CL1、CL2碼樹比較大,碼字長度普遍比較長,大部分很可能接近于中間值,那么這個時候PK的方法可能就體現出優勢了。不得不說,對一個算法或者數據結構的優化程度,簡直完全取決于程序員對那個東西細節的理解的深度。當我仔細研究了ZIP壓縮算法的過程之后,我對PK這種深夜埋頭冥思苦想的大牛佩服得五體投地。

到此為止,ZIP壓縮算法的結果已經完畢。這個算法命名為Deflate算法。總結一下其編碼流程為:

?

8、Deflate壓縮數據格式

ZIP的格式實際上就是Deflate壓縮碼流外面套了一層文件相關的信息,這里先介紹Deflate壓縮碼流格式。其格式為:

Header:3個比特,第一個比特如果是1,表示此部分為最后一個壓縮數據塊;否則表示這是.ZIP文件的某個中間壓縮數據塊,但后面還有其他數據塊。這是ZIP中使用分塊壓縮的標志之一;第2、3比特表示3個選擇:壓縮數據中沒有使用Huffman、使用靜態Huffman、使用動態Huffman,這是對LZ77編碼后的literal/length/distance進行進一步編碼的標志。我們前面分析的都是動態Huffman,其實Deflate也支持靜態Huffman編碼,靜態Huffman編碼原理更為簡單,無需記錄碼表(因為PK自己定義了一個固定的碼表),但壓縮率不高,所以大多數情況下都是動態Huffman。

HLIT:5比特,記錄literal/length碼樹中碼長序列(CL1)個數的一個變量。后面CL1個數等于HLIT+257(因為至少有0-255總共256個literal,還有一個256表示解碼結束,但length的個數不定)。

HDIST:5比特,記錄distance碼樹中碼長序列(CL2)個數的一個變量。后面CL2個數等于HDIST+1。哪怕沒有1個重復字符串,distance都為0也是一個CL。

HCLEN:4比特,記錄Huffman碼表3中碼長序列(CCL)個數的一個變量。后面CCL個數等于HCLEN+4。PK認為CCL個數不會低于4個,即使對于整個文件只有1個字符的情況。

接下來是3比特編碼的CCL,一共HCLEN+4個,用以構造Huffman碼表3;

接下來是對CL1(碼長)序列經過游程編碼(SQ1:縮短的整數序列)后,并對SQ1繼續用Huffman編碼后的比特流。包含HLIT+257個CL1,其解碼碼表為Huffman碼表3,用以構造Huffman碼表1;

接下來是對CL2(碼長)序列經過游程編碼(SQ2:縮短的整數序列)后,并對SQ2繼續用Huffman編碼后的比特流。包含HDIST+1個CL2,其解碼碼表為Huffman碼表3,用于構造Huffman碼表2;

總之,上面的數據都是為了構造LZ解碼需要的2個Huffman碼表。

接下來才是經過Huffman編碼的壓縮數據,解碼碼表為Huffman碼表1和碼表2。
最后是數據塊結束標志,即literal/length這個碼表輸入符號位256的編碼比特。
對倒數第1、2內容塊進行解碼時,首先利用Huffman碼表1進行解碼,如果解碼所得整數位于0-255之間,表示literal未匹配字符,接下來仍然利用Huffman碼表1解碼;如果位于257-285之間,表示length匹配長度,之后需要利用Huffman碼表2進行解碼得到distance偏移距離;如果等于256,表示數據塊解碼結束。

?

9、ZIP文件格式解析

?上面各節對ZIP的原理進行了分析,這一節我們來看一個實際的例子,為了更好地描述,我們盡量把這個例子舉得簡單一些。下面是我隨便從一本書拷貝出來的一段較短的待壓縮的英文文本數據:

As mentioned above,there are many kinds of wireless systems other than cellular.

這段英文文本長度為80字節。經過ZIP壓縮后,其內容如下:

可以看到,第1、2字節就是PK??粗趺幢仍倪€長,這怎么叫壓縮?實際上,這里面大部分內容是ZIP的文件標記開銷,真正壓縮的內容(也就是我們前面提到的Deflate數據,劃線部分都是ZIP文件開銷)其實肯定要比原文短(否則ZIP不會啟用壓縮),我們這個例子是個短文本,但對于更長的文本而言,那ZIP文件總體長度肯定是要短于原始文本的。上面的這個ZIP文件,可以看到好幾個以PK開頭的區域,也就是不同顏色的劃線區域,這些其實都是ZIP文件本身的開銷。

所以,我們首先來看一看ZIP的格式,其格式定義為:

[local file header 1]
[file data 1]
[data descriptor 1]
……….
[local file header n]
[file data n]
[data descriptor n]
[archive decryption header]?
[archive extra data record]?
[central directory]
[zip64 end of central directory record]
[zip64 end of central directory locator]?
[end of central directory record]
local file header+file data+data descriptor這是一段ZIP壓縮數據,在一個ZIP文件里,至少有一段,至多那就不好說了,假如你要壓縮的文件一共有10個,那這個地方至少會有10段,ZIP對每個文件進行了獨立壓縮,RAR在此進行了改進,將多個文件聯合起來進行壓縮,提高了壓縮率。local file header的格式如下:

可見,起始的4個字節就是0x50(P)、0x4B(K)、0x03、0x04,因為是低字節優先,所以Signature=0x03044B50.接下來的內容按照上面的格式解析,十分簡單,這個區域在上面ZIP數據的那個圖里面是紅色劃線區域,之后則是壓縮后的Deflate數據。在文件的尾部,還有ZIP尾部數據,上面這個例子包含了central directory和end of central directory record,一般這兩部分也是必須的。central directory以0x50、0x4B、0x01、0x02開頭;end of central directory record以0x50、0x4B、0x05、0x06開頭,其含義比較簡單,分別對應于上面ZIP數據那個圖的藍色和綠色部分,下面是兩者的格式:

end of central directory record格式:

這幾張圖是我從網上找的,寫得比較清晰。對于其中的含義,解釋起來也比較簡單,我分析的結果如下:注意ZIP采用的低字節優先,在一個字節里面低位優先,需要反過來看。

Local File Header: (38B,304b)
00001010110100101100000000100000 (signature)
0000000000010100 (version:20)
0000000000000000 (generalBitFlag)
0000000000001000 (compressionMethod:8)
0100110110001110 (lastModTime:19854)
0100010100100101 (lastModDate:17701)
01010100101011010100001100111100 (CRC32)
00000000000000000000000001001000 (compressedSize:72)
00000000000000000000000001010000 (uncompressedSize:80)
0000000000001000 (filenameLength:8)
0000000000000000 (extraFieldLength:0)
0010101010100110110011100010111001110100001011100001111000101110 (fileName:Test.txt)
?(extraField)

Central File Header: (54B,432b)
00001010110100101000000001000000 (signature)
0000000000010100 (versionMadeBy:20)
0000000000010100 (versionNeeded:20)
0000000000000000 (generalBitFlag)
0000000000001000 (compressionMethod:8)
0100110110001110 (lastModTime:19854)
0100010100100101 (lastModDate:17701)
01010100101011010100001100111100 (CRC32)
00000000000000000000000001001000 (compressedSize:72)
00000000000000000000000001010000 (uncompressedSize:80)
0000000000001000 (filenameLength:8)
0000000000000000 (extraFieldLength:0)
0000000000000000 (fileCommenLength:0)
0000000000000000 (diskNumberStart)
0000000000000001 (internalFileAttr)
10000001100000000000000000100000 (externalFileAttr)
00000000000000000000000000000000 (relativeOffsetLocalHeader)
0010101010100110110011100010111001110100001011100001111000101110 (fileName:Test.txt)
?(extraField)
?(fileComment)

end of Central Directory Record: (22B,176b)
00001010110100101010000001100000 (signature)
0000000000000000 (numberOfThisDisk:0)
0000000000000000 (numberDiskCentralDirectory:0)
0000000000000001 (EntriesCentralDirectDisk:1)
0000000000000001 (EntriesCentralDirect:1)
00000000000000000000000000110110 (sizeCentralDirectory:54)
00000000000000000000000001101110 (offsetStartCentralDirectory:110)
0000000000000000 (fileCommentLength:0)
?(fileComment)

Local File Header Length:304
Central File Header Length:432
End Central Directory Record Length:176

可見,開銷總的長度為38+54+22=114字節,整個文件長度為186字節,因此Deflate壓縮數據長度為72字節(576比特)。盡管這里看起來只是從80字節壓縮到72字節,那是因為這是一段短文本,重復字符串出現較少,但如果文本較長,那壓縮率就會增加,這里只是舉個例子。

下面對其中的關鍵部分,也就是Deflate壓縮數據進行解析。

?

10,Deflate解碼過程實例分析

我們按照ZIP格式把Deflate壓縮數據(72字節)提取出來,如下(每行8字節):

1010100001010011100010111011000000000001000001000011000010100010
1000101110101010011110110000000001100011101110000011100010100101
0101001111001100000010001101001010010010000101101010101100001101
1011110100011111100011101111111001110010011101110110011100010101
0010110100010100101100110001100100000100110111101101111000011101
0010001001100110111001000010011001101010101000110110000001110101
0100011010010011100010110111001000111101101001011100101010010111
0111000011111000011110000011010111001011011111111100100010001001
1010001100001110000010101010111101101010100101111101011111100000

Deflate格式除了上面的介紹,也可以參考RFC1951,解析如下:

Header:101, 第一個比特是1,表示此部分為最后一個壓縮數據塊;后面的兩個比特01表示采用動態哈夫曼、靜態哈夫曼、或者沒有編碼的標志,01表示采用動態Huffman;在RFC1951里面是這么說明的:

00 – no compression

01 – compressed with fixed Huffman codes

10 – compressed with dynamic Huffman codes

11 – reserved (error)

注意,這里需要按照低比特在先的方式去看,否則會誤以為是靜態Huffman。

接下來:
HLIT:01000,記錄literal/length碼樹中碼長序列個數的一個變量,表示HLIT=2(低位在前),說明后面存在HLIT + 257=259個CL1,CL1即0-258被編碼后的長度,其中0-255表示Literal,256表示無效符號,257、258分別表示Length=3、4(length從3開始)。因此,這里實際上只出現了兩種重復字符串的長度,即3和4?;仡欉@個圖可以更清楚:

繼續:
HDIST:01010,記錄distance碼樹中碼長序列個數的一個變量,表示HDIST=10,說明后面存在HDIST+1=11個CL2,CL2即Distance Code=0-10被編碼的長度。

繼續:

HCLEN:0111,記錄Huffman碼樹3中碼長序列個數的一個變量,表示HCLEN=14(1110二進制),即說明緊接著跟著HCLEN+4=18個CCL,前面已經分析過,CCL記錄了一個Huffman碼表,這個碼表可以用一個碼長序列表示,根據這個碼長序列可以得到碼表。于是接下來我們把后面的18*3=54個比特拷貝出來,上面的碼流目前解析為下面的結果:

101(Header) 01000(HLIT) 01010(HDIST) 0111(HCLEN)?
000 101 110 110 000 000 000 010 000 010 000 110 000 101 000 101 000 101 (CCL)
110101010011110110000000001100011101110000011100010100101
0101001111001100000010001101001010010010000101101010101100001101
1011110100011111100011101111111001110010011101110110011100010101
0010110100010100101100110001100100000100110111101101111000011101
0010001001100110111001000010011001101010101000110110000001110101
0100011010010011100010110111001000111101101001011100101010010111
0111000011111000011110000011010111001011011111111100100010001001
1010001100001110000010101010111101101010100101111101011111100000

標準的CCL長度為19(回憶一下:CCL范圍為0-18,按照整數大小排序記錄各自的碼字長度),因此最后一個補0。得到序列:

000 101 110 110 000 000 000 010 000 010 000 110 000 101 000 101 000 101 000

其長度分別為(低位在前):
0、5、3、3、0、0、0、2、0、2、0、3、0、5、0、5、0、5、0
前面已經分析過,這個CCL序列實際上是經過一次置換操作得到的,需要進行相反的置換,置換后為:

3、5、5、5、3、2、2、0、0、0、0、0、0、0、0、0、0、5、3
這個就是對應于0-18的碼字長度序列。
根據Deflate樹的構造方式,得到下面的碼表(Huffman碼表3):

00????? <–>?? 5
01????? <–>?? 6
100???? <–>? 0
101???? <–>? 4
110???? <–>? 18
11100?? <–>1
11101?? <–>2
11110?? <–>3
11111?? <–>17

接下來就是CL1序列,按照前面的指示,一共有259個,分別對應于literal/length:0-258對應的碼字長度序列,我們隊跟著CCL后面的比特按照上面獲得的碼表進行逐步解碼,在解碼之前,實際上并不知道CL1的比特流長度有多少,需要根據259這個數字來判定,解完了259個整數,表明解析CL1完畢:

101(Header) 01000(HLIT) 01010(HDIST) 0111(HCLEN)?
000 101 110 110 000 000 000 010 000 010 000 110 000 101 000 101 000 101 (CCL)

110(18)1010100(7比特,記錄連續的11-138個0,此處一共0010101b=21,即記錄21+11=32個0)

11110(3)110(18)0000000(7比特,記錄連續的11-138個0,此處為全0,即記錄0+11=11個0)

01(6)100(0)01(6)110(18)1110000(7比特,記錄連續的11-138個0,此處為111b=7,即記錄7+11=18個0)

01(6)110(18)0010100(7比特,記錄連續的11-138個0,此處為10100b=20,即記錄20+11=31個0)

101(4)01(6)01(6)00(5)11110(3)01(6)100(0)00(5)00(5)100(0)01(6)101(4)

00(5)101(4)00(5)100(0)100(0)00(5)101(4)101(4)01(6)01(6)01(6)100(0)

00(5)110(18)1101111(7比特,記錄連續的11-138個0,此處為1111011b=123,即記錄123+11=134個0)

統計一下,上面已經解了32+11+18+31+134+30=256個數了,因為總共259個,還差三個:

01(6)00(5)01(6)

好了,CL1比特流解析完畢了,得到的CL1碼長序列為:

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0?3?0 0 0 0 0 0 0?
0 0 0 0 6 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0?
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4 6 6 5 3 6 0 5 5 0 6 4 5 4 5 0 0 5 4 4 6 6 6?
0 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0?
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0?
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0?
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0?6 5 6

總共259個,每行40個。根據這個序列,同樣按照Deflate樹構造方法,得到literal/length碼表(Huffman碼表1)為:

000???? –> (System.Char)(看前面的CL1序列,空格對應的ASCII為0x20=32,碼字長度3,即上面序列中第一個3)
001???? –>e(System.Char)
0100??? –>a(System.Char)
0101??? –>l(System.Char)
0110??? –>n(System.Char)
0111??? –>s(System.Char)
1000??? –>t(System.Char)
10010?? –>d(System.Char)
10011?? –>h(System.Char)
10100?? –>i(System.Char)
10101?? –>m(System.Char)
10110?? –>o(System.Char)
10111?? –>r(System.Char)
11000?? –>y(System.Char)
11001?? –>3(System.Int32)(看前面的CL1序列,對應257,碼字長度5)
110100? –>,(System.Char)
110101? –>.(System.Char)
110110? –>A(System.Char)
110111? –>b(System.Char)
111000? –>c(System.Char)
111001? –>f(System.Char)
111010? –>k(System.Char)
111011? –>u(System.Char)
111100? –>v(System.Char)
111101? –>w(System.Char)
111110? –>-1(System.Int32)(看前面的CL1序列,對應256,碼字長度6)
111111? –>4(System.Int32)(看前面的CL1序列,對應258,碼字長度6)

可以看出,碼表里存在兩個重復字符串長度3和4,當解碼結果為-1(上面進行了處理,即256),或者說遇到111110的時候,表示Deflate碼流結束。

按照同樣的道理,對CL2序列進行解析,前面已經知道HDIST=10,即有11個CL2整數序列:

11111(17)000(3比特,記錄連續的3-10個0,此處為0,即記錄3個0)

11101(2)11111(17)100(3比特,記錄連續的3-10個0,此處為001b=1,即記錄4個0)

11100(1)100(0)11101(2)

已經結束,總共11個。

于是CL2序列為:

0 0 0 2 0 0 0 0 1 0 2

分別記錄的是distance碼為0-10的碼字長度,根據下面的對應關系,需要進行擴展:

比如,第1個碼長2記錄的是Code=3的長度,即Distance=4對應的碼字為:

10????? –>4(System.Int32)

第1個碼長1記錄的是Code=8的長度(碼字為0,擴展三位000-111),即Distance=17-24對應的碼字為(注意,低比特優先):

0 000??? –>17(System.Int32)
0 100??? –>18(System.Int32)
0 010??? –>19(System.Int32)
0 110??? –>20(System.Int32)
0 001??? –>21(System.Int32)
0 101??? –>22(System.Int32)
0 011??? –>23(System.Int32)
0 111??? –>24(System.Int32)

注意,擴展的時候還是低比特優先。

最后1個碼長2記錄的是Code=10的長度(其實是碼字:11,擴展四位0000-1111),即Distance=33-48對應的碼字為:

11 0000? –>33(System.Int32)
11 1000? –>34(System.Int32)
11 0100? –>35(System.Int32)
11 1100? –>36(System.Int32)
11 0010? –>37(System.Int32)
11 1010? –>38(System.Int32)
11 0110? –>39(System.Int32)
11 1110? –>40(System.Int32)
11 0001? –>41(System.Int32)
11 1001? –>42(System.Int32)
11 0101? –>43(System.Int32)
11 1101? –>44(System.Int32)
11 0011? –>45(System.Int32)
11 1011? –>46(System.Int32)
11 0111? –>47(System.Int32)
11 1111? –>48(System.Int32)

至此為止,Huffman碼表1、Huffman碼表2已經還原出來,接下來是對LZ壓縮所得到的literal、distance、length進行解碼,目前剩余的比特流如下,先按照Huffman碼表1解碼,如果解碼結果是長度(>256),則接下來按照Huffman碼表2解碼,逐步解碼即可:

[As ]:110110(A)0111(s)000(空格)

[mentioned ]:10101(m)001(e)0110(n)1000(t)10100(i)10110(o)0110(n)001(e)10010(d)000(空格)

[above,]:0100(a)110111(b)10110(o)111100(v)001(e)110100(,)

[there ]:1000(t)10011(h)001(e)10111(r)001(e)000(空格)

[are ]:0100(a)11001(長度3,表示下一個需要用Huffman解碼)10(Distance=4,即重復字符串為re空格)

[many ]:10101(m)0100(a)0110(n)11000(y)000(空格)

[kinds ]:111010(k)10100(i)0110(n)10010(d)0111(s)000(空格)

[of ]:10110(o)111001(f)000(空格)

[wireless ]:111101(w)10100(i)10111(r)001(e)0101(l)001(e)0111(s)0111(s)000(空格)

[systems o]:0111(s)11000(y)0111(s)1000(t)001(e)10101(m)11001(長度指示=3,接下來根據distance解碼)0110(Distance=20,即重復字符串為s o)

[ther ]:111111(長度指示=4,接下來根據distance解碼)111001(Distance=42,即重復字符串為ther)000(空格)

[than ]:1000(t)10011(h)0100(a)0110(n)000(空格)

[cellular.]:111000(c)001(e)0101(l)0101(l)111011(u)0101(l)0100(a)10111(r)110101(.)

[256,結束標志]111110(結束標志)0000(字節補齊的0)

于是解壓縮結果為:

As mentioned above,there are many kinds of wireless systems other than cellular.

再來回顧我們的解碼過程:

譯碼過程:
1、根據HCLEN得到截尾信息,并參照固定置換表,根據CCL比特流得到CCL整數序列;
2、根據CCL整數序列構造出等價于CCL的二級Huffman碼表3;
3、根據二級Huffman碼表3對CL1、CL2比特流進行解碼,得到SQ1整數序列,SQ2整數序列;
4、根據SQ1整數序列,SQ2整數序列,利用游程編碼規則得到等價的CL1整數序列、CL2整數序列;
5、根據CL1整數序列、CL2整數序列分別構造兩個一級Huffman碼表:literal/length碼表、distance碼表;
6、根據兩個一級Huffman碼表對后面的LZ壓縮數據進行解碼得到literal/length/distance流;
7、根據literal/length/distance流按照LZ規則進行解碼。

Deflate碼流長度總共為72字節=576比特,其中:

3比特Header;

5比特HLIT;

5比特HDIST;

4比特HCLEN;

54比特CCL序列碼流;

133比特CL1序列碼流;

34比特CL2序列碼流;

338比特LZ壓縮后的literal/length/distance碼流。

?

11、ZIP的其它說明

上面各個環節已經詳細分析了ZIP壓縮的過程以及解碼流程,通過對一個實例的解壓縮過程分析,可以徹底地掌握ZIP壓縮和解壓縮的原理和過程。還有一些情況需要說明:

(1)上面的算法復雜度主要在于壓縮一端,因為需要統計literal/length/distance,創建動態Huffman碼表,相反解壓只需要還原碼表后,逐比特解析即可,這也是壓縮軟件的一個典型特點,解壓速度遠快于壓縮速度。

(2)上面我們分析了動態Huffman,對于LZ壓縮后的literal/length/distance,也可以采用靜態Huffman編碼,這主要取決于ZIP在壓縮中看哪種方式更節省空間,靜態Huffman編碼不需要記錄碼表,因為這個碼表是固定的,在RFC1951里面也有說明。對于literal/length碼表來說,需要對0-285進行編碼,其碼表為:

對于Distance來說,需要對Code=0-29的數進行編碼,則直接采用5比特表示。Distance和動態Huffman一樣,在此基礎上進行擴展。

(3)ZIP中使用的LZ77算法是一種改進的LZ77。主要區別有兩點:

1)標準LZ77在找到重復字符串時輸出三元組(length, distance, 下一個未匹配的字符)(有興趣可以關注LZ77那篇論文);Deflate在找到重復字符串時僅輸出雙元組(length, distance)。
2)標準LZ77使用”貪婪“的方式解析,尋找的都是最長匹配字符串。Deflate中不完全如此。David Salomon的書里給了一個例子:

對于上面這個例子,標準LZ77在滑動窗口中查找最長匹配字符串,找到的是”the”與前面的there的前三個字符匹配,這種貪婪解析方式邏輯簡單,但編碼效率不一定最高。Deflate則不急于輸出,跳過t繼續往后查看,發現”th ne”這5個字符存在重復字符串,因此,Deflate算法會選擇將t作為未匹配字符輸出,而對后面的匹配字符串用(length, distance)編碼輸出。顯然,這樣就提高了壓縮效率,因為標準的LZ77找到的重復字符串長度為3,而Deflate找到的是5。換句話說,Deflate算法并不是簡單的尋找最長匹配后輸出,而是會權衡幾種可行的編碼方式,用其中最高效的方式輸出。

?

12、總結

本篇博文對ZIP中使用的壓縮算法進行了詳細分析,從一個簡單地例子出發,一步步地分析了PK設計Deflate算法的思路。最后,通過一個實際例子,分析了其解壓縮流程??偟膩砜?#xff0c;ZIP的核心在于如何對LZ壓縮后的literal、length、distance進行Huffman編碼,以及如何以最小空間記錄Huffman碼表。整個過程充滿了對數據結構尤其是樹的深入優化利用。按照上面的分析,如果要對ZIP進行進一步改進,可以考慮的地方也有不少,典型的有:

(1)擴大LZ編碼的滑動窗口的大小;

(2)將Huffman編碼改進為算術編碼等壓縮率更高的方法,畢竟,Huffman的碼字長度必須為整數,這就從理論上限制了它的壓縮率只能接近于理論極限,但難以達到。我記得在JPEG圖像編碼領域,以前的JPEG采用了DCT變換編碼+Huffman的方式,現在JPEG2000將其改為小波變換+算數編碼,所以數據壓縮也可以嘗試類似的思路;

(3)將多個文件進行合并壓縮,ZIP中,不同的文件壓縮過程沒有關系,獨立進行,如果將它們合并起來一起進行壓縮,壓縮率可以得到進一步提高。

?

描述分析有誤的地方,敬請指正。針對數據壓縮相關的話題,后續會對HBase列壓縮等等進行分析,看看ZIP這種文件壓縮和HBase這種數據庫數據壓縮的區別和聯系。

?

原文:?cnblog(esingchan)

總結

以上是生活随笔為你收集整理的ZIP压缩算法详细分析及解压实例的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。

狠狠激情中文字幕 | 精品1区2区3区 | 日韩精品视频在线观看网址 | 精品视频免费在线 | 国产一级片久久 | 婷婷久久网 | 美女视频黄频大全免费 | 国产成人精品免高潮在线观看 | 欧美一区二区三区在线视频观看 | 天天躁日日躁狠狠躁 | 人人草天天草 | 欧美日韩国产一二三区 | 亚洲黄色av | 国产成人精品一区在线 | 国产免费久久 | 午夜国产福利在线 | 日韩毛片在线免费观看 | 久久精品一二三 | 中文字幕有码在线 | 欧美日韩久 | 深夜精品福利 | 欧美在线91| 国产又粗又硬又爽的视频 | 日韩av免费在线看 | 亚洲三级黄色 | 日韩大片在线观看 | 国产一区影院 | 国产成人久久精品77777综合 | www亚洲视频 | 国产黄a三级三级三级三级三级 | 久久1电影院| 国内精品久久久久 | 97在线免费视频 | 久久66热这里只有精品 | 久久天天躁狠狠躁亚洲综合公司 | 久草精品国产 | 99热精品久久 | 色综合久久网 | 伊人天天操| 激情网五月天 | 丁香花在线视频观看免费 | av网站有哪些 | 91麻豆看国产在线紧急地址 | 成人av中文字幕在线观看 | 国产福利一区二区在线 | 色婷婷激情电影 | 美女网站色免费 | 国产成人精品一区二区三区福利 | 国产区免费在线 | 99精品久久久久久久久久综合 | 狠狠狠综合 | 亚洲天堂社区 | 午夜三级在线 | 日本中文字幕高清 | 人人爽人人av | 在线导航av | 久久国产欧美日韩精品 | 久99久精品视频免费观看 | 麻豆国产视频 | 日本中文不卡 | 日韩大片在线看 | 99热最新精品 | 日韩av成人在线 | 福利视频一二区 | 欧美在线视频日韩 | 超碰公开在线 | 久久久99精品免费观看 | 国产精品免费人成网站 | 亚洲精品国产电影 | 国产色网站 | 男女啪啪免费网站 | 久久人人爽人人人人片 | 久久国产精品99久久人人澡 | 国产美女视频网站 | av高清一区二区三区 | 久久久久成人精品免费播放动漫 | 免费在线观看av片 | 久久精品中文字幕少妇 | 免费网站在线观看成人 | 在线观看免费av网站 | av久久久| 婷婷综合影院 | 在线看国产日韩 | 99精品国产高清在线观看 | 久久久婷 | 97电影在线| 一区二区三区精品在线 | 亚洲天天看 | 久久五月婷婷丁香社区 | 日日操操操 | 亚洲免费a| 亚洲精品黄色在线观看 | 日韩在线视频精品 | 日韩精品中文字幕一区二区 | 欧美一区二区三区在线观看 | 欧美精品小视频 | 国产精品美女免费看 | 欧美激情精品久久久久久变态 | 麻花豆传媒mv在线观看 | 亚洲精品美女免费 | 亚洲日本va午夜在线电影 | 粉嫩av一区二区三区入口 | 久久精品日本啪啪涩涩 | 日韩欧三级 | 四虎天堂| 欧美精品久久久久久久久老牛影院 | 中文字幕在线一二 | 日韩免费高清 | 欧美色图亚洲图片 | 在线观看成人国产 | 午夜在线免费观看视频 | 精品三级av | 欧美日韩视频免费看 | 国产黄色电影 | 射九九| 91探花国产综合在线精品 | 亚洲精品女人久久久 | 午夜精品久久久久久久99热影院 | 91高清在线 | 久久免费视频网站 | 国产高清av免费在线观看 | 日韩黄色免费在线观看 | 在线观看成人国产 | 91激情视频在线观看 | 亚洲国产成人在线 | av软件在线观看 | 日韩有色| 国产精品18久久久久久不卡孕妇 | 天天天天综合 | 亚洲午夜久久久影院 | 亚洲www天堂com | 国产91精品看黄网站 | 欧美日韩在线观看一区 | 黄色片视频在线观看 | 午夜精品一区二区三区视频免费看 | 久久久人 | 毛片区 | 国产在线观看91 | 久久免费a | 国产一区观看 | 天天干天天射天天操 | 免费看国产视频 | 久久免费片| 久久国色夜色精品国产 | 69国产精品视频 | 91完整版 | 夜夜干天天操 | 五月激情姐姐 | 久久这里只有精品久久 | 少妇bbw搡bbbb搡bbb | 久久免费视频这里只有精品 | 五月婷婷色综合 | 在线观看免费观看在线91 | av东方在线| 久久久亚洲影院 | 亚洲一区二区三区毛片 | 国产午夜一级毛片 | 国产一级片播放 | 日韩三级不卡 | 人人干人人干人人干 | 国产精品igao视频网网址 | 蜜桃视频日本 | 国产精品欧美一区二区三区不卡 | 色五丁香 | 亚洲精品综合一二三区在线观看 | 免费毛片一区二区三区久久久 | 福利视频第一页 | 99久久婷婷国产精品综合 | 精品国产1区二区 | 五月婷婷毛片 | av一级免费 | 精品欧美一区二区在线观看 | 久久久久久草 | 亚洲三级在线播放 | 国产精品涩涩屋www在线观看 | 久草色在线观看 | 国产99久久九九精品 | 欧美a级一区二区 | 久久综合电影 | 亚洲资源在线 | 国产资源av | 婷婷5月色 | 91av中文字幕 | 亚洲国产欧美在线看片xxoo | 丁香五婷 | 亚洲三区在线 | 亚洲国产资源 | 日本电影久久 | 国产在线国产 | 国产日韩精品一区二区三区 | 91免费网址 | 成在人线av| 国色综合 | 日韩激情av在线 | 综合五月 | 国产成人精品一区二区三区在线 | 天天天综合网 | 99精品一区二区 | 色资源中文字幕 | 狠狠狠色| 在线观看完整版免费 | 精品不卡av | 国产又粗又猛又色又黄视频 | 婷婷色中文字幕 | 国产精品 欧美 日韩 | 免费十分钟 | 日日夜夜干| 免费午夜网站 | 国产视频精品久久 | 日韩在线视频精品 | 三级免费黄色 | 天天拍天天色 | 91精品国自产在线观看 | 中文字幕在线观看av | 综合网伊人 | 国产视频午夜 | 国产美女主播精品一区二区三区 | 美女网站视频一区 | 成人永久免费 | 亚洲传媒在线 | 天天草天天 | 在线视频1卡二卡三卡 | 国产黄色av影视 | 在线观看va| 最新av在线网站 | 91精品视频网站 | 丝袜精品视频 | 免费三级av | 久久免费播放视频 | 国产色网站 | 999一区二区三区 | 欧美久久久久久久久中文字幕 | 国产999久久久| 黄色毛片视频免费观看中文 | 精品亚洲免费 | 免费看的黄色片 | 精品v亚洲v欧美v高清v | 成人av资源网 | 色综合久久综合网 | 欧美一区二区在线免费看 | 午夜视频免费在线观看 | 五月天综合激情 | 在线视频 亚洲 | 中文字幕乱偷在线 | 日日操操 | 久久久久久美女 | 中文字幕在线观看完整版 | 在线看国产视频 | 黄色网址国产 | 天天操天天干天天爽 | 亚洲一区网站 | 日韩中文字幕视频在线观看 | a√资源在线 | 精品亚洲男同gayvideo网站 | 免费看黄视频 | 91禁看片 | 精品久久久亚洲 | 久久经典国产 | 18国产精品白浆在线观看免费 | 国产一区二区精品91 | 97精品超碰一区二区三区 | 黄色三级免费 | 久久亚洲福利 | 国产亚洲精品久久久久久网站 | 亚洲涩涩涩涩涩涩 | 日韩av电影中文字幕 | 国产一区二区在线影院 | 五月婷婷在线综合 | 欧美亚洲久久 | 婷婷色在线观看 | 日韩手机在线观看 | 亚洲久草网 | 国产伦精品一区二区三区无广告 | 亚洲美女视频在线 | av免费观看高清 | 国产高清在线免费观看 | 97久久久免费福利网址 | 亚洲最新精品 | 91你懂的 | 色网站在线看 | 亚洲精品裸体 | 国产一区二区在线免费播放 | www.777奇米 | 波多野结衣在线播放视频 | 色综合亚洲精品激情狠狠 | 四虎国产精品成人免费影视 | 国产成人久久精品77777综合 | 麻豆影视在线免费观看 | 天天插天天操天天干 | 国产免费又爽又刺激在线观看 | 欧洲精品亚洲精品 | 日韩成人中文字幕 | 久久久这里有精品 | 91电影福利 | 久草剧场| 91视频在线观看下载 | 精品久久影院 | 黄色精品网站 | 欧美精品久久久久久久久老牛影院 | 亚洲国产美女精品久久久久∴ | 亚洲国产剧情 | 久久久受www免费人成 | 国产精品成人一区二区三区 | 91久久奴性调教 | 97精品国产97久久久久久春色 | 看片一区二区三区 | 人人澡人人草 | 国产精品久久在线观看 | 综合天天 | 日韩视频一区二区在线 | 亚洲综合精品在线 | 免费网站在线观看成人 | 91麻豆网站 | 一区二区激情视频 | 一区二区欧美在线观看 | 在线不卡的av| 中文字幕乱码日本亚洲一区二区 | 日韩欧三级 | 久久亚洲综合色 | 国产99精品在线观看 | 西西大胆免费视频 | 成全在线视频免费观看 | 欧美日韩视频一区二区 | 欧美亚洲国产日韩 | 欧美精品国产精品 | 97色se| 日韩1级片 | 久久综合网色—综合色88 | 日韩一区二区三区免费电影 | 国产黄色视 | 99久久精品国 | 亚洲国产精品久久久久婷婷884 | 波多野结衣精品 | 青草视频免费观看 | 美女精品久久久 | 成人免费视频播放 | 在线免费观看黄色小说 | 狠狠干天天色 | 激情综合色综合久久综合 | 成人在线一区二区 | 韩国av不卡 | 国产在线视频资源 | 一区二区三区高清 | 日韩区在线观看 | 亚洲经典中文字幕 | 91精品视频网站 | 天天躁日日躁狠狠 | 欧美性脚交| 久久免费av电影 | 久亚洲| 天天干夜夜想 | 97在线观看免费视频 | 久久免费视频网 | 久99久精品视频免费观看 | 日日干夜夜干 | 久久久久亚洲精品男人的天堂 | 精品国内| 国产精品观看在线亚洲人成网 | 在线观看欧美成人 | 成人黄大片 | 三级在线视频观看 | 久草观看 | 国产区高清在线 | 成人在线免费视频 | 日本精品小视频 | 亚洲精品在线一区二区 | 精品一二三四五区 | 成人久久毛片 | 久久精品专区 | av观看免费在线 | 看国产黄色大片 | 又黄又爽免费视频 | 日韩一级电影网站 | 在线成人免费av | 精品天堂av| 色婷婷婷 | 国产精品h在线观看 | 午夜av不卡 | 国产精品99久久久久 | 久久精品导航 | www亚洲视频 | 91九色成人| 国产色小视频 | 中文字幕在线视频一区 | 免费在线观看成人av | 日本三级久久 | 国产成人精品在线 | 国产精品美女久久久久久久 | 久久久亚洲电影 | 一级片免费观看视频 | 中文字幕成人在线 | 成人在线观看免费 | 精品国产三级a∨在线欧美 免费一级片在线观看 | 国产丝袜 | 91九色在线观看视频 | 国产成人精品久 | 日本中文字幕在线电影 | 天天艹天天操 | 在线观看免费av网站 | 6080yy午夜一二三区久久 | 国产中文字幕在线 | 日日干夜夜爱 | 亚洲乱码精品 | 五月婷婷在线综合 | 在线一二区| 国产精品视频地址 | 91久草视频 | 日韩高清在线不卡 | 99国内精品久久久久久久 | 国产高清视频免费在线观看 | 色偷偷97| 91视频一8mav| 国产成人一区二区三区在线观看 | 视频一区亚洲 | 国产午夜麻豆影院在线观看 | 99免费在线播放99久久免费 | 国产一级免费片 | 91中文在线 | 四虎免费在线观看 | 丁香视频全集免费观看 | 在线之家免费在线观看电影 | 天天草天天干天天 | 中文字幕在线观看免费高清完整版 | 99精品在线免费在线观看 | 久久久久免费视频 | 国产亚洲精品久久久久久无几年桃 | 国产一区二区久久久 | 婷婷在线精品视频 | 国产一区二区三精品久久久无广告 | 亚洲伊人网在线观看 | 91免费高清 | 少妇性bbb搡bbb爽爽爽欧美 | 色婷婷激情电影 | 久久99热这里只有精品国产 | 国产中文字幕免费 | www.91成人 | 激情黄色av| 日韩中文字| 亚洲一区二区三区四区在线视频 | 成人动图 | 三级av在线免费观看 | 国产伦理久久精品久久久久_ | 波多在线视频 | 国产成人一级 | 欧美日韩一区二区三区不卡 | 久久综合婷婷国产二区高清 | av 一区二区三区四区 | 久久婷婷网 | 亚洲精品久久激情国产片 | 成人久久电影 | 国产亚洲小视频 | 国产最新视频在线 | 亚洲精品日韩一区二区电影 | 美腿丝袜一区二区三区 | 最近日韩中文字幕中文 | 久热色超碰 | 亚洲日本在线视频观看 | 中文av资源站 | 福利视频区 | 91黄色在线看| 国产亚洲va综合人人澡精品 | 日韩视频免费 | 97超碰成人在线 | 日韩在线免费视频 | 二区三区av| 中文字幕在线免费观看视频 | 六月色婷婷 | 久精品视频免费观看2 | 91超碰在线播放 | 国产一区视频在线观看免费 | 国产精品久久久久久模特 | 久久精品99国产精品亚洲最刺激 | 激情伊人 | 综合色综合 | 伊人天天狠天天添日日拍 | 国产成人一区二区三区影院在线 | 久久精品3 | 久久精品永久免费 | 日韩美精品视频 | 欧美污在线观看 | 中文av字幕在线观看 | 国产专区视频 | 国产亚洲观看 | 丁香六月国产 | 亚洲视频精品在线 | 免费毛片aaaaaa| 人人dvd| 综合天天色 | 欧美精品九九99久久 | 黄污网| 久久高视频 | 日韩av三区 | 国内亚洲精品 | 免费看三级网站 | 国产日韩视频在线播放 | 亚洲专区免费观看 | 一级成人免费视频 | 日韩二区在线观看 | 国产a网站| 一区二区三区日韩精品 | 久久精品日产第一区二区三区乱码 | 成人黄色大片在线免费观看 | 全久久久久久久久久久电影 | 三级黄色免费片 | 菠萝菠萝在线精品视频 | 一区二区在线影院 | 18国产精品白浆在线观看免费 | 成人av手机在线 | 中文字幕在线免费观看 | 国产美女无遮挡永久免费 | 国产在线传媒 | 色射爱| 成人久久18免费网站 | 国产综合福利在线 | 四虎成人精品在永久免费 | 97在线视 | 国产又粗又猛又黄又爽的视频 | 在线观看亚洲精品 | 久久久免费 | 精品视频亚洲 | 免费看黄在线观看 | 日韩欧美高清在线观看 | 久久香蕉电影 | 精品久久久一区二区 | 欧美日韩视频一区二区三区 | 日本在线观看一区二区三区 | 色综合人人 | 中文字幕在线播出 | 黄色在线视频网址 | www.com.日本一级 | av成人亚洲 | 国产中文在线播放 | 麻豆高清免费国产一区 | 日日夜夜综合 | 久久久夜色 | 国产中文 | 92精品国产成人观看免费 | 国内精品久久久久久中文字幕 | 日韩中文字幕免费看 | 日本不卡一区二区 | 免费看日韩片 | 91看片在线看片 | 国产v在线 | 男女啪啪免费网站 | 日韩欧美一区二区三区视频 | 91免费网站在线观看 | 日韩av影视在线 | 人人干人人添 | 高清精品视频 | 精品久久久久久亚洲综合网站 | 国产伦理一区二区三区 | 超碰在线公开免费 | 久久免费视频在线观看6 | 国产精品久久久久久久午夜 | 日韩欧美精品在线观看视频 | 欧美a级在线免费观看 | 亚洲精品久 | 夜夜操夜夜干 | 国产精品久久久久四虎 | 99精品国产aⅴ | 日韩在线精品 | 久久夜色精品国产亚洲aⅴ 91chinesexxx | 超碰免费公开 | 欧美一级片在线 | 成年人天堂com| 成人影视片 | 国内外成人免费在线视频 | 婷婷在线不卡 | 亚洲国产成人在线播放 | 在线国产精品视频 | 高清国产午夜精品久久久久久 | 99一区二区三区 | 亚洲精品高清视频 | 午夜精品久久久久久久99无限制 | se视频网址 | 国内揄拍国产精品 | 免费观看日韩av | 中文字幕一区二区三区乱码不卡 | 天天激情综合网 | www.夜夜| 97爱爱爱 | 911香蕉视频 | 天天干天天操天天干 | 一区二区精品久久 | 国产第一福利 | 亚洲日本色 | 麻豆视频一区 | 日本精品视频一区二区 | 五月天九九| 久久黄色a级片 | 亚洲黄色成人av | 在线观看亚洲精品视频 | 久久久久成人精品免费播放动漫 | 久久综合九色综合欧美就去吻 | 免费热情视频 | av永久网址 | 久草久草视频 | 99国产一区 | 久久激情五月激情 | 国产精品久久久久久av | 日韩欧美一区二区三区视频 | 色在线观看网站 | 成人久久久精品国产乱码一区二区 | 亚洲国产小视频在线观看 | 国产精品久久久久久久久久久免费 | 日日婷婷夜日日天干 | 蜜臀久久99精品久久久久久网站 | 少妇性色午夜淫片aaaze | 天天干.com| 欧美日韩免费在线观看视频 | 99亚洲精品视频 | 天天操天天操天天操天天操天天操天天操 | 国产淫片 | 国产黄色观看 | 在线观看 亚洲 | 亚洲欧美日韩不卡 | 麻豆一区在线观看 | 成人av在线网 | 97视频在线观看成人 | 中文字幕在线电影 | 亚洲国产午夜视频 | 国产小视频在线免费观看 | 三级毛片视频 | 免费在线观看黄网站 | 人人草网站 | 免费在线色视频 | www.亚洲精品在线 | 97精品国产91久久久久久 | 狠狠色丁香婷婷综合视频 | 亚洲无毛专区 | 天天爱天天射天天干天天 | av观看免费在线 | 欧美精品一区二区在线播放 | 九九精品毛片 | 在线视频成人 | 国产一区二区精品久久 | 天天综合视频在线观看 | 成人wwwxxx视频| 高清视频一区 | 国产免费久久 | 我要看黄色一级片 | 00av视频| 日韩免| 天天爱天天草 | 五月天国产| 亚洲特级片| 成人在线免费看视频 | 丁香花在线视频观看免费 | 久久视奸 | 国产999视频在线观看 | 精品久久免费看 | 国产91对白在线 | 亚洲视频 中文字幕 | 免费成人黄色av | 在线观看91久久久久久 | 香蕉网站在线观看 | 欧美日本啪啪无遮挡网站 | 特级西西www44高清大胆图片 | 久草视频在线免费看 | 日韩视频一区二区三区 | 免费性网站 | 精品国内| av日韩精品| 97福利视频 | 免费在线观看日韩 | 亚洲,国产成人av | 免费亚洲成人 | 在线看国产视频 | 亚洲国产成人精品电影在线观看 | 亚洲精品永久免费视频 | 精品国产一二三 | 麻豆一区在线观看 | 狠狠干狠狠艹 | 99久久这里只有精品 | 免费福利片2019潦草影视午夜 | 久久久久欧美精品 | 免费在线观看一区二区三区 | 精品视频亚洲 | 人人澡人人添人人爽一区二区 | 久久嗨| 亚洲国产免费看 | 国产午夜精品一区二区三区在线观看 | 免费在线观看污网站 | 激情开心色| 超碰成人免费电影 | 欧美综合在线视频 | 成人毛片久久 | 精品国产免费一区二区三区五区 | 一区三区视频在线观看 | 国产精彩视频一区 | 精品国产伦一区二区三区免费 | 国产一区91 | 日韩精品在线免费播放 | 天天插日日操 | 久久草视频 | 精品一区 在线 | 日韩欧美高清一区二区 | 午夜精品福利一区二区 | 国产麻豆视频在线观看 | 丁香六月婷婷 | 天天操天天操天天操天天 | 中文字幕 国产视频 | 九色精品免费永久在线 | 香蕉精品视频在线观看 | 永久免费精品视频网站 | 午夜999 | 亚洲精品免费在线 | 正在播放五月婷婷狠狠干 | 99久久精品国产一区二区三区 | 久久不射电影院 | 成人h在线| 日韩精品在线一区 | 五月天激情综合网 | 日韩欧美99 | 中文字幕在线观看视频免费 | 又黄又爽的免费高潮视频 | 日韩不卡高清 | 99色99| 最新免费av在线 | 婷婷av综合| 欧美激情视频在线观看免费 | 最新中文字幕在线播放 | 国产一区免费 | 国产成人黄色片 | 亚洲人久久 | 97人人爽 | 国产女教师精品久久av | 91日韩在线视频 | 亚洲激情在线观看 | 国产精品福利久久久 | 久久久综合电影 | 国产成人一区二区啪在线观看 | 精品日韩在线一区 | 久免费视频 | 中文字幕一区二区三区在线观看 | 91在线中文字幕 | 欧美伦理一区二区三区 | 亚洲精品白浆高清久久久久久 | 亚洲精欧美一区二区精品 | 一区二区三区四区精品 | 精品一区中文字幕 | 国产精品2018| 成人一区不卡 | 成人午夜久久 | 丰满少妇一级片 | 99在线精品免费视频九九视 | 热re99久久精品国产99热 | 久久久麻豆精品一区二区 | wwwwww黄| 天堂在线一区二区三区 | 丝袜美腿在线 | 99久热 | 久久这里有精品 | 亚洲精品国产电影 | 韩国精品福利一区二区三区 | 97香蕉久久国产在线观看 | 成年在线观看 | 999男人的天堂 | 亚洲一区二区麻豆 | 精品一二三四视频 | 狠狠操操| 日日夜夜天天久久 | 亚洲人在线7777777精品 | 亚洲精品乱码久久久久久按摩 | 亚洲欧美日韩国产一区二区三区 | 精品视频 | 97色婷婷| 日本三级在线观看中文字 | 日韩免费网址 | 亚洲精品456在线播放第一页 | 欧美韩国日本在线观看 | 五月婷婷在线综合 | 麻豆视频免费网站 | 一区二区av | 久久久久久久久久亚洲精品 | 最新黄色av网址 | 久久精品成人 | 欧美日韩国产精品久久 | 日日夜夜人人天天 | 婷婷色综合 | 丁香5月婷婷久久 | 九色精品免费永久在线 | 91色偷偷| 九月婷婷综合网 | 在线精品国产 | 麻豆视频在线免费看 | 久久久国产一区二区三区 | 最近中文字幕免费视频 | 日本精品视频免费 | 久久成人毛片 | 我要色综合天天 | 亚洲网站在线看 | 久久久久综合 | 日韩网站免费观看 | 欧美性受极品xxxx喷水 | 久久黄色精品视频 | 久久在草 | 婷婷久久一区 | 天天干天天拍天天操 | 超碰在线人人艹 | 婷五月天激情 | 91视频传媒 | 亚洲精品国偷自产在线91正片 | 成人在线视频免费看 | 国产粉嫩在线观看 | 96av在线| 中文免费观看 | 亚洲无吗av | 天天摸天天弄 | 玖玖玖影院 | 亚洲精品婷婷 | 麻豆超碰 | 日韩av电影手机在线观看 | 国产馆在线播放 | 麻豆传媒视频在线免费观看 | 99久久精品免费看国产 | 国产精品中文久久久久久久 | 日韩美女高潮 | 中文字幕在线观看第一区 | 黄色国产大片 | 国产一区二区三区网站 | 欧美色图亚洲图片 | 九九视频在线播放 | 欧美999| 国产精品自拍在线 | 国产精品精品视频 | 久久久久伊人 | av在线直接看 | 天天操月月操 | 精品国产伦一区二区三区观看体验 | 婷婷六月天丁香 | 亚洲精品88欧美一区二区 | 日韩中文字幕国产精品 | 色婷婷亚洲精品 | 国产伦精品一区二区三区高清 | 国产视频 亚洲视频 | 免费看三片 | 日韩高清dvd | 激情久久综合网 | 最近中文字幕视频网 | 国产高清免费在线播放 | 欧美日韩综合在线 | 久久综合中文字幕 | 最近免费中文字幕大全高清10 | 97操碰 | 色综合久久久久久久 | 69国产精品视频 | 天天干天天玩天天操 | 亚洲高清视频在线观看 | 亚洲成人黄色av | 亚洲国产人午在线一二区 | 国产精品99爱 | 91插插插网站 | 天天射日 | 97精品国产97久久久久久春色 | 二区三区毛片 | 日韩免费一级a毛片在线播放一级 | 成人三级网站在线观看 | 免费中文字幕在线观看 | 日韩欧美视频在线 | 91精品久久久久久久久久久久久 | 成人久久18免费网站麻豆 | 一级黄色大片在线观看 | 久久精品小视频 | 五月天久久精品 | 91麻豆精品国产91 | 在线免费观看视频a | 久久精品国产一区二区电影 | 91麻豆网| 日韩视频免费看 | 国产精品久久嫩一区二区免费 | 久久久免费精品 | 99视频免费播放 | 国产精品video | 日韩理论在线播放 | 日韩视频在线播放 | 丁香在线观看完整电影视频 | 亚洲精品资源 | 91视频在线自拍 | 99九九99九九九视频精品 | 色亚洲网 | 亚洲成人精品久久 | 91在线区| 欧美激情va永久在线播放 | 精品国产一二三四区 | 日本精品视频免费观看 | 天天操天天干天天综合网 | 国产精品免费观看国产网曝瓜 | 亚洲激情视频在线 | 中文在线资源 | 亚洲国产黄色片 | 天天色天天射天天干 | 久久久久在线观看 | 国产视频日韩 | 免费的黄色的网站 | 亚洲精品婷婷 | 日韩视频专区 | 国产精品久久久久久久久久99 | 国产成人精品av在线 | 新版资源中文在线观看 | 天天色综合久久 | 91麻豆精品国产91久久久无限制版 | 999抗病毒口服液 | 综合中文字幕 | 久久99久久99精品免费看小说 | 免费人成在线观看网站 | 成人av网站在线观看 | 国产一区在线播放 | 午夜免费福利视频 | 亚洲片在线观看 | 日本黄色免费大片 | 久久综合导航 | 中文字幕中文字幕在线中文字幕三区 | 欧美日韩一区二区在线观看 | 91av观看| 久久艹人人 | 国产成人综 | 久久综合九色综合欧美就去吻 | 国产美女精品人人做人人爽 | 在线亚洲人成电影网站色www | 日本一区二区三区免费看 | 最新色站| 免费日韩 精品中文字幕视频在线 | 精品久操 | 国产一区久久 | 成年人免费av | 久草视频中文在线 | 丝袜美腿av | 99久久99久久精品国产片果冰 | 亚洲精品合集 | 欧美日韩视频一区二区三区 | 在线小视频你懂的 | 久久婷婷精品 | 在线免费视频一区 | 亚洲国产精品999 | 久久久久久国产精品999 | 亚洲国产网站 | 91自拍91 | 欧美午夜性生活 | 91av看片| 国产中文字幕在线视频 | 天天爱天天舔 | 国产精品久久一区二区无卡 | www五月 | 在线精品亚洲 | 视频91| 91精品国产自产91精品 | 成人黄色资源 | 亚洲精品欧美专区 | 欧美亚洲一区二区在线 | 亚洲五月激情 | 久久精品专区 | 色先锋资源网 | 国产日韩高清在线 | 免费网站色 | 天天操天天摸天天射 | 久久免费观看少妇a级毛片 久久久久成人免费 | 欧美另类xxx | 狠狠网站| 国产一级不卡视频 | 91九色精品女同系列 | 波多野结衣在线播放一区 | 久久网站最新地址 | 久久久黄色av | 国产精品色视频 | 久久在线一区 | 天天干人人干 | 中文字幕在线影院 | 天天操天天干天天插 | 中文字幕在线视频国产 | 丁香婷婷射 | 国产98色在线 | 日韩 | 成全免费观看视频 | 日本不卡一区二区三区在线观看 | 人人澡人人澡人人 | 午夜精品区 | 五月天亚洲综合小说网 | 色婷婷综合成人av | 成人av资源在线 | 色网站免费在线看 | 午夜美女视频 | 色在线高清 | 亚洲国产999 | 亚洲精品欧美专区 | 欧美一区二视频在线免费观看 | 亚洲综合射 | 在线观看中文字幕dvd播放 | 色综合天天综合网国产成人网 | 日韩av一区二区在线播放 | 国产不卡视频在线 | 在线观看www. | 99re国产视频 | 在线天堂中文在线资源网 | 国产精品中文 | 久久国产麻豆 | 韩日电影在线免费看 | 欧美韩国日本在线 | 欧美日韩中文国产 |