高精度算法模板
一、高精度乘法
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
struct bign{
int d[10000];
int len;
bign(){
memset(d,0,sizeof(d));
len=0;
}
};
bign change(string str){
bign a;
a.len=str.size();
for(int i=0;i<a.len;i++){
a.d[i]=str[a.len-i-1]-'0';
}
return a;
}
bign mul(bign a,bign b){
bign c;
for(int i=0;i<a.len;i++){
for(int j=0;j<b.len;j++){
c.d[i+j]+=a.d[i]*b.d[j];
c.d[i+j+1]+=c.d[i+j]/10;
c.d[i+j]%=10;
}
}
c.len=a.len+b.len+1;
while(c.d[c.len-1]==0&&c.len>1) c.len--;
return c;
}
int main(){
string str1,str2;
cin>>str1>>str2;
bign a=change(str1);
bign b=change(str2);
bign c=mul(a,b);
for(int i=c.len-1;i>=0;i--){
cout<<c.d[i];
}
return 0;
}
二、高精度减法
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
struct bign{
int d[10000];
int len;
bign(){
memset(d,0,sizeof(d));
len=0;
}
};
bign change(string str){
bign a;
a.len=str.size();
for(int i=0;i<a.len;i++){
a.d[i]=str[a.len-i-1]-'0';
}
return a;
}
bool compare(bign a,bign b){
if(a.len>b.len) return true;
else if(a.len<b.len) return false;
else{
for(int i=a.len-1;i>=0;i--){
if(a.d[i]>b.d[i]) return true;
if(a.d[i]<b.d[i]) return false;
}
return true;
}
}
bign sub(bign a,bign b){
bign c;
if(compare(a,b)){
for(int i=0;i<a.len;i++){
c.d[c.len++]=a.d[i]-b.d[i];
if(c.d[c.len-1]<0){
a.d[i+1]--;
c.d[c.len-1]+=10;
}
}
while(c.d[c.len-1]==0&&c.len>1) c.len--;
return c;
}
else{
cout<<'-';
for(int i=0;i<b.len;i++){
c.d[c.len++]=b.d[i]-a.d[i];
if(c.d[c.len-1]<0){
b.d[i+1]--;
c.d[c.len-1]+=10;
}
}
while(c.d[c.len-1]==0&&c.len>1) c.len--;
return c;
}
}
int main(){
string str1,str2;
cin>>str1>>str2;
bign a=change(str1);
bign b=change(str2);
bign c=sub(a,b);
for(int i=c.len-1;i>=0;i--){
cout<<c.d[i];
}
return 0;
}
三、高精度加法
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
struct bign{
int d[1000];
int len;
bign(){
memset(d,0,sizeof(d));
len=0;
}
};
bign change(string str){
bign a;
a.len = str.size();
for(int i=0;i<a.len;i++){
a.d[i]=str[a.len-i-1] - '0';
}
return a;
}
bign add(bign a,bign b){
bign c;
int temp,carry=0;
for(int i=0;i<a.len || i<b.len;i++){
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;
}
int main(){
string str1,str2;
cin>>str1>>str2;
bign a=change(str1);
bign b=change(str2);
bign c=add(a,b);
for(int i=c.len-1;i>=0;i--){
cout<<c.d[i];
}
return 0;
}