C程序训练：大数相乘与阶乘的计算

两个大数相乘，我们可以利用小学生列竖式做乘法的方法编写程序即可。例如，计算123*23，可以按以下步骤做：

answer = 0；
temp=123*3 =369
answer = answer + temp
temp = 123 * 20 = 2460
answer = answer + temp

按该方法编写程序，需要设计：

计算两个数的加，完成上面的3和4步。
完成被乘数与一位数的乘法，完成上面的第2步和第4步。注意低位的0，比如第4步，123和乘数的十位2相乘，结果还需乘10。如果有百位，计算结果需要乘以100。在实现时，可以预先留出0的位数，计算的结果放在前面。

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#define DIGIT 500
void reserved(char * str, int len) // 逆置
{
for (int m = 0, n = len - 1; m < n; m++, n--) {
char temp = str[m];
str[m] = str[n];
str[n] = temp;
}
}
// 加法 a=a+b, 返回a的长度。
int add(char *a, char *b, int len_a, int len_b) {
int len, i;
len = (len_a>=len_b)?len_a:len_b;
for(i=0; i<len; i++){
a[i] = a[i] + b[i];
a[i+1] += a[i] / 10;
a[i] %= 10;
}
if(a[i]) i++;
return i;
}

int multiply(char *a, char *b, char *answer, int len1, int len2)
{
char temp[DIGIT]={0};//每一轮乘法的值
int anslen = 0;
int k;
for (int i = 0; i <len2; i++) {
memset(temp,0,DIGIT);
k = i;
for(int j=0; j<len1; j++,k++) {
temp[k] += a[j] * b[i];
temp[k+1] = temp[k] / 10;
temp[k] %= 10;
}
if(temp[k]) k++;
anslen = add(answer, temp, anslen, k);//加法运算
}
reserved(answer, anslen);
return anslen;
}

int main()
{
char num1[DIGIT]={0}, num2[DIGIT]={0}, answer[DIGIT]={0};
int i,len1,len2,len;
scanf("%s", num1);
scanf("%s", num2);
len1 = strlen(num1);
i=len1;
while(i) {
num1[--i] -= '0';
if(num1[i]<0 || num1[i]>9) {
printf("输入错。\n");
exit(1);
}
}
len2 = strlen(num2);
i=len2;
while(i) {
num2[--i] -= '0';
if(num2[i]<0 || num2[i]>9) {
printf("输入错。\n");
exit(2);
}
}
reserved(num1, len1);
reserved(num2, len2);
len = multiply(num1, num2, answer, len1, len2);
for(i=0; i<len; i++)
printf("%d", answer[i]);
printf("\n");
return 0;
}

对上面的算法进行优化，将add和multiply两个函数合并在一起，代码更紧凑。

优化后的multiply代码如下。

复制代码

int multiply(char *a, char *b, char *answer, int len1, int len2) {
  int anslen = 0;  
  for (int i = 0; i <len2; i++) {
    for (int j = 0; j <len1; j++) {
      answer[i+j] += a[j] * b[i];
      answer[i+j+1] += answer[i+j]/10;
      answer[i+j] %= 10;
    }    
  }
  for(anslen=len2+len1;!answer[anslen]; anslen--);
  reserved(answer, ++anslen);
  return anslen;
}

同样，按此思路，我们可以计算出任何数的阶乘。

程序代码如下：

复制代码

#include <stdio.h>
#define DIGIT 5000

int multiply(int *answer,int n) 
{
  int length = 1;
  int k,i,j,c=0;  
  for (i = 2; i <= n; i++) {
    for(k=0; k<length; k++) {
       answer[k] = answer[k] * i + c;
       c = answer[k] / 10;    
       answer[k] %= 10;        
    }
    while(c) {
       answer[length++]=c % 10;        
       c /= 10;
    }        
  }
  return length;
}

int main() {
  int answer[DIGIT]={1};
  int n,len;
  scanf("%d",&n);
  len = multiply(answer,n);
  for(int i=len-1; i>=0; i--) 
    printf("%d", answer[i]);
  printf("\n");
  return 0;
}

为了让程序看起来更加容易，对程序进行改进、优化，代码如下。

复制代码

#include <stdio.h>
#define DIGIT 5000
int multiply(int answer[],int len, int n) 
{  
  int i,c=0;  
  for (i = 0; i < len; i++) {
       answer[i] = answer[i] * n + c;
     c = answer[i] / 10;    
     answer[i] %= 10;        
    }
    while(c){
     answer[i++] = c%10;  
     c = c / 10;         
  }
  return i;
}

int main() {
  int answer[DIGIT]={1};
  int n,len=1;
  scanf("%d",&n);
  for(int i=2; i<=n; i++) {
    len = multiply(answer,len,i);
  }
  for(int i=len-1; i>=0; i--) 
    printf("%d", answer[i]);
  printf("\n");
  return 0;
}

参考文献：

$1$ 李红卫，李秉璋. C程序设计与训练（第四版） $M$ ，大连，大连理工大学出版社，2003.

$2$ https://pan.baidu.com/s/17ZXphwqySNIsIgcGtYMjvg?pwd=lhwc