[算法系列之八]大数问题(高精度运算)

【大数相加】
【代码一】

/*********************************
*   日期:2015-01-28
*   作者:SJF0115
*   题目: 大数加法(高精度加法)
*   博客:
**********************************/
#include <iostream>
using namespace std;
 
string AddString(string num1,string num2){
    int len1 = num1.length();
    int len2 = num2.length();
    // 容错处理
    if(len1 <= 0){
        return num2;
    }//if
    if(len2 <= 0){
        return num1;
    }
    string result;
    int i = len1-1,j = len2-1;
    int a,b,sum,carry = 0;
    // 倒序相加
    while(i >= 0 || j >= 0 || carry > 0){
        a = i >= 0 ? num1[i] - '0' : 0;
        b = j >= 0 ? num2[j] - '0' : 0;
        // 按位相加并加上进位
        sum = a + b + carry;
        // 进位
        carry = sum / 10;
        result.insert(result.begin(),sum % 10 + '0');
        --i;
        --j;
    }//while
    return result;
}
 
 
int main(){
    string num1("72");
    string num2("874");
    string result = AddString(num1,num2);
    // 输出
    cout<<result<<endl;
    return 0;
}

 

【代码二】

/*********************************
*   日期:2015-01-28
*   作者:SJF0115
*   题目: 大数加法(高精度加法)
*   博客:
**********************************/
#include <iostream>
using namespace std;
 
string AddString(string num1,string num2){
    int len1 = num1.length();
    int len2 = num2.length();
    int max = len1 > len2 ? (len1 + 2) : (len2 + 2);
    char* result = new char[max];
    int index = 0;
    int i = len1 - 1,j = len2 - 1;
    int numA,numB,carry = 0,sum = 0;
    //加法
    while(i >= 0 || j >= 0 || carry > 0){
        numA = i >= 0 ? (num1[i] - '0') : 0;
        numB = j >= 0 ? (num2[j] - '0') : 0;
        sum = numA + numB + carry;
        carry = sum / 10;
        result[index++] = sum % 10 + '0';
        --i;
        --j;
    }//while
    result[index] = '\0';
    //反转
    for(i = 0,j = index - 1;i < j;++i,--j){
        char temp = result[i];
        result[i] = result[j];
        result[j] = temp;
    }//for
    return string(result);
}
 
 
int main(){
    string num1("872");
    string num2("874");
    string result = AddString(num1,num2);
    // 输出
    cout<<result<<endl;
    return 0;
}

 

【代码三】

#include<stdio.h>  
#include<string.h>  
  
char a[10001],b[10001],sum[10002];  
  
int BigIntegerAdd(){  
    //两个数的长度  
    int lena = strlen(a);  
    int lenb = strlen(b);  
    //进位标记  
    int c = 0;  
    int i,j,k;  
    //初始化数组sum  
    memset(sum,0,10001);  
    //倒序相加  
    k = 0;  
    for(i = lena-1,j = lenb-1;i >= 0 && j >= 0;i--,j--){  
        sum[k] = a[i] + b[j] - '0' + c;  
        c = 0;  
        //如果相加大于等于10  
        if(sum[k] > '9'){  
            sum[k] -= 10;  
            c = 1;  
        }  
        k++;  
    }  
    //a > b  
    while(i >= 0){  
        sum[k] = a[i] + c;  
        c = 0;  
        //如果相加大于等于10  
        if(sum[k] > '9'){  
            sum[k] -= 10;  
            c = 1;  
        }  
        k++;  
        i--;  
    }  
    //b > a  
    while(j >= 0){  
        sum[k] = b[j] + c;  
        c = 0;  
        //如果相加大于等于10  
        if(sum[k] > '9'){  
            sum[k] -= 10;  
            c = 1;  
        }  
        k++;  
        j--;  
    }  
    //如果最后两个数相加有进位的情况  
    //例如:67 + 51:5+6有进位  
    if(c == 1){  
        sum[k++] = '1';  
    }  
    //翻转  
    char temp;  
    for(i = 0,j = k-1;i < j;i++,j--){  
        temp = sum[i];  
        sum[i] = sum[j];  
        sum[j] = temp;  
    }  
    return 0;  
}  
  
int main(){  
    while(scanf("%s %s",a,b) != EOF){  
        BigIntegerAdd();  
        puts(sum);  
    }  
}  

 

【大数相减】

【代码一】

string MinusString(string num1, string num2) {
    int len1 = num1.length();
    int len2 = num2.length();
    // 相等
    if(num1 == num2){
        return "0";
    }//if
    // 正负
    bool positive = true;
    if(len1 < len2 || (len1 == len2 && num1 < num2)){
        positive = false;
        // 交换使之num1 > num2
        string tmp = num1;
        num1 = num2;
        num2 = tmp;
        int temp = len1;
        len1 = len2;
        len2 = temp;
    }//if
    string result;
    int i = len1 - 1,j = len2 - 1;
    int a,b,sum,carray = 0;
    // 从低位到高位对位做减法
    while(i >= 0 || j >= 0){
        a = i >= 0 ? num1[i] - '0' : 0;
        b = j >= 0 ? num2[j] - '0' : 0;
        sum = a - b + carray;
        carray = 0;
        // 不够减
        if(sum < 0){
            sum += 10;
            carray = -1;
        }//if
        result.insert(result.begin(),sum + '0');
        --i;
        --j;
    }//while
    // 删除前导0
    string::iterator it = result.begin();
    while(it != result.end() && *it == '0'){
        ++it;
    }//while
    result.erase(result.begin(),it);
    return positive ? result : "-"+result;
}

 

【代码二】

/* 
1.比较减数与被减数的长度,确定正负号 
2.模拟减法运算 
    (1)a[i] == b[j]  
    (2)a[i] < b[j] //借位 
    (3)a[i] > b[j]  
 */  
#include<stdio.h>  
#include<string.h>  
  
char a[10001],b[10001],sum[10001],temp[10001];  
  
int BigIntegerMinus(){  
    //两个数的长度  
    int lena = strlen(a);  
    int lenb = strlen(b);  
    //借位标记  
    int c = 0;  
    int i,j,k,positive = 1;  
    //判断正负号  
    if(lena < lenb || (lena == lenb && strcmp(a,b) < 0) ){  
        strcpy(temp,a);  
        strcpy(a,b);  
        strcpy(b,temp);  
        positive = -1;//负号  
        lena = strlen(a);  
        lenb = strlen(b);  
    }  
    //倒序相减   
    k = 0;    
    for(i = lena-1,j = lenb-1;i >= 0 && j >= 0;i--,j--){    
        sum[k] = a[i] - b[j] + '0' + c;    
        c = 0;    
        //如果相减小于0    
        if(sum[k] < '0'){    
            sum[k] += 10;    
            c = -1;    
        }    
        k++;    
    }    
    //a > b    
    while(i >= 0){    
        sum[k] = a[i] + c;    
        c = 0;    
        //如果相减小于0    
        if(sum[k] < '0'){    
            sum[k] += 10;    
            c = -1;    
        }    
        k++;    
        i--;    
    }    
    //b > a    
    while(j >= 0){    
        sum[k] = b[j] + c;    
        c = 0;    
        //如果相减小于0    
        if(sum[k] < '0'){    
            sum[k] += 10;    
            c = -1;    
        }    
        k++;    
        j--;    
    }  
    int index = k - 1;  
    //检索高位,有可能相减为0  
    while(sum[index] == '0' && index > 0){  
        index--;  
    }  
    //正号  
    if(positive == 1){  
        sum[index+1] = '\0';  
    }  
    //负号  
    else{  
        sum[++index] = '-';  
        sum[index+1] = '\0';  
    }  
    //翻转  
    char t;  
    for(i = 0,j = index;i < j;i++,j--){  
        t = sum[i];  
        sum[i] = sum[j];  
        sum[j] = t;  
    }  
    return 0;  
}  
  
int main(){  
    while(scanf("%s %s",a,b) != EOF){  
        BigIntegerMinus();  
        puts(sum);  
    }  
}  

 

【大数乘法】

【代码一】

/*********************************
*   日期:2015-01-28
*   作者:SJF0115
*   题目: 高精度乘法(大数乘法)
*   博客:
**********************************/
#include <iostream>
#include <cstring>
using namespace std;
 
string MultiplyString(string num1, string num2) {
    int len1 = num1.length();
    int len2 = num2.length();
    // 容错处理
    if(len1 <= 0 || len2 <= 0){
        return "";
    }//if
    int sum = 0;
    int len3 = len1 + len2;
    char result[len3];
    memset(result,'0',sizeof(result[0])*(len3+1));
    for(int i = len1 - 1,m = 0;i >= 0;--i,++m){
        for(int j = len2 - 1,n = 0;j >= 0;--j,++n){
            sum = (num1[i] - '0') * (num2[j] - '0') + result[m+n] - '0';
            result[m+n] = sum % 10 + '0';
            // 进位
            result[m+n+1] += sum / 10;
        }//for
    }//for
    //去掉前导0 确定乘积的位数
    while(result[len3] == '0' && len3 > 0){
        --len3;
    }//while
    //注意:加'\0'
    result[len3+1] = '\0';
    //翻转
    int temp;
    for(int i = 0,j = len3;i < j;++i,--j){
        temp = result[i];
        result[i] = result[j];
        result[j] = temp;
    }//for
    return string(result);
}
 
int main(){
    string num1("2");
    string num2("123");
    string result = MultiplyString(num1,num2);
    // 输出
    cout<<result<<endl;
    return 0;
}<span style="color:#ff0000;">
</span>

 

【代码二】

#include<stdio.h>  
#include<string.h>  
  
char a[10001],b[10001],c[20002];  
  
int BigIntegerMulti()  
{  
    int i,j,lena,lenb,lenc;  
    //初始化  
    memset(c,'0',20002);  
    //两个数的长度  
    lena = strlen(a);  
    lenb = strlen(b);  
    //乘积最大位数  
    lenc = lena + lenb - 1;  
    //翻转    
    char temp;    
    for(i = 0,j = lena-1;i < j;i++,j--){    
        temp = a[i];    
        a[i] = a[j];    
        a[j] = temp;    
    }    
    for(i = 0,j = lenb-1;i < j;i++,j--){    
        temp = b[i];    
        b[i] = b[j];    
        b[j] = temp;    
    }   
    //倒序相乘  
    /* 
          123 
        *  54 
        ------ 
          492 
         615 
        ----- 
         6642  
    */  
    int t;  
    for(i = 0;i < lena;i++){  
        for(j = 0;j < lenb;j++){  
            /* c[i+j] = (a[i] - '0') * (b[j] - '0') + c[i+j];  
               9 * 9 + '0' 超出char范围越界所以设置一个int临时变量t    */  
            t = (a[i] - '0') * (b[j] - '0') + c[i+j] - '0';  
            //进位  
            c[i+j+1] = t / 10 + c[i+j+1];  
            c[i+j] = t % 10 + '0';  
        }  
    }  
    //确定乘积的位数  
    while(c[lenc] == '0' && lenc >0){  
        lenc--;  
    }  
    //注意:加'\0'  
    c[lenc+1] = '\0';  
    //翻转  
    for(i = 0,j = lenc;i < j;i++,j--){    
        temp = c[i];    
        c[i] = c[j];    
        c[j] = temp;    
    }   
    return 0;  
}  
  
int main(){  
    while(scanf("%s %s",a,b) != EOF){  
        BigIntegerMulti();  
        puts(c);  
    }  
    return 0;  
}  

 

【代码三】

/*********************************
*   日期:2015-01-29
*   作者:SJF0115
*   题目: Karatsuba快速相乘算法
*   博客:
**********************************/
#include <iostream>
#include <string>
#include <algorithm>
using namespace std;
 
// given two unequal sized bit strings, converts them to
// same length by adding leading 0s in the smaller string. Returns the
// the new length
int MakeSameLen(string& num1,string& num2){
    int len1 = num1.length();
    int len2 = num2.length();
    if(len1 < len2){
        for(int i = 0;i < len2 - len1;++i){
            num1 = "0" + num1;
        }//for
        return len2;
    }//if
    else{
        for(int i = 0;i < len1 - len2;++i){
            num2 = "0" + num2;
        }//for
        return len1;
    }//else
}
// big number minus function
string MinusString(string num1, string num2) {
    int len1 = num1.length();
    int len2 = num2.length();
    // 相等
    if(num1 == num2){
        return "0";
    }//if
    // 正负
    bool positive = true;
    if(len1 < len2 || (len1 == len2 && num1 < num2)){
        positive = false;
        // 交换使之num1 > num2
        string tmp = num1;
        num1 = num2;
        num2 = tmp;
        int temp = len1;
        len1 = len2;
        len2 = temp;
    }//if
    string result;
    int i = len1 - 1,j = len2 - 1;
    int a,b,sum,carray = 0;
    // 从低位到高位对位做减法
    while(i >= 0 || j >= 0){
        a = i >= 0 ? num1[i] - '0' : 0;
        b = j >= 0 ? num2[j] - '0' : 0;
        sum = a - b + carray;
        carray = 0;
        // 不够减
        if(sum < 0){
            sum += 10;
            carray = -1;
        }//if
        result.insert(result.begin(),sum + '0');
        --i;
        --j;
    }//while
    // 删除前导0
    string::iterator it = result.begin();
    while(it != result.end() && *it == '0'){
        ++it;
    }//while
    result.erase(result.begin(),it);
    return positive ? result : "-"+result;
}
// big number add function
string AddString(string num1,string num2){
    int len1 = num1.length();
    int len2 = num2.length();
    // 容错处理
    if(len1 <= 0){
        return num2;
    }//if
    if(len2 <= 0){
        return num1;
    }
    string result;
    int i = len1-1,j = len2-1;
    int a,b,sum,carry = 0;
    // 倒序相加
    while(i >= 0 || j >= 0 || carry > 0){
        a = i >= 0 ? num1[i] - '0' : 0;
        b = j >= 0 ? num2[j] - '0' : 0;
        // 按位相加并加上进位
        sum = a + b + carry;
        // 进位
        carry = sum / 10;
        result.insert(result.begin(),sum % 10 + '0');
        --i;
        --j;
    }//while
    return result;
}
// 移位
string ShiftString(string num,int len){
    if(num == "0"){
        return num;
    }//if
    for(int i = 0;i < len;++i){
        num += "0";
    }//for
    return num;
}
// Karatsuba快速相乘算法
string KaratsubaMultiply(string num1, string num2) {
    int len = MakeSameLen(num1,num2);
    if(len == 0){
        return 0;
    }//if
    // all digit are one
    if(len == 1){
        return to_string((num1[0] - '0')*(num2[0] - '0'));
    }//if
    int mid = len / 2;
    // Find the first half and second half of first string.
    string x1 = num1.substr(0,mid);
    string x0 = num1.substr(mid,len - mid);
    // Find the first half and second half of second string
    string y1 = num2.substr(0,mid);
    string y0 = num2.substr(mid,len - mid);
    // Recursively computer
    string z0 = KaratsubaMultiply(x0,y0);
    string z1 = KaratsubaMultiply(AddString(x1,x0),AddString(y1,y0));
    string z2 = KaratsubaMultiply(x1,y1);
    // (z2*10^(2*m))+((z1-z2-z0)*10^(m))+(z0)
    // z2*10^(2*m)
    string r1 = ShiftString(z2,2*(len - mid));
    // (z1-z2-z0)*10^(m)
    string r2 = ShiftString(MinusString(MinusString(z1,z2),z0),len - mid);
    return  AddString(AddString(r1,r2),z0);
}
 
 
int main(){
    string num1("12345001");
    string num2("1006789");
    string result = KaratsubaMultiply(num1,num2);
    // 输出
    cout<<result<<endl;
    return 0;
}

 

本文原文地址:https://blog.csdn.net/SunnyYoona/article/details/43226505

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
1 2/3维图像分割工具箱 2 PSORT粒子群优化工具箱 3 matlab计量工具箱Lesage 4 MatCont7p1 5 matlab模糊逻辑工具箱函数 6 医学图像处理工具箱 7 人工蜂群工具箱 8 MPT3安装包 9 drEEM toolbox 10 DOMFluor Toolbox v1.7 11 Matlab数学建模工具箱 12 马尔可夫决策过程(MDP)工具箱MDPtoolbox 13 国立SVM工具箱 14 模式识别与机器学习工具箱 15 ttsbox1.1语音合成工具箱 16 分数阶傅里叶变换的程序FRFT 17 魔方模拟器与规划求解 18 隐马尔可夫模型工具箱 HMM 19 图理论工具箱GrTheory 20 自由曲线拟合工具箱ezyfit 21 分形维数计算工具箱FracLab 2.2 22 For-Each 23 PlotPub 24 Sheffield大学最新遗传算法工具箱 25 Camera Calibration 像机标定工具箱 26 Qhull(二维三维三角分解、泰森图)凸包工具箱 2019版 27 jplv7 28 MatlabFns 29 张量工具箱Tensor Toolbox 30 海洋要素计算工具箱seawater 31 地图工具箱m_map 32 othercolor配色工具包 33 Matlab数学建模工具箱 34 元胞自动机 35 量子波函数演示工具箱 36 图像局域特征匹配工具箱 37 图像分割graphcut工具箱 38 NSGA-II工具箱 39 chinamap中国地图数据工具箱(大陆地区) 40 2D GaussFit高斯拟合工具箱 41 dijkstra最小成本路径算法 42 多维数据快速矩阵乘法 43 约束粒子群优化算法 44 脑MRI肿瘤的检测与分类 45 Matlab数值分析算法程序 46 matlab车牌识别完整程序 47 机器人工具箱robot-10.3.1 48 cvx凸优化处理工具箱 49 hctsa时间序列分析工具箱 50 神经科学工具箱Psychtoolbox-3-PTB 51 地震数据处理工具CREWES1990版 52 经济最优化工具箱CompEcon 53 基于约束的重构分析工具箱Cobratoolbox 54 Schwarz-Christoffel Toolbox 55 Gibbs-SeaWater (GSW)海洋学工具箱 56 光声仿真工具箱K-Wave-toolbox-1.2.1 57 语音处理工具箱Sap-Voicebox 58 贝叶斯网工具箱Bayes Net Toolbox(BNT) 59 计算机视觉工具箱VFfeat-0.9.21 60 全向相机校准工具箱OCamCalib_v3.0 61 心理物理学数据分析工具箱Palamedes1_10_3 62 生理学研究工具箱EEGLAB 63 磁共振成像处理工具箱CONN 18b 64 matlab 复杂网络工具箱 65 聚类分析工具箱FuzzyClusteringToolbox 66 遗传规划matlab工具箱 67 粒子群优化工具箱 68 数字图像处理工具箱DIPUM Toolbax V1.1.3 69 遗传算法工具箱 70 鱼群算法工具箱OptimizedAFSAr 71 蚁群算法工具箱 72 matlab优化工具箱 73 数据包络分析工具箱 74 图像分割质量评估工具包 75 相关向量机工具箱 76 音频处理工具箱 77 nurbs工具箱 78 Nurbs-surface工具箱 79 grabit数据提取工具箱 80 量子信息工具箱QLib 81 DYNAMO工具箱 82 NEDC循环的整车油耗量 83 PlotHub工具箱 84 MvCAT_Ver02.01 85 Regularization Tools Version 4.1 86 MatrixVB 4.5(含注册) 87 空间几何工具箱 matGeom-1.2.2 88 大数计算工具箱 VariablePrecisionIntegers 89 晶体织构分析工具包 mtex-5.7.0 90 Minimal Paths 2工具箱 91 Matlab数学建模工具箱
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值