a[i] -= '0';
for(i=n;i >=0;i--)
b[i] -='0';
c[m+n+2] ='/0';
carry =0;
for(i=m+n;i>=0;i--)
{
sum=carry;
if((j=(i-m))<0)
j=0;
for(;j<=i&& j <=n;j++)
sum += a[i-j]b[j];
c[i+1] = sum %10 + '0'; /*算出保留的数字*/
carry = sum/10;
}
if((c[0]=carry+'0')=='0') /* if no carry*/
c[0] = '/040'; /* space */
}
int main()
{
char a[MAXLENGTH],b[MAXLENGTH],c[MAXLENGTH*2];
puts("Input multiplier");
gets(a);
puts("Input multiplier");
gets(b);
compute(a,b,c);
puts("Answer:");
puts(c);
getchar();
}
效率分析:用以上算法计算m位整数乘以n位整数,需要先进行m*n次乘法,再进行约m+n次加法运算和m+n次取模运算(实为整数除法) 。把这个程序稍加修改,让它自己生产乘数和被乘数,然后计算机随机的7200为整数互乘 。
经过改进,此算法效率可以提高约9倍 。
注意到以下事实:8216547*96785 将两数从个位起,每3位分为节,列出乘法表,将斜线间的数字相加:
8 216 547
96 785
8
216
547
*
768
20736
52512
96
6250
169560
429395
785
768
20736
52512
6250
169560
429395
768
27016
222072
429395
将表中最后一行进行如下处理:从个位数开始,每一个方格里只保留三个数字,超出1000的部分进位到前一个方格里:
768
27016
222072
429395
768+27=795
27016+222=27238
222072+429=222501
留395进429
795
238
501
395
所以8216547*96785 = 795238501395
也就是说我们在计算生成这个二维表时,不必一位一位的乘,而可以三位三位的乘;在累加时也是满1000进位 。这样,我们计算m位整数乘以n位整数,只需要进行m*n/9次乘法运算,再进行约(m+n)/3次加法运算和(m+n)/3次去摸运算 。总体看来,效率是前一种算法的9倍 。
有人可能会想:既然能用三位三位的乘,为什么不能4位4位的乘,甚至5位 。听我解来:本算法在累加表中斜线间的数字时,如果用无符号长整数(范围0至~4294967295)作为累加变量,在最不利的情况下(两个乘数的所有数字均为9),能够累加约4294967295/(999*999)=4300次,也就是能够准确计算任意两个约不超过12900(每次累加的结果“值”三位,故4300*3=12900)位的整数相乘 。如果4位4位地乘,在最不利的情况下,能过累加月4294967295/(9999*9999)=43次,仅能够确保任意两个不超过172位的整数相乘,没什么实用价值,更不要说5位了 。
[cpp] view plain copy
#include <stdio.h>
#include <string.h>
#include <conio.h>
#include <stdlib.h>
#include <time.h>
#define N 7200 //做72xx位的整数乘法
int max(int a,int b,int c)
{
int d = (a >b)?a:b;
return (d>c)?d:c;
}
int initarray(int a[])
{
int q,p,i;
q = N + random(100);
if(q%3 ==0)
p =q/3;
else
p =q/3+1;
for(i=0;i <p;i++)
a[i] = random(1000);
if(q%3 ==0)
a[0] = 100+random(900);
if(q%3 == 2)
a[0] = 10 + random(90);
if(q%3 == 1)
a[0] = 1 + random(9);
return p;
}
void write(int a[],int l)
{
int i;
char string[10];
for(i=1;i<l;i++)
{
itoa(a[i],string,10);
【电脑如何输入乘法符号,在笔记本电脑中做表格时里边做函数公式的乘号怎样打上去?】 if(strlen(string)==1)
fprintf(fp,"00");
if(strlen(string)==2)
fprintf(fp,"0");
fprintf(fp,"%s",string);
if((i+1)%25 == 0)
fprintf(fp,"/n");
}
fprintf(fp,"/n");
fprintf(fp,"/n");
}
FILE * fp;
int main()
{
int a[5000]={0},b[5000]={0},k[10001]={0};
unsigned long c,d,e;//申明作累加用的无符号长整数变量
int i,j,la,lb,ma,mi,p,q,t;
randomize();//初始化随机数
la = initarray(a);//被乘数
lb = initarray(b);//乘数
if(la < lb)//如果被乘数长度小于乘数,则交换被乘数与乘数
{
p = (lb > la)?lb:la;
for(q=0;q<p;q++)
t=a[q],a[q]=b[q],b[q]=t;
t =la,la=lb,lb =t;
}
c=d=0;
for(i=la+lb-2;i>=0;i--)//累加斜线间的数,i位横纵坐标之和
{
c=d;//将前一位的进位标志存入累加变量C
推荐阅读
- 有没有大神啊,推荐个玩 方舟生存进化 顶配的 电脑参数和价格啊?
- 怎么辨别真假蜂王浆 如何辨别蜂王浆的真假?
- 我想配一套七千块的电脑.玩方舟用的.自己感觉 i512490然后想上个3070在来个一t的硬盘?
- 发泡如何清洗 发泡清洗方法
- 社区团购好做吗?做社区团购该如何经营?
- 房间蟑螂如何彻底清除掉 房间蟑螂如何彻底清除
- 如何清理纱窗上的毛毛
- 抖音流量是如何变现的 抖音视频浏览量可以变现吗
- 大棚樱桃花期温度管理 大棚樱桃花期温度如何管理
- 辣烧排骨怎么做 排骨应该如何腌制