#include <stdlib.h>
#include <stdio.h>
struct node;
typedef struct node* PtrToNode;
struct node {
long long n;
PtrToNode next;
};
PtrToNode factorial(int n, PtrToNode init);
void multiply(long long n, PtrToNode init, long long carry);
int main() {
int n;
while (1) {
scanf("%d", &n);
if (n > 0) {
break;
} else if (n == 0){
printf("1\n");
return 0;
} else {
printf("Retry.\n");
}
}
PtrToNode init = malloc(sizeof(struct node));
init->n = 1;
init->next = NULL;
PtrToNode head = factorial(n, init);
PtrToNode tail = reverse(head);
PtrToNode backup = tail;
printf("%lld", tail->n);
tail = tail->next;
while(tail != NULL) {
printf("%04lld", tail->n);
tail = tail->next;
}
PtrToNode t;
while (backup != NULL) {
t = backup;
backup = backup->next;
free(t);
}
}
PtrToNode factorial(int n, PtrToNode init) {
for (int i = 1; i <= n; i++)
multiply(i, init, 0);
return init;
}
void multiply(long long n, PtrToNode init, long long carry) {
long long temp = n * init->n + carry;
init->n = temp % 10000;
carry = temp / 10000;
if (init->next != NULL) {
multiply(n, init->next, carry);
} else if (carry > 0) {
init->next = malloc(sizeof(struct node));
init->next->n = carry;
init->next->next = NULL;
} else {
return;
}
}
这是我计算10000阶乘的函数,我已经计算出与在线数据相比的正确答案。但是,当我将n的范围限制为[0.999]时,我甚至无法计算2000因子。为什么当我改变n&#39;范围为[0,9999],那么它可以计算2000!甚至10000!?
答案 0 :(得分:10)
将数字群集限制为三位数的问题是,将三位数字乘以1,000以上的值可能会产生一个数字四位数。
虽然这不会立即产生问题,但是当您将值传递到计算链时,错误将会累积。 2000年的结果超过5000位!进位溢出肯定会在结果中产生可见的错误。这发生在计算2000年的第1258步!。
这不会发生四位数群集和10,000,因为进位可能“溢出”到下一个数字的唯一位置是最后一个群集的计算,long long
有足够的空间容纳这没有溢出。