大整数的乘法

#include<stdio.h>
#include<string.h>

/* 一个数组元素表示4 个十进制位，即数组是万进制的*/
#define BIGINT_RADIX 10000 //超过1万进1，类比10进制，一位就是一个数字，超过10进1
#define RADIX_LEN 4 //每个元素表示4位
#define MAX_LEN (200/RADIX_LEN+1) /* 整数的最大位数，之所以要加1，是因为当遇到不能整除的时候，除法是向下取整的，必须多要一位才能存储；当然，这里两个不超过200位的整数相乘，200刚好能被4整除，可以不需要加1*/
char a[MAX_LEN * RADIX_LEN], b[MAX_LEN * RADIX_LEN]; //采用字符串形式输入被乘数和乘数
int x[MAX_LEN], y[MAX_LEN], z[MAX_LEN * 2];//整数数组形式表示的被乘数，乘数，积

/**
* @brief 打印大整数.
* @param[in] x 大整数，用数组表示，低位在低地址
* @param[in] n 数组x 的长度
* @return 无
*/
void bigint_print(const int x[], const int n) {
	int i;
	int start_output = 0; /* 用于跳过前导0 */
	for (i = n - 1; i >= 0; --i) {
		if (start_output) { /* 如果多余的0 已经都跳过，则输出*/
			printf("%04d", x[i]);//结果的低位，一个整数元素表示结果的4位，不足补0
		} else if (x[i] > 0) {
			printf("%d", x[i]);//结果的最高1-4位，不需要输出前导0
			start_output = 1; /* 碰到第一个非0 的值，就说明多余的0 已经都跳过*/
		}
	}
	if(!start_output) printf("0"); /* 当x 全为0 时*/
}

/**
* @brief 将输入的字符串转化为大整数.
* @param[in] s 输入的字符串
* @param[out] x 大整数，用数组表示，低位在低地址
* @return 无
*/
void bigint_input(const char s[], int x[]) {
	int i, j = 0;
	const int len = strlen(s);
	for (i = 0; i < MAX_LEN; i++) x[i] = 0;
	for (i = len; i > 0; i -= RADIX_LEN) { 
		int temp = 0;
		int k;
		const int low = i-RADIX_LEN > 0 ? i-RADIX_LEN : 0;
		for (k = low; k < i; k++) {//从右边起（低位开始），每4位表示成1个数组元素
			temp = temp * 10 + s[k] - '0';
		}
		x[j++] = temp;
	}
}

/**
* @brief 大整数乘法.
* @param[in] x x
* @param[in] y y
* @param[out] z z=x*y
* @return 无
*/
void bigint_mul(const int x[], const int y[], int z[]) {
	int i, j;
	for (i = 0; i < MAX_LEN * 2; i++) z[i] = 0;
	for (i = 0; i < MAX_LEN; i++) {
		for (j = 0; j < MAX_LEN; j++) { /* 用y[i]，去乘以x 的各位*/
			z[i + j] += y[i] * x[j]; /* 两数第i, j 位相乘，累加到结果的第i+j 位*/
			if (z[i + j] >= BIGINT_RADIX) { /* 看是否要进位*/
				z[i + j + 1] += z[i + j] / BIGINT_RADIX; /* 进位*/
				z[i + j] %= BIGINT_RADIX;
			}
		}
	}
}

int main() {
	scanf("%s%s", a, b);
	bigint_input(a, x);
	bigint_input(b, y);
	bigint_mul(x, y, z);
	bigint_print(z, MAX_LEN * 2);
	printf("\n");
	return 0;
}

例子如下图：