大数的阶乘

#include<stdio.h>
#include<string.h>
/* 一个数组元素表示4 个十进制位,即数组是万进制的*/
#define BIGINT_RADIX 10000
#define RADIX_LEN 4
/* 10000! 有35660 位*/
#define MAX_LEN (35660/RADIX_LEN+1) /* 整数的最大位数*/
int x[MAX_LEN + 1];
/**
* @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]);
		} else if (x[i] > 0) {
			printf("%d", x[i]);
			start_output = 1; /* 碰到第一个非0 的值,就说明多余的0 已经都跳过*/
		}
	}
	if(!start_output) printf("0"); /* 当x 全为0 时*/
}
/**
* @brief 大整数乘法, x = x*y.
* @param[inout] x x
* @param[in] y y
* @return 无
*/
void bigint_mul(int x[], const int y) {
	int i;
	int c = 0; /* 保存进位*/
	for (i = 0; i < MAX_LEN; i++) { /* 用y,去乘以x 的各位*/
		const int tmp = x[i] * y + c;
		x[i] = tmp % BIGINT_RADIX;
		c = tmp / BIGINT_RADIX;
	}
}
/**
* @brief 计算n 的阶乘
* @param[in] n
* @param[out] x 存放结果
* @return 无
*/
void bigint_factorial(int n, int x[]) {
	int i;
	memset(x, 0, sizeof(int) * (MAX_LEN + 1));
	x[0] = 1;
	for (i = 2; i <= n; i++) {
		bigint_mul(x, i);
	}
}
int main() {
	int n;
	while (scanf("%d", &n) != EOF) {//计算10000以内数字的阶乘
		bigint_factorial(n, x);
		bigint_print(x, MAX_LEN + 1);
		printf("\n");
	}
	return 0;
}

©️2020 CSDN 皮肤主题: 大白 设计师:CSDN官方博客 返回首页