高精度 <加减乘除>

// by spli
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<iostream>
using namespace std;

const int N=10010;
const int p=4;
const int base=10000;

struct node{
    char s[N];
    int a[N];
    void cl(){memset(a,0,sizeof(a));}
    void update(char *t){
        memset(a,0,sizeof(a));
        int len=strlen(t);
        reverse(t,t+len);
        a[0]=(p+len-1)/p;
        for(int i=0,c=0,w;i<len;++i,w*=10){
            if(i%p==0) c++,w=1;
            a[c]+=w*(t[i]-'0');
        }
        //for(int i=1;i<=a[0];++i) cout<<a[i];cout<<endl;
    }
    void pushback(int x){a[++a[0]]=x;}
    void rev(){reverse(a+1,a+1+a[0]);}
    bool operator < (const node &b)const{   //reverse后的比较
        if(a[0]!=b.a[0]) return a[0]<b.a[0];
        for(int i=a[0];i>=1;--i){
            if(a[i]!=b.a[i]) return a[i]<b.a[i];
        }
        return 0;
    }
    node operator + (const node &b)const{
        node ret;
        ret.cl();
        ret.a[0]=max(a[0],b.a[0]);
        for(int i=1;i<=ret.a[0];++i){
            ret.a[i]+=a[i]+b.a[i];
            ret.a[i+1]+=ret.a[i]/base;
            ret.a[i]%=base;
        }
        if(ret.a[ret.a[0]+1]) ret.a[0]++;
        //reverse(ret.a+1,ret.a+1+ret.a[0]);
        return ret;
    }
    node operator - (const node &b)const{
        node ret;
        for(int i=1;i<=a[0];++i) ret.a[i]=a[i];
        ret.a[0]=a[0];
        for(int i=1;i<=ret.a[0];++i){
            ret.a[i]-=b.a[i];
            if(ret.a[i]<0) ret.a[i]+=base,ret.a[i+1]-=1;
        }
        while(ret.a[0]>0&&ret.a[ret.a[0]]==0) ret.a[0]--;
        return ret;
    }
    node operator * (const node &b)const{
        node ret;
        ret.cl();
        ret.a[0]=a[0]+b.a[0];
        for(int i=1,k=0;i<=a[0];++i,++k)
            for(int j=1;j<=b.a[0];++j){
                ret.a[j+k]+=a[i]*b.a[j];
                ret.a[j+k+1]+=ret.a[j+k]/base;
                ret.a[j+k]%=base;
            }
        if(ret.a[ret.a[0]+1]) ret.a[0]++;
        while(ret.a[0]>0&&!ret.a[ret.a[0]]) ret.a[0]--;
        return ret;
    }
    node operator / (const node &b)const{
        node ret,k;
        ret.cl();k.cl();
        for(int i=a[0];i>=1;--i){
            k.pushback(a[i]);
            k.rev();
            while(!(k<b)) k=k-b,ret.a[i]++;
            k.rev();
        }
        ret.a[0]=a[0];
        while(ret.a[0]>0&&!ret.a[ret.a[0]]) ret.a[0]--;
        return ret;
    }
    void print(int fx){
        printf("%d",fx*a[a[0]]);
        for(int i=a[0]-1;i>=1;--i) printf("%0*d",p,a[i]);
        puts("");
    }
}x[3];

void add(node u,node v){
        node ans=u+v;
        ans.print(1);
    }
void sub(node u,node v){
        node ans;
        if(v<u) ans=u-v,ans.print(1);
        else ans=v-u,ans.print(-1);
    }
void mul(node u,node v){
        node ans=u*v;
        ans.print(1);
    }
void div(node u,node v){
        node ans;
        if(v<u) ans=u/v;
        else ans.cl();
        ans.print(1);
}

int main(){
    scanf("%s%s",x[1].s+1,x[2].s+1);
    x[1].update(x[1].s+1);
    x[2].update(x[2].s+1);
    add(x[1],x[2]);
    sub(x[1],x[2]);
    mul(x[1],x[2]);
    div(x[1],x[2]);
    return 0;
}

/*
19987787889 98776665
*
ans=1974327028402810185
*/

/*
19987787889 98776665
/
ans=202
*/










评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值