高精模板

来源：互联网发布：淘宝店进货渠道有哪些编辑：程序博客网时间：2024/06/02 20:02

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

#define M 1000#define P 10000struct 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