大整数的运算
大整数的存储:
用数组存储即可:
如将整数235813存储到数组中,则有d[0]=3,d[1]=1,d[2]=8,d[3]=5,d[4]=3,d[5]=2,即整数的高位存储在数组的高位,整数的低位存储在数组的低位。不反过来存储的原因是,在进行运算的时候都是从整数的低位到高位进行枚举,顺位存储和这种思维相合。但是这也会产生一个需要注意问题:把整数按字符串%s读入的时候,实际上是逆位存储,机str[0]=‘2’,str[1]=‘3’,…,str[5]=‘3’,因此在读入之后需要在另存为至d[]数组的时候反转一下。而为了方便随时获取大整数的长度,一般都会定义一个int型变量len来记录其长度,并和d[]数组组合成结构体:
大整数的存储
struct bign{
int d[1000];
int len;
bign(){
memset(d,0,sizeof(d));
len=0;
}
}
大整数的输入:
bign change(char str[]){
bign a;
a.len=strlen(str);
for(int i=0;i<a.len;i++){
a.d[i]=str[a.len-i-1]-'0';
}
return a;
}
大整数的比较:
int compare(bign a,bign b){
if(a.len>b.len)return 1;
else if(a.len<b.len)return -1;
else{
for(int i=a.len-1;i>=0;i--){
if(a.d[i]>b.d[i])return 1;
else if(a.d[i]<b.d[i])return -1;
}
return 0;
}
}
大整数的四则运算:
这边只介绍高精度加法,高精度减法,低精度乘法,低精度除法。
高精度加法:(前提是两个数都是正数,如果有一方是负数的话,可以在转化数组时把负号去掉,然后转化成高精度的减法;如果两个都是负数的话,可以去掉负号后用高精度的加法,最后再把负号加上)
bign add(bign a,bign b){
bign c;
int carry=0; //carry是进位
for(int i=0;i<a.len||i<b.len;i++){
int temp=a.d[i]+b.d[i]+carry;
c.d[c.len++]=temp%10;
carry=temp/10;
}
if(carry!=0)c.d[c.len++]=carry;
return c;
}
高精度减法:(使用高精度减法之前要比较两个数的大小,当被减数小于减数时,要将两个数交换然后输出符号后再使用sub函数)
bign sub(bign a,bign b){
bign c;
for(int i=0;i<a.len||i<b.len;i++){
if(a.d[i]<b.d[i]){
a.d[i+1]--;
a.d[i]+=10;
}
c.d[c.len++]=a.d[i]-b.d[i];
}
while(c.len-1>=1&&c.d[c.len-1]==0)c.len--;
return c;
}
高精度与低精度的乘法:指的是bign类型和int类型的乘法
步骤:取bign的某位与int型整体相乘,再与进位相加,所得结果的个位数作为该位结果,高位部分作为新进位
如果a和b中存在负数,需要先记录下来,然后去它们的绝对值带入函数
bign multi(bign a,int b){
bign c;
int carry=0;
for(int i=0;i<a.len;i++){
int temp=a.d[i]*b+carry;
c.d[c.len++]=temp%10;
carry=temp/10;
}
while(carry!=0){
c.d[c.len++]=carry%10;
carry/=10;
}
return c;
}
高精度与低精度的除法:指的是bign类型和int型的除法
步骤:上一步的余数乘以10加上该步上的位,得到该步临时的被除数,将其与除数比较:如果不够除,则该位的商为0;如果够除,则商即为对应的商,余数即为对应的余数。最后一步要注意高位后面可能有多余的0,要忽视他们,但也要保住结果至少有一位
bign divide(bign a,int b,int &r){ //r为余数
bign c;
c.len=a.len;
for(int i=a.len-1;i>=0;i--){
r=r*10+a.d[i];
if(r<d)c.d[i]=0;
else{
c.d[i]=r/d;
r=r%d;
}
}
while(c.len-1>=1&&c.d[c.len-1]==0)c.len--;
return c;
}