本發明涉及一種消息完整性認證方法，尤其是涉及一種基于全同態mac的消息完整性認證方法。
背景技術：
現有的消息完整性認證方法是通過構造消息認證碼(messageauthenticationcode,mac)來對消息進行認證，其原理是在原始消息發送前，發送方通過散列函數計算得到原始消息的摘要值mac1，在雙方(發送方和接收方)共享的會話密鑰下，發送方由摘要值mac1獲得消息認證碼，然后將消息認證碼和原始消息一起發送；接收方收到消息認證碼和消息后，先用會話密鑰對消息認證碼進行解密，得到原始消息的摘要值mac1，再利用與發送方計算原始消息的摘要值mac1相同的散列函數計算所收到消息的摘要值mac2，如果mac1和mac2相等，則原始消息沒有被篡改，可認證消息的完整性。
但是，現有的消息完整性認證方法均采用的是原始消息在信道中傳輸，使消息存在被泄露的威脅，導致攻擊者可以自由篡改或創建信道中信息。bellare等基于消息認證碼提出了xor-mac認證方法，雖然在一定程度上降低了攻擊者成功偽造的概率，但是其所傳輸的仍然是原始消息，仍然存在消息泄露的風險，安全性能不夠高。
技術實現要素：
本發明所要解決的技術問題是提供一種將消息以密文的形式在信道中傳輸，能使發送方和接收方基于密文完成消息的認證，避免了消息被泄露的風險，安全性較高的基于全同態mac的消息完整性認證方法。
本發明解決上述技術問題所采用的技術方案為：一種基于全同態mac的消息完整性認證方法，包括以下步驟：
①將原始消息記為m，m∈{0，1}，利用全同態加密技術對m進行處理，得到密文c，c為128位二進制數據；
②采用md5算法對密文c進行處理，得到128位摘要值mac1；
③發送方將密文c和128位摘要值mac1通過信道進行傳輸，發送給接收方；
④接收方接收到密文和128位摘要值mac1，將接收方接收到的密文記為c1；
⑤采用與所述的步驟②相同的md5算法對密文c1進行處理，得到128位摘要值mac2；
⑥對128位摘要值mac1和128位摘要值mac2進行比較，如果兩者一致，則消息完整，驗證通過，如果兩者不一致，則消息已被篡改，驗證不通過。
所述的步驟②中采用md5算法對密文c進行處理的具體過程為：
a.在密文c的最低位即右側先用一個1補位，然后在補位后的1的右側用若干個0進行補位，得到中間消息，中間消息的長度為512的n倍，n為大于等于1的整數，將該中間消息從左邊第一位開始，以512位為單位進行分組，得到n個單元消息，將第k個單元消息記為nk，k＝1，2，…，n；將每個單元消息從左邊第一位開始，以32位為單位再次進行分組，每個單元消息被分為16個子單元，其中第k個單元消息nk中第j個子單元用mk-j表示，j＝1，2，…，16；
b.將md5算法中的四個32位的鏈接變量記為a、b、c、d，對a、b、c、d進行初始化，使a、b、c、d的初始值為標準幻數；
c.設置四個中間變量o、p、q、r，并對四個中間變量進行初始化：o＝a，p＝b，q＝c，r＝d；
d.依次采用第1個單元消息n1的16個子單元m1-1-m1-16至第n個單元消息nn的16個子單元mn-1-mn-16對四個中間變量o、p、q、r進行n次更新，每次更新中對四個中間變量o、p、q、r進行4輪循環更新，每一輪循環更新包括16步運算：
d-1.設定每次更新時，對四個中間變量o、p、q、r進行4輪循環更新的四個非線性函數為：
f(x,y,z)＝(x&y)|((～x)&z)(1)
g(x,y,z)＝(x&y)|(y&(～z))(2)
其中，&為與運算符號，|為或運算符號，～為非運算符號，⊕為異或運算符號。
d-2.設定第k次更新，對四個中間變量o、p、q、r進行第1輪循環更新時，與非線性函數(1)～(4)有關的操作運算，其中q為待更新變量，wk-j和tk-j分別為對四個中間變量o、p、q、r第k次更新時設定的已知參數，ff(q,x,y,z,mk-j,wk-j,tk-j)表示將x+(q+(f(x,y,z)+mk-j+tk-j)<<<wk-j)的值賦予q，即更新q的值，<<<wk-j表示循環左移wk-j位；在對o進行更新時，q＝o，x＝p，y＝q，z＝r；在對p進行更新時，q＝p，x＝q，y＝r，z＝o；在對q進行更新時，q＝q，x＝r，y＝o，z＝p；在對r進行更新時，q＝r，x＝o，y＝p，z＝q；
d-3.按照o、r、q、p、o、r、q、p、o、r、q、p、o、r、q、p的順序，對四個中間變量o、p、q、r進行第k次第1輪循環更新，結果如下：
第1步運算：o＝ff(o,p,q,r,mk-1,7,0xd76aa478)
第2步運算：r＝ff(r,o,p,q,mk-2,12,0xe8c7b756)
第3步運算：q＝ff(q,r,o,p,mk-3,17,0x242070db)
第4步運算：p＝ff(p,q,r,o,mk-4,22,0xc1bdceee)
第5步運算：o＝ff(o,p,q,r,mk-5,7,0xf57c0faf)
第6步運算：r＝ff(r,o,p,q,mk-6,12,0x4787c62a)
第7步運算：q＝ff(q,r,o,p,mk-7,17,0xa8304613)
第8步運算：p＝ff(p,q,r,o,mk-8,22,0xfd469501)
第9步運算：o＝ff(o,p,q,r,mk-9,7,0x698098d8)
第10步運算：r＝ff(r,o,p,q,mk-10,12,0x8b44f7af)
第11步運算：q＝ff(q,r,o,p,mk-11,17,0xffff5bb1)
第12步運算：p＝ff(p,q,r,o,mk-12,22,0x895cd7be)
第13步運算：o＝ff(o,p,q,r,mk-13,7,0x6b901122)
第14步運算：r＝ff(r,o,p,q,mk-14,12,0xfd987193)
第15步運算：q＝ff(q,r,o,p,mk-15,17,0xa679438e)
第16步運算：p＝ff(p,q,r,o,mk-16,22,0x49b40821)
其中，第k次第1輪循環更新的第1-16步運算中，o、p、q、r如果還沒有開始更新，采用其初始值進行計算，如果已更新，則采用其當前運算前被更新后的最新值進行計算；
d-4.按照o、r、q、p、o、r、q、p、o、r、q、p、o、r、q、p的順序，對四個中間變量o、p、q、r進行第k次第2輪循環更新，結果如下：
第1步運算：o＝gg(o,p,q,r,mk-2,5,0xf61e2562)
第2步運算：r＝gg(r,o,p,q,mk-7,9,0xc040b340)
第3步運算：q＝gg(q,r,o,p,mk-12,14,0x265e5a51)
第4步運算：p＝gg(p,q,r,o,mk-1,20,0xe9b6c7aa)
第5步運算：o＝gg(o,p,q,r,mk-6,5,0xd62f105d)
第6步運算：r＝gg(r,o,p,q,mk-11,9,0x02441453)
第7步運算：q＝gg(q,r,o,p,mk-16,14,0xd8a1e681)
第8步運算：p＝gg(p,q,r,o,mk-5,20,0xe7d3fbc8)
第9步運算：o＝gg(o,p,q,r,mk-10,5,0x21e1cde6)
第10步運算：r＝gg(r,o,p,q,mk-159,0xc33707d6)
第11步運算：q＝gg(q,r,o,p,mk-414,0xf4d50d87)
第12步運算：p＝gg(p,q,r,o,mk-9,20,0x455a14ed)
第13步運算：o＝gg(o,p,q,r,mk-14,5,0xa9e3e905)
第14步運算：r＝gg(r,o,p,q,mk-3,9,0xfcefa3f8)
第15步運算：q＝gg(q,r,o,p,mk-8,14,0x676f02d9)
第16步運算：p＝gg(p,q,r,o,mk-13,20,0x8d2a4c8a)
其中，第k次第2輪循環更新的第1-16步運算中，o、p、q、r采用其當前運算前被更新后的最新值進行計算；
d-5.按照o、r、q、p、o、r、q、p、o、r、q、p、o、r、q、p的順序，對四個中間變量o、p、q、r進行第k次第3輪循環更新，結果如下：
第1步運算：o＝hh(o,p,q,r,mk-6,4,0xfffa3942)
第2步運算：r＝hh(r,o,p,q,mk-9,11,0x8771f681)
第3步運算：q＝hh(q,r,o,p,mk-12,16,0x6d9d6122)
第4步運算：p＝hh(p,q,r,o,mk-15,23,0xfde5380c)
第5步運算：o＝hh(o,p,q,r,mk-2,4,0xa4beea44)
第6步運算：r＝hh(r,o,p,q,mk-5,11,0x4bdecfa9)
第7步運算：q＝hh(q,r,o,p,mk-8,16,0xf6bb4b60)
第8步運算：p＝hh(p,q,r,o,mk-11,23,0xbebfbc70)
第9步運算：o＝hh(o,p,q,r,mk-14,4,0x289b7ec6)
第10步運算：r＝hh(r,o,p,q,mk-1,11,0xeaa127fa)
第11步運算：q＝hh(q,r,o,p,mk-4,16,0xd4ef3085)
第12步運算：p＝hh(p,q,r,o,mk-7,23,0x04881d05)
第13步運算：o＝hh(o,p,q,r,mk-10,4,0xd9d4d039)
第14步運算：r＝hh(r,o,p,q,mk-13,11,0xe6db99e5)
第15步運算：q＝hh(q,r,o,p,mk-16,16,0x1fa27cf8)
第16步運算：p＝hh(p,q,r,o,mk-3,23,0xc4ac5665)
其中，第k次第3輪循環更新的第1-16步運算中，o、p、q、r采用其當前運算前被更新后的最新值進行計算；
d-6.按照o、r、q、p、o、r、q、p、o、r、q、p、o、r、q、p的順序，對四個中間變量o、p、q、r進行第k次第4輪循環更新，結果如下：
第1步運算：o＝ii(o,p,q,r,mk-1,6,0xf4292244)
第2步運算：r＝ii(r,o,p,q,mk-8,10,0x432aff97)
第3步運算：q＝ii(q,r,o,p,mk-15,15,0xab9423a7)
第4步運算：p＝ii(p,q,r,o,mk-6,21,0xfc93a039)
第5步運算：o＝ii(o,p,q,r,mk-13,6,0x655b59c3)
第6步運算：r＝ii(r,o,p,q,mk-4,10,0x8f0ccc92)
第7步運算：q＝ii(q,r,o,p,mk-11,15,0xffeff47d)
第8步運算：p＝ii(p,q,r,o,mk-2,21,0x85845dd1)
第9步運算：o＝ii(o,p,q,r,mk-9,6,0x6fa87e4f)
第10步運算：r＝ii(r,o,p,q,mk-16,10,0xfe2ce6e0)
第11步運算：q＝ii(q,r,o,p,mk-7,15,0xa3014314)
第12步運算：p＝ii(p,q,r,o,mk-14,21,0x4e0811a1)
第13步運算：o＝ii(o,p,q,r,mk-5,6,0xf7537e82)
第14步運算：r＝ii(r,o,p,q,mk-12,10,0xbd3af235)
第15步運算：q＝ii(q,r,o,p,mk-3,15,0x2ad7d2bb)
第16步運算：p＝ii(p,q,r,o,mk-10,21,0xeb86d391)
其中，第k次第4輪循環更新的第1-16步運算中，o、p、q、r采用其當前運算前被更新后的最新值進行計算；
d-7.按照步驟d-3至步驟d-6依次完成四個中間變量o、p、q、r的n次更新后，得到第n次4輪循環更新后中間變量o、p、q、r的最終更新值，將第n次4輪循環更新后o、p、q、r的最終更新值與其初始值分別相加，得到第n次4輪循環更新后的32位的鏈接變量數據源an、bn、cn、dn，an為第n次4輪循環更新后o的最新值與其初始值的和，bn為第n次4輪循環更新后p的最新值與其初始值的和，cn為第n次4輪循環更新后q的最新值與其初始值的和，dn為第n次4輪循環更新后r的最新值與其初始值的和。
d-8.將an、bn、cn、dn拼接，得到128位摘要值mac1，mac1＝dncnbnan，其中an為最低位，dn為最高位。采用該md5算法對密文c進行處理，對密文c的任何改動，哪怕改動很小，所計算出的128位摘要值mac2與未改動的密文c的128位摘要值mac1之間會存在很大的區別，同時對于一個已知的原始消息和其mac值，想要找到一個具有相同的mac值的原始消息(即偽造消息)是非常困難，并且采用該算法由密文c計算獲得128位摘要值mac1這個過程比較簡單。
與現有技術相比，本發明的優點在于通過將原始消息記為m，m∈{0，1}，利用全同態加密技術對m進行處理，得到密文c，c為128位二進制數據，采用md5算法對密文c進行處理，得到128位摘要值mac1，發送方將密文c和128位摘要值mac1通過信道進行傳輸，發送給接收方，接收方接收到密文和128位摘要值mac1，將接收方接收到的密文記為c1，采用md5算法對密文c1進行處理，得到128位摘要值mac2，對128位摘要值mac1和128位摘要值mac2進行比較，如果兩者一致，則消息完整，驗證通過，如果兩者不一致，則消息已被篡改，驗證不通過，本發明利用全同態允許對密文執行特定的操作后將其輸出進行解密，解密所得的結果與對明文(原始消息)進行相同操作的結果相等的特性，先對原始消息(即明文)進行了全同態加密，將消息以密文的形式在信道中傳輸，再對其構造消息認證碼來完成消息的認證，由此能使發送方和接收方基于密文完成消息的認證，認證過程是基于加密后的密文進行處理的，避免了原始消息在傳輸信道中的傳輸，可以有效地解決消息被泄露的風險，安全性較高。
具體實施方式
以下結合實施例對本發明作進一步詳細描述。
實施例一：一種基于全同態mac的消息完整性認證方法，包括以下步驟：
①將原始消息記為m，m∈{0，1}，利用全同態加密技術對m進行處理，得到密文c，c為128位二進制數據；
②采用md5算法對密文c進行處理，得到128位摘要值mac1；
③發送方將密文c和128位摘要值mac1通過信道進行傳輸，發送給接收方；
④接收方接收到密文和128位摘要值mac1，將接收方接收到的密文記為c1；
⑤采用與步驟②相同的md5算法對密文c1進行處理，得到128位摘要值mac2；
⑥對128位摘要值mac1和128位摘要值mac2進行比較，如果兩者一致，則消息完整，驗證通過，如果兩者不一致，則消息已被篡改，驗證不通過。
實施例二：一種基于全同態mac的消息完整性認證方法，包括以下步驟：
①將原始消息記為m，m∈{0，1}，利用全同態加密技術對m進行處理，得到密文c，c為128位二進制數據；
②采用md5算法對密文c進行處理，得到128位摘要值mac1；
③發送方將密文c和128位摘要值mac1通過信道進行傳輸，發送給接收方；
④接收方接收到密文和128位摘要值mac1，將接收方接收到的密文記為c1；
⑤采用與步驟②相同的md5算法對密文c1進行處理，得到128位摘要值mac2；
⑥對128位摘要值mac1和128位摘要值mac2進行比較，如果兩者一致，則消息完整，驗證通過，如果兩者不一致，則消息已被篡改，驗證不通過。
本實施例中，步驟②中采用md5算法對密文c進行處理的具體過程為：
a.在密文c的最低位即右側先用一個1補位，然后在補位后的1的右側用若干個0進行補位，得到中間消息，中間消息的長度為512的n倍，n為大于等于1的整數，將該中間消息從左邊第一位開始，以512位為單位進行分組，得到n個單元消息，將第k個單元消息記為nk，k＝1，2，…，n；將每個單元消息從左邊第一位開始，以32位為單位再次進行分組，每個單元消息被分為16個子單元，其中第k個單元消息nk中第j個子單元用mk-j表示，j＝1，2，…，16；
b.將md5算法中的四個32位的鏈接變量記為a、b、c、d，對a、b、c、d進行初始化，使a、b、c、d的初始值為標準幻數；
c.設置四個中間變量o、p、q、r，并對四個中間變量進行初始化：o＝a，p＝b，q＝c，r＝d；
d.依次采用第1個單元消息n1的16個子單元m1-1-m1-16至第n個單元消息nn的16個子單元mn-1-mn-16對四個中間變量o、p、q、r進行n次更新，每次更新中對四個中間變量o、p、q、r進行4輪循環更新，每一輪循環更新包括16步運算：
d-1.設定每次更新時，對四個中間變量o、p、q、r進行4輪循環更新的四個非線性函數為：
f(x,y,z)＝(x&y)|((～x)&z)(1)
g(x,y,z)＝(x&y)|(y&(～z))(2)
其中，&為與運算符號，|為或運算符號，～為非運算符號，
d-2.設定第k次更新，對四個中間變量o、p、q、r進行第1輪循環更新時，與非線性函數(1)～(4)有關的操作運算，其中q為待更新變量，wk-j和tk-j分別為對四個中間變量o、p、q、r第k次更新時設定的已知參數，ff(q,x,y,z,mk-j,wk-j,tk-j)表示將x+(q+(f(x,y,z)+mk-j+tk-j)<<<wk-j)的值賦予q，即更新q的值，<<<wk-j表示循環左移wk-j位；在對o進行更新時，q＝o，x＝p，y＝q，z＝r；在對p進行更新時，q＝p，x＝q，y＝r，z＝o；在對q進行更新時，q＝q，x＝r，y＝o，z＝p；在對r進行更新時，q＝r，x＝o，y＝p，z＝q；
d-3.按照o、r、q、p、o、r、q、p、o、r、q、p、o、r、q、p的順序，對四個中間變量o、p、q、r進行第k次第1輪循環更新，結果如下：
第1步運算：o＝ff(o,p,q,r,mk-1,7,0xd76aa478)
第2步運算：r＝ff(r,o,p,q,mk-2,12,0xe8c7b756)
第3步運算：q＝ff(q,r,o,p,mk-3,17,0x242070db)
第4步運算：p＝ff(p,q,r,o,mk-4,22,0xc1bdceee)
第5步運算：o＝ff(o,p,q,r,mk-5,7,0xf57c0faf)
第6步運算：r＝ff(r,o,p,q,mk-6,12,0x4787c62a)
第7步運算：q＝ff(q,r,o,p,mk-7,17,0xa8304613)
第8步運算：p＝ff(p,q,r,o,mk-8,22,0xfd469501)
第9步運算：o＝ff(o,p,q,r,mk-9,7,0x698098d8)
第10步運算：r＝ff(r,o,p,q,mk-10,12,0x8b44f7af)
第11步運算：q＝ff(q,r,o,p,mk-11,17,0xffff5bb1)
第12步運算：p＝ff(p,q,r,o,mk-12,22,0x895cd7be)
第13步運算：o＝ff(o,p,q,r,mk-13,7,0x6b901122)
第14步運算：r＝ff(r,o,p,q,mk-14,12,0xfd987193)
第15步運算：q＝ff(q,r,o,p,mk-15,17,0xa679438e)
第16步運算：p＝ff(p,q,r,o,mk-16,22,0x49b40821)
其中，第k次第1輪循環更新的第1-16步運算中，o、p、q、r如果還沒有開始更新，采用其初始值進行計算，如果已更新，則采用其當前運算前被更新后的最新值進行計算；
d-4.按照o、r、q、p、o、r、q、p、o、r、q、p、o、r、q、p的順序，對四個中間變量o、p、q、r進行第k次第2輪循環更新，結果如下：
第1步運算：o＝gg(o,p,q,r,mk-2,5,0xf61e2562)
第2步運算：r＝gg(r,o,p,q,mk-7,9,0xc040b340)
第3步運算：q＝gg(q,r,o,p,mk-12,14,0x265e5a51)
第4步運算：p＝gg(p,q,r,o,mk-1,20,0xe9b6c7aa)
第5步運算：o＝gg(o,p,q,r,mk-6,5,0xd62f105d)
第6步運算：r＝gg(r,o,p,q,mk-11,9,0x02441453)
第7步運算：q＝gg(q,r,o,p,mk-16,14,0xd8a1e681)
第8步運算：p＝gg(p,q,r,o,mk-5,20,0xe7d3fbc8)
第9步運算：o＝gg(o,p,q,r,mk-10,5,0x21e1cde6)
第10步運算：r＝gg(r,o,p,q,mk-159,0xc33707d6)
第11步運算：q＝gg(q,r,o,p,mk-414,0xf4d50d87)
第12步運算：p＝gg(p,q,r,o,mk-9,20,0x455a14ed)
第13步運算：o＝gg(o,p,q,r,mk-14,5,0xa9e3e905)
第14步運算：r＝gg(r,o,p,q,mk-3,9,0xfcefa3f8)
第15步運算：q＝gg(q,r,o,p,mk-8,14,0x676f02d9)
第16步運算：p＝gg(p,q,r,o,mk-13,20,0x8d2a4c8a)
其中，第k次第2輪循環更新的第1-16步運算中，o、p、q、r采用其當前運算前被更新后的最新值進行計算；
d-5.按照o、r、q、p、o、r、q、p、o、r、q、p、o、r、q、p的順序，對四個中間變量o、p、q、r進行第k次第3輪循環更新，結果如下：
第1步運算：o＝hh(o,p,q,r,mk-6,4,0xfffa3942)
第2步運算：r＝hh(r,o,p,q,mk-9,11,0x8771f681)
第3步運算：q＝hh(q,r,o,p,mk-12,16,0x6d9d6122)
第4步運算：p＝hh(p,q,r,o,mk-15,23,0xfde5380c)
第5步運算：o＝hh(o,p,q,r,mk-2,4,0xa4beea44)
第6步運算：r＝hh(r,o,p,q,mk-5,11,0x4bdecfa9)
第7步運算：q＝hh(q,r,o,p,mk-8,16,0xf6bb4b60)
第8步運算：p＝hh(p,q,r,o,mk-11,23,0xbebfbc70)
第9步運算：o＝hh(o,p,q,r,mk-14,4,0x289b7ec6)
第10步運算：r＝hh(r,o,p,q,mk-1,11,0xeaa127fa)
第11步運算：q＝hh(q,r,o,p,mk-4,16,0xd4ef3085)
第12步運算：p＝hh(p,q,r,o,mk-7,23,0x04881d05)
第13步運算：o＝hh(o,p,q,r,mk-10,4,0xd9d4d039)
第14步運算：r＝hh(r,o,p,q,mk-13,11,0xe6db99e5)
第15步運算：q＝hh(q,r,o,p,mk-16,16,0x1fa27cf8)
第16步運算：p＝hh(p,q,r,o,mk-3,23,0xc4ac5665)
其中，第k次第3輪循環更新的第1-16步運算中，o、p、q、r采用其當前運算前被更新后的最新值進行計算；
d-6.按照o、r、q、p、o、r、q、p、o、r、q、p、o、r、q、p的順序，對四個中間變量o、p、q、r進行第k次第4輪循環更新，結果如下：
第1步運算：o＝ii(o,p,q,r,mk-1,6,0xf4292244)
第2步運算：r＝ii(r,o,p,q,mk-8,10,0x432aff97)
第3步運算：q＝ii(q,r,o,p,mk-15,15,0xab9423a7)
第4步運算：p＝ii(p,q,r,o,mk-6,21,0xfc93a039)
第5步運算：o＝ii(o,p,q,r,mk-13,6,0x655b59c3)
第6步運算：r＝ii(r,o,p,q,mk-4,10,0x8f0ccc92)
第7步運算：q＝ii(q,r,o,p,mk-11,15,0xffeff47d)
第8步運算：p＝ii(p,q,r,o,mk-2,21,0x85845dd1)
第9步運算：o＝ii(o,p,q,r,mk-9,6,0x6fa87e4f)
第10步運算：r＝ii(r,o,p,q,mk-16,10,0xfe2ce6e0)
第11步運算：q＝ii(q,r,o,p,mk-7,15,0xa3014314)
第12步運算：p＝ii(p,q,r,o,mk-14,21,0x4e0811a1)
第13步運算：o＝ii(o,p,q,r,mk-5,6,0xf7537e82)
第14步運算：r＝ii(r,o,p,q,mk-12,10,0xbd3af235)
第15步運算：q＝ii(q,r,o,p,mk-3,15,0x2ad7d2bb)
第16步運算：p＝ii(p,q,r,o,mk-10,21,0xeb86d391)
其中，第k次第4輪循環更新的第1-16步運算中，o、p、q、r采用其當前運算前被更新后的最新值進行計算；
d-7.按照步驟d-3至步驟d-6依次完成四個中間變量o、p、q、r的n次更新后，得到第n次4輪循環更新后中間變量o、p、q、r的最終更新值，將第n次4輪循環更新后o、p、q、r的最終更新值與其初始值分別相加，得到第n次4輪循環更新后的32位的鏈接變量數據源an、bn、cn、dn，an為第n次4輪循環更新后o的最新值與其初始值的和，bn為第n次4輪循環更新后p的最新值與其初始值的和，cn為第n次4輪循環更新后q的最新值與其初始值的和，dn為第n次4輪循環更新后r的最新值與其初始值的和；
d-8.將an、bn、cn、dn拼接，得到128位摘要值mac1，mac1＝dncnbnan，其中an為最低位，dn為最高位。