生成AES S盒的实现和解析

ome…

于 2024-06-06 08:11:59 发布

阅读量908

点赞数 5

文章标签：算法密码学

本文链接：https://blog.csdn.net/2401_83140188/article/details/139486279

版权

在高级加密标准（AES）中，S盒（Substitution box）是一个非线性替换函数，用于提高加密的复杂性和安全性。本文将通过一个C程序，详细介绍如何生成AES S盒，并对代码进行解析，以帮助读者理解其实现原理。

代码实现

以下是生成AES S盒的C程序：

#include <stdio.h>

// 0x11B为AES中求逆要模的不可约多项式 x^8 + x^4 + x^3 + x + 1
// 0x63为AES仿射变换要加的列向量 
void gen_tabs(int sbx_tab[]) {
    int pow_tab[256]; // 存储2^i在GF(2^8)下的值
    int log_tab[256]; // 存储GF(2^8)下的2^i对应的系数 
    int mid_tab[256]; // 存储GF(2^8)上256个元素的乘法逆元 
    int b[8] = {0xf1, 0xe3, 0xc7, 0x8f, 0x1f, 0x3e, 0x7c, 0xf8}; 
    // 存放仿射变换的矩阵 
    int i, j, k, p;
    for (i = 0, p = 1; i < 256; i++) {
        pow_tab[i] = p; 
        log_tab[p] = i; 
        // 计算下一个2的i次方，如果p的最高位为1，就需要异或一个固定值0x11b
        p = p ^ (p << 1) ^ (p & 0x80 ? 0x11b : 0);
    }
    for (i = 0; i < 256; i++) {
        // 查表法求逆元 
        mid_tab[i] = (i ? pow_tab[255 - log_tab[i]] : 0);
    }
    // 生成S盒
    for (i = 0; i < 256; i++) {
        int t = 0, m = 0, mid = 0, tab = 0;
        // 对每个输入值i的乘法逆元，与8个字节b[j]进行运算
        for (j = 0; j < 8; j++) {
            m = mid = (b[j] & mid_tab[i]); 
            // 有限域上的乘法
            for (k = 0; k < 8; k++) {
                int n = (mid >> 1); // 将mid右移1位
                // 如果m不等于(n<<1)，说明m的最高位为1，需要异或一个固定值0x1b
                if (m != (n << 1))
                    t++; // 统计位移操作次数
                mid = n; 
                m = mid; 
            }
            // t存放仿射变换矩阵的行与逆元相与后结果中1的个数 
            if (t % 2 > 0) { 
                int temp = 1;
                // 左移j位，相当于对应的系数为1，其余系数为0
                for (k = 0; k < j; k++)
                    temp <<= 1;
                tab += temp; 
            }
            t = 0; 
        }
        sbx_tab[i] = tab ^ 0x63; // 仿射变换 
    }
}

int main() {
    int sbx_tab[256]; 
    gen_tabs(sbx_tab); // 生成S盒
    printf("\n"); 
    // 输出S盒内容
    for (int i = 0; i < 256; i++) {
        printf("%02x ", sbx_tab[i]);
        if ((i+1) % 16 == 0) printf("\n"); // 每输出16个字节换行
    }
    return 0;
}

代码解析

1. 初始化

首先，代码定义了一些变量和常量：

pow_tab[256] 和 log_tab[256]：用于存储有限域GF(2^8)中的指数和对数表。
mid_tab[256]：用于存储GF(2^8)上的乘法逆元。
b[8]：存储仿射变换的矩阵。

2. 生成指数表和对数表

代码通过以下循环生成2的指数和对数表：

for (i = 0, p = 1; i < 256; i++) {
    pow_tab[i] = p; 
    log_tab[p] = i; 
    // 计算下一个2的i次方，如果p的最高位为1，就需要异或一个固定值0x11b
    p = p ^ (p << 1) ^ (p & 0x80 ? 0x11b : 0);
}

3. 计算乘法逆元

代码通过查表法计算GF(2^8)上的每个元素的乘法逆元：

for (i = 0; i < 256; i++) {
    mid_tab[i] = (i ? pow_tab[255 - log_tab[i]] : 0);
}

4. 生成S盒

最后，通过以下代码生成S盒：

for (i = 0; i < 256; i++) {
    int t = 0, m = 0, mid = 0, tab = 0;
    // 对每个输入值i的乘法逆元，与8个字节b[j]进行运算
    for (j = 0; j < 8; j++) {
        m = mid = (b[j] & mid_tab[i]); 
        // 有限域上的乘法
        for (k = 0; k < 8; k++) {
            int n = (mid >> 1); // 将mid右移1位
            // 如果m不等于(n<<1)，说明m的最高位为1，需要异或一个固定值0x1b
            if (m != (n << 1))
                t++; // 统计位移操作次数
            mid = n; 
            m = mid; 
        }
        // t存放仿射变换矩阵的行与逆元相与后结果中1的个数 
        if (t % 2 > 0) { 
            int temp = 1;
            // 左移j位，相当于对应的系数为1，其余系数为0
            for (k = 0; k < j; k++)
                temp <<= 1;
            tab += temp; 
        }
        t = 0; 
    }
    sbx_tab[i] = tab ^ 0x63; // 仿射变换 
}

5. 输出S盒内容

在main函数中，调用gen_tabs函数生成S盒，并将其输出：

int main() {
    int sbx_tab[256]; 
    gen_tabs(sbx_tab); // 生成S盒
    printf("\n"); 
    // 输出S盒内容
    for (int i = 0; i < 256; i++) {
        printf("%02x ", sbx_tab[i]);
        if ((i+1) % 16 == 0) printf("\n"); // 每输出16个字节换行
    }
    return 0;
}