从零构建SM4国密算法：仅用stdio.h的C语言实战指南
密码学作为数字世界的基石，其核心算法实现往往被神秘化。本文将打破这种认知壁垒，带你用最基础的C语言工具——仅需stdio.h头文件，从零开始构建中国商用密码标准SM4算法。不同于依赖现成库的快速实现，我们将深入每个字节的处理细节，理解密码算法背后的设计哲学。
1. 环境准备与基础概念
在开始编码之前，我们需要明确几个关键概念。SM4是一种分组密码算法，采用128位分组长度和128位密钥长度。整个算法包含32轮非线性迭代结构，通过反复的混淆和扩散操作实现数据加密。
所需开发环境：
- 任意C语言编译器（GCC、Clang、MSVC等）
- 文本编辑器或IDE
- 终端或命令行界面
提示：虽然示例代码在Windows下的Dev-CPP测试通过，但标准C代码具有跨平台特性，可在Linux/macOS等系统直接编译
基础数据类型定义是整个项目的基石，我们将使用以下类型别名简化代码：
#include <stdio.h> #define u8 unsigned char // 8位无符号整数 #define u32 unsigned long // 32位无符号整数2. SM4核心组件实现
2.1 S盒变换实现
S盒是SM4算法的非线性核心，实现字节级的替换操作。其本质是一个256字节的查找表，每个输入字节被映射为另一个字节。这种非线性特性为算法提供了混淆能力。
const u8 Sbox[256] = { 0xd6,0x90,0xe9,0xfe,0xcc,0xe1,0x3d,0xb7,0x16,0xb6,0x14,0xc2,0x28,0xfb,0x2c,0x05, // ... 完整S盒数据见文末完整代码 };实现S盒查询函数时，需要注意大端序处理：
u32 functionB(u32 b) { u8 a[4]; a[0] = b >> 24; // 最高位字节 a[1] = b >> 16; a[2] = b >> 8; a[3] = b; // 最低位字节 return (Sbox[a[0]] << 24) | (Sbox[a[1]] << 16) | (Sbox[a[2]] << 8) | Sbox[a[3]]; }2.2 循环移位与线性变换
循环移位操作是密码算法中实现扩散效果的关键技术。SM4采用32位字的循环左移：
u32 loopLeft(u32 a, short length) { return (a << length) | (a >> (32 - length)); }基于循环移位，我们可以实现算法中两个关键的线性变换L和L'：
| 变换类型 | 公式 | 作用 |
|---|---|---|
| L | B⊕(B<<<2)⊕(B<<<10)⊕(B<<<18)⊕(B<<<24) | 加密流程中的线性扩散 |
| L' | B⊕(B<<<13)⊕(B<<<23) | 密钥扩展中的线性变换 |
对应代码实现：
u32 functionL1(u32 a) { return a ^ loopLeft(a, 2) ^ loopLeft(a, 10) ^ loopLeft(a, 18) ^ loopLeft(a, 24); } u32 functionL2(u32 a) { return a ^ loopLeft(a, 13) ^ loopLeft(a, 23); }3. 密钥扩展算法剖析
SM4的密钥扩展过程将128位初始密钥扩展为32个轮密钥（每个32位）。这个过程同样采用32轮迭代结构，但使用不同的线性变换。
3.1 初始化阶段
首先定义算法所需的常量和固定参数：
const u32 FK[4] = { 0xa3b1bac6, 0x56aa3350, 0x677d9197, 0xb27022dc }; const u32 CK[32] = { 0x00070e15, 0x1c232a31, 0x383f464d, 0x545b6269, // ... 完整CK数组见文末 };密钥扩展第一步将初始密钥与FK进行异或：
void extendFirst(u32 MK[], u32 K[]) { for(int i = 0; i < 4; i++) { K[i] = MK[i] ^ FK[i]; } }3.2 轮密钥生成
密钥扩展的核心迭代过程生成32个轮密钥：
void extendSecond(u32 RK[], u32 K[]) { for(short i = 0; i < 32; i++) { u32 temp = K[(i+1)%4] ^ K[(i+2)%4] ^ K[(i+3)%4] ^ CK[i]; K[(i+4)%4] = K[i%4] ^ functionT(temp, 2); // 使用L'变换 RK[i] = K[(i+4)%4]; } }合成变换T根据模式选择不同的线性变换：
u32 functionT(u32 a, short mode) { return mode == 1 ? functionL1(functionB(a)) // 加密使用L : functionL2(functionB(a)); // 密钥扩展使用L' }4. 加密与解密流程实现
4.1 加密算法核心
SM4加密采用32轮Feistel结构，每轮使用一个轮密钥：
void iterate32(u32 X[], u32 RK[]) { for(short i = 0; i < 32; i++) { u32 temp = X[(i+1)%4] ^ X[(i+2)%4] ^ X[(i+3)%4] ^ RK[i]; X[(i+4)%4] = X[i%4] ^ functionT(temp, 1); // 使用L变换 } }加密完成后需要进行反序操作：
void reverse(u32 X[], u32 Y[]) { for(short i = 0; i < 4; i++) { Y[i] = X[3 - i]; } } void encryptSM4(u32 X[], u32 RK[], u32 Y[]) { iterate32(X, RK); reverse(X, Y); }4.2 解密算法实现
SM4算法的对合特性使得解密过程与加密几乎相同，只需逆序使用轮密钥：
void decryptSM4(u32 X[], u32 RK[], u32 Y[]) { u32 reverseRK[32]; for(short i = 0; i < 32; i++) { reverseRK[i] = RK[31-i]; // 密钥逆序 } iterate32(X, reverseRK); reverse(X, Y); }5. 完整测试案例
为验证算法正确性，我们使用标准测试向量：
int main(void) { u32 X[4] = {0x01234567, 0x89abcdef, 0xfedcba98, 0x76543210}; // 明文 u32 MK[4] = {0x01234567, 0x89abcdef, 0xfedcba98, 0x76543210}; // 密钥 u32 RK[32], K[4], Y[4]; printf("************** 密钥扩展过程 **************\n"); getRK(MK, K, RK); printf("************** 加密测试 **************\n"); encryptSM4(X, RK, Y); printf("密文：%08x %08x %08x %08x\n", Y[0], Y[1], Y[2], Y[3]); printf("************** 解密测试 **************\n"); decryptSM4(Y, RK, X); printf("明文：%08x %08x %08x %08x\n", X[0], X[1], X[2], X[3]); return 0; }预期输出结果应满足：
- 加密密文：681edf34 d206965e 86b3e94f 536e4246
- 解密结果应与原始明文一致
6. 性能优化与内存安全
虽然我们的实现注重教学清晰度，但在实际应用中还需考虑：
性能优化方向：
- 将S盒查询展开为32位预计算表
- 使用宏代替函数调用减少开销
- 循环展开技术加速迭代过程
内存安全注意事项：
- 敏感数据(如密钥)使用后应立即清零
- 避免密钥数据被交换到磁盘
- 考虑加入抗侧信道攻击措施
void secureZero(u32 *arr, size_t len) { volatile u32 *p = arr; while(len--) *p++ = 0; }通过这150行左右的精炼代码，我们完整实现了SM4国密算法的所有核心功能。这种从零开始的实现方式不仅帮助理解算法本质，也为嵌入式等受限环境提供了可移植的解决方案。
