高精度算法模板

转自九野的博客
因为偶尔要用,转过来备用。。

//1
#define MAXN 720
struct HP
{
    int len,s[MAXN];
    HP()
    {
        memset(s,0,sizeof(s));
        len=1;
    }
    HP operator =(const char *num)
    {
        len=strlen(num);
        for(int i=0;i<len;i++) s[i]=num[len-i-1]-'0';
    }
    HP operator =(int num)
    {
        char s[MAXN];
        sprintf(s,"%d",num);
        *this=s;
        return *this;
    }
    HP(int num) { *this=num;}
    HP(const char*num) {*this=num;}
    string str()const
    {
        string res="";
        for(int i=0;i<len;i++) res=(char)(s[i]+'0')+res;
        if(res=="") res="0";
        return res;
    }
    HP operator +(const HP& b) const
    {
        HP c;
        c.len=0;
        for(int i=0,g=0;g||i<max(len,b.len);i++)
        {
            int x=g;
            if(i<len) x+=s[i];
            if(i<b.len) x+=b.s[i];
            c.s[c.len++]=x%10;
            g=x/10;
        }
        return c;
    }
    void clean()
    {
        while(len > 1 && !s[len-1]) len--;
    }
    HP operator += (const HP& b)
    {
        *this = *this + b;
        return *this;
    }
};

istream& operator >>(istream &in, HP& x)
{
    string s;
    in >> s;
    x = s.c_str();
    return in;
}

ostream& operator <<(ostream &out, const HP& x)
{
    out << x.str();
    return out;
}
//2
struct BigNum{
    char str[maxn];     //数字反着存方便运算
    int len;

    BigNum():len(0){
        memset(str, 0, sizeof str);
    }

    BigNum operator+(BigNum &num)
    {
        BigNum ret;
        for(int i=0, g=0; g || i<max(len, num.len); i++){
            int ts = g;
            if(i<len) ts+=(int)str[i]-'0';
            if(i<num.len) ts+=(int)num.str[i]-'0';
            g = ts/10;
            ret.str[ret.len++] = ts%10+'0';
        }
        ret.str[ret.len] = '\0';
        return ret;
    }

    BigNum operator+(int n)
    {
        BigNum num, ret;
        sprintf(num.str, "%d", n);
        num.len = strlen(num.str);
        for(int i=0; i<num.len/2; i++){
            swap(num.str[i], num.str[num.len-i-1]);
        }
        ret = num+*this;
        return ret;
    }

    BigNum operator=(char c[]){
        BigNum ret;
        ret.len = strlen(c);
        for(int i=ret.len-1; i>=0; i--){
            ret.str[i] = c[i];
        }
        ret.str[len] = '\0';
        return ret;
    }

    BigNum operator=(int num){
        sprintf(this->str, "%d", num);
        this->len = strlen(this->str);
        for(int i=0; i<this->len/2; i++){
            swap(this->str[i], this->str[this->len-i-1]);
        }
        return *this;
    }
};

ostream& operator <<(ostream &out, const BigNum& x){
    for(int i=x.len-1; i>=0; i--){
        out << x.str[i];
    }
    return out;
}

istream& operator >>(istream &in, BigNum& x){
    in >> x.str;;
    x.len = strlen(x.str);
    for(int i=0; i<x.len/2; i++){
        swap(x.str[i], x.str[x.len-i-1]);
    }
    return in;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
高精度算法是指能够处理超出计算机基本数据类型范围的数字运算问题的算法。常见的高精度算法有大整数加减乘除、高精度开方、高精度取模等。 以下是acwing的高精度算法模板: 1. 大整数加法 C++ 代码: vector<int> add(vector<int> &A, vector<int> &B) { vector<int> C; int t = 0; for (int i = 0; i < A.size() || i < B.size(); i++) { if (i < A.size()) t += A[i]; if (i < B.size()) t += B[i]; C.push_back(t % 10); t /= 10; } if (t) C.push_back(1); return C; } 2. 大整数减法 C++ 代码: bool cmp(vector<int> &A, vector<int> &B) { if (A.size() != B.size()) return A.size() > B.size(); for (int i = A.size() - 1; i >= 0; i--) { if (A[i] != B[i]) return A[i] > B[i]; } return true; } vector<int> sub(vector<int> &A, vector<int> &B) { vector<int> C; for (int i = 0, t = 0; i < A.size(); i++) { t = A[i] - t; if (i < B.size()) t -= B[i]; C.push_back((t + 10) % 10); if (t < 0) t = 1; else t = 0; } while (C.size() > 1 && C.back() == 0) C.pop_back(); return C; } 3. 大整数乘法 C++ 代码: vector<int> mul(vector<int> &A, int b) { vector<int> C; int t = 0; for (int i = 0; i < A.size() || t; i++) { if (i < A.size()) t += A[i] * b; C.push_back(t % 10); t /= 10; } return C; } vector<int> mul(vector<int> &A, vector<int> &B) { vector<int> C(A.size() + B.size(), 0); for (int i = 0; i < A.size(); i++) { int t = 0; for (int j = 0; j < B.size() || t; j++) { if (j < B.size()) t += A[i] * B[j]; t += C[i + j]; C[i + j] = t % 10; t /= 10; } } while (C.size() > 1 && C.back() == 0) C.pop_back(); return C; } 4. 大整数除法 C++ 代码: int cmp(vector<int> &A, vector<int> &B) { if (A.size() != B.size()) return A.size() < B.size() ? -1 : 1; for (int i = A.size() - 1; i >= 0; i--) { if (A[i] != B[i]) return A[i] < B[i] ? -1 : 1; } return 0; } vector<int> div(vector<int> &A, int b, int &r) { vector<int> C; r = 0; for (int i = A.size() - 1; i >= 0; i--) { r = r * 10 + A[i]; C.push_back(r / b); r %= b; } reverse(C.begin(), C.end()); while (C.size() > 1 && C.back() == 0) C.pop_back(); return C; } vector<int> div(vector<int> &A, vector<int> &B) { vector<int> C; for (int i = A.size() - 1; i >= 0; i--) { C.insert(C.begin(), A[i]); while (cmp(C, B) >= 0) { vector<int> t = sub(C, B); C = t; } } while (C.size() > 1 && C.back() == 0) C.pop_back(); return C; } 5. 高精度开方 C++ 代码: int cmp(vector<int> &A, vector<int> &B) { if (A.size() != B.size()) return A.size() < B.size() ? -1 : 1; for (int i = A.size() - 1; i >= 0; i--) { if (A[i] != B[i]) return A[i] < B[i] ? -1 : 1; } return 0; } vector<int> sqrt(vector<int> &A) { vector<int> C; if (A.size() % 2 == 1) A.push_back(0); for (int i = A.size() - 2; i >= 0; i -= 2) { int res = 0; for (int j = 9; j >= 0; j--) { vector<int> t = mul(C, 20); t.push_back(j * j); if (cmp(t, A) <= 0) { res = j; C.push_back(j); break; } } } reverse(C.begin(), C.end()); while (C.size() > 1 && C.back() == 0) C.pop_back(); return C; } 6. 高精度取模 C++ 代码: int mod(vector<int> &A, int b) { int r = 0; for (int i = A.size() - 1; i >= 0; i--) { r = r * 10 + A[i]; r %= b; } return r; } vector<int> mod(vector<int> &A, vector<int> &B) { vector<int> C; for (int i = A.size() - 1; i >= 0; i--) { C.insert(C.begin(), A[i]); while (cmp(C, B) >= 0) { vector<int> t = sub(C, B); C = t; } } return C; }
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值