高精模板

高精加和乘很简单,主要是高精减和高精除,均是模拟竖式运算,同时以10000进制,高精除里面用了二分.
关键的一点就是高精乘除里面最后的判前导0部分.

#define M 1000
#define P 10000
struct BigInt{
    int num[M];
    int len;
    BigInt(){
        memset(num,0,sizeof(num));
        len=0;  
    }
    void rd(){
        char str[M];
        scanf("%s",str);
        int n=strlen(str);
        len=0;
        for(int i=n-1;i>=0;i-=4){
            int p=0;
            for(int j=max(0,i-3);j<=i;j++)p=(p<<3)+(p<<1)+(str[j]^48);
            num[len++]=p;
        }
    }
    bool operator <(const BigInt &tmp)const{
        if(len!=tmp.len)return len<tmp.len;
        for(int i=len-1;i>=0;--i)
            if(num[i]!=tmp.num[i])return num[i]<tmp.num[i];
        return 0;
    }
    BigInt operator +(const BigInt &tmp)const{
        BigInt C;
        C.len=max(len,tmp.len);
        for(int i=0;i<C.len;i++){
            C.num[i]+=num[i]+tmp.num[i];
            if(C.num[i]>=P)C.num[i]-=P,C.num[i+1]++;
        }
        if(C.num[C.len]>0)C.len++;
        return C;
    }
    BigInt operator *(const BigInt &tmp)const{
        BigInt C;
        C.len=len+tmp.len;
        for(int i=0;i<len;i++)
            for(int j=0;j<tmp.len;j++){
                C.num[i+j]+=num[i]*tmp.num[j];
                if(C.num[i+j]>=P)C.num[i+j+1]+=C.num[i+j]/P,C.num[i+j]%=P;
            }
        while(C.len>1&&C.num[C.len-1]==0)C.len--;
        return C;
    }
    BigInt operator -(const BigInt &tmp)const{
        BigInt C;
        for(int i=0;i<len;i++){
            C.num[i]+=num[i]-tmp.num[i];
            C.num[i+1]+=C.num[i]/P;
            C.num[i]%=P;
            if(C.num[i]<0)C.num[i]+=P,C.num[i+1]--;
        }
        C.len=len;
        while(C.len>1&&C.num[C.len-1]==0)C.len--;
        return C;
    }
    BigInt operator /(const BigInt &tmp)const{
        BigInt C,ans;

        for(int i=0;i<len;i++)C.num[i]=num[i];
        C.len=len;

        for(int i=len-tmp.len;i>=0;i--){
            BigInt L,R,res;
            L.len=R.len=i+1;
            L.num[i]=0,R.num[i]=P-1;
            while(L.num[i]<=R.num[i]){
                BigInt mid;
                mid.len=i+1;
                mid.num[i]=L.num[i]+R.num[i]>>1;
                if(C<mid*tmp)R.num[i]=mid.num[i]-1;
                else res=mid,L.num[i]=mid.num[i]+1;
            }
            C=C-res*tmp;
            ans=ans+res;
        }

        while(ans.len>1&&ans.num[ans.len-1]==0)ans.len--;
        return ans;
    }
    void pf(){
        printf("%d",num[len-1]);
        for(int i=len-2;i>=0;--i)printf("%04d",num[i]);
        putchar('\n');
    }
};
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值