例68 大整数乘法
问题描述
求两个不超过200位的非负整数的积。
输入
有两行,每行是一个不超过200位的非负整数,没有多余的前导0。
输出
一行,即相乘后的结果。结果里不能有多余的前导0,即如果结果是342,那么就不能输出为0342。
输入样例
12345678900
98765432100
输出样例
1219326311126352690000
(1)编程思路。
将大整数用字符串保存,编写函数void mul(char *a,char *b,char *c)实现大整数c=a*b。
在函数中用int x[201]和int y[201]分别存放两个乘数,用int z[401]来存放积。计算的中间结果也都存在数组z中。数组z长度取401 是因为两个200 位的数相乘,积最多会有400 位。x[0], y[0], z[0]都表示个位。
计算的过程基本上和小学生列竖式做乘法相同。为编程方便,不急于处理进位,而将进位问题留待最后统一处理。
在乘法过程中,数组x的第i 位和y的第j 位相乘所得的数,一定是要累加到z的第i+j 位上。这里下标i, j 都是从右往左,从0 开始数。
(2)源程序。
#include <stdio.h>
#include <string.h>
void mul(char *a,char *b,char *c)
{
int len1=strlen(a),len2=strlen(b);
int x[201],y[201],z[401];
int len=len1+len2;
memset(x,0,sizeof(x));
memset(y,0,sizeof(y));
memset(z,0,sizeof(z));
int i,j;
for (i=len1-1;i>=0;i--)
x[len1-1-i]=a[i]-'0';
for (i=len2-1;i>=0;i--)
y[len2-1-i]=b[i]-'0';
for (i=0;i<len1;i++)
{
for (j=0;j<len2;j++)
{
z[i+j]+=x[i]*y[j];
}
}
for (i=0;i<len;i++)
{
if (z[i]>=10)
{
z[i+1]+=z[i]/10;
z[i]=z[i]%10;
}
}
while (len>0 && z[len-1]==0) // 去前置0
len--;
if (len==0) // a*b=0时特判
{
c[0]='0';
c[1]='\0';
return ;
}
for (i=0;i<len;i++)
c[i]=z[len-1-i]+'0';
c[len]='\0';
}
int main()
{
char s1[201],s2[201],ans[401];
scanf("%s%s",s1,s2);
mul(s1,s2,ans);
printf("%s\n",ans);
return 0;
}
习题68
68-1 大整数除法
问题描述
求两个大的正整数相除的商。
输入
第1行是被除数,第2行是除数。每个数均不超过100位。
输出
一行,相应的商的整数部分
输入样例
2376
24
输出样例
99
(1)编程思路。
将大整数用字符串保存,编写函数void div(char *a,char *b,char *q)实现q=a/b。
实现大整数除法的基本的思想是反复做减法,从被除数里最多能减去多少个除数,商就是多少。但是不能一个一个地减除数,这样既慢又不好处理。
实际处理方法应该这样,将除数扩大10倍、100倍、…或10k倍,使得除数和被除数等长,之后再进行减法操作,依次减除数的10k倍、…、除数的100倍、除数的10倍、除数本身。每次记下够减的次数,所得结果就是商。
下面以231568除以328为例进行说明。
先将除数328扩大1000倍,即后面加上3个0,使得除数328000与被除数231568等长,231568减去328000,不够减,因此记下减的次数为0;之后,去掉除数后面的1个0,得到32800(即除数扩大100倍),231568减去32800,够减7次,余下 1968,同样记下减的次数为7;再去掉除数后面的1个0,得到3280(即除数扩大10倍),1968减3280,不够减,同样记下减的次数为0;再去掉除数后面的1个0,此时就是除数本身328,1968减328,够减6次,余下0,同样记下减的次数为6,结束运算。依次记下的减的次数为0706,去掉前导的0,商就是706。
(2)源程序。
#include <stdio.h>
#include <string.h>
void sub(char *a,char *b,char *c)
{
int len1=strlen(a),len2=strlen(b);
int x[111],y[111],z[111];
int len=len1;
memset(x,0,sizeof(x));
memset(y,0,sizeof(y));
memset(z,0,sizeof(z));
int i,j;
for (i=len1-1,j=0;i>=0;i--,j++)
x[j]=a[i]-'0';
for (i=len2-1,j=0;i>=0;i--,j++)
y[j]=b[i]-'0';
int cf=0;
for (i=0;i<len;i++)
{
z[i]=x[i]-y[i]+cf;
if (z[i]<0)
{
z[i]+=10;
cf=-1;
}
else
cf=0;
}
while (len>0 && z[len-1]==0) // 去前置0
len--;
if (len==0) // a-b=0时特判
{
c[0]='0';
c[1]='\0';
return ;
}
for (i=0;i<len;i++)
c[i]=z[len-1-i]+'0';
c[len]='\0';
}
int bigger(char *a,char *b)
{
if (strlen(a)>strlen(b)) return 1;
if (strlen(a)<strlen(b)) return 0;
if (strcmp(a,b)>=0) return 1;
else return 0;
}
void div(char *a,char *b,char *q)
{
if (bigger(a,b)==0) // 被除数a小于除数b,商为0
{
q[0]='0';
q[1]='\0';
return ;
}
int len1=strlen(a);
int len2=strlen(b);
int i;
for (i=len2;i<len1;i++) // 将b的末尾添加0到与a等长
b[i]='0';
b[i]='\0';
for (i=0;i<=len1-len2;i++)
{
int times=0;
while (bigger(a,b))
{
times++;
sub(a,b,a);
}
q[i]=times+'0';
b[strlen(b)-1]='\0';
}
q[i]='\0';
for (i=0;q[i]=='0';i++);
strcpy(q,&q[i]);
}
int main()
{
char a[115],b[115],q[115];
scanf("%s", a);
scanf("%s", b);
div(a,b,q);
printf("%s\n",q);
return 0;
}
68-2 阶乘之和
问题描述
用高精度计算出 S = 1!+2!+3!+⋯+n!(n≤50)。
其中“!”表示阶乘,例如:5! = 5×4×3×2×1。
输入格式
一个正整数n。
输出格式
一个正整数 S,表示计算结果。
输入样例
3
输出样例
9
(1)编程思路。
编写函数void bigNumMul(char a[],int b,char c[])实现一个大整数a乘以一个int型整数,得到乘积为大整数c。
编写函数void bigNumAdd(char a[],char b[],char c[])实现大整数c=a+b。
(2)源程序。
#include <stdio.h>
#include <string.h>
#define MAX_LEN 201
void bigNumAdd(char a[],char b[],char c[])
{
int len1=strlen(a),len2=strlen(b);
int x[MAX_LEN],y[MAX_LEN],z[MAX_LEN];
int len=len1>len2?len1:len2;
memset(x,0,sizeof(x));
memset(y,0,sizeof(y));
memset(z,0,sizeof(z));
int i;
for (i=len1-1;i>=0;i--)
x[len1-1-i]=a[i]-'0';
for (i=len2-1;i>=0;i--)
y[len2-1-i]=b[i]-'0';
int cf=0;
for (i=0;i<len;i++)
{
z[i]=(x[i]+y[i]+cf)%10;
cf=(x[i]+y[i]+cf)/10;
}
z[len++]=cf;
while (len>0 && z[len-1]==0) // 去前置0
len--;
if (len==0) // a+b=0时特判
{
c[0]='0';
c[1]='\0';
return ;
}
for (i=0;i<len;i++)
c[i]=z[len-1-i]+'0';
c[len]='\0';
}
void bigNumMul(char a[],int b,char c[]) // C=a*b
{
int num[MAX_LEN]={0},result[MAX_LEN+10]={0};
// 将a中存储的字符串形式的整数转换到num中去,
// num[0]对应于个位、num[1]对应于十位、……
int len1 = strlen(a);
int i,j;
for (i=len1-1,j=0;i>=0; i--)
num[j++] = a[i] - '0';
int len2 =0;
int t=b;
do {
len2++;
t=t/10;
} while (t!=0);
for (i=0;i < len1; i++)
result[i] = num[i]*b;
// 统一处理进位问题
int len=len1+len2;
for (i = 0; i < len; i++)
{
if (result[i] >= 10)
{
result[i+1] += result[i] / 10;
result[i] %= 10;
}
}
while (len>0 && result[len-1]==0) // 去前置0
len--;
if (len==0) // a*b=0时特判
{
c[0]='0';
c[1]='\0';
return ;
}
for (i=0;i<len;i++)
c[i]=result[len-1-i]+'0';
c[len]='\0';
}
int main()
{
char a[MAX_LEN],s[MAX_LEN];
int n,i;
scanf("%d",&n);
a[0]='1', a[1]='\0';
s[0]='1'; s[1]='\0';
for (i=2;i<=n;i++)
{
bigNumMul(a,i,a);
bigNumAdd(s,a,s);
}
printf("%s\n",s);
return 0;
}
68-3 Hillwer编码
问题描述
Hillwer编码的转换规则如下:对于每一条原码S,保证仅由26个大写字母组成。将每个字母后移R位,得到中转码S1(当S=‘XYZ’,R=2时,S1=‘ZAB’。即变成当前字母后R个字母,超过‘Z’则从‘A’开始)。接着,将中转码进行“符转数”操作,将S1每一位的ACS码(即ASCLL码)相乘,得到数串Q。转换后的编码即为Q。
输入
第1行,读入n,R。第2~n+1行,每行一条编码S。
输出
共n*2行,奇数行,每行一条中转码S1; 偶数行,每行一条转换后的编码Q。
输入样例
2 6
HELLOWORLD
LETUSGO
输出样例
NKRRUCUXRJ
10167740864629920000
RKZAYMU
20957073637500
(1)编程思路1。
将S1中每一个字符的ASCLL码相乘得到的数串Q超出了长整数的表示范围,因此需要进行高精度乘法运算。
编写函数void bigNumMul(char a[],int b,char c[])实现一个大整数a乘以一个int型整数,得到乘积为大整数c。
(2)源程序1。
#include <stdio.h>
#include <string.h>
void bigNumMul(char a[],int b,char c[]) // C=a*b
{
int num[1300]={0},result[1300+10]={0};
// 将a中存储的字符串形式的整数转换到num中去,
// num[0]对应于个位、num[1]对应于十位、……
int len1 = strlen(a);
int i,j;
for (i=len1-1,j=0;i>=0; i--)
num[j++] = a[i] - '0';
int len2 =0;
int t=b;
do {
len2++;
t=t/10;
} while (t!=0);
for (i=0;i < len1; i++)
result[i] = num[i]*b;
// 统一处理进位问题
int len=len1+len2;
for (i = 0; i < len; i++)
{
if (result[i] >= 10)
{
result[i+1] += result[i] / 10;
result[i] %= 10;
}
}
while (len>0 && result[len-1]==0) // 去前置0
len--;
if (len==0) // a*b=0时特判
{
c[0]='0';
c[1]='\0';
return ;
}
for (i=0;i<len;i++)
c[i]=result[len-1-i]+'0';
c[len]='\0';
}
int main()
{
int n,r;
scanf("%d%d",&n,&r);
r=r%26;
char s[600],q[1300];
while(n--)
{
scanf("%s",s);
q[0]='1';
q[1]='\0';
int i;
for (i=0;s[i]!='\0';i++)
{
if (s[i]+r<='Z') s[i]=s[i]+r;
else s[i]=s[i]+r-'Z'+'A'-1;
bigNumMul(q,(int)s[i],q);
}
printf("%s\n",s);
printf("%s\n",q);
}
return 0;
}
(3)编程思路2。
在上面的源程序1中,实现大整数乘法时,保存大整数的数组num中,每个数组元素保存1位数字,这样既浪费存储空间,又耗费计算时间。实际上,一个数组元素中可以保存大整数的6位数字,因为它与一个ASCII码(最多为3位整数)相乘,不会超过int型整数的表数范围。下面的源程序2就是按这个想法编写的。
(4)源程序2。
#include <stdio.h>
#include <string.h>
#define MOD 1000000
void work(char s[])
{
int q[205];
memset(q,0,sizeof(q));
int len=1;
q[1]=1;
int i,j;
for (i=0;s[i]!='\0';i++)
{
int t=(char)s[i];
for (j=1;j<=len;j++)
q[j]=q[j]*t;
int cf;
for (j=1;j<=len;j++)
{
cf=q[j]/MOD;
q[j]=q[j]%MOD;
q[j+1]+=cf;
}
if (q[len+1]!=0)
{
len++;
}
}
printf("%d",q[len]);
for (i=len-1;i>=1;i--)
printf("%06d",q[i]);
printf("\n");
}
int main()
{
int n,r;
scanf("%d%d",&n,&r);
while (n--)
{
char s[605];
scanf("%s",s);
int i;
for (i=0;s[i]!='\0';i++)
{
s[i]=(s[i]-'A'+r)%26+'A';
}
printf("%s\n",s);
work(s);
}
return 0;
}