高精度模板

本模板资瓷正数范围内的以下操作,不资瓷负数(其实加一个符号判定就可以了):

打钩√表示直接资瓷,打叉x表示不直接资瓷,详见括号。

高精度 + 整数 x(最naive的操作,懒得打,也许以后会补上?)
高精度 + 高精度 √

高精度 - 整数 √
高精度 - 高精度 √

高精度 * 整数 √
高精度 * 高精度 √

高精度 / 整数 √
高精度 / 高精度 √(BASE调成10后复杂度才是对的)

高精度 % 整数 x (直接用乘法和减法做,复杂度应该也是对的)
高精度 % 高精度 √(用完除法模板之后的this就是余数)

UPD (2018.06.13) : 感谢 DntcryBecthlev 的指正

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define BASE 10000
#define L 20005
using namespace std;
int p;
char s[10*L];
struct bigint
{
    int num[L], len;
    bigint(int x = 0)
    {
        memset(num,0,sizeof(num));
        len = 1;
        num[0] = x;
    }
    bigint operator + (bigint b)
    {
        bigint c;
        c.len = max(b.len, len);
        for(int i = 0; i < c.len; i++)
        {
            c.num[i] += num[i] + b.num[i];
            c.num[i+1] = c.num[i] / BASE;
            c.num[i] %= BASE;
        }
        if(c.num[c.len])c.len++;
        return c;
    }
    bigint operator - (bigint b)
    {
        bigint c;
        c.len = max(len, b.len);
        for(int i = 0; i < c.len; i++)
        {
            c.num[i] += num[i] - b.num[i];
            if(c.num[i] < 0)
            {
                c.num[i] += BASE;
                c.num[i+1]--;
            }
        }
        while(!c.num[c.len-1] && c.len > 1)c.len--;
        return c;
    } 
    bigint operator - (int b)
    {
        bigint c;
        c.len = len;
        num[0] -= b;
        for(int i = 0; i < len; i++)
        {
            num[i+1] += num[i] / BASE;
            num[i] %= BASE;
            if(num[i] < 0)num[i] += BASE, num[i+1]--;
            c.num[i] = num[i];
        }
        while(!c.num[c.len-1] && c.len > 1)c.len--;
        return c;
    }
    bigint operator * (bigint b)
    {
        bigint c;
        c.len = len + b.len;
        for(int i = 0; i < len; i++)
        {
            for(int j = 0; j < b.len; j++)
            {
                c.num[i+j] += num[i] * b.num[j];
                c.num[i+j+1] += c.num[i+j] / BASE;
                c.num[i+j] %= BASE;
            }
        }
        if(!c.num[c.len-1] && c.len > 1)c.len--;
        return c;
    }
    bigint operator * (int b)
    {
        bigint c;
        for(int i = 0; i < len; i++)
            c.num[i] = num[i] * b;
        for(int i = 0; i < len; i++)
        {
            c.num[i+1] += c.num[i] / BASE;
            c.num[i] %= BASE;
        }
        c.len = len;
        while(c.num[c.len])c.len++;
        return c;
    }
    bool substract(bigint b, int pos)
    {
        if(len < b.len - pos)return false;
        else if(len == b.len-pos)
            for(int i = len-1; i>=0; i--)
                if(num[i] < b.num[i+pos])return false;
                else if(num[i] > b.num[i+pos])break;
        for(int i = 0; i < len; i++)
        {
            num[i] -= b.num[i+pos];
            if(num[i] < 0)
            {
                num[i] += BASE;
                num[i+1] --;
            }
        }
        while(!num[len-1] && c.len > 1)len--;
        return true;
    }
    bigint operator / (bigint b)
    {
        bigint c;
        if(len < b.len)return c;
        int k = len - b.len;
        c.len = k + 1;  
        for(int i = len-1; i>=0; i--)
        {
            if(i>=k)b.num[i] = b.num[i-k];
            else b.num[i] = 0;
        }
        b.len = len;
        for(int i = 0; i <= k; i++)
            while(this->substract(b,i)) c.num[k-i]++;
        for(int i = 0; i < c.len; i++)
        {
            c.num[i+1] += c.num[i] / BASE;
            c.num[i] %= BASE;
        }
        while(!c.num[c.len-1] && c.len > 0) c.len--; 
        return c;
    }
    bigint operator / (int b)
    {
        bigint c;
        int tmp = 0;
        for(int i = len-1; i>=0; i--)
        {
            tmp = tmp * BASE + num[i];
            c.num[i] = tmp / b;
            tmp %= b;
        }
        for(c.len = len; !c.num[c.len-1] && c.len > 1; c.len--);
        return c;
    }
    bool scan()
    {
        int n = -1;
        char ch = getchar();
        while(ch < '0' || ch > '9') if(ch == EOF)return false; else ch = getchar();
        while(ch >= '0' && ch <= '9') s[++n] = ch - '0', ch = getchar();
        len = 0;
        for(int i = n; i >= 0; i-=4)
        {
            num[len] += s[i];
            if(i>=1)num[len] += s[i-1] * 10;
            if(i>=2)num[len] += s[i-2] * 100;
            if(i>=3)num[len] += s[i-3] * 1000;
            ++len;
        }
        return true;
    }
    void clr()
    {
        memset(num,0,sizeof(num));
    }
    void print()
    {
        printf("%d",num[len-1]);
        for(int i = len-2; i>=0; i--)
            printf("%04d",num[i]);
        printf("\n");
    }
};
int main()
{ 
    printf("by ziqian\n");
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值