ACM模板——大数

最新推荐文章于 2020-01-11 13:04:35 发布
最新推荐文章于 2020-01-11 13:04:35 发布
阅读量69 收藏
点赞数
原文链接：http://www.cnblogs.com/Asurudo/p/10456022.html
版权 
class HP
{
    public : int len,s[maxn]; HP() {(*this) = 0;};
    HP(int inte) {(*this)=inte;}; HP(const char*str) {(*this)=str;};
    friend ostream& operator << (ostream &cout,const HP &x);
    HP operator = (int inte); HP operator = (const char*str);
    HP operator * (const HP &b);HP operator + (const HP &b);
    HP operator - (const HP &b);HP operator / (const HP &b);
    HP operator % (const HP &b);int Compare(const HP &b);
    bool operator < (const HP &b);
};

ostream& operator << (ostream &cout,const HP &x)
{for(int i = x.len;i >= 1;i --) cout<<x.s[i];return cout;}

HP HP::operator = (const char *str)
{
    len = strlen(str);
    for(int i = 1;i <= len;i ++) s[i] = str[len-i]-'0';
    return (*this);
}

HP HP::operator = (int inte)
{
    if(inte==0) {len = 1;s[1] = 0;return (*this);};
    for(len = 0;inte > 0;) {s[++len] = inte%10;inte /= 10;};
    return (*this);
}

HP HP::operator * (const HP&b)
{
    int i,j;HP c;c.len = len+b.len;
    for(i = 1;i <= c.len;i ++) c.s[i] = 0;
    for(i = 1;i <=len;i ++) for(j = 1;j <=b.len;j ++) c.s[i+j-1]+=s[i]*b.s[j];
    for(i = 1;i < c.len;i ++) {c.s[i+1]+=c.s[i]/10;c.s[i]%=10;}
    while(c.s[i]) {c.s[i+1]=c.s[i]/10;c.s[i]%=10;i ++;}
    while(i>1&&!c.s[i]) i--;c.len = i;
    return c;
}

HP HP::operator+(const HP &b)
{
    int i;HP c;c.s[1] = 0;
    for(i = 1;i <=len || i<=b.len || c.s[i];i ++)
    {
        if(i<=len) c.s[i]+=s[i];
        if(i<=b.len) c.s[i]+=b.s[i];
        c.s[i+1]=c.s[i]/10;c.s[i]%=10;
    }
    c.len = i-1;if(c.len==0) c.len = 1;
    return c;
}

HP HP::operator-(const HP&b)
{
    int i, j;HP c;
    for(i = 1,j = 0;i <= len;i ++)
    {
        c.s[i] = s[i]-j;if(i<=b.len) c.s[i]-=b.s[i];
        if(c.s[i]<0){j = 1;c.s[i]+=10;}else j = 0;
    }
    c.len = len;while(c.len>1&&!c.s[c.len]) c.len--;
    return c;
}

int HP::Compare(const HP &y)
{
    if(len>y.len) return 1;
    if(len<y.len) return -1;
    int i = len;
    while((i>1)&&(s[i]==y.s[i])) i--;
    return s[i]-y.s[i];
}

bool HP::operator < (const HP &y)
{
    if(this->Compare(y)>=0)    return false;
    return true;
}

HP HP::operator / (const HP&b)
{
    int i,j;HP d(0),c;
    for(i = len;i > 0;i --)
    {
        if(!(d.len==1 && d.s[1]==0))
            {for(j = d.len;j > 0;j --) d.s[j+1]=d.s[j];++d.len;}
        d.s[1] = s[i]; c.s[i] = 0;
        while((j = d.Compare(b))>=0)
            {d=d-b;c.s[i]++;if(j==0) break;}
    }
    c.len = len;while((c.len>1)&&(c.s[c.len]==0)) c.len--;
    return c;
}

HP HP::operator%(const HP&b)
{
    int i,j;HP d(0);
    for(i = len;i > 0;i --)
    {
        if(!(d.len==1 && d.s[1]==0))
        {for(j = d.len;j > 0;j --) d.s[j+1]=d.s[j];++d.len;}
    d.s[1] = s[i];
    while((j = d.Compare(b))>=0) {d = d-b;if(j==0)break;}
    }
    return d;
}

int main()
{
    HP a[6] {234324234,6,9,11,2,99};
    sort(a,a+6);
    _for(i,0,6)
         cout << a[i] << " ";
    cout << endl;
    return 0;
}
整数

 

#define MAXL 50000 
typedef long long ll;
const int bit=1000;
struct bign
{
    int s[MAXL],len;
    bign()
    {
        memset(s,0,sizeof s);
        len=0;
    }
    void clean()
    {
        while(len>1&&!s[len-1]) len--;
    }
    void read()
    {
        char num[MAXL];
        scanf("%s",num);
        *this=num;
    }
    void print()
    {
        printf("%d",s[len-1]);
        for(int i=len-2; i>=0; i--)
            printf("%03d",s[i]);
        putchar(10);
    }
    bign operator=(const int &num)
    {
        int tmp=num;
        len=0;
        while(tmp)
        {
            s[len++]=tmp%bit;
            tmp/=bit;
        }
        return *this;
    }
    bign operator=(const char *num)
    {
        int up=strlen(num),g=1;
        len=0;
        for(int i=up-1; i>=0; i--)
        {
            s[len]+=(num[i]-'0')*g;
            if(g==bit/10) len++;
            g=((g!=bit/10)?g*10:1);
        }
        len++;
        clean();
        return *this;
    }
    bign operator+(const bign &b)const
    {
        bign c;
        int up=max(len,b.len);
        for(int i=0,g=0; i<up||g; i++)
        {
            if(i<len) g+=s[i];
            if(i<b.len) g+=b.s[i];
            c.s[c.len++]=g%bit;
            g/=bit;
        }
        c.clean();
        return c;
    }
    bign operator-(const bign &b)const
    {
        bign c;
        bool g=0;
        for(int i=0; i<len; i++)
        {
            int x=s[i]-g;
            if(i<b.len) x-=b.s[i];
            if(x>=0) g=0;
            else x+=bit,g=1;
            c.s[c.len++]=x;
        }
        c.clean();
        return c;
    }
    bign operator*(const bign &b)const
    {
        bign c;
        c.len=len+b.len;
        for(int i=0; i<len; i++)
            for(int j=0; j<b.len; j++)
                c.s[i+j]+=s[i]*b.s[j];
        for(int i=0; i<c.len; i++)
            c.s[i+1]+=c.s[i]/bit,c.s[i]%=bit;
        c.clean();
        return c;
    }
    bign operator*(const int &b)const
    {
        bign c;
        c.len=len+32;
        for(int i=0; i<c.len; i++)
        {
            c.s[i]+=s[i]*b;
            c.s[i+1]+=c.s[i]/bit;
            c.s[i]%=bit;
        }
        c.clean();
        return c;
    }
    bign operator/(const bign &b)const
    {
        bign c,f;
        c.len=len;
        for(int i=len-1; i>=0; i--)
        {
            f=f*bit;
            f.s[0]=s[i];
            while(f>=b)
                f=f-b,c.s[i]++;
        }
        c.clean();
        return c;
    }
    bign operator/(const int b)const
    {
        bign c;
        c.len=len;
        int f=0;
        for(int i=len-1; i>=0; i--)
        {
            f=f*bit+s[i];
            if(f>=b)
                c.s[i]=f/b,f%=b;
        }
        c.clean();
        return c;
    }
    bign operator%(const bign &b)const
    {
        bign f;
        for(int i=len-1; i>=0; i--)
        {
            f=f*bit;
            f.s[0]=s[i];
            while(f>=b) f=f-b;
        }
        return f;
    }
    int operator%(const int &b)const
    {
        int f=0;
        for(int i=len-1; i>=0; i--)
        {
            f=f*bit+s[i];
            if(f>=b) f%=b;
        }
        return f;
    }
    bool operator<(const bign &b)const
    {
        if(len!=b.len) return len<b.len;
        for(int i=len-1; i>=0; i--)
            if(s[i]!=b.s[i]) return s[i]<b.s[i];
        return false;
    }
    bool operator>=(const bign &b)const
    {
        return !(*this<b);
    }
    bool operator>(const bign &b)const
    {
        if(len!=b.len) return len>b.len;
        for(int i=len-1; i>=0; i--)
            if(s[i]!=b.s[i]) return s[i]>b.s[i];
        return false;
    }
    bool operator<=(const bign &b)const
    {
        return !(*this>b);
    }
    bool operator==(const bign &b)const
    {
        if(len!=b.len) return false;
        for(int i=0; i<len; i++)
            if(s[i]!=b.s[i]) return false;
        return true;
    }
};
bign bign_sqrt(bign n)
{
    bign l,r,mid,one,ans;
    r=n;
    l.s[l.len++]=1;
    one.s[one.len++]=1;
    ans=(l+r)/2;
    while(l<r)
    {
        mid=(l+r)/2;
        if(mid*mid>n) r=mid;
        else l=mid+one,ans=mid;
    }
    return ans;
}
整数2

 

const int fprec = 100;

HP zero = 0;
class FS
{
    public:
    FS(); void SetZero();
    FS(int inte) {(*this) = inte;}
    FS(char *s) {(*this)=s;}
    FS operator=(char *s); FS operator=(int inte);
    FS operator+(FS b); FS operator-(FS b);
    FS operator*(FS b);FS operator/(FS b);
    friend ostream& operator << (ostream &cout,FS x);
    int sign,prec;
    HP num;
};

void FS::SetZero() {sign=1;num=0;prec=0;}
FS::FS() {SetZero();}

ostream& operator<<(ostream &cout,FS x)
{
    if(x.sign<0) cout << "-";
    int i,k,low=1;
    for(i = x.num.len;i > x.prec;i --) cout << x.num.s[i];
    if(x.num.len<=x.prec) cout << "0";
    if(x.num.Compare(zero)==0) {cout << ".0";return cout;}
    k = i;
    while(k>0 && x.num.s[k]==0) k--;
    if(k==0) {cout << ".0";return cout;}
    cout << ".";
    if(x.num.len<x.prec)for(int j = 0;j < x.prec-x.num.len;j ++)cout << "0";
    while(i>=low) cout << x.num.s[i--];
    return cout;
}

FS FS::operator=(int inte)
{
    prec = 0;
    if(inte>=0) {sign = 1;num = inte;}
    else {sign = -1;num = -inte;}
    return (*this);
}

FS FS::operator=(char *s)
{
    int p,i,j,l;
    SetZero();
    if(s[0]=='-') {sign = -1;s ++;};
    if(s[0]=='+') {sign = 1;s ++;};
    l = strlen(s);
    for(p = 0;p < l;p ++) if(s[p]=='.') break;
    if(p==l) prec = 0;else prec = l-1-p;
    for(i = l-1,j = 0;i >= 0;i --) if(s[i]!='.') num.s[++j] = s[i]-'0';
    while(j>1 && num.s[j]==0) --j; num.len = j;
    return (*this); 
}

void LShift(FS &a,int sl)
{
    a.prec += sl; a.num.len+=sl; int i;
    for(i = a.num.len;i > sl;i --) a.num.s[i] = a.num.s[i-sl];
    while(i>0) a.num.s[i--] = 0;
}

void RShift(FS &a,int sl)
{
    a.prec-=sl; a.num.len -= sl;int i;
    for(i = 1;i <= a.num.len;i ++) a.num.s[i] = a.num.s[i+sl];
}

FS FS::operator+(FS b)
{
    FS c;
    if(prec > b.prec) LShift(b,prec-b.prec);else
    if(prec < b.prec) LShift((*this),b.prec-prec);
    if(sign==b.sign)
    {
        c.sign = sign;c.prec = prec;c.num = num+b.num;
        if(c.num.Compare(zero)==0) c.SetZero();
    }
    else
    {
        c.prec = prec;
        if(num.Compare(b.num)==0) c.SetZero();
        else if(num.Compare(b.num)>0) {c.sign = sign;c.num = num-b.num;}
        else if(num.Compare(b.num)<0) {c.sign = b.sign;c.num = b.num-num;}
        if(c.num.Compare(zero)==0) c.SetZero();
    }
    if(c.prec > fprec) RShift(c,c.prec-fprec);
    return c;
}

FS FS::operator-(FS b)
{
    b.sign = -b.sign;
    FS c = (*this) + b;
    b.sign = -b.sign;
    return c;
}

FS FS::operator*(FS b)
{
    FS c;
    c.sign = sign*b.sign;
    c.prec = prec+b.prec;
    c.num = num*b.num;
    if(c.num.Compare(zero)==0) c.SetZero();
    if(c.prec > fprec) RShift(c,c.prec-fprec);
    return c;
}

FS FS::operator/(FS b)
{
    FS c,d;
    d = (*this);LShift(d,fprec);
    c.sign = d.sign*b.sign;
    c.prec = d.prec;
    LShift(d,b.prec);
    c.num = d.num/b.num;
    if(c.prec > fprec) RShift(c,c.prec-fprec);
    return c;
}
浮点数

 

# define FOR(i, a, b) for (int i = a; i <= b; ++i)
# define _FOR(i, a, b) for (int i = a; i >= b; --i)

using namespace std;

const int NR = 1100;

struct BigInt {
    static const int M = 1100, Dight = 4, Base = 10000;
    // M表示数组最大长度,P表示压几位,Base = 10 ^ P.
    static char str[M << 2];
    int num[M];
    BigInt()
    {
        memset(num, 0, sizeof(num));
        num[0] = 1;
    }

    void read()
    {
        scanf("%s", str);
        num[0] = 0;
        int sz = strlen(str);
        for (int i = sz - 1; i >= 0; i -= Dight) {
            num[++num[0]] = 0;
            FOR(j, max(0, i - Dight + 1), i)
                num[num[0]] = (num[num[0]] << 3) + (num[num[0]] << 1) + (str[j] ^ 48);
        }
    }

    void print()
    {
        printf("%d", num[num[0]]);
        _FOR(i, num[0] - 1, 1) printf("%0*d", Dight, num[i]);
        puts("");
    }

    bool operator < (const BigInt &cmp) const
    {
        if (num[0] != cmp.num[0]) return num[0] < cmp.num[0];
        _FOR(i, num[0], 1)
            if (num[i] != cmp.num[i]) return num[i] < cmp.num[i];
        return false;
    }
    bool operator > (const BigInt &cmp) const {return cmp < *this;}
    bool operator <= (const BigInt &cmp) const {return !(cmp < *this);}
    // bool operator != (const BigInt &cmp) const {return cmp < *this || *this < cmp;}
    // bool operator == (const BigInt &cmp) const {return !(cmp < *this || *this < cmp);}

    BigInt operator + (const BigInt &A) const
    {
        BigInt B;
        B.num[0] = max(A.num[0], num[0]);
        FOR(i, 1, B.num[0]){
            B.num[i] += num[i] + A.num[i];
            if (B.num[i] >= Base){
                B.num[i] -= Base;
                B.num[i + 1]++;
            }
        }
        while (B.num[B.num[0] + 1]) ++B.num[0];
        return B;
    }

    BigInt operator - (const BigInt &A) const
    {
        BigInt B;
        B.num[0] = max(A.num[0], num[0]);
        FOR(i, 1, B.num[0]){
            B.num[i] += num[i] - A.num[i];
            if (B.num[i] < 0){
                B.num[i] += Base;
                B.num[i + 1]--;
            }
        }
        while (B.num[0] > 1 && !B.num[B.num[0]]) --B.num[0];
        return B;
    }

    BigInt operator * (const int &x) const
    {
        BigInt B;
        B.num[0] = num[0];
        FOR(i, 1, num[0]){
            B.num[i] += num[i] * x;
            B.num[i + 1] += B.num[i] / Base;
            B.num[i] %= Base;
        }
        while (B.num[B.num[0] + 1]) ++B.num[0];
        return B;
    }

    BigInt operator * (const BigInt &A) const
    {
        BigInt B;
        B.num[0] = A.num[0] + num[0] - 1;
        FOR(i, 1, num[0])
            FOR(j, 1, A.num[0]){
                B.num[i + j - 1] += num[i] * A.num[j];
                B.num[i + j] += B.num[i + j - 1] / Base;
                B.num[i + j - 1] %= Base;
            }
        while (B.num[B.num[0] + 1]) ++B.num[0];
        return B;
    }

    BigInt operator / (const int &x) const
    {
        BigInt B = *this;
        _FOR(i, num[0], 2){
            B.num[i - 1] += B.num[i] % x * Base;
            B.num[i] /= x;
        }
        B.num[1] /= x;
        while (B.num[0] > 1 && !B.num[B.num[0]]) --B.num[0];
        return B;
    }

    BigInt operator / (const BigInt &A) const
    {
        stack<BigInt> s1, s2;
        BigInt Power = A, K;
        K.num[0] = 1; K.num[1] = 1;
        while (Power <= *this) {
            s1.push(Power); s2.push(K);
            Power = Power + Power;  K = K + K;
        }
        BigInt C = *this, B;
        while (!s1.empty()) {
            if (s1.top() <= C) {
                C = C - s1.top();
                B = B + s2.top();
            }
            s1.pop(); s2.pop();
        }
        return B;
    }
};
char BigInt::str[M << 2];
整数3

 

转载于:https://www.cnblogs.com/Asurudo/p/10456022.html

weixin_30888413
关注
ACM大数的高精度模板一套，来自zju
zck921031的专栏
06-26 3275
大数(只能处理正数) //功能挺全的，以后再作说明吧。 #include #include   #define DIGIT       4 #define DEPTH     10000 #define MAX     100 typedef int bignum_t[MAX+1];   int read(bignum_t a,istream& is=cin){
acm大数高精度模板
01-19
转载自他人博客 分享使用
ACM大数模板(c/c++)
07-20
ACM中的常用高精度模板，内容：大数加法，大数乘小数，大数乘大数，大数除法。
ACM大数模版ACM大数模版
07-29
ACM大数模版ACM大数模版ACM大数模版ACM大数模版ACM大数模版ACM大数模版ACM大数模版ACM大数模版
ACM模板——大数相加
结构之美 算法之道
04-09 2531
代码如下： #include #include using namespace std; int main() { string sum(string s1,string s2); string s1,s2; cin>>s1>>s2; cout<<sum(s1,s2)<<endl; return 0; } string sum(string s1,st
常用函数模板——ACMer者的精选.doc
06-10
在ACMer的编程竞赛生涯中，掌握一些常用的函数模板能够极大地提高解题效率。这些模板涵盖了从基础的数值计算到复杂的算法实现，是每一位热衷于算法竞赛的程序员的必备工具。以下是一些重要的函数模板及其详细说明： ...
ACM模板-清华大学
05-17
精度计算——大数阶乘 - **功能**：计算一个较大整数的阶乘。 - **输入参数**：`int n` —— 需要计算阶乘的整数。 - **输出**：返回阶乘结果的位数，并打印出阶乘的具体数值。 - **实现思路**：由于阶乘的结果可能...
史上最全ACM算法模板集——223页免费PDF
3. **大数模板，字符读入**：这部分内容可能涉及到处理大整数的算法和技巧，以及如何高效地从输入读取字符数据，这是ACM竞赛中经常遇到的问题，特别是在处理大规模数据时。 4. **数论算法**：这部分包含了一系列...
史上最全ACM算法模板集——223页精华PDF
该模板集旨在为ACM竞赛提供标准代码库，涵盖了一般函数、STL、重要公式与定理、大数处理、数论算法、图论算法等多个方面的内容，对于ACM竞赛参与者来说极具参考价值。" 在ACM算法模板集中，我们可以找到以下几个...
ACM比赛模板
08-27
精度计算——大数阶乘** 大数阶乘是指当数字很大时，普通的整型或浮点型数据类型无法存储其结果。因此，需要设计一种特殊的数据结构来处理这种大数的阶乘运算。通常会使用数组来存储每一位的数字，并通过特殊的...
ACM大数常见处理模板
09-24
帮助初学者处理acm中的常见大数问题！ 其中一些模板值得参考！！！
ACM大数四则运算模板函数
12-20
ACM大数四则运算模板函数 包含加减乘除四则运算 比较常用的算法模板
大数模板(C++)大数加法、大数乘法、大数除法
03-29
用C++写的重载的大数模板 大数加法、大数乘法、大数除法、大数减法 带有注释
c++ 大数类 大数模板
weixin_30462049的博客
04-17 256
转自：http://blog.csdn.net/hackbuteer1/article/details/6595881 这个模板 用c++类写的, 有点晕(麻烦) 自己有添加了几个模板 快速幂取余模板 const int MAXN=10000; typedef long long ll; ll fun(ll n,ll x) { ll ans=1; ...
ACM之大数模板整理
wqq的博客
01-11 508
目录：一.c++版本1.大数运算2.大数进制转换3.卡特兰数二.java大数 一.c++版本 1.大数运算 typedef long long ll; const ll m=1e8;//大数数组每个元素为8位 struct Bigint{ ll s[50];int l;//l为长度 void print(){ printf("%lld",s[l]); for(int i=l-1;i&g...
大整数类模板
weixin_30367169的博客
06-01 60
　　为了方便写，加了个符号变量neg...重载了各种比较符号,支持加减乘除取模以及模快速幂...快速幂不知道怎么用位运算,效率会偏低.. 1 #include<iostream> 2 #include<cstdio> 3 #include<string> 4 #include<vector> 5 #include&l...
大数模板
某科学的程序员
08-05 476
数字不会超大的情况下可以用 #include #include #include #include #include #define rep(i,j,k) for(int i=(int)j;i<=(int)k;i++) using namespace std; struct huge{ #define N_huge 850 #define base 100000000
C++大数模板
编程资料大全
07-10 98
分别使用C++中的运算符重载的方法来实现大数之间的数学运算，包括加法、减法、乘法、除法、n次方、取模、大小比较、赋值以及输入流、输出流的重载。。 并且使用这个大数模板，顺利AC了HDOJ上的1134这个题目的Catalan数计数问题。。http://acm.hdu.edu.cn/showproblem.php?pid=1134 大数模板的代码如下： #include&lt;iostream&g...
ACM大数的处理
zddream的专栏
05-04 1110
数据精度以及定义： 类型                      字节数                 范围                  2进制位数           输入输出 int                            4                 -2 e9~2e9                 2^31             %d  %o  %
ACM竞赛大数运算模板：加减乘除
"ACM大数四则运算模板函数提供了处理大整数加法、减法的基本算法，适用于解决计算竞赛（如ACM/ICPC）中的大数运算问题。这些模板函数通过字符串来表示大整数，并进行相应的数学操作。" 在编程竞赛或算法设计中，...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值