|
3
|
1 #include <stdio.h> |
|
|
2 |
|
|
3 unsigned E[16] = { 1, 0xb, 9, 0xc, 0xd, 6, 0xf, 3, 0xe, 8, 7, 4, 0xa, 2, 5, 0 }; |
|
|
4 unsigned Ei[16]; |
|
|
5 unsigned R[16] = { 7, 0xc, 0xb, 0xd, 0xe, 4, 9, 0xf, 6, 3, 8, 0xa, 2, 5, 1, 0 }; |
|
|
6 unsigned cir[8][8] = { |
|
|
7 {1, 1, 4, 1, 8, 5, 2, 9 }, |
|
|
8 }; |
|
|
9 |
|
|
10 |
|
|
11 unsigned gf_mul(unsigned a, unsigned b) |
|
|
12 { |
|
|
13 unsigned r; |
|
|
14 |
|
|
15 r = 0; |
|
|
16 while (a) { |
|
|
17 if (a & 1) r ^= b; |
|
|
18 a >>= 1; |
|
|
19 b = (b << 1) ^ (b & 0x80 ? 0x11d : 0x00); |
|
|
20 } |
|
|
21 return r; |
|
|
22 } |
|
|
23 |
|
|
24 unsigned sbox(unsigned x) |
|
|
25 { |
|
|
26 unsigned a, b, w; |
|
|
27 |
|
|
28 a = x >> 4; |
|
|
29 b = x & 15; |
|
|
30 |
|
|
31 a = E[a]; b = Ei[b]; |
|
|
32 w = a ^ b; w = R[w]; |
|
|
33 a = E[a ^ w]; b = Ei[b ^ w]; |
|
|
34 |
|
|
35 |
|
|
36 return (a << 4) | b; |
|
|
37 } |
|
|
38 |
|
|
39 int main(void) |
|
|
40 { |
|
|
41 unsigned x, y; |
|
|
42 |
|
|
43 for (x = 0; x < 16; x++) Ei[E[x]] = x; |
|
|
44 |
|
|
45 // for (x = 0; x < 16; x++) printf("%2x ", sbox(x)); |
|
|
46 for (y = 1; y < 8; y++) { |
|
|
47 for (x = 0; x < 8; x++) { |
|
|
48 cir[y][x] = cir[y-1][(x-1)&7]; |
|
|
49 } |
|
|
50 } |
|
|
51 |
|
|
52 /* |
|
|
53 printf("\n"); |
|
|
54 for (y = 0; y < 8; y++) { |
|
|
55 for (x = 0; x < 8; x++) printf("%2d ", cir[y][x]); |
|
|
56 printf("\n"); |
|
|
57 } |
|
|
58 */ |
|
|
59 |
|
|
60 for (y = 0; y < 8; y++) { |
|
|
61 printf("static const ulong64 sbox%d[] = {\n", y); |
|
|
62 for (x = 0; x < 256; ) { |
|
|
63 printf("CONST64(0x%02x%02x%02x%02x%02x%02x%02x%02x)", |
|
|
64 gf_mul(sbox(x), cir[y][0]), |
|
|
65 gf_mul(sbox(x), cir[y][1]), |
|
|
66 gf_mul(sbox(x), cir[y][2]), |
|
|
67 gf_mul(sbox(x), cir[y][3]), |
|
|
68 gf_mul(sbox(x), cir[y][4]), |
|
|
69 gf_mul(sbox(x), cir[y][5]), |
|
|
70 gf_mul(sbox(x), cir[y][6]), |
|
|
71 gf_mul(sbox(x), cir[y][7])); |
|
|
72 if (x < 255) printf(", "); |
|
|
73 if (!(++x & 3)) printf("\n"); |
|
|
74 } |
|
|
75 printf("};\n\n"); |
|
|
76 } |
|
|
77 |
|
|
78 printf("static const ulong64 cont[] = {\n"); |
|
|
79 for (y = 0; y <= 10; y++) { |
|
|
80 printf("CONST64(0x"); |
|
|
81 for (x = 0; x < 8; x++) { |
|
|
82 printf("%02x", sbox((8*y + x)&255)); |
|
|
83 } |
|
|
84 printf("),\n"); |
|
|
85 } |
|
|
86 printf("};\n\n"); |
|
|
87 return 0; |
|
|
88 |
|
|
89 } |
|
|
90 |
|
|
91 |
