C实现大数据计算(+ - * / %)

最新推荐文章于 2023-02-22 15:21:24 发布

sauphy

最新推荐文章于 2023-02-22 15:21:24 发布

阅读量528

点赞数

分类专栏： C

C 专栏收录该内容

13 篇文章

订阅专栏

<span style="font-family: Arial, Helvetica, sans-serif;">#include<iostream></span>

#include<cstring>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<vector>
#include<cmath>
#include<stdexcept>
using namespace std;

class integer
{
    friend istream& operator>>(istream& is,integer&);
    friend ostream& operator<<(ostream& os,const integer&);//
    friend bool operator>(const integer &,const integer&);//
    friend bool operator<(const integer &,const integer&);//
    friend bool operator>=(const integer &,const integer&);//
    friend bool operator<=(const integer &,const integer&);//
    friend bool operator==(const integer &,const integer&);//
    friend bool operator!=(const integer &,const integer&);//
    friend integer operator+(const integer&,const integer&);//
    friend integer operator-(const integer&,const integer&);//
    friend integer operator*(const integer&,const integer&);//
    friend integer operator/(const integer&,const integer&);//
    friend integer operator%(const integer&,const integer&);//
    friend integer abs(const integer&);//
    friend integer operator-(const integer& num);//
public:
    integer& operator=(const integer &);//
    integer& operator++();
    integer operator++(int);
    integer& operator--();
    integer operator--(int);
    integer& operator+=(const integer&);
    integer& operator-=(const integer&);
    integer& operator*=(const integer&);
    integer& operator/=(const integer&);
    integer& operator%=(const integer&);
    integer(const char* s="0");//
    integer(const long long int& num);//
    integer(const integer& num);//
    virtual ~integer();//

private:
    struct _int
    {
public:
        _int(const char*str="0");
        _int(const long long&);
        _int(const _int&);
        void show()const;
        ~_int();
        _int(size_t,short);
        static _int add(const _int&,const _int&);
        static _int sub(const _int&,const _int&);
        static _int mul(const _int&,const _int&);
        static _int div(const _int&,const _int&);
        static _int mod(const _int&,const _int&);
        static void intcpy(_int&,const _int&);
        static bool big(const _int&,const _int&);
        static bool low(const _int&,const _int&);
        static bool equ(const _int&,const _int&);
        short *p;
        size_t len;
    }*p;
    bool sign;
};

integer::_int::_int(const char* str):p(0),len(1)
{
    char* t_str=(char*)str;
    if(strlen(str)==0)
    {
        len=1;
        p=new short[len];
        p[0]=0;
        return;
    }
    while(*t_str=='0')t_str++;
    if(t_str==&str[strlen(str)]&&*(t_str-1)=='0')
    {
        len=1;
        p=new short[len];
        p[0]=0;
    }
    else
    {
        int t_len=strlen(t_str),t_sum=0;
        len=t_len%4==0?t_len/4:t_len/4+1;
        int index=0;
        p=new short[len];
        string s(t_str);
        reverse(s.begin(),s.end());
        string::iterator it=s.begin();
        for(it=s.begin(); it<s.end()-4; it+=4)
        {
            reverse(it,it+4);
        }
        reverse(it,s.end());
        for(int i=0; i<t_len; ++i)
        {
            t_sum=t_sum*10+s.c_str()[i]-'0';
            if((i!=0&&(i+1)%4==0)||i==t_len-1)
            {
                p[index++]=t_sum;
                t_sum=0;
            }
        }
    }
}
integer::_int::_int(const long long& num):p(0),len(1)
{
    int t_len=0;
    long long temp(num);
    while(temp>0)temp/=10,t_len++;
    len=t_len%4==0?t_len/4:t_len/4+1;
    p=new short [len];
    temp=num;
    for(int i=0; i<(int)len; ++i)
    {
        p[i]=temp%10000;
        temp/=10000;
    }
}

integer::_int::_int(const _int&other):p(0),len(1)
{
    this->len=other.len;
    p=new short [len];
    for(size_t i=0; i<len; i++)
    {
        p[i]=other.p[i];
    }
}
integer::_int::_int(size_t l,short v):p(0),len(1)
{
    len=l;
    p=new short[l];
    for(int i=0; i<(int)l; i++)p[i]=v;
}

integer::_int integer::_int:: add(const _int& a,const _int& other)
{
    int max_len,min_len;
    short *p1=0,*p2=0;
    if(other.len>=a.len)
    {
        p1=a.p,p2=other.p,max_len=other.len,min_len=a.len;
    }
    else
    {
        p1=other.p,p2=a.p,max_len=a.len,min_len=other.len;
    }
    _int sum(max_len+1,0);
    int t_num=0;
    for(int i=0; i<max_len; i++)
    {
        if(i<min_len)t_num=p1[i];
        else t_num=0;
        sum.p[i]+=t_num+p2[i];
        if(sum.p[i]>=10000)
        {
            sum.p[i]-=10000;
            sum.p[i+1]++;
        }
    }
    if(sum.p[max_len]==0)sum.len--;
    return sum;
}
integer::_int integer::_int:: sub(const _int& a,const _int& b)
{
    int max_len,min_len;
    short *p1=0,*p2=0;
    max_len=a.len,min_len=b.len;
    p1=a.p,p2=b.p;
    _int sub(max_len+1,0);
    int t_num=0;
    for(int i=0; i<max_len; ++i)
    {
        if(i<min_len)t_num=p2[i];
        else t_num=0;
        sub.p[i]+=p1[i]-t_num;
        if(sub.p[i]<0)
        {
            sub.p[i+1]--;
            sub.p[i]+=10000;
        }
    }
    sub.len--;
    while(sub.len>1&&sub.p[sub.len-1]==0)sub.len--;
    return sub;
}
integer::_int integer::_int:: mul(const _int& a,const _int& other)
{
    _int result(a.len+other.len,0);;
    for(int i=0; i<(int)result.len; i++)
        result.p[i]=0;
    int index=0;
    for(int i=0; i<(int)a.len; i++)
    {
        index=i;
        for(int j=0; j<(int)other.len; j++)
        {
            int mul=a.p[i]*other.p[j];
            result.p[index]+=mul%10000;
            if(result.p[index]>=10000)result.p[index+1]++,result.p[index]-=10000;
            result.p[++index]+=mul/10000;
            if(result.p[index]>=10000)result.p[index+1]++,result.p[index]-=10000;
        }
    }
    for(int i=result.len-1; i>0; i--)
    {
        if(result.p[i]==0)
            result.len--;
        else break;
    }
    return result;
}

integer::_int integer::_int:: div(const _int& a,const _int& b)
{
    vector<_int>v;
    v.reserve(10000);
    if(low(a,b))return "0";
    _int res("0");
    _int zero("0");
    v.push_back(_int("1"));
    _int base(a);
    _int sub_num;
    while(!low(base,b))
    {
        intcpy(sub_num,mul(v[v.size()-1],b));
        if(!low(base,sub_num))
        {
            ;
            intcpy(base,sub(base,sub_num));
            intcpy((res),add(res,v[v.size()-1]));
            v.push_back(mul(v[v.size()-1],2));
        }
        else if(equ(v[v.size()-1],"1"))
        {
            break;
        }
        else
            v.erase(v.end()-1);
    }
    return res;
}

integer::_int integer::_int:: mod(const _int& a,const _int& b)
{
    vector<_int>v;
    v.reserve(10000);
    if(low(a,b))return "0";
    _int res("0");
    _int zero("0");
    v.push_back(_int("1"));
    _int base(a);
    _int sub_num;
    while(!low(base,b))
    {
        intcpy(sub_num,mul(v[v.size()-1],b));
        if(!low(base,sub_num))
        {
            ;
            intcpy(base,sub(base,sub_num));
            intcpy((res),add(res,v[v.size()-1]));
            v.push_back(mul(v[v.size()-1],2));
        }
        else if(equ(v[v.size()-1],"1"))
        {
            break;
        }
        else
            v.erase(v.end()-1);
    }
    return base;
}
void integer::_int:: intcpy(_int& dest,const _int&src)
{
    if(&dest!=&src)
    {
        delete []dest.p;
        dest.p=0;
        dest.p=new short[src.len];
        for(int i=0; i<(int)src.len; i++)
            dest.p[i]=src.p[i];
        dest.len=src.len;
    }
}

bool integer::_int:: big(const _int& a,const _int& b)
{
    if(a.len>b.len)return true;
    if(a.len<b.len)return false;
    for(int i=a.len-1; i>=0; i--)
    {
        if(a.p[i]<b.p[i])return false;
    }
    if(a.p[0]==b.p[0])return false;
    return true;
}
bool integer::_int:: low(const _int& a,const _int& b)
{
    if(a.len<b.len)return true;
    if(a.len>b.len)return false;
    for(int i=a.len-1; i>=0; i--)
    {
        if(a.p[i]<b.p[i])return true;
    }
    if(a.p[0]==b.p[0])return false;
    return false;
}
bool integer::_int:: equ(const _int& a,const _int& b)
{
    return (!big(a,b))&&(!low(a,b));
}
void integer::_int::show()const
{
    int l=this->len;
    while(l-1>-0&&0==p[l-1])l--;
    printf("%d",p[l-1]);
    for(int i=l-2; i>=0; --i)
    {
        printf("%.4d",p[i]);
    }
}

integer::_int::~_int()
{
    delete []p;
}


integer::integer(const char* s):p(0),sign(0)
{
    if(s[0]=='-')p=new _int(s+1),sign=1;
    else p=new _int(s);
}
integer::integer(const long long int& n):p(0),sign(0)
{
    if(n<0)p=new _int(-n),sign=1;
    else p=new _int(n);
}
integer::integer(const integer& n):p(0),sign(0)
{
    sign=n.sign;
    p=new _int;
    integer::_int::intcpy(*p,*n.p);
}
istream& operator>>(istream& is,integer& num)
{
    string s;
    is>>s;
    num=integer(s.c_str());
    return is;
}
ostream& operator<<(ostream& os,const integer&num)
{
    if(abs(num)=="0")
    {
        printf("0");
        return os;
    }
    if(num.sign==1)printf("-");
    num.p->show();
    return os;
}

bool operator>(const integer &a,const integer&b)
{
    if(a.sign^b.sign)
    {
        if(a.sign)return false;
        else return true;
    }
    else
    {
        return integer::_int:: big(*a.p,*b.p);
    }
}
bool operator<(const integer &a,const integer&b)
{
    if(a.sign^b.sign)
    {
        if(a.sign)return true;
        else return false;
    }
    else
    {
        return integer::_int:: low(*a.p,*b.p);
    }
}
bool operator>=(const integer &a,const integer&b)
{
    return a>b||a==b;
}
bool operator<=(const integer &a,const integer&b)
{
    return a<b||a==b;
}
bool operator==(const integer &a,const integer&b)
{
    return !(a<b||b<a);
}
bool operator!=(const integer &a,const integer&b)
{
    return !(a==b);
}
integer operator+(const integer&a,const integer&b)
{
    if(a.sign^b.sign)
    {
        if(a.sign)
        {
            return abs(b)-abs(a);
        }
        else
        {
            return abs(a)-abs(b);
        }
    }
    else
    {
        integer t;
        t.sign=a.sign;
        integer::_int::intcpy(*t.p,integer::_int::add(*a.p,*b.p));
        // cout<<t<<endl;
        return t;
    }
}
integer operator-(const integer&a,const integer&b)
{
    if(a.sign^b.sign)
    {
        if(a.sign)
        {
            integer t=abs(a)+abs(b);
            t.sign=a.sign;
            return t;
        }
        else
        {
            return abs(a)+abs(b);
        }
    }
    else
    {
        if(a.sign==1)
        {
            return -(abs(a)-abs(b));
        }
        else
        {
            integer t;
            if(abs(a)>=abs(b))
            {

                t.sign=0;
                integer::_int::intcpy(*t.p,integer::_int::sub(*a.p,*b.p));
            }
            else
            {
                t.sign=1;
                integer::_int::intcpy(*t.p,integer::_int::sub(*b.p,*a.p));
            }
            return t;
        }
    }
}
integer operator*(const integer&a,const integer&b)
{
    integer t;
    t.sign=a.sign^b.sign;
    integer::_int::intcpy(*t.p,integer::_int::mul(*a.p,*b.p));
    return t;
}
integer operator/(const integer&a,const integer&b)
{
    integer t;
    if(abs(b)=="0")throw;
    t.sign=a.sign^b.sign;
    integer::_int::intcpy(*t.p,integer::_int::div(*a.p,*b.p));
    return t;
}
integer operator%(const integer&a,const integer&b)
{
    if(abs(b)=="0")throw ;
    if(abs(a)<abs(b))return a;
    integer t;
    t.sign=a.sign^b.sign;
    integer::_int::intcpy(*t.p,integer::_int::mod(*a.p,*b.p));
    return t;
}
integer abs(const integer & a)
{
    integer t(a);
    t.sign=0;
    return t;
}
integer operator -(const integer&a)
{
    integer t(a);
    t.sign^=1;
    if(abs(a)=="0")t.sign=0;
    return t;
}
integer& integer::operator =(const integer&a)
{
    if(this!=&a)
    {
        sign=a.sign;
        integer::_int::intcpy(*p,*a.p);
    }
    return *this;
}

integer& integer::operator++()
{
    return *this+=1;
}
integer integer:: operator++(int)
{
    integer t=*this;
    *this+=1;
    return t;
}
integer& integer::operator--()
{
    return *this-=1;
}
integer integer::operator--(int)
{
    integer t=*this;
    *this-=1;
    return t;
}
integer& integer::operator+=(const integer&a)
{
    return *this=*this+a;
}
integer& integer::operator-=(const integer&a)
{
    return *this=*this-a;
}
integer& integer::operator*=(const integer&a)
{
    return *this=*this*a;
}
integer& integer::operator/=(const integer&a)
{
    return *this=*this/a;
}
integer& integer::operator%=(const integer&a)
{
    return *this=*this%a;
}
integer::~integer()
{
    delete p;
}

int main()
{
    integer a,b;
    cout<<"input [a,b]"<<endl;
    while(cin>>a>>b,!cin.eof())
    {
        cout<<"a+b= "<<a+b<<endl;
        cout<<"a-b= "<<a-b<<endl;
        cout<<"a*b= "<<a*b<<endl;
        cout<<"a/b= "<<a/b<<endl;
        cout<<"a%b= "<<a%b<<endl;
        //cout<<"a++= "<<a++<<endl;
        //cout<<"a--= "<<a--<<endl;
        //cout<<"++b= "<<++b<<endl;
        //cout<<"--b= "<<--b<<endl;

        //cout<<"a>b= "<<(a>b)<<endl;
        //cout<<"a<b= " <<(a<b)<<endl;
        //cout<<"a==b = "<<(a==b)<<endl;
    }
    return 0;
}

原文链接：http://blog.youkuaiyun.com/ultimater/article/details/8764230