API tools faq
paste
Advertisement
Guest User

Untitled

a guest
Dec 22nd, 2011
204
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C 25.10 KB | None | 0 0
raw download clone embed print report
  1. /*
  2.  *****************************************************************************
  3.  *
  4.  * File:    rigndael.c
  5.  *
  6.  * Purpose: rijndael - An implementation of the Rijndael cipher.
  7.  *
  8.  * Copyright (C) 2000, 2001 Rafael R. Sevilla <sevillar@team.ph.inter.net>
  9.  *
  10.  * Currently maintained by brian d foy, <bdfoy@cpan.org>
  11.  *
  12.  *  License (GNU Public License):
  13.  *
  14.  *  This program is free software; you can redistribute it and/or
  15.  *  modify it under the terms of the GNU General Public License
  16.  *  as published by the Free Software Foundation; either version 2
  17.  *  of the License, or (at your option) any later version.
  18.  *
  19.  *  This program is distributed in the hope that it will be useful,
  20.  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
  21.  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  22.  *  GNU General Public License for more details.
  23.  *
  24.  *  You should have received a copy of the GNU General Public License
  25.  *  along with this program; if not, write to the Free Software
  26.  *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307
  27.  *  USA
  28.  *
  29.  *****************************************************************************
  30. */
  31. #include "rijndael.h"
  32. #include <stdlib.h>
  33. #include <string.h>
  34. #include <stdio.h>
  35.  
  36. /* These tables combine both the S-boxes and the mixcolumn transformation, so
  37.    that we can perform a round's encryption or by means of four table lookups
  38.    and four XOR's per column of state.  They were generated by the
  39.    makertbls.pl script. */
  40. uint32_t dtbl[] = {
  41.     0xa56363c6, 0x847c7cf8, 0x997777ee, 0x8d7b7bf6,
  42.     0x0df2f2ff, 0xbd6b6bd6, 0xb16f6fde, 0x54c5c591,
  43.     0x50303060, 0x03010102, 0xa96767ce, 0x7d2b2b56,
  44.     0x19fefee7, 0x62d7d7b5, 0xe6abab4d, 0x9a7676ec,
  45.     0x45caca8f, 0x9d82821f, 0x40c9c989, 0x877d7dfa,
  46.     0x15fafaef, 0xeb5959b2, 0xc947478e, 0x0bf0f0fb,
  47.     0xecadad41, 0x67d4d4b3, 0xfda2a25f, 0xeaafaf45,
  48.     0xbf9c9c23, 0xf7a4a453, 0x967272e4, 0x5bc0c09b,
  49.     0xc2b7b775, 0x1cfdfde1, 0xae93933d, 0x6a26264c,
  50.     0x5a36366c, 0x413f3f7e, 0x02f7f7f5, 0x4fcccc83,
  51.     0x5c343468, 0xf4a5a551, 0x34e5e5d1, 0x08f1f1f9,
  52.     0x937171e2, 0x73d8d8ab, 0x53313162, 0x3f15152a,
  53.     0x0c040408, 0x52c7c795, 0x65232346, 0x5ec3c39d,
  54.     0x28181830, 0xa1969637, 0x0f05050a, 0xb59a9a2f,
  55.     0x0907070e, 0x36121224, 0x9b80801b, 0x3de2e2df,
  56.     0x26ebebcd, 0x6927274e, 0xcdb2b27f, 0x9f7575ea,
  57.     0x1b090912, 0x9e83831d, 0x742c2c58, 0x2e1a1a34,
  58.     0x2d1b1b36, 0xb26e6edc, 0xee5a5ab4, 0xfba0a05b,
  59.     0xf65252a4, 0x4d3b3b76, 0x61d6d6b7, 0xceb3b37d,
  60.     0x7b292952, 0x3ee3e3dd, 0x712f2f5e, 0x97848413,
  61.     0xf55353a6, 0x68d1d1b9, 0x00000000, 0x2cededc1,
  62.     0x60202040, 0x1ffcfce3, 0xc8b1b179, 0xed5b5bb6,
  63.     0xbe6a6ad4, 0x46cbcb8d, 0xd9bebe67, 0x4b393972,
  64.     0xde4a4a94, 0xd44c4c98, 0xe85858b0, 0x4acfcf85,
  65.     0x6bd0d0bb, 0x2aefefc5, 0xe5aaaa4f, 0x16fbfbed,
  66.     0xc5434386, 0xd74d4d9a, 0x55333366, 0x94858511,
  67.     0xcf45458a, 0x10f9f9e9, 0x06020204, 0x817f7ffe,
  68.     0xf05050a0, 0x443c3c78, 0xba9f9f25, 0xe3a8a84b,
  69.     0xf35151a2, 0xfea3a35d, 0xc0404080, 0x8a8f8f05,
  70.     0xad92923f, 0xbc9d9d21, 0x48383870, 0x04f5f5f1,
  71.     0xdfbcbc63, 0xc1b6b677, 0x75dadaaf, 0x63212142,
  72.     0x30101020, 0x1affffe5, 0x0ef3f3fd, 0x6dd2d2bf,
  73.     0x4ccdcd81, 0x140c0c18, 0x35131326, 0x2fececc3,
  74.     0xe15f5fbe, 0xa2979735, 0xcc444488, 0x3917172e,
  75.     0x57c4c493, 0xf2a7a755, 0x827e7efc, 0x473d3d7a,
  76.     0xac6464c8, 0xe75d5dba, 0x2b191932, 0x957373e6,
  77.     0xa06060c0, 0x98818119, 0xd14f4f9e, 0x7fdcdca3,
  78.     0x66222244, 0x7e2a2a54, 0xab90903b, 0x8388880b,
  79.     0xca46468c, 0x29eeeec7, 0xd3b8b86b, 0x3c141428,
  80.     0x79dedea7, 0xe25e5ebc, 0x1d0b0b16, 0x76dbdbad,
  81.     0x3be0e0db, 0x56323264, 0x4e3a3a74, 0x1e0a0a14,
  82.     0xdb494992, 0x0a06060c, 0x6c242448, 0xe45c5cb8,
  83.     0x5dc2c29f, 0x6ed3d3bd, 0xefacac43, 0xa66262c4,
  84.     0xa8919139, 0xa4959531, 0x37e4e4d3, 0x8b7979f2,
  85.     0x32e7e7d5, 0x43c8c88b, 0x5937376e, 0xb76d6dda,
  86.     0x8c8d8d01, 0x64d5d5b1, 0xd24e4e9c, 0xe0a9a949,
  87.     0xb46c6cd8, 0xfa5656ac, 0x07f4f4f3, 0x25eaeacf,
  88.     0xaf6565ca, 0x8e7a7af4, 0xe9aeae47, 0x18080810,
  89.     0xd5baba6f, 0x887878f0, 0x6f25254a, 0x722e2e5c,
  90.     0x241c1c38, 0xf1a6a657, 0xc7b4b473, 0x51c6c697,
  91.     0x23e8e8cb, 0x7cdddda1, 0x9c7474e8, 0x211f1f3e,
  92.     0xdd4b4b96, 0xdcbdbd61, 0x868b8b0d, 0x858a8a0f,
  93.     0x907070e0, 0x423e3e7c, 0xc4b5b571, 0xaa6666cc,
  94.     0xd8484890, 0x05030306, 0x01f6f6f7, 0x120e0e1c,
  95.     0xa36161c2, 0x5f35356a, 0xf95757ae, 0xd0b9b969,
  96.     0x91868617, 0x58c1c199, 0x271d1d3a, 0xb99e9e27,
  97.     0x38e1e1d9, 0x13f8f8eb, 0xb398982b, 0x33111122,
  98.     0xbb6969d2, 0x70d9d9a9, 0x898e8e07, 0xa7949433,
  99.     0xb69b9b2d, 0x221e1e3c, 0x92878715, 0x20e9e9c9,
  100.     0x49cece87, 0xff5555aa, 0x78282850, 0x7adfdfa5,
  101.     0x8f8c8c03, 0xf8a1a159, 0x80898909, 0x170d0d1a,
  102.     0xdabfbf65, 0x31e6e6d7, 0xc6424284, 0xb86868d0,
  103.     0xc3414182, 0xb0999929, 0x772d2d5a, 0x110f0f1e,
  104.     0xcbb0b07b, 0xfc5454a8, 0xd6bbbb6d, 0x3a16162c,
  105. };
  106.  
  107. uint32_t itbl[] = {
  108.     0x50a7f451, 0x5365417e, 0xc3a4171a, 0x965e273a,
  109.     0xcb6bab3b, 0xf1459d1f, 0xab58faac, 0x9303e34b,
  110.     0x55fa3020, 0xf66d76ad, 0x9176cc88, 0x254c02f5,
  111.     0xfcd7e54f, 0xd7cb2ac5, 0x80443526, 0x8fa362b5,
  112.     0x495ab1de, 0x671bba25, 0x980eea45, 0xe1c0fe5d,
  113.     0x02752fc3, 0x12f04c81, 0xa397468d, 0xc6f9d36b,
  114.     0xe75f8f03, 0x959c9215, 0xeb7a6dbf, 0xda595295,
  115.     0x2d83bed4, 0xd3217458, 0x2969e049, 0x44c8c98e,
  116.     0x6a89c275, 0x78798ef4, 0x6b3e5899, 0xdd71b927,
  117.     0xb64fe1be, 0x17ad88f0, 0x66ac20c9, 0xb43ace7d,
  118.     0x184adf63, 0x82311ae5, 0x60335197, 0x457f5362,
  119.     0xe07764b1, 0x84ae6bbb, 0x1ca081fe, 0x942b08f9,
  120.     0x58684870, 0x19fd458f, 0x876cde94, 0xb7f87b52,
  121.     0x23d373ab, 0xe2024b72, 0x578f1fe3, 0x2aab5566,
  122.     0x0728ebb2, 0x03c2b52f, 0x9a7bc586, 0xa50837d3,
  123.     0xf2872830, 0xb2a5bf23, 0xba6a0302, 0x5c8216ed,
  124.     0x2b1ccf8a, 0x92b479a7, 0xf0f207f3, 0xa1e2694e,
  125.     0xcdf4da65, 0xd5be0506, 0x1f6234d1, 0x8afea6c4,
  126.     0x9d532e34, 0xa055f3a2, 0x32e18a05, 0x75ebf6a4,
  127.     0x39ec830b, 0xaaef6040, 0x069f715e, 0x51106ebd,
  128.     0xf98a213e, 0x3d06dd96, 0xae053edd, 0x46bde64d,
  129.     0xb58d5491, 0x055dc471, 0x6fd40604, 0xff155060,
  130.     0x24fb9819, 0x97e9bdd6, 0xcc434089, 0x779ed967,
  131.     0xbd42e8b0, 0x888b8907, 0x385b19e7, 0xdbeec879,
  132.     0x470a7ca1, 0xe90f427c, 0xc91e84f8, 0x00000000,
  133.     0x83868009, 0x48ed2b32, 0xac70111e, 0x4e725a6c,
  134.     0xfbff0efd, 0x5638850f, 0x1ed5ae3d, 0x27392d36,
  135.     0x64d90f0a, 0x21a65c68, 0xd1545b9b, 0x3a2e3624,
  136.     0xb1670a0c, 0x0fe75793, 0xd296eeb4, 0x9e919b1b,
  137.     0x4fc5c080, 0xa220dc61, 0x694b775a, 0x161a121c,
  138.     0x0aba93e2, 0xe52aa0c0, 0x43e0223c, 0x1d171b12,
  139.     0x0b0d090e, 0xadc78bf2, 0xb9a8b62d, 0xc8a91e14,
  140.     0x8519f157, 0x4c0775af, 0xbbdd99ee, 0xfd607fa3,
  141.     0x9f2601f7, 0xbcf5725c, 0xc53b6644, 0x347efb5b,
  142.     0x7629438b, 0xdcc623cb, 0x68fcedb6, 0x63f1e4b8,
  143.     0xcadc31d7, 0x10856342, 0x40229713, 0x2011c684,
  144.     0x7d244a85, 0xf83dbbd2, 0x1132f9ae, 0x6da129c7,
  145.     0x4b2f9e1d, 0xf330b2dc, 0xec52860d, 0xd0e3c177,
  146.     0x6c16b32b, 0x99b970a9, 0xfa489411, 0x2264e947,
  147.     0xc48cfca8, 0x1a3ff0a0, 0xd82c7d56, 0xef903322,
  148.     0xc74e4987, 0xc1d138d9, 0xfea2ca8c, 0x360bd498,
  149.     0xcf81f5a6, 0x28de7aa5, 0x268eb7da, 0xa4bfad3f,
  150.     0xe49d3a2c, 0x0d927850, 0x9bcc5f6a, 0x62467e54,
  151.     0xc2138df6, 0xe8b8d890, 0x5ef7392e, 0xf5afc382,
  152.     0xbe805d9f, 0x7c93d069, 0xa92dd56f, 0xb31225cf,
  153.     0x3b99acc8, 0xa77d1810, 0x6e639ce8, 0x7bbb3bdb,
  154.     0x097826cd, 0xf418596e, 0x01b79aec, 0xa89a4f83,
  155.     0x656e95e6, 0x7ee6ffaa, 0x08cfbc21, 0xe6e815ef,
  156.     0xd99be7ba, 0xce366f4a, 0xd4099fea, 0xd67cb029,
  157.     0xafb2a431, 0x31233f2a, 0x3094a5c6, 0xc066a235,
  158.     0x37bc4e74, 0xa6ca82fc, 0xb0d090e0, 0x15d8a733,
  159.     0x4a9804f1, 0xf7daec41, 0x0e50cd7f, 0x2ff69117,
  160.     0x8dd64d76, 0x4db0ef43, 0x544daacc, 0xdf0496e4,
  161.     0xe3b5d19e, 0x1b886a4c, 0xb81f2cc1, 0x7f516546,
  162.     0x04ea5e9d, 0x5d358c01, 0x737487fa, 0x2e410bfb,
  163.     0x5a1d67b3, 0x52d2db92, 0x335610e9, 0x1347d66d,
  164.     0x8c61d79a, 0x7a0ca137, 0x8e14f859, 0x893c13eb,
  165.     0xee27a9ce, 0x35c961b7, 0xede51ce1, 0x3cb1477a,
  166.     0x59dfd29c, 0x3f73f255, 0x79ce1418, 0xbf37c773,
  167.     0xeacdf753, 0x5baafd5f, 0x146f3ddf, 0x86db4478,
  168.     0x81f3afca, 0x3ec468b9, 0x2c342438, 0x5f40a3c2,
  169.     0x72c31d16, 0x0c25e2bc, 0x8b493c28, 0x41950dff,
  170.     0x7101a839, 0xdeb30c08, 0x9ce4b4d8, 0x90c15664,
  171.     0x6184cb7b, 0x70b632d5, 0x745c6c48, 0x4257b8d0,
  172. };
  173.  
  174.  
  175. /* Needed only for the key schedule and for final rounds */
  176. uint8_t sbox[256] = {
  177.     99, 124, 119, 123, 242, 107, 111, 197,  48,   1, 103,  43, 254, 215, 171,
  178.     118, 202, 130, 201, 125, 250,  89,  71, 240, 173, 212, 162, 175, 156, 164,
  179.     114, 192, 183, 253, 147,  38,  54,  63, 247, 204,  52, 165, 229, 241, 113,
  180.     216,  49,  21,  4, 199,  35, 195,  24, 150,   5, 154,   7,  18, 128, 226,
  181.     235,  39, 178, 117,  9, 131,  44,  26,  27, 110,  90, 160,  82,  59, 214,
  182.     179,  41, 227,  47, 132, 83, 209,   0, 237,  32, 252, 177,  91, 106, 203,
  183.     190,  57,  74,  76,  88, 207, 208, 239, 170, 251,  67,  77,  51, 133,  69,
  184.     249,   2, 127,  80,  60, 159, 168, 81, 163,  64, 143, 146, 157,  56, 245,
  185.     188, 182, 218,  33,  16, 255, 243, 210, 205,  12,  19, 236,  95, 151,  68,
  186.     23, 196, 167, 126,  61, 100,  93,  25, 115, 96, 129,  79, 220,  34,  42,
  187.     144, 136,  70, 238, 184,  20, 222,  94,  11, 219, 224,  50,  58,  10,  73,
  188.     6,  36,  92, 194, 211, 172,  98, 145, 149, 228, 121, 231, 200,  55, 109,
  189.     141, 213,  78, 169, 108,  86, 244, 234, 101, 122, 174,   8, 186, 120,  37,
  190.     46,  28, 166, 180, 198, 232, 221, 116,  31,  75, 189, 139, 138, 112,  62,
  191.     181, 102,  72,   3, 246,  14,  97,  53,  87, 185, 134, 193,  29, 158, 225,
  192.     248, 152,  17, 105, 217, 142, 148, 155,  30, 135, 233, 206,  85,  40, 223,
  193.     140, 161, 137,  13, 191, 230,  66, 104,  65, 153,  45,  15, 176,  84, 187,
  194.     22,
  195. };
  196.  
  197. uint8_t isbox[256] = {
  198.     82,   9, 106, 213,  48,  54, 165,  56, 191,  64, 163, 158, 129, 243, 215,
  199.     251, 124, 227,  57, 130, 155,  47, 255, 135,  52, 142,  67,  68, 196, 222,
  200.     233, 203, 84, 123, 148,  50, 166, 194,  35,  61, 238,  76, 149,  11,  66,
  201.     250, 195,  78,  8,  46, 161, 102,  40, 217,  36, 178, 118,  91, 162,  73,
  202.     109, 139, 209,  37, 114, 248, 246, 100, 134, 104, 152,  22, 212, 164,  92,
  203.     204,  93, 101, 182, 146, 108, 112,  72,  80, 253, 237, 185, 218,  94,  21,
  204.     70,  87, 167, 141, 157, 132, 144, 216, 171,   0, 140, 188, 211,  10, 247,
  205.     228,  88,   5, 184, 179,  69,   6, 208,  44,  30, 143, 202,  63,  15,   2,
  206.     193, 175, 189,   3,   1,  19, 138, 107,  58, 145,  17,  65,  79, 103, 220,
  207.     234, 151, 242, 207, 206, 240, 180, 230, 115, 150, 172, 116,  34, 231, 173,
  208.     53, 133, 226, 249,  55, 232,  28, 117, 223, 110,  71, 241,  26, 113,  29,
  209.     41, 197, 137, 111, 183,  98,  14, 170,  24, 190,  27, 252,  86,  62,  75,
  210.     198, 210, 121,  32, 154, 219, 192, 254, 120, 205,  90, 244,  31, 221, 168,
  211.     51, 136,   7, 199,  49, 177,  18,  16,  89,  39, 128, 236,  95,  96,  81,
  212.     127, 169,  25, 181,  74,  13,  45, 229, 122, 159, 147, 201, 156, 239, 160,
  213.     224,  59,  77, 174,  42, 245, 176, 200, 235, 187,  60, 131,  83, 153,  97,
  214.     23,  43,   4, 126, 186, 119, 214,  38, 225, 105,  20,  99,  85,  33,  12,
  215.     125,
  216. };
  217.  
  218. /* Used only by the key schedule */
  219. uint8_t Logtable[256] = {
  220.     0,   0,  25,   1,  50,   2,  26, 198,  75, 199,  27, 104,  51, 238, 223,  3,
  221.     100,   4, 224,  14,  52, 141, 129, 239,  76, 113,   8, 200, 248, 105,  28,
  222.     193, 125, 194,  29, 181, 249, 185,  39, 106,  77, 228, 166, 114, 154, 201,
  223.     9, 120, 101,  47, 138,   5,  33,  15, 225,  36,  18, 240, 130,  69,  53,
  224.     147, 218, 142, 150, 143, 219, 189,  54, 208, 206, 148,  19,  92, 210, 241,
  225.     64,  70, 131,  56, 102, 221, 253,  48, 191,   6, 139,  98, 179,  37, 226,
  226.     152,  34, 136, 145,  16, 126, 110,  72, 195, 163, 182,  30,  66,  58, 107,
  227.     40,  84, 250, 133,  61, 186, 43, 121,  10,  21, 155, 159,  94, 202,  78,
  228.     212, 172, 229, 243, 115, 167,  87, 175,  88, 168,  80, 244, 234, 214, 116,
  229.     79, 174, 233, 213, 231, 230, 173, 232, 44, 215, 117, 122, 235,  22,  11,
  230.     245,  89, 203,  95, 176, 156, 169,  81, 160, 127,  12, 246, 111,  23, 196,
  231.     73, 236, 216,  67,  31,  45, 164, 118, 123, 183, 204, 187,  62,  90, 251,
  232.     96, 177, 134,  59,  82, 161, 108, 170,  85,  41, 157, 151, 178, 135, 144,
  233.     97, 190, 220, 252, 188, 149, 207, 205,  55,  63,  91, 209, 83,  57, 132,
  234.     60,  65, 162, 109,  71,  20,  42, 158,  93,  86, 242, 211, 171, 68,  17,
  235.     146, 217,  35,  32,  46, 137, 180, 124, 184,  38, 119, 153, 227, 165, 103,
  236.     74, 237, 222, 197,  49, 254,  24,  13,  99, 140, 128, 192, 247, 112,   7,
  237. };
  238.  
  239. uint8_t Alogtable[256] = {
  240.     1,   3,   5,  15,  17,  51,  85, 255,  26,  46, 114, 150, 161, 248,  19,
  241.     53, 95, 225,  56,  72, 216, 115, 149, 164, 247,   2,   6,  10,  30,  34,
  242.     102, 170, 229,  52,  92, 228,  55,  89, 235,  38, 106, 190, 217, 112, 144,
  243.     171, 230,  49,  83, 245,   4,  12,  20,  60,  68, 204,  79, 209, 104, 184,
  244.     211, 110, 178, 205,  76, 212, 103, 169, 224,  59,  77, 215,  98, 166, 241,
  245.     8,  24,  40, 120, 136, 131, 158, 185, 208, 107, 189, 220, 127, 129, 152,
  246.     179, 206,  73, 219, 118, 154, 181, 196,  87, 249,  16,  48,  80, 240,  11,
  247.     29,  39, 105, 187, 214,  97, 163, 254,  25,  43, 125, 135, 146, 173, 236,
  248.     47, 113, 147, 174, 233,  32,  96, 160, 251,  22,  58,  78, 210, 109, 183,
  249.     194,  93, 231,  50,  86, 250,  21,  63,  65, 195,  94, 226,  61,  71, 201,
  250.     64, 192,  91, 237,  44, 116, 156, 191, 218, 117, 159, 186, 213, 100, 172,
  251.     239,  42, 126, 130, 157, 188, 223, 122, 142, 137, 128, 155, 182, 193,  88,
  252.     232,  35, 101, 175, 234,  37, 111, 177, 200,  67, 197,  84, 252,  31,  33,
  253.     99, 165, 244,   7,   9,  27,  45, 119, 153, 176, 203,  70, 202,  69, 207,
  254.     74, 222, 121, 139, 134, 145, 168, 227,  62,  66, 198,  81, 243,  14,  18,
  255.     54,  90, 238,  41, 123, 141, 140, 143, 138, 133, 148, 167, 242,  13,  23,
  256.     57,  75, 221, 124, 132, 151, 162, 253,  28,  36, 108, 180, 199,  82, 246, 1,
  257. };
  258.  
  259. #define ROTBYTE(x) (((x) >> 8) | (((x) & 0xff) << 24))
  260. #define ROTRBYTE(x) (((x) << 8) | (((x) >> 24) & 0xff))
  261. #define SUBBYTE(x, box) (((box)[((x) & 0xff)]) | \
  262.         ((box)[(((x) >> 8) & 0xff)] << 8) | \
  263.         ((box)[(((x) >> 16) & 0xff)] << 16) | \
  264.         ((box)[(((x) >> 24) & 0xff)] << 24))
  265.  
  266. static uint8_t
  267. xtime(uint8_t a)
  268. {
  269.     uint8_t b;
  270.  
  271.     b = (a & 0x80) ? 0x1b : 0;
  272.     a<<=1;
  273.     a^=b;
  274.     return(a);
  275. }
  276.  
  277. static uint8_t
  278. mul(uint8_t a, uint8_t b)
  279. {
  280.     if (a && b) return Alogtable[(Logtable[a] + Logtable[b])%255];
  281.     else return 0;
  282. }
  283.  
  284. static void
  285. inv_mix_column(uint32_t *a, uint32_t *b)
  286. {
  287.     uint8_t c[4][4];
  288.     int i, j;
  289.  
  290.     for(j = 0; j < 4; j++) {
  291.         for(i = 0; i < 4; i++) {
  292.             c[j][i] = mul(0xe, (a[j] >> i*8) & 0xff)
  293.                 ^ mul(0xb, (a[j] >> ((i+1)%4)*8) & 0xff)
  294.                 ^ mul(0xd, (a[j] >> ((i+2)%4)*8) & 0xff)
  295.                 ^ mul(0x9, (a[j] >> ((i+3)%4)*8) & 0xff);
  296.         }
  297.     }
  298.     for(i = 0; i < 4; i++) {
  299.         b[i] = 0;
  300.         for(j = 0; j < 4; j++)
  301.             b[i] |= c[i][j] << (j*8);
  302.     }
  303. }
  304.  
  305. void
  306. rijndael_setup(RIJNDAEL_context *ctx, size_t keysize, const uint8_t *key)
  307. {
  308.     int nk, nr, i, lastkey;
  309.     uint32_t temp, rcon;
  310.  
  311.     /* Truncate keysizes to the valid key sizes provided by Rijndael */
  312.     if (keysize >= 32) {
  313.         nk = 8;
  314.         nr = 14;
  315.     } else if (keysize >= 24) {
  316.         nk = 6;
  317.         nr = 12;
  318.     } else {            /* must be 16 or more */
  319.         nk = 4;
  320.         nr = 10;
  321.     }
  322.  
  323.     lastkey = (RIJNDAEL_BLOCKSIZE/4) * (nr + 1);
  324.     ctx->nrounds = nr;
  325.     rcon = 1;
  326.     for (i=0; i<nk; i++) {
  327.         ctx->keys[i] = key[i*4] + (key[i*4+1]<<8) + (key[i*4+2]<<16) +
  328.             (key[i*4+3]<<24);
  329.     }
  330.  
  331.     for (i=nk; i<lastkey; i++) {
  332.         temp = ctx->keys[i-1];
  333.         if (i % nk == 0) {
  334.             temp = SUBBYTE(ROTBYTE(temp), sbox) ^ rcon;
  335.             rcon = (uint32_t)xtime((uint8_t)rcon&0xff);
  336.         } else if (nk > 6 && (i%nk) == 4) {
  337.             temp = SUBBYTE(temp, sbox);
  338.         }
  339.         ctx->keys[i] = ctx->keys[i-nk] ^ temp;
  340.     }
  341.     /* Generate the inverse keys */
  342.     for (i=0; i<4; i++) {
  343.         ctx->ikeys[i] = ctx->keys[i];
  344.         ctx->ikeys[lastkey-4 + i] = ctx->keys[lastkey-4 + i];
  345.     }
  346.     for (i=4; i<lastkey-4; i+=4)
  347.         inv_mix_column(&(ctx->keys[i]), &(ctx->ikeys[i]));
  348. }
  349.  
  350. /* Key addition that also packs every byte in the key to a word rep. */
  351. static void
  352. key_addition_8to32(const uint8_t *txt, uint32_t *keys, uint32_t *out)
  353. {
  354.     const uint8_t *ptr;
  355.     int i, j;
  356.     uint32_t val;
  357.  
  358.     ptr = txt;
  359.     for (i=0; i<4; i++) {
  360.         val = 0;
  361.         for (j=0; j<4; j++)
  362.             val |= (*ptr++ << 8*j);
  363.         out[i] = keys[i]^val;
  364.     }
  365. }
  366.  
  367. static void
  368. key_addition32(const uint32_t *txt, uint32_t *keys, uint32_t *out)
  369. {
  370.     int i;
  371.  
  372.     for (i=0; i<4; i++)
  373.         out[i] = keys[i] ^ txt[i];
  374. }
  375.  
  376. static void
  377. key_addition32to8(const uint32_t *txt, uint32_t *keys, uint8_t *out)
  378. {
  379.     uint8_t *ptr;
  380.     int i, j;
  381.     uint32_t val;
  382.  
  383.     ptr = out;
  384.     for (i=0; i<4; i++) {
  385.         val = txt[i] ^ keys[i];
  386.         for (j=0; j<4; j++)
  387.             *ptr++ = (val >> 8*j) & 0xff;
  388.     }
  389. }
  390.  
  391. static int idx[4][4] = {
  392.     { 0, 1, 2, 3 },
  393.     { 1, 2, 3, 0 },
  394.     { 2, 3, 0, 1 },
  395.     { 3, 0, 1, 2 } };
  396.  
  397. void
  398. rijndael_encrypt(RIJNDAEL_context *ctx,
  399.         const uint8_t *plaintext,
  400.         uint8_t *ciphertext)
  401. {
  402.     int r, j;
  403.     uint32_t wtxt[4], t[4];     /* working ciphertext */
  404.     uint32_t e;
  405.  
  406.     key_addition_8to32(plaintext, &(ctx->keys[0]), wtxt);
  407.     for (r=1; r<ctx->nrounds; r++) {
  408.         for (j=0; j<4; j++) {
  409.             t[j] = dtbl[wtxt[j] & 0xff] ^
  410.                 ROTRBYTE(dtbl[(wtxt[idx[1][j]] >> 8) & 0xff]^
  411.                         ROTRBYTE(dtbl[(wtxt[idx[2][j]] >> 16) & 0xff] ^
  412.                             ROTRBYTE(dtbl[(wtxt[idx[3][j]] >> 24) & 0xff])));
  413.         }
  414.         key_addition32(t, &(ctx->keys[r*4]), wtxt);
  415.     }
  416.     /* last round is special: there is no mixcolumn, so we can't use the big
  417.        tables. */
  418.     for (j=0; j<4; j++) {
  419.         e = wtxt[j] & 0xff;
  420.         e |= (wtxt[idx[1][j]]) & (0xff  << 8 );
  421.         e |= (wtxt[idx[2][j]]) & (0xff  << 16);
  422.         e |= (wtxt[idx[3][j]]) & (0xffU << 24);
  423.         t[j] = e;
  424.     }
  425.     for (j=0; j<4; j++)
  426.         t[j] = SUBBYTE(t[j], sbox);
  427.     key_addition32to8(t, &(ctx->keys[4*ctx->nrounds]), ciphertext);
  428. }
  429.  
  430. static int iidx[4][4] = {
  431.     { 0, 1, 2, 3 },
  432.     { 3, 0, 1, 2 },
  433.     { 2, 3, 0, 1 },
  434.     { 1, 2, 3, 0 } };
  435.  
  436. void
  437. rijndael_decrypt(RIJNDAEL_context *ctx,
  438.         const uint8_t *ciphertext,
  439.         uint8_t *plaintext)
  440. {
  441.     int r, j;
  442.     uint32_t wtxt[4], t[4];     /* working ciphertext */
  443.     uint32_t e;
  444.  
  445.     key_addition_8to32(ciphertext, &(ctx->ikeys[4*ctx->nrounds]), wtxt);
  446.     for (r=ctx->nrounds-1; r> 0;  r--) {
  447.         for (j=0; j<4; j++) {
  448.             t[j] = itbl[wtxt[j] & 0xff] ^
  449.                 ROTRBYTE(itbl[(wtxt[iidx[1][j]] >> 8) & 0xff]^
  450.                         ROTRBYTE(itbl[(wtxt[iidx[2][j]] >> 16) & 0xff] ^
  451.                             ROTRBYTE(itbl[(wtxt[iidx[3][j]] >> 24) & 0xff])));
  452.         }
  453.         key_addition32(t, &(ctx->ikeys[r*4]), wtxt);
  454.     }
  455.     /* last round is special: there is no mixcolumn, so we can't use the big
  456.        tables. */
  457.     for (j=0; j<4; j++) {
  458.         e = wtxt[j] & 0xff;
  459.         e |= (wtxt[iidx[1][j]]) & (0xff  << 8);
  460.         e |= (wtxt[iidx[2][j]]) & (0xff  << 16);
  461.         e |= (wtxt[iidx[3][j]]) & (0xffU << 24);
  462.         t[j] = e;
  463.     }
  464.     for (j=0; j<4; j++)
  465.         t[j] = SUBBYTE(t[j], isbox);
  466.     key_addition32to8(t, &(ctx->ikeys[0]), plaintext);
  467. }
  468.  
  469. void
  470. block_encrypt(RIJNDAEL_context *ctx, uint8_t *input, int inputlen,
  471.         uint8_t *output, uint8_t *iv)
  472. {
  473.     int i, j, nblocks, carry_flg;
  474.     uint8_t block[RIJNDAEL_BLOCKSIZE], block2[RIJNDAEL_BLOCKSIZE];//, oldptxt;
  475.  
  476.     nblocks = inputlen / RIJNDAEL_BLOCKSIZE;
  477.  
  478.     switch (ctx->mode) {
  479.         case MODE_ECB:      /* electronic code book */
  480.             for (i = 0; i<nblocks; i++) {
  481.                 rijndael_encrypt(ctx, &input[RIJNDAEL_BLOCKSIZE*i],
  482.                         &output[RIJNDAEL_BLOCKSIZE*i]);
  483.             }
  484.             break;
  485.         case MODE_CBC:      /* Cipher block chaining */
  486.             /* set initial value */
  487.             memcpy(block, iv, RIJNDAEL_BLOCKSIZE);
  488.             for (i=0; i< nblocks; i++) {
  489.                 for (j=0; j<RIJNDAEL_BLOCKSIZE; j++)
  490.                     block[j] ^= input[i*RIJNDAEL_BLOCKSIZE + j] & 0xff;
  491.                 rijndael_encrypt(ctx, block, block);
  492.                 memcpy(&output[RIJNDAEL_BLOCKSIZE*i], block, RIJNDAEL_BLOCKSIZE);
  493.             }
  494.             break;
  495.         case MODE_CFB:      /* 128-bit cipher feedback */
  496.             memcpy(block, iv, RIJNDAEL_BLOCKSIZE);
  497.             for (i=0; i<nblocks; i++) {
  498.                 rijndael_encrypt(ctx, block, block);
  499.                 for (j=0; j<RIJNDAEL_BLOCKSIZE; j++)
  500.                     block[j] ^= input[i*RIJNDAEL_BLOCKSIZE + j];
  501.                 memcpy(&output[RIJNDAEL_BLOCKSIZE*i], block, RIJNDAEL_BLOCKSIZE);
  502.             }
  503.             break;
  504.         case MODE_OFB:      /* 128-bit output feedback */
  505.             memcpy(block, iv, RIJNDAEL_BLOCKSIZE);
  506.             for (i=0; i<nblocks; i++) {
  507.                 rijndael_encrypt(ctx, block, block);
  508.                 for (j=0; j<RIJNDAEL_BLOCKSIZE; j++) {
  509.                     output[RIJNDAEL_BLOCKSIZE*i + j] = block[j] ^
  510.                         input[RIJNDAEL_BLOCKSIZE*i + j];
  511.                 }
  512.             }
  513.             break;
  514.         case MODE_CTR:      /* counter */
  515.             memcpy(block, iv, RIJNDAEL_BLOCKSIZE);
  516.             for (i=0; i<nblocks; i++) {
  517.                 rijndael_encrypt(ctx, block, block2);
  518.                 for (j=0; j<RIJNDAEL_BLOCKSIZE; j++) {
  519.                     output[RIJNDAEL_BLOCKSIZE*i + j] = block2[j] ^
  520.                         input[RIJNDAEL_BLOCKSIZE*i + j];
  521.                 }
  522.                 block[RIJNDAEL_BLOCKSIZE-1]++;
  523.                 carry_flg = block[RIJNDAEL_BLOCKSIZE-1] != 0 ? 0 : 1;
  524.                 for (j=RIJNDAEL_BLOCKSIZE-2; j>=0; j--) {
  525.                     if (carry_flg) {
  526.                         block[j]++;
  527.                         carry_flg = block[j] != 0 ? 0 : 1;
  528.                     } else
  529.                         break;
  530.                 }
  531.             }
  532.             break;
  533.         default:
  534.             break;
  535.     }
  536. }
  537.  
  538. void
  539. block_decrypt(RIJNDAEL_context *ctx, uint8_t *input, int inputlen,
  540.         uint8_t *output, uint8_t *iv)
  541. {
  542.     int i, j, nblocks, carry_flg;
  543.     uint8_t block[RIJNDAEL_BLOCKSIZE], block2[RIJNDAEL_BLOCKSIZE];
  544.  
  545.     nblocks = inputlen / RIJNDAEL_BLOCKSIZE;
  546.     switch (ctx->mode) {
  547.         case MODE_ECB:
  548.             for (i = 0; i<nblocks; i++) {
  549.                 rijndael_decrypt(ctx, &input[RIJNDAEL_BLOCKSIZE*i],
  550.                         &output[RIJNDAEL_BLOCKSIZE*i]);
  551.             }
  552.             break;
  553.         case MODE_CBC:
  554.             /* first block */
  555.             rijndael_decrypt(ctx, input, block);
  556.             /* XOR the block with the IV to get the output */
  557.             for (i=0; i<RIJNDAEL_BLOCKSIZE; i++)
  558.                 output[i] = block[i] ^ iv[i];
  559.             for (i=1; i<nblocks; i++) {
  560.                 rijndael_decrypt(ctx, &input[i*RIJNDAEL_BLOCKSIZE], block);
  561.                 for (j=0; j<RIJNDAEL_BLOCKSIZE; j++) {
  562.                     output[i*RIJNDAEL_BLOCKSIZE + j] = block[j] ^
  563.                         input[(i-1)*RIJNDAEL_BLOCKSIZE + j];
  564.                 }
  565.             }
  566.             break;
  567.         case MODE_CFB:      /* 128-bit cipher feedback */
  568.             memcpy(block, iv, RIJNDAEL_BLOCKSIZE);
  569.             for (i=0; i<nblocks; i++) {
  570.                 rijndael_encrypt(ctx, block, block); /* ENCRYPT is right! */
  571.                 for (j=0; j<RIJNDAEL_BLOCKSIZE; j++) {
  572.                     output[RIJNDAEL_BLOCKSIZE*i + j] = block[j] ^
  573.                         input[RIJNDAEL_BLOCKSIZE*i + j];
  574.                 }
  575.                 memcpy(block, &input[RIJNDAEL_BLOCKSIZE*i], RIJNDAEL_BLOCKSIZE);
  576.             }
  577.             break;
  578.         case MODE_OFB:      /* 128-bit output feedback */
  579.             /* this is exactly the same as encryption in OFB...in fact you can use
  580.                the encryption in OFB mode to decrypt! */
  581.             memcpy(block, iv, RIJNDAEL_BLOCKSIZE);
  582.             for (i=0; i<nblocks; i++) {
  583.                 rijndael_encrypt(ctx, block, block);
  584.                 for (j=0; j<RIJNDAEL_BLOCKSIZE; j++) {
  585.                     output[RIJNDAEL_BLOCKSIZE*i + j] = block[j] ^
  586.                         input[RIJNDAEL_BLOCKSIZE*i + j];
  587.                 }
  588.             }
  589.             break;
  590.         case MODE_CTR:      /* counter */
  591.             memcpy(block, iv, RIJNDAEL_BLOCKSIZE);
  592.             for (i=0; i<nblocks; i++) {
  593.                 rijndael_encrypt(ctx, block, block2);
  594.                 for (j=0; j<RIJNDAEL_BLOCKSIZE; j++) {
  595.                     output[RIJNDAEL_BLOCKSIZE*i + j] = block2[j] ^
  596.                         input[RIJNDAEL_BLOCKSIZE*i + j];
  597.                 }
  598.                 block[RIJNDAEL_BLOCKSIZE-1]++;
  599.                 carry_flg = block[RIJNDAEL_BLOCKSIZE-1] != 0 ? 0 : 1;
  600.                 for (j=RIJNDAEL_BLOCKSIZE-2; j>=0; j--) {
  601.                     if (carry_flg) {
  602.                         block[j]++;
  603.                         carry_flg = block[j] != 0 ? 0 : 1;
  604.                     } else
  605.                         break;
  606.                 }
  607.             }
  608.             break;
  609.         default:
  610.             break;
  611.     }
  612. }
  613.  
  614. /***EOF***/
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement
Public Pastes
Advertisement
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.  OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!