大整数乘法(运算符重载) (25分)/高精算法

BigInteger类表示不超过1000位的无符号大整数。试重载>><<*,以支持无符号大整数的输入、输出与乘法。

重载面向BigInteger类对象的运算符:

>>
<<
*

裁判测试程序样例:

#include <iostream>
#include <string>
using namespace std;

/* 请在这里填写答案 */

int main() 
{
    BigInteger a, b, c;
    cin >> a >> b;
    c = a * b;
    cout << a << "*" << b << "=" << c << endl;
    return 0;
}

输入样例:

123456789
987654321

输出样例:

123456789*987654321=121932631112635269

 代码实现:

没有满分的。。。

class BigInteger{
    public:
        BigInteger(){}
        friend ostream& operator<<(ostream &ou,BigInteger &a);
        friend istream& operator>>(istream &in,BigInteger &a);
        friend BigInteger operator*(BigInteger &a,BigInteger &b);
        string aa;
        int na[1002],re[2005],l;
};
ostream& operator<<(ostream &ou,BigInteger &a){
    for(int i=a.l-1;i>=0;i--)
        ou<<a.re[i];
    return ou;
}
istream& operator>>(istream &in,BigInteger &a){
    cin>>a.aa;
//    cout<<a.aa<<endl;
    a.l=a.aa.length();
    for(int i=0;i<a.l;i++)
        a.re[i]=a.na[i]=(a.aa[a.l-1-i]-'0');
    return in;
}
BigInteger operator*(BigInteger &a,BigInteger &b){
    BigInteger c;
    c.l=a.l+b.l;
    for(int i=0;i<a.l;i++)
        for(int j=0;j<b.l;j++)
            c.re[i+j]+=a.na[i]*b.na[j];
    while(c.re[c.l]==0)
        c.l--;
    for(int i=0;i<c.l;i++){
        c.re[i+1]+=c.re[i]/10;
        c.re[i]=c.re[i]%10;
    }
    if(c.re[c.l]!=0) c.l++;
    return c;    
}

满分(借鉴别人的来纠错)

class BigInteger
{
private:
    string str;
public:
    BigInteger() {}
    BigInteger(string s)
    {
        str=s;
    }
    friend ostream&operator<<(ostream &os, BigInteger &bigint)
    {
        os<<bigint.str;
        return os;
    }
    friend istream&operator>>(istream &is,BigInteger &bigint)
    {
        is>>bigint.str;
        return is;
    }
    friend BigInteger operator*( BigInteger &bigint1, BigInteger &bigint2);
};
string reserve(string s)
{
    int l=s.length();
    int j=l-1;
    char temp;
    for(int i=0; i<l/2; i++,j--)
    {
        temp=s[i];
        s[i]=s[j];
        s[j]=temp;
    }
    return s;
}
 
BigInteger operator*( BigInteger &bigint1, BigInteger &bigint2)
{
    string a=reserve(bigint1.str);
    string b=reserve(bigint2.str);
    if(bigint1.str[0]=='0' || bigint2.str[0]=='0')
    {
        string st="0";
        return BigInteger(st);
    }
    int m,n;
    int *c=new int[a.length()+b.length()+2]();
    for(m=0; m<a.length(); m++)
    {
        for(n=0; n<b.length(); n++)
        {
            c[m+n]+=(a[m]-'0')*(b[n]-'0');
        }
    }
 
    m--;
    n--;
    int carry=0;
    char *ch=new char[a.length()+b.length()+5];
    int p;
    for( p=0; p<=m+n; p++)
    {
        if(c[p]>9)
        {
            ch[p]=(c[p]+carry)%10+'0';
            carry=(c[p]+carry)/10;
        }
        else
        {
            if(carry+c[p]>9)
            {
                ch[p]=(carry+c[p])%10+'0';
                carry=(carry+c[p])/10;
            }
            else
            {
                ch[p]=carry+c[p]+'0';
                carry=0;
            }
        }
 
    }
    if(carry>0)
    {
        ch[p]=carry+'0';
        ch[++p]='\0';
    }
    else
        ch[p]='\0';
    string st(ch);
    string s =reserve(st);
    return BigInteger(s);
 
}

#include <iostream> #include <cstdio> #include <algorithm> #include <cstring> #include <cmath> using namespace std; const int MAXN = 85, Mod = 10000; //高精压缩大法好 int n, m; int ar[MAXN]; struct HP { int p[505], len; HP() { memset(p, 0, sizeof p); len = 0; } //这是构造函数,用于直接创建一个高精度变量 void print() { printf("%d", p[len]); for (int i = len - 1; i > 0; i--) { if (p[i] == 0) { printf("0000"); continue; } for (int k = 10; k * p[i] < Mod; k *= 10) printf("0"); printf("%d", p[i]); } } //四压缩的出 } f[MAXN][MAXN], base[MAXN], ans; HP operator + (const HP &a, const HP &b) { HP c; c.len = max(a.len, b.len); int x = 0; for (int i = 1; i <= c.len; i++) { c.p[i] = a.p[i] + b.p[i] + x; x = c.p[i] / Mod; c.p[i] %= Mod; } if (x > 0) c.p[++c.len] = x; return c; } //高精+高精 HP operator * (const HP &a, const int &b) { HP c; c.len = a.len; int x = 0; for (int i = 1; i <= c.len; i++) { c.p[i] = a.p[i] * b + x; x = c.p[i] / Mod; c.p[i] %= Mod; } while (x > 0) c.p[++c.len] = x % Mod, x /= Mod; return c; } //高精*单精 HP max(const HP &a, const HP &b) { if (a.len > b.len) return a; else if (a.len < b.len) return b; for (int i = a.len; i > 0; i--) if (a.p[i] > b.p[i]) return a; else if (a.p[i] < b.p[i]) return b; return a; } //比较取最大值 void BaseTwo() { base[0].p[1] = 1, base[0].len = 1; for (int i = 1; i <= m + 2; i++){ //这里是m! m! m! 我TM写成n调了n年... base[i] = base[i - 1] * 2; } } //预处理出2的幂 int main(void) { scanf("%d%d", &n, &m); BaseTwo(); while (n--) { memset(f, 0, sizeof f); for (int i = 1; i <= m; i++) scanf("%d", &ar[i]); for (int i = 1; i <= m; i++) for (int j = m; j >= i; j--) { //因为终值是小区间,DP自然就从大区间开始 f[i][j] = max(f[i][j], f[i - 1][j] + base[m - j + i - 1] * ar[i - 1]); f[i][j] = max(f[i][j], f[i][j + 1] + base[m - j + i - 1] * ar[j + 1]); } //用结构体重载运算符写起来比较自然 HP Max; for (int i = 1; i <= m; i++) Max = max(Max, f[i][i] + base[m] * ar[i]); ans = ans + Max; //记录到总答案中 } ans.print(); //出 return 0; } 润色一下
最新发布
06-16
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值