C++无符号整数四则运算 压位高精度

C++无符号整数四则运算 压位高精度

储存：

利用一个数组储存，每个数组位代表按照10的整数倍数进制（所谓压位）存储一个整数，输出时指定固定位，补前导0。数组的第一位（下标为0）储存数字长度，从低位向高位存储。

加法：

直接从低位到高位按位相加，加完后看是否进位。结果位数最多增加1，注意更新位数。

减法：

同理从低位向高位按位相减，如果是负的结果就借位。结果要去除前导0。

乘法：

结论：结果的第i+j-1位是两乘数的第i位和第j位相乘的结果。注意和加法一样可能溢出一位。

除法：

模拟手算除法，从高位到低位每次取一个数，取完之后加在上次取出的数的末尾，直到这个数大于除数为止，用这个数反复减除数，直到结果小于除数，减的次数就是这一位取得数对应的位的商。

代码：

/*
 压位整数高精度
 无符号
 */
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

const int power = 1;// 每次运算的位数为10的power次方
const int base = 10;// 每一位的基数，即10的power次方，用于压位

// 压位的时候只需改power和base，如压万位高精，则power=4，base=10000

const int MAXL = 1001;// 数组的长度

char a[MAXL],b[MAXL];// 进行运算的两个数，字符串形式

struct num{
    int a[MAXL];// 储存
    num(){
        memset(a, 0, sizeof(a));
    }

    // 将输入的字符串转换为大数形式储存
    num(char *s){
        memset(a, 0, sizeof(a));
        int len = strlen(s);
        a[0] = (len + power - 1) / power;// 数的长度，向上取整算法
        // 处理字符串
        for (int i=0, t=0, w; i < len ;w *= 10, ++i)
        {
            if (i % power == 0) { w = 1, ++t; }
            a[t] += w * (s[i]-'0');
        }
    }
    void add(int k){if(k || a[0]) a[++a[0]] = k;} // 在末尾添加一个数，除法时用，假如该数为0（长度0）且要添0，则不可添加
    void re(){reverse(a+1,a+a[0]+1);} // 把数反过来，除法时用
    void print(){
        // 先打印最高位，后面还有压位修改
        printf("%d", a[a[0]]);
        // 后续各位有前导0
        for(int i = a[0]-1;i > 0;--i)
            printf("%0*d", power, a[i]);
        printf("\n");
    }
} p,q,ans;

bool operator < (const num &p, const num &q){
    if(p.a[0] < q.a[0]) return true;
    if(p.a[0] < q.a[0]) return false;
    for(int i = p.a[0];i > 0;--i){
        if(p.a[i] != q.a[i]) return p.a[i] < q.a[i];
    }
    return false;
}

num operator + (const num &p, const num &q){
    num c;
    c.a[0] = max(p.a[0], q.a[0]);
    for(int i = 1;i <= c.a[0];++i){
        c.a[i] += p.a[i] + q.a[i];
        c.a[i+1] += c.a[i] / base;
        c.a[i] %= base;
    }
    
    if(c.a[c.a[0] + 1]) ++c.a[0];
    
    return c;
}

num operator - (const num &p, const num &q){
    num c = p;
    for(int i = 1;i <= c.a[0];++i){
        c.a[i] -= q.a[i];
        if(c.a[i] < 0) {c.a[i] += base;--c.a[i+1];}
    }
    // 去除前导0,数字0的长度视为0，换句话说用长度0来表示数字0
    while(c.a[0] > 0 && !c.a[c.a[0]]) --c.a[0];
    return c;
}

num operator * (const num &p, const num &q){ // 第i位和第j位相乘得到的是结果的第i+j-1位
    num c;
    c.a[0] = p.a[0] + q.a[0] - 1;
    for(int i = 1;i <= p.a[0];++i)
        for(int j = 1;j <= q.a[0];++j){
            c.a[i+j-1] = p.a[i] * q.a[j];
            c.a[i+j] += c.a[i+j-1] / base;
            c.a[i+j-1] %= base;
    }
    // 和加法一样更新一下位数，这个也只会溢出一位
    if(c.a[c.a[0]+1]) ++c.a[0];
    return c;
}

num operator / (const num &p, const num &q){
    num x,y;// y是每次的待除数，x是每次的商
    // 从最高位开始取数，添加到y的末尾
    for(int i = p.a[0];i >= 1;--i){
        y.add(p.a[i]);// 此时是高位在前
        y.re(); // 变为低位在前，从而可以计算
        while(!(y < q)) // 一直减到剩下的数小于除数为止，其减的次数就是该位的商，如果一开始y<q的话则不会做任何操作，这样刚好代表这一位是0
            y = y - q, ++x.a[i];
        y.re(); // 方便下次计算，改回高位在前
    }
    x.a[0] = p.a[0];
    while(x.a[0] > 0 && !x.a[x.a[0]]) --x.a[0];// 和减法一样的方法除去前导0
    return x;
}

int main(){
    scanf("%s", a);
    scanf("%s", b);
    
    reverse(a, a + strlen(a));
    reverse(b, b + strlen(b));
    
    p = num(a);
    q = num(b);
    
    (p + q).print();
    
    (p - q).print();
    
    (p * q).print();
    
    (p / q).print();
}