高精度

From kuangbin

支持加法，乘法，int char*对 BigInt 初始化

#define mst(name, value) memset(name,value,sizeof(name))
struct BigInt {
    const static int mod=1e4;
    const static int blo=4;
    int a[600], len;
    BigInt() { mst(a, 0), len=1; }
    BigInt(int v)
    {
        mst(a, 0);
        len=0;
        do
        {
            a[len++]=v%mod;
            v/=mod;
        } while (v);
    }
    BigInt(const char s[])
    {
        mst(a, 0);
        int slen=strlen(s);
        len=slen/blo;
        if (slen%blo) len++;
        int index=0;
        for (int i=slen-1; i>=0; i-=blo)
        {
            int t=0, k=i-blo+1;
            if (k<0) k=0;
            for (int j=k; j<=i; ++j)
                t=t*10+s[j]-'0';
            a[index++]=t;
        }
    }
    BigInt operator+(const BigInt &b) const
    {
        BigInt res;
        res.len=max(len, b.len);
        for (int i=0; i<=res.len; ++i)
            res.a[i]=0;
        for (int i=0; i<res.len; ++i)
            res.a[i]+=((i<len) ? a[i] : 0)+((i<b.len) ? b.a[i] : 0),
            res.a[i+1]+=res.a[i]/mod, res.a[i]%=mod;
        if (res.a[res.len]>0) res.len++;
        return res;
    }
    BigInt operator*(const BigInt &b) const
    {
        BigInt res;
        for (int i=0; i<len; ++i)
        {
            int up=0;
            for (int j=0; j<b.len; ++j)
            {
                int temp=a[i]*b.a[j]+res.a[i+j]+up;
                res.a[i+j]=temp%mod, up=temp/mod;
            }
            if (up!=0) res.a[i+b.len]=up;
        }
        res.len=len+b.len;
        while (res.a[res.len-1]==0 && res.len>1) res.len--;
        return res;
    }
    void output()
    {
        printf("%d", a[len-1]);
        for (int i=len-2; i>=0; --i)
            printf("%04d", a[i]);
        puts("");
    }
};