高精度算法

高精度加法

例题：P1601【洛谷】

A+B Problem（高精）

题目描述

高精度加法，相当于 a+b problem，不用考虑负数。

输入格式

分两行输入。$a,b\le 10^{500}$。

输出格式

输出只有一行，代表 $a+b$ 的值。

样例 #1

样例输入 #1

1
1

样例输出 #1

2

样例 #2

样例输入 #2

1001
9099

样例输出 #2

10100

提示

$20\%$ 的测试数据，$0\le a,b\le 10^9$；

$40\%$ 的测试数据，$0\le a,b\le 10^{18}$。

#include<bits/stdc++.h>
using namespace std;
int main()
{
    string s1,s2;
    cin>>s1>>s2;
    int a[505],b[505],c[505];
    memset(c,0,sizeof(c));
    int s1len=s1.size(),s2len=s2.size();
    //倒序存储
    for(int i=0;i<s1len;i++) a[s1len-i]=s1[i]-'0';
    for(int i=0;i<s2len;i++) b[s2len-i]=s2[i]-'0';
    int idx=1,x=0;
    while(idx<=s1len||idx<=s2len)
    {
        c[idx]=(a[idx]+b[idx]+x)%10;
        //x表示进位的数
        x=(a[idx]+b[idx]+x)/10;
        idx++;
    }
    if(x) c[idx]=x;
    //考虑前导0情况
    while(!c[idx]) idx--;
    for(int i=idx;i>0;i--) cout<<c[i];
    //考虑结果为0情况
    if(idx<=0) cout<<0<<endl;
    return 0;
}

高精度减法

例题：P2142【洛谷】

高精度减法

题目描述

高精度减法。

输入格式

两个整数 $a,b$（第二个可能比第一个大）。

输出格式

结果（是负数要输出负号）。

样例 #1

样例输入 #1

2
1

样例输出 #1

1

提示

• $20\%$ 数据 $a,b$ 在 long long 范围内；
• $100\%$ 数据 $0<a,b\le 10^{10086}$。

#include<bits/stdc++.h>
using namespace std;
int main()
{
    string s1,s2;
    cin>>s1>>s2;
    int a[10100],b[10100],c[10100];
    memset(c,0,sizeof(c));
    int s1len=s1.size(),s2len=s2.size();
    //判断，令s1总是大于等于s2
    if((s1len<s2len)||(s1len==s2len&&s1<s2))
    {
        swap(s1,s2);
        swap(s1len,s2len);
        cout<<"-";
    }
    //倒序存储
    for(int i=0;i<s1len;i++) a[s1len-i]=s1[i]-'0';
    for(int i=0;i<s2len;i++) b[s2len-i]=s2[i]-'0';
    int idx=1;
    while(idx<=s1len)
    {
        //借位操作
        if(a[idx]<b[idx])
        {
            a[idx+1]--;
            a[idx]+=10;
        }
        c[idx]=a[idx]-b[idx];
        idx++;
    }
    //处理前导0
    while(!c[idx]) idx--;
    for(int i=idx;i>0;i--) cout<<c[i];
    //处理结果为0
    if(idx<=0) cout<<0<<endl;
    return 0;
}

高精度乘法

例题：P1303【洛谷】

A*B Problem

题目背景

高精度乘法模板题。

题目描述

给出两个非负整数，求它们的乘积。

输入格式

输入共两行，每行一个非负整数。

输出格式

输出一个非负整数表示乘积。

样例 #1

样例输入 #1

1 
2

样例输出 #1

2

提示

每个非负整数不超过 $10^{2000}$。

#include<bits/stdc++.h>
using namespace std;
//注意：高精乘数组要开两倍，因为结果的长度为s1len+s2len
int c[4000006];
int main()
{
    string s1,s2;
    cin>>s1>>s2;
    int a[2005],b[2005];
    int s1len=s1.size(),s2len=s2.size();
    //倒序存储
    for(int i=0;i<s1len;i++) a[s1len-i]=s1[i]-'0';
    for(int i=0;i<s2len;i++) b[s2len-i]=s2[i]-'0';
    for(int i=1;i<=s1len;i++)
    {
        //进位
        int x=0;
        for(int j=1;j<=s2len;j++)
        {
            int k=c[i+j-1];
            c[i+j-1]=(a[i]*b[j]+x+k)%10;
            x=(a[i]*b[j]+x+k)/10;
        }
        if(x) c[i+s2len]=x;
    }
    int idx=s1len+s2len;
    //处理前导0
    while(!c[idx]) idx--;
    for(int i=idx;i>0;i--) cout<<c[i];
    //处理结果为0
    if(idx<=0) cout<<0<<endl;
    return 0;
}

高精度除以低精度

例题：P1480【洛谷】

A/B Problem

题目描述

输入两个整数 $a,b$，输出它们的商。

输入格式

两行，第一行是被除数，第二行是除数。

输出格式

一行，商的整数部分。

样例 #1

样例输入 #1

10
2

样例输出 #1

5

提示

$0\le a\le 10^{5000}$，$1\le b\le 10^9$。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int main()
{
    string s;
    //后面这个可能会爆，要开long long
    ll b;
    cin>>s>>b;
    int a[5005],c[5005];
    memset(c,0,sizeof(c));
    int len=s.size();
    //不用倒序存储
    for(int i=0;i<len;i++) a[i]=s[i]-'0';
    int blen=to_string(b).size();
    ll x;
    (blen>1)?x=stoi(s.substr(0,blen-1)):x=0;
    for(int i=blen-1;i<len;i++)
    {
        ll temp=x*10+a[i];
        if(temp<b)
        {
            x=temp;
            continue;
        }
        c[i]=temp/b;
        x=temp-b*c[i];
    }
    int idx=blen-1;
    while(!c[idx]) idx++;
    for(int i=idx;i<len;i++) cout<<c[i];
    if(idx>=len) cout<<0<<endl;
    return 0;
}

高精度除以高精度

#include<iostream>
#include<cstring>
using namespace std;
int a[101],b[101],c[101],i;
//输入函数
void init(int a[]){
    string s;
    cin>>s;
    a[0]=s.length();
    for(i=1;i<=a[0];i++)
        a[i]=s[a[0]-i]-'0';//减法倒序存储
}
//输出函数
void print(int a[]){
    int i;
    if(a[0]==0){
        cout<<0<<endl;
        return;
    }
    for(i=a[0];i>0;i--)
        cout<<a[i];
    cout<<endl;
    return;  //函数执行完毕回到主程序
}
//比较函数
int compare(int a[],int b[]){
    int i;
     if(a[0]>b[0])
         return 1;
     if(a[0]<b[0])
         return-1;
     for(i=a[0];i>0;i--){ //如果两数位数相等，则按位比大小
         if(a[i]>b[i])
             return 1;
          if(a[i]<b[i])
              return -1;  //按位比较若该位数相同，则判断下一位
     }
    return 0;//如果返回0则表示两数相等
}
//减法模拟除法
void jian(int a[],int b[]){
    int flag,i;
    flag=compare(a,b);
    if(flag==0){
        a[0]=0;
        return;
    }
    if(flag==1){
        for(i=1;i<=a[0];i++){
            if(a[i]<b[i]){
                a[i+1]--;
                a[i]=a[i]+10;
            }
            a[i]-=b[i];
        }
        //处理前导0
        while(a[0]>0&&a[a[0]]==0)
            a[0]--;
        return;
    }
}
//复制数组
//低位补0，为了方便减法
void numcpy(int p[],int q[],int det){
    for(int i=1;i<=p[0];i++)
        q[i+det-1]=p[i];
    q[0]=p[0]+det-1;
    /*
    for(int i=q[0];i>0;i--)
        cout<<q[i];
    cout<<endl;
    打印复制后的数字，方便理解算法，此算法主要采用低位补0做减法
    */
}
//除法计算
void chugao(int a[],int b[],int c[]){
    int i,tmp[101];
    c[0]=a[0]-b[0]+1;   //商的位数不超过被除数的位数－除数的位数+1
    for(i=c[0];i>0;i--){  //每次循环确定某位商的的值，从高位开始
        memset(tmp,0,sizeof(tmp));
        numcpy(b,tmp,i);
        while(compare(a,tmp)>=0){
            c[i]++;
            jian(a,tmp);
        }
    }
    //处理前导0
    while(c[0]>0&&c[c[0]]==0)
        c[0]--;
    return;
}
//主函数
int main(){
    init(a);
    init(b);
    chugao(a,b,c);
    print(c);
    print(a);
    return 0;
}