DCU_APP/GD32E508_DCU_APP/sys/aes128.cpp

/**
 * @file aes128.cpp
 * @brief AES-128 ECB 加解密实现（移植自电机 BOOT 工程 sys/aes128.cpp）。
 */
#include "sys/aes128.h"

/* S盒 */
static const uint8_t sbox[256] = {
    0x63,0x7c,0x77,0x7b,0xf2,0x6b,0x6f,0xc5,0x30,0x01,0x67,0x2b,0xfe,0xd7,0xab,0x76,
    0xca,0x82,0xc9,0x7d,0xfa,0x59,0x47,0xf0,0xad,0xd4,0xa2,0xaf,0x9c,0xa4,0x72,0xc0,
    0xb7,0xfd,0x93,0x26,0x36,0x3f,0xf7,0xcc,0x34,0xa5,0xe5,0xf1,0x71,0xd8,0x31,0x15,
    0x04,0xc7,0x23,0xc3,0x18,0x96,0x05,0x9a,0x07,0x12,0x80,0xe2,0xeb,0x27,0xb2,0x75,
    0x09,0x83,0x2c,0x1a,0x1b,0x6e,0x5a,0xa0,0x52,0x3b,0xd6,0xb3,0x29,0xe3,0x2f,0x84,
    0x53,0xd1,0x00,0xed,0x20,0xfc,0xb1,0x5b,0x6a,0xcb,0xbe,0x39,0x4a,0x4c,0x58,0xcf,
    0xd0,0xef,0xaa,0xfb,0x43,0x4d,0x33,0x85,0x45,0xf9,0x02,0x7f,0x50,0x3c,0x9f,0xa8,
    0x51,0xa3,0x40,0x8f,0x92,0x9d,0x38,0xf5,0xbc,0xb6,0xda,0x21,0x10,0xff,0xf3,0xd2,
    0xcd,0x0c,0x13,0xec,0x5f,0x97,0x44,0x17,0xc4,0xa7,0x7e,0x3d,0x64,0x5d,0x19,0x73,
    0x60,0x81,0x4f,0xdc,0x22,0x2a,0x90,0x88,0x46,0xee,0xb8,0x14,0xde,0x5e,0x0b,0xdb,
    0xe0,0x32,0x3a,0x0a,0x49,0x06,0x24,0x5c,0xc2,0xd3,0xac,0x62,0x91,0x95,0xe4,0x79,
    0xe7,0xc8,0x37,0x6d,0x8d,0xd5,0x4e,0xa9,0x6c,0x56,0xf4,0xea,0x65,0x7a,0xae,0x08,
    0xba,0x78,0x25,0x2e,0x1c,0xa6,0xb4,0xc6,0xe8,0xdd,0x74,0x1f,0x4b,0xbd,0x8b,0x8a,
    0x70,0x3e,0xb5,0x66,0x48,0x03,0xf6,0x0e,0x61,0x35,0x57,0xb9,0x86,0xc1,0x1d,0x9e,
    0xe1,0xf8,0x98,0x11,0x69,0xd9,0x8e,0x94,0x9b,0x1e,0x87,0xe9,0xce,0x55,0x28,0xdf,
    0x8c,0xa1,0x89,0x0d,0xbf,0xe6,0x42,0x68,0x41,0x99,0x2d,0x0f,0xb0,0x54,0xbb,0x16
};

/* 逆S盒 */
static const uint8_t inv_sbox[256] = {
    0x52,0x09,0x6a,0xd5,0x30,0x36,0xa5,0x38,0xbf,0x40,0xa3,0x9e,0x81,0xf3,0xd7,0xfb,
    0x7c,0xe3,0x39,0x82,0x9b,0x2f,0xff,0x87,0x34,0x8e,0x43,0x44,0xc4,0xde,0xe9,0xcb,
    0x54,0x7b,0x94,0x32,0xa6,0xc2,0x23,0x3d,0xee,0x4c,0x95,0x0b,0x42,0xfa,0xc3,0x4e,
    0x08,0x2e,0xa1,0x66,0x28,0xd9,0x24,0xb2,0x76,0x5b,0xa2,0x49,0x6d,0x8b,0xd1,0x25,
    0x72,0xf8,0xf6,0x64,0x86,0x68,0x98,0x16,0xd4,0xa4,0x5c,0xcc,0x5d,0x65,0xb6,0x92,
    0x6c,0x70,0x48,0x50,0xfd,0xed,0xb9,0xda,0x5e,0x15,0x46,0x57,0xa7,0x8d,0x9d,0x84,
    0x90,0xd8,0xab,0x00,0x8c,0xbc,0xd3,0x0a,0xf7,0xe4,0x58,0x05,0xb8,0xb3,0x45,0x06,
    0xd0,0x2c,0x1e,0x8f,0xca,0x3f,0x0f,0x02,0xc1,0xaf,0xbd,0x03,0x01,0x13,0x8a,0x6b,
    0x3a,0x91,0x11,0x41,0x4f,0x67,0xdc,0xea,0x97,0xf2,0xcf,0xce,0xf0,0xb4,0xe6,0x73,
    0x96,0xac,0x74,0x22,0xe7,0xad,0x35,0x85,0xe2,0xf9,0x37,0xe8,0x1c,0x75,0xdf,0x6e,
    0x47,0xf1,0x1a,0x71,0x1d,0x29,0xc5,0x89,0x6f,0xb7,0x62,0x0e,0xaa,0x18,0xbe,0x1b,
    0xfc,0x56,0x3e,0x4b,0xc6,0xd2,0x79,0x20,0x9a,0xdb,0xc0,0xfe,0x78,0xcd,0x5a,0xf4,
    0x1f,0xdd,0xa8,0x33,0x88,0x07,0xc7,0x31,0xb1,0x12,0x10,0x59,0x27,0x80,0xec,0x5f,
    0x60,0x51,0x7f,0xa9,0x19,0xb5,0x4a,0x0d,0x2d,0xe5,0x7a,0x9f,0x93,0xc9,0x9c,0xef,
    0xa0,0xe0,0x3b,0x4d,0xae,0x2a,0xf5,0xb0,0xc8,0xeb,0xbb,0x3c,0x83,0x53,0x99,0x61,
    0x17,0x2b,0x04,0x7e,0xba,0x77,0xd6,0x26,0xe1,0x69,0x14,0x63,0x55,0x21,0x0c,0x7d
};

/* 轮常量 */
static const uint8_t rcon[10] = {0x01,0x02,0x04,0x08,0x10,0x20,0x40,0x80,0x1b,0x36};

static uint8_t xtime(uint8_t x) { return static_cast<uint8_t>((x << 1) ^ (((x >> 7) & 1) * 0x1b)); }

static uint8_t gf_mul(uint8_t a, uint8_t b) {
    uint8_t p = 0;
    while (b) {
        if (b & 1) p ^= a;
        const uint8_t hi = a & 0x80;
        a <<= 1;
        if (hi) a ^= 0x1b;
        b >>= 1;
    }
    return p;
}

static void add_round_key(uint8_t *state, const uint8_t *rk) {
    for (uint8_t i = 0; i < 16; i++) state[i] ^= rk[i];
}

static void sub_bytes(uint8_t *state) {
    for (uint8_t i = 0; i < 16; i++) state[i] = sbox[state[i]];
}

static void inv_sub_bytes(uint8_t *state) {
    for (uint8_t i = 0; i < 16; i++) state[i] = inv_sbox[state[i]];
}

static void shift_rows(uint8_t *s) {
    uint8_t t;
    t=s[1]; s[1]=s[5]; s[5]=s[9]; s[9]=s[13]; s[13]=t;
    t=s[2]; s[2]=s[10]; s[10]=t; t=s[6]; s[6]=s[14]; s[14]=t;
    t=s[3]; s[3]=s[15]; s[15]=s[11]; s[11]=s[7]; s[7]=t;
}

static void inv_shift_rows(uint8_t *s) {
    uint8_t t;
    t=s[13]; s[13]=s[9]; s[9]=s[5]; s[5]=s[1]; s[1]=t;
    t=s[2]; s[2]=s[10]; s[10]=t; t=s[6]; s[6]=s[14]; s[14]=t;
    t=s[7]; s[7]=s[11]; s[11]=s[15]; s[15]=s[3]; s[3]=t;
}

static void mix_columns(uint8_t *s) {
    for (uint8_t c = 0; c < 4; c++) {
        const uint8_t i = c * 4;
        const uint8_t a0=s[i], a1=s[i+1], a2=s[i+2], a3=s[i+3];
        s[i]   = xtime(a0) ^ xtime(a1) ^ a1 ^ a2 ^ a3;
        s[i+1] = a0 ^ xtime(a1) ^ xtime(a2) ^ a2 ^ a3;
        s[i+2] = a0 ^ a1 ^ xtime(a2) ^ xtime(a3) ^ a3;
        s[i+3] = xtime(a0) ^ a0 ^ a1 ^ a2 ^ xtime(a3);
    }
}

static void inv_mix_columns(uint8_t *s) {
    for (uint8_t c = 0; c < 4; c++) {
        const uint8_t i = c * 4;
        const uint8_t a0=s[i], a1=s[i+1], a2=s[i+2], a3=s[i+3];
        s[i]   = gf_mul(0x0e,a0)^gf_mul(0x0b,a1)^gf_mul(0x0d,a2)^gf_mul(0x09,a3);
        s[i+1] = gf_mul(0x09,a0)^gf_mul(0x0e,a1)^gf_mul(0x0b,a2)^gf_mul(0x0d,a3);
        s[i+2] = gf_mul(0x0d,a0)^gf_mul(0x09,a1)^gf_mul(0x0e,a2)^gf_mul(0x0b,a3);
        s[i+3] = gf_mul(0x0b,a0)^gf_mul(0x0d,a1)^gf_mul(0x09,a2)^gf_mul(0x0e,a3);
    }
}

static void key_expansion(aes128_context *ctx) {
    memcpy(ctx->round_keys[0], ctx->key, 16);
    for (uint8_t i = 1; i <= 10; i++) {
        uint8_t temp[4];
        temp[0] = ctx->round_keys[i-1][12];
        temp[1] = ctx->round_keys[i-1][13];
        temp[2] = ctx->round_keys[i-1][14];
        temp[3] = ctx->round_keys[i-1][15];
        const uint8_t t = temp[0];
        temp[0]=temp[1]; temp[1]=temp[2]; temp[2]=temp[3]; temp[3]=t;
        temp[0]=sbox[temp[0]]; temp[1]=sbox[temp[1]];
        temp[2]=sbox[temp[2]]; temp[3]=sbox[temp[3]];
        temp[0] ^= rcon[i-1];
        ctx->round_keys[i][0] = ctx->round_keys[i-1][0] ^ temp[0];
        ctx->round_keys[i][1] = ctx->round_keys[i-1][1] ^ temp[1];
        ctx->round_keys[i][2] = ctx->round_keys[i-1][2] ^ temp[2];
        ctx->round_keys[i][3] = ctx->round_keys[i-1][3] ^ temp[3];
        for (uint8_t j = 1; j < 4; j++) {
            ctx->round_keys[i][j*4+0] = ctx->round_keys[i][(j-1)*4+0] ^ ctx->round_keys[i-1][j*4+0];
            ctx->round_keys[i][j*4+1] = ctx->round_keys[i][(j-1)*4+1] ^ ctx->round_keys[i-1][j*4+1];
            ctx->round_keys[i][j*4+2] = ctx->round_keys[i][(j-1)*4+2] ^ ctx->round_keys[i-1][j*4+2];
            ctx->round_keys[i][j*4+3] = ctx->round_keys[i][(j-1)*4+3] ^ ctx->round_keys[i-1][j*4+3];
        }
    }
}

void aes128_set_key(aes128_context *ctx, const uint8_t *key) {
    memcpy(ctx->key, key, AES128_KEY_LEN);
    key_expansion(ctx);
}

void aes128_encrypt_block(aes128_context *ctx, const uint8_t *input, uint8_t *output) {
    memcpy(ctx->state, input, AES_BLOCK_SIZE);
    add_round_key(ctx->state, ctx->round_keys[0]);
    for (uint8_t r = 1; r <= 9; r++) {
        sub_bytes(ctx->state);
        shift_rows(ctx->state);
        mix_columns(ctx->state);
        add_round_key(ctx->state, ctx->round_keys[r]);
    }
    sub_bytes(ctx->state);
    shift_rows(ctx->state);
    add_round_key(ctx->state, ctx->round_keys[10]);
    memcpy(output, ctx->state, AES_BLOCK_SIZE);
}

void aes128_decrypt_block(aes128_context *ctx, const uint8_t *input, uint8_t *output) {
    memcpy(ctx->state, input, AES_BLOCK_SIZE);
    add_round_key(ctx->state, ctx->round_keys[10]);
    for (uint8_t r = 9; r >= 1; r--) {
        inv_shift_rows(ctx->state);
        inv_sub_bytes(ctx->state);
        add_round_key(ctx->state, ctx->round_keys[r]);
        inv_mix_columns(ctx->state);
    }
    inv_shift_rows(ctx->state);
    inv_sub_bytes(ctx->state);
    add_round_key(ctx->state, ctx->round_keys[0]);
    memcpy(output, ctx->state, AES_BLOCK_SIZE);
}

void aes128_ecb_encrypt_48bytes(aes128_context *ctx, const uint8_t *input, uint8_t *output) {
    aes128_encrypt_block(ctx, input,      output);
    aes128_encrypt_block(ctx, input + 16, output + 16);
    aes128_encrypt_block(ctx, input + 32, output + 32);
}

void aes128_ecb_decrypt_48bytes(aes128_context *ctx, const uint8_t *input, uint8_t *output) {
    aes128_decrypt_block(ctx, input,      output);
    aes128_decrypt_block(ctx, input + 16, output + 16);
    aes128_decrypt_block(ctx, input + 32, output + 32);
}