加入收藏 | 设为首页 | 会员中心 | 我要投稿 李大同 (https://www.lidatong.com.cn/)- 科技、建站、经验、云计算、5G、大数据,站长网!
当前位置: 首页 > 大数据 > 正文

大数的阶乘

发布时间:2020-12-14 02:14:40 所属栏目:大数据 来源:网络整理
导读:#includestdio.h#includestring.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 打印
#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,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;
}

(编辑:李大同)

【声明】本站内容均来自网络,其相关言论仅代表作者个人观点,不代表本站立场。若无意侵犯到您的权利,请及时与联系站长删除相关内容!

    推荐文章
      热点阅读