Twofish

Twofish הוא צופן בלוקים סימטרי שפותח ב-1998 בחברת Counterpane Labs על ידי צוות קריפטוגרפים בראשות ברוס שנייר, על מנת לשמש כמועמד לתקן ההצפנה המתקדם. דורג שלישי בחמשת המועמדים המובילים והפסיד לריינדל. האלגוריתם פועל על בלוקים בגודל 128 סיביות וטווח מפתח מ-128 עד 256 סיביות והוא ממוטב למעבדי 32 סיביות. בעת התחרות לתקן הוכרז על ידי NIST כמו לגבי יתר המועמדים המובילים שלא התגלו בו חולשות ולא ידוע על התקפה יעילה נגדו. הצופן נחשב כממשיכו של Blowfish וגם הוא חופשי לשימוש ואינו מוגן בזכויות יוצרים או פטנט. הוא בנוי בסגנון רשת פייסטל בשישה-עשר סבבים וכולל; פונקציה חד-חד-ערכית ועל המיוצגת על ידי ארבע תיבות החלפה (S-box) לא ליניאריות תלויות מפתח, בגודל 8x8 סיביות. הכפלה במטריצת MDS קיצור של maximum distance separable, שהיא פונקציה ליניארית קבועה 4x4 מעל שדה סופי F 2 8 {\displaystyle F_{2^{8}}} וצמצום בפולינום פרימיטיבי קבוע. התמרת פסבדו-הדמר PHT. פונקציית ערבוב הפועלת על זוגות משתנים בני 32 סיביות כדלהלן; בהינתן a , b {\displaystyle a,b} הפלט הוא a ′ = ( a + b )  mod  2 32 {\displaystyle a'=(a+b){\mbox{ mod }}2^{32}} ו- b ′ = ( a + 2 b )  mod  2 32 {\displaystyle b'=(a+2b){\mbox{ mod }}2^{32}} . הזזה מעגלית בסיביות (bitwise rotation), כלומר לאחר כל הזזה סיבית הגלישה הנפלטת מצד אחד מוחזרת מהצד השני. הלבנה שהיא XOR עם 128 סיביות מהמפתח לפני הסבב הראשון ואחרי הסבב האחרון. תהליך הרחבת מפתח המתואר להלן. יישום ממוטב של האלגוריתם במחשב 32 סיביות צורך כ-18 מחזורי שעון לבית או 1820 מחזורי שעון במעבד 8 סיביות. Twofish ניתן ליישום בחומרה בעלות של 14,000 שערים. תהליך הרחבת המפתח וכן סבב ההצפנה מתחשב במגוון אפשרויות כדי לאפשר איזון בין ביצועים לביטחון בהתאם לצורך. ההתקפה הטובה ביותר הידועה כנגד האלגוריתם היא התקפת גלוי נבחר עם 2 22.5 {\displaystyle 2^{22.5}} טקסטים ו- 2 51 {\displaystyle 2^{51}} ניסיונות בחמשה סבבים בלבד. עם מפתח בגודל 128 סיביות יישום בתוכנה של Twofish איטי מעט בהשוואה לריינדל ואילו עם מפתח 256 סיביות קיים יתרון קל ל-Twofish. אלגוריתם Twofish פותח במיוחד לפי דרישות תקן ההצפנה של המכון הלאומי לתקנים וטכנולוגיה של ארצות הברית: אורך הבלוק נקבע ל-128 סיביות. אורכי מפתח נקבעו ל-128, 192 או 256 סיביות. יעילות הן בחומרה והן בתוכנה. גמישות; יישום במגוון אפשרויות על פלטפורמות שונות ועם אורכי מפתח שונים וכן התאמה לשימוש כצופן זרם, פונקציית גיבוב או קוד אימות מסרים. מבנה פשוט המסייע ביישום ובניתוח האלגוריתם. מהירות הצפנה של כ-500 מחזורי שעון בגרסה הממוטבת. פעולות אלגבריות פשוטות שאינן מערימות קשיים ביישום בתוכנה או בחומרה. ביטחון; ההתקפה הטובה ביותר האפשרית כנגד האלגוריתם במלוא הסבבים תהיה בסיבוכיות של לא פחות מ- 2 80 {\displaystyle 2^{80}} . קלט: 128 סיביות טקסט קריא המיוצגים על ידי 16 בתים p 0 , p 1 , . . . , p 15 {\displaystyle p_{0},p_{1},...,p_{15}} , מפתח M {\displaystyle M} בגודל 128, 192 או 256 סיביות. פלט: 16 בתי צופן c 0 , c 1 , . . . , c 15 {\displaystyle c_{0},c_{1},...,c_{15}} . 1. הקלט מפוצל לארבע מילים, P 0 , P 1 , P 2 , P 3 {\displaystyle P_{0},P_{1},P_{2},P_{3}} של 32 סיביות כל אחת לפי סדר בתים קטן (little-endian) כלומר הספרה הפחות משמעותית מוצבת בבית הראשון שבכתובת הזיכרון. הנוסחה היא: P i = ∑ j = 0 3 p ( 4 i + j ) ⋅ 2 8 j ,           i = 0 , . . . , 3 {\displaystyle \textstyle P_{i}=\sum _{j=0}^{3}p_{(4i+j)}\cdot 2^{8_{j}},\ \ \ \ \ i=0,...,3} (הערה: במעבד עם סדר בתים קטן (כמו בארכיטקטורת x86), בשפות C ו-C++‎ אין צורך בפעולת חישוב כלשהי בהמרה זו ובהמרה האחרונה בסעיף 6, כיוון שניתן להתייחס למצביע הקלט כאל כמצביע לשלם-32 במקום לבתים או להשתמש במבנה Union, כך שסיביות הקלט והפלט נותרות במקומן). 2. שלב הלבנה; מבצעים XOR עם 4 מילות מהמפתח המורחב K {\displaystyle K} (הכנת המפתח בהמשך): R 0 , i = P i ⊕ K i ,           i = 0 , . . . , 3 {\displaystyle R_{0,i}=P_{i}\oplus K_{i},\ \ \ \ \ i=0,...,3} 3. בכל סבב מ-16 הסבבים, שתי המילים הראשונות R 0 , R 1 {\displaystyle R_{0},R_{1}} משמשות כקלט לפונקציה F {\displaystyle F} (להלן) שמקבלת כפרמטר גם את מספר הסבב, המילה השלישית R 2 {\displaystyle R_{2}} מחוברת ב-XOR עם התוצאה הראשונה של F {\displaystyle F} ולאחר מכן מסובבת בהזזה מעגלית סיבית אחת ימינה. המילה הרביעית R 3 {\displaystyle R_{3}} מסובבת בהזזה מעגלית שמאלה סיבית אחת ואז מחוברת ב-XOR עם התוצאה השנייה של F {\displaystyle F} . לאחר מכן שני החצאים מחליפים מקומות וחוזר חלילה. באופן פורמלי עבור r = 0 {\displaystyle r=0} עד 15 {\displaystyle 15} מבצעים: 3.1 ( F r , 0 , F r , 1 ) = F ( R r , 0 , R r , 1 , r ) {\displaystyle (F_{r,0},F_{r,1})=F(R_{r,0},R_{r,1},r)} 3.2 R r + 1 , 0 = ROR 1 ( R r , 2 ⊕ F r , 0 ) {\displaystyle R_{r+1,0}={\mbox{ROR}}_{1}(R_{r,2}\oplus F_{r,0})} 3.3 R r + 1 , 1 = ROL 1 ( R r , 3 ) ⊕ F r , 1 {\displaystyle R_{r+1,1}={\mbox{ROL}}_{1}(R_{r,3})\oplus F_{r,1}} 3.4 R r + 1 , 2 = R r , 0 {\displaystyle R_{r+1,2}=R_{r,0}} 3.5 R r + 1 , 3 = R r , 1 {\displaystyle R_{r+1,3}=R_{r,1}} כאן ROR מייצג Rotate Right שהיא פונקציית הזזה מעגלית לימין במספר סיביות לפי המציין התחתי ו-ROL מייצג Rotate Left. 4. מבצעים החלפה נוספת בין שני חצאי הקלט (למעשה מבטלים את ההחלפה מהסבב האחרון) 5. מבצעים הלבנה נוספת עם ארבע מילות מפתח מורחב נוספות. C i = R 16 , ( i + 2 )  mod  4 ⊕ K i + 4 ,         i = 0 , . . . , 3 {\displaystyle C_{i}=R_{16,(i+2){\mbox{ mod }}4}\oplus K_{i+4},\ \ \ \ i=0,...,3} 6. התוצאה מומרת בחזרה ל-16 בתים לפי אותו סדר בתים: c i = ⌊ C ⌊ i / 4 ⌋ 2 8 ( i  mod  4 ) ⌋  mod  2 8 ,         i = 0 , . . . , 15 {\displaystyle c_{i}=\left\lfloor {\frac {C_{\left\lfloor i/4\right\rfloor }}{2^{8(i{\mbox{ mod }}4)}}}\right\rfloor {\mbox{ mod }}2^{8},\ \ \ \ i=0,...,15} הפונקציה F {\displaystyle F} מקבלת שלושה פרמטרים: שתי מילים R 0 , R 1 {\displaystyle R_{0},R_{1}} ומספר הסבב r {\displaystyle r} שמהווה אינדקס לערך המתאים במפתח המורחב. R 0 {\displaystyle R_{0}} משמש קלט לפונקציית g {\displaystyle g} שמתוארת להלן המחזירה את T 0 {\displaystyle T_{0}} . ו- R 1 {\displaystyle R_{1}} מסובב בהזזה מעגלית 8 סיביות שמאלה ואז משמש גם הוא קלט לפונקציה g {\displaystyle g} שמחזירה את T 1 {\displaystyle T_{1}} . התוצאות משולבות יחד באמצעות פונקציית PHT בה מוסיפים שתי מילים מתאימות מהמפתח המורחב לפי r {\displaystyle r} , הפלט הוא הזוג ( F 0 , F 1 ) {\displaystyle (F_{0},F_{1})} ובניסוח פורמלי: T 0 = g ( R 0 ) {\displaystyle T_{0}=g(R_{0})} T 1 = g ( ROL 8 ( R 1 ) ) {\displaystyle T_{1}=g({\mbox{ROL}}_{8}(R_{1}))} F 0 = ( T 0 + T 1 + K 2 r + 8 )  mod  2 32 {\displaystyle F_{0}=(T_{0}+T_{1}+K_{2r+8}){\mbox{ mod }}2^{32}} F 1 = ( T 0 + 2 T 1 + K 2 r + 9 )  mod  2 32 {\displaystyle F_{1}=(T_{0}+2T_{1}+K_{2r+9}){\mbox{ mod }}2^{32}} הפונקציה g {\displaystyle g} מקבלת כפרמטר מילה אחת X {\displaystyle X} בגודל 32 סיביות, מפצלת אותה לארבעה בתים, כל אחד מהם משמש אינדקס לתיבת ההחלפה המתאימה שמחזירה בית אחר. מייצגים את התוצאות כווקטור באורך 4 מעל השדה F 2 8 {\displaystyle F_{2^{8}}} ומכפילים במטריצת MDS בגודל 4x4 (ראו תרשים). התוצאה מומרת בחזרה למילה אחת: MDS = ( 01 EF 5B 5B 5B EF EF 01 EF 5B 01 EF EF 01 EF 5B ) {\displaystyle {\mbox{MDS}}={\begin{pmatrix}{\mbox{01 EF 5B 5B}}\\{\mbox{5B EF EF 01}}\\{\mbox{EF 5B 01 EF}}\\{\mbox{EF 01 EF 5B}}\end{pmatrix}}} x i = ⌊ X / 2 8 i ⌋  mod  2 8           i = 0 , . . . , 3 {\displaystyle x_{i}=\lfloor X/2^{8i}\rfloor {\mbox{ mod }}2^{8}\ \ \ \ \ i=0,...,3} y i = s i [ x i ]           i = 0 , . . . , 3 {\displaystyle y_{i}=s_{i}[x_{i}]\ \ \ \ \ i=0,...,3} ( z 0 z 1 z 2 z 3 ) = ( ⋅ ⋯ ⋅ ⋮ MDS ⋮ ⋅ ⋯ ⋅ ) ⋅ ( y 0 y 1 y 2 y 3 ) {\displaystyle {\begin{pmatrix}z_{0}\\z_{1}\\z_{2}\\z_{3}\end{pmatrix}}={\begin{pmatrix}\cdot &\cdots &\cdot \\\vdots &{\mbox{MDS}}&\vdots \\\cdot &\cdots &\cdot \end{pmatrix}}\cdot {\begin{pmatrix}y_{0}\\y_{1}\\y_{2}\\y_{3}\end{pmatrix}}} כאשר s i {\displaystyle s_{i}} היא תיבת החלפה ו- Z {\displaystyle Z} הוא התוצאה של g {\displaystyle g} . לצורך המרה בין אלמנטים בשדה F 2 8 {\displaystyle F_{2^{8}}} לבין בתי הקלט משתמשים בפולינום הפרימיטיבי v ( x ) = x 8 + x 6 + x 5 + x 3 + 1 {\displaystyle v(x)=x^{8}+x^{6}+x^{5}+x^{3}+1} ממעלה 8. האלמנט a = ∑ i = 0 7 a i x i {\displaystyle a=\textstyle \sum _{i=0}^{7}a_{i}x^{i}} כאשר מקדמיו בינאריים, מתאים לבית ∑ i = 0 7 a i 2 i {\displaystyle \textstyle \sum _{i=0}^{7}a_{i}2^{i}} . זהו מיפוי טבעי במובן שפעולות חיבור הן בעצם XOR בלבד. לצורך יעילות ביישום הן בתוכנה והן בחומרה מיישמים את הכפל במטריצה באמצעות סדרת פעולות XOR ו-Shift תוך שימוש בטכניקה שנקראת LFSR, כדלהלן: z 0 = y 0     ⊕     ( y 1     ⊕     L1 ( y 1 )     ⊕     L2 ( y 1 ) )     ⊕     ( ( y 2     ⊕     L2 ( y 2 ) )     ⊕     ( ( y 3     ⊕     L2 ( y 3 ) ) {\displaystyle z_{0}=y_{0}\ \ \oplus \ \ (y_{1}\ \ \oplus \ \ {\mbox{L1}}(y_{1})\ \ \oplus \ \ {\mbox{L2}}(y_{1}))\ \ \oplus \ \ ((y_{2}\ \ \oplus \ \ {\mbox{L2}}(y_{2}))\ \ \oplus \ \ ((y_{3}\ \ \oplus \ \ {\mbox{L2}}(y_{3}))} z 1 = ( y 0     ⊕     L2 ( y 0 ) )     ⊕     ( y 1     ⊕     L1 ( y 1 )     ⊕     L2 ( y 1 ) )     ⊕     ( y 2     ⊕     L1 ( y 2 )     ⊕     L2 ( y 2 ) )     ⊕     y 3 {\displaystyle z_{1}=(y_{0}\ \ \oplus \ \ {\mbox{L2}}(y_{0}))\ \ \oplus \ \ (y_{1}\ \ \oplus \ \ {\mbox{L1}}(y_{1})\ \ \oplus \ \ {\mbox{L2}}(y_{1}))\ \ \oplus \ \ (y_{2}\ \ \oplus \ \ {\mbox{L1}}(y_{2})\ \ \oplus \ \ {\mbox{L2}}(y_{2}))\ \ \oplus \ \ y_{3}} z 2 = ( y 0     ⊕     L1 ( y 0 )     ⊕     L2 ( y 0 ) )     ⊕     ( ( y 1     ⊕     L2 ( y 1 ) )     ⊕     y 2     ⊕     ( y 3     ⊕     L1 ( y 3 )     ⊕     L2 ( y 3 ) ) {\displaystyle z_{2}=(y_{0}\ \ \oplus \ \ {\mbox{L1}}(y_{0})\ \ \oplus \ \ {\mbox{L2}}(y_{0}))\ \ \oplus \ \ ((y_{1}\ \ \oplus \ \ {\mbox{L2}}(y_{1}))\ \ \oplus \ \ y_{2}\ \ \oplus \ \ (y_{3}\ \ \oplus \ \ {\mbox{L1}}(y_{3})\ \ \oplus \ \ {\mbox{L2}}(y_{3}))} z 3 = ( y 0     ⊕     L1 ( y 0 )     ⊕     L2 ( y 1 ) )     ⊕     y 1     ⊕     ( y 2     ⊕     L1 ( y 2 )     ⊕     L2 ( y 2 ) )     ⊕     ( ( y 3     ⊕     L2 ( y 3 ) ) {\displaystyle z_{3}=(y_{0}\ \ \oplus \ \ {\mbox{L1}}(y_{0})\ \ \oplus \ \ {\mbox{L2}}(y_{1}))\ \ \oplus \ \ y_{1}\ \ \oplus \ \ (y_{2}\ \ \oplus \ \ {\mbox{L1}}(y_{2})\ \ \oplus \ \ {\mbox{L2}}(y_{2}))\ \ \oplus \ \ ((y_{3}\ \ \oplus \ \ {\mbox{L2}}(y_{3}))} והפונקציות L1, L2 הן אוגרי זיזה (LFSR) הראשון מבוסס על הסיבית הראשונה והשני על הסיביות הראשונה והשנייה והן מוגדרות כך: L1 ( x ) = ( x ≫ 1 )     ⊕     ( x     AND     0x01     ?     0xB4     :     0 ) {\displaystyle {\mbox{L1}}(x)=(x\gg 1)\ \ \oplus \ \ (x\ \ {\mbox{AND}}\ \ {\mbox{0x01}}\ \ ?\ \ {\mbox{0xB4}}\ \ :\ \ 0)} L2 ( x ) = ( x ≫ 2 )     ⊕     ( x     AND     0x02     ?     0xB4     :     0 )     ⊕     ( x     AND     0x01     ?     0x5A     :     0 ) {\displaystyle {\mbox{L2}}(x)=(x\gg 2)\ \ \oplus \ \ (x\ \ {\mbox{AND}}\ \ {\mbox{0x02}}\ \ ?\ \ {\mbox{0xB4}}\ \ :\ \ 0)\ \ \oplus \ \ (x\ \ {\mbox{AND}}\ \ {\mbox{0x01}}\ \ ?\ \ {\mbox{0x5A}}\ \ :\ \ 0)} (הסימן ≫ {\displaystyle \gg } מייצג הזזה לימין והסימן AND שנקרא וגם בודק האם הסיבית המתאימה דולקת, 0x02 היא סיבית LSB השנייה). התוצאה היא המרה של ערכי z i {\displaystyle z_{i}} למילה בגודל 32 סיביות בסדר בתים קטן לפי הנוסחה Z = ∑ i = 0 3 z i ⋅ 2 8 i {\displaystyle \textstyle Z=\sum _{i=0}^{3}z_{i}\cdot 2^{8i}} (במעבד עם סדר בתים קטן, התוצאה מתקבלת במבנה הרצוי ואין צורך בהמרה כלשהי). המפתח המורחב מכיל מערך של 40 מילים K 0 , . . . , K 39 {\displaystyle K_{0},...,K_{39}} וארבע תיבות החלפה מהפונקציה g {\displaystyle g} . המפתח כאמור יכול להיות N = 128 {\displaystyle N=128} , N = 192 {\displaystyle N=192} או N = 256 {\displaystyle N=256} סיביות. לצורך הכנת המפתח המורחב יהי k = N / 64 {\displaystyle k=N/64} , מכינים מהמפתח מערך M {\displaystyle M} בגודל 8 k {\displaystyle 8k} בתים m 0 , . . . , m 8 k − 1 {\displaystyle m_{0},...,m_{8k-1}} (מרפדים באפסים אם נדרש), ממירים ל- 2 k {\displaystyle 2k} מילים בגודל 32 סיביות לפי סדר בתים little-endian. ואז ממירים את M {\displaystyle M} לשני וקטורים (זוגי ואי זוגי) כלומר M e = ( M 0 , M 2 , . . . , M 2 k − 2 ) {\displaystyle M_{e}=(M_{0},M_{2},...,M_{2k-2})} וכן M o = ( M 1 , M 3 , . . . , M 2 k − 1 ) {\displaystyle M_{o}=(M_{1},M_{3},...,M_{2k-1})} . וקטור נוסף בגודל k {\displaystyle k} מיוצר מהמפתח בקבוצות של 8 סיביות, אותם מפרשים כווקטור מעל F 2 8 {\displaystyle F_{2^{8}}} ומכפילים במטריצה 4x8 של קוד RS (ריד-סולומון) כמתואר בתרשים להלן. פעולת הכפל במטריצה נראית כך: RS = ( 01 A4 55 87 5A 58 DB 9E A4 56 82 F3 1E C6 68 E5 02 A1 FC C1 47 AE 3D 19 A4 55 87 5A 58 DB 9E 03 ) {\displaystyle {\text{RS}}={\begin{pmatrix}{\text{01 A4 55 87 5A 58 DB 9E}}\\{\text{A4 56 82 F3 1E C6 68 E5}}\\{\text{02 A1 FC C1 47 AE 3D 19}}\\{\text{A4 55 87 5A 58 DB 9E 03}}\end{pmatrix}}} ( s i , 0 s i , 1 s i , 2 s i , 3 ) = ( ⋅ ⋯ ⋅ ⋮ RS ⋮ ⋅ ⋯ ⋅ ) ⋅ ( m 8 i m 8 i + 1 m 8 i + 2 m 8 i + 3 m 8 i + 4 m 8 i + 5 m 8 i + 6 m 8 i + 7 ) {\displaystyle {\begin{pmatrix}s_{i,0}\\s_{i,1}\\s_{i,2}\\s_{i,3}\end{pmatrix}}={\begin{pmatrix}\cdot &\cdots &\cdot \\\vdots &{\mbox{RS}}&\vdots \\\cdot &\cdots &\cdot \end{pmatrix}}\cdot {\begin{pmatrix}m_{8i}\\m_{8i+1}\\m_{8i+2}\\m_{8i+3}\\m_{8i+4}\\m_{8i+5}\\m_{8i+6}\\m_{8i+7}\end{pmatrix}}} ואז מחשבים עבור i = 0 , . . . , k − 1 {\displaystyle i=0,...,k-1} בסדר הפוך S = ( S k − 1 , S k − 2 , . . . , S 0 ) {\displaystyle S=(S_{k-1},S_{k-2},...,S_{0})} : S i = ∑ j = 0 3 s i , j ⋅ 2 8 j {\displaystyle S_{i}=\sum _{j=0}^{3}s_{i,j}\cdot 2^{8j}} להכפלה במטריצה בשיטת ריד-סולומון, משתמשים בפולינום פרימיטיבי w ( x ) = x 8 + x 6 + x 3 + x 2 + 1 {\displaystyle w(x)=x^{8}+x^{6}+x^{3}+x^{2}+1} . שלושת הווקטורים M e , M o , S {\displaystyle M_{e},M_{o},S} הם הבסיס לתהליך הרחבת המפתח. הפונקציה המשמשת להרחבת המפתח מקבלת פרמטר X {\displaystyle X} בגודל 32 סיביות ורשימה L = ( L 0 , . . . , L k − 1 ) {\displaystyle L=(L_{0},...,L_{k-1})} של מילות 32 סיביות באורך k {\displaystyle k} ומפיקה מילה אחת. הפונקציה פועלת באופן איטרטיבי k {\displaystyle k} שלבים, בכל שלב 4 בתים עוברים בתיבות החלפה קבועות ומחוברים ב-XOR עם בתים מהרשימה L {\displaystyle L} ולסיום מועברים שוב בתיבות ההחלפה ומוכפלים במטריצת MDS. כדלהלן: תחילה מפצלים לבתים כאשר i = 0 , . . . , k − 1 {\displaystyle i=0,...,k-1} ו- j = 0 , . . . , 3 {\displaystyle j=0,...,3} : l u , j = ⌊ L i / 2 8 j ⌋  mod  2 8 {\displaystyle l_{u,j}=\left\lfloor L_{i}/2^{8j}\right\rfloor {\mbox{ mod }}2^{8}} x j = ⌊ X / 2 8 j ⌋  mod  2 8 {\displaystyle x_{j}=\left\lfloor X/2^{8j}\right\rfloor {\mbox{ mod }}2^{8}} בשלב החלפה ו-XOR תחילה מציבים y k , j = x j {\displaystyle y_{k,j}=x_{j}} כאשר j = 0 , . . . , 3 {\displaystyle j=0,...,3} ואז מחשבים לפי ערך k {\displaystyle k} : כאשר q 0 , q 1 {\displaystyle q_{0},q_{1}} הן תמורות קבועות של ערכים בגודל 8 סיביות המוגדרות להלן. הווקטור y i {\displaystyle y_{i}} מוכפל במטריצת MDS כמו בפונקציה g {\displaystyle g} לעיל: ( z 0 z 1 z 2 z 3 ) = ( ⋅ ⋯ ⋅ ⋮ MDS ⋮ ⋅ ⋯ ⋅ ) ⋅ ( y 0 y 1 y 2 y 3 ) {\displaystyle {\begin{pmatrix}z_{0}\\z_{1}\\z_{2}\\z_{3}\end{pmatrix}}={\begin{pmatrix}\cdot &\cdots &\cdot \\\vdots &{\mbox{MDS}}&\vdots \\\cdot &\cdots &\cdot \end{pmatrix}}\cdot {\begin{pmatrix}y_{0}\\y_{1}\\y_{2}\\y_{3}\end{pmatrix}}} Z = ∑ i = 0 3 z i ⋅ 2 8 i {\displaystyle Z=\sum _{i=0}^{3}z_{i}\cdot 2^{8i}} כאמור תיבות ההחלפה תלויות במפתח, ליתר דיוק בווקטור S {\displaystyle S} שהוכן מהמפתח לעיל. דהיינו ארבע התיבות מחושבות מתוך x i , y i {\displaystyle x_{i},y_{i}} עבור i = 0 , . . . , 3 {\displaystyle i=0,...,3} בפונקציה h {\displaystyle h} כאשר הרשימה L {\displaystyle L} היא S {\displaystyle S} בתמציתיות: g ( X ) = h ( X , S ) {\displaystyle g(X)=h(X,S)} . עם השלמת הגדרת הפונקציה h {\displaystyle h} , אפשר להשלים את תהליך הרחבת המפתח כדלהלן: ρ = 2 24 + 2 16 + 2 8 + 1 {\displaystyle \rho =2^{24}+2^{16}+2^{8}+1} A i = h ( 2 i ρ , M e ) {\displaystyle A_{i}=h(2i\rho ,M_{e})} B i = ROL ( h ( ( 2 i + 1 ) ρ , M o ) , 8 ) {\displaystyle B_{i}={\mbox{ROL}}(h((2i+1)\rho ,M_{o}),8)} K 2 i = ( A i + B i )  mod  2 32 {\displaystyle K_{2i}=(A_{i}+B_{i}){\mbox{ mod }}2^{32}} K 2 i + 1 = ROL ( ( A i + 2 B i )  mod  2 32 , 9 ) {\displaystyle K_{2i+1}={\mbox{ROL}}((A_{i}+2B_{i}){\mbox{ mod }}2^{32},9)} הקבוע ρ {\displaystyle \rho } משמש כאן להכפלת בתים, התכונה שלו היא שהמילה i ρ {\displaystyle i\rho } מורכבת מארבעה בתים שווים בערכם ל- i {\displaystyle i} . שימו לב שנעשה שימוש בהזזה מעגלית, בווקטורים M e , M o {\displaystyle M_{e},M_{o}} ובפונקציית PHT. התמורות הקבועות q 0 {\displaystyle q_{0}} ו- q 1 {\displaystyle q_{1}} מחושבות מתוך 4 תמורות שונות, עבור קלט x {\displaystyle x} מחשבים את הפלט y {\displaystyle y} על ידי: a 0 , b 0 = ⌊ x / 16 ⌋ , x  mod  16 {\displaystyle a_{0},b_{0}=\left\lfloor x/16\right\rfloor ,x{\mbox{ mod }}16} a 1 = a 0 ⊕ b 0 {\displaystyle a_{1}=a_{0}\oplus b_{0}} b 1 = a 0 ⊕ ROR 4 ( b 0 , 1 ) ⊕ 8 a 0  mod  16 {\displaystyle b_{1}=a_{0}\oplus {\mbox{ROR}}_{4}(b_{0},1)\oplus 8a_{0}{\mbox{ mod }}16} a 2 , b 2 = t 0 [ a 1 ] , t 1 [ b 1 ] {\displaystyle a_{2},b_{2}=t_{0}[a_{1}],t_{1}[b_{1}]} a 3 = a 2 ⊕ b 2 {\displaystyle a_{3}=a_{2}\oplus b_{2}} b 3 = a 2 ⊕ ROR 4 ( b 2 , 1 ) ⊕ 8 a 2  mod  16 {\displaystyle b_{3}=a_{2}\oplus {\mbox{ROR}}_{4}(b_{2},1)\oplus 8a_{2}{\mbox{ mod }}16} a 4 , b 4 = t 2 [ a 3 ] , t 3 [ b 3 ] {\displaystyle a_{4},b_{4}=t_{2}[a_{3}],t_{3}[b_{3}]} y = 16 b 4 + a 4 {\displaystyle y=16b_{4}+a_{4}} תחילה מפצלים את הבית לשני ניבלים, שעוברים שלב ערבוב עם פונקציה חד-חד-ערכית ועל, כל ניבל עובר בנוסף בתיבת החלפה של 4 סיביות קבועה משלו ואז מאחדים את הניבלים. תיבות ההחלפה 4-סיביות הן: צופן Twofish משתמש בתהליך הכנה מעט שונה עבור פענוח. השוני הוא בסדר ההפוך של המפתח המורחב K {\displaystyle K} . בגרסאות הממוטבות ביישום האלגוריתם בפועל, נמנעים בדרך כלל מפעולות אריתמטיות מסובכות, אותם מחליפים בטבלאות החלפה פשוטות בהם הערכים קודדו מראש בהתאם. ובמקרה של יישום בחומרה מוגבלת כמו כרטיס חכם, מפרקים פעולות לפקודות בסיסיות ביותר כמו Shift ו-XOR, לעיתים על חשבון יעילות. AES Threefish ברוס שנייר, Twofish, Schneier on Security