大数乘法的深入讨论二(分治)

《大数乘法的深入讨论一》中介绍了最基本的模拟手工算法,它的时间复杂度为O(n2),接下来将介绍一种更高效的算法——分治算法。

(原理&思路写得不太好,留个坑,以后有空了再改,但是代码一定是可以用的,亲自调过的)

算法的理论:

该算法的来源——高斯技巧(Gauss Trick):
两个复数的乘法 涉及到4个乘法:


通过下式,可以将4个乘法化简到三个:

应用到大数乘法中:
将n位二进制整数X和Y都分为2段,每段的长为n/2

由此得到:
X=A2^{n\over2}+B
,
Y=C2^{n\over2}+D

XY=(A2^{n\over2}+B)(C2^{n\over2}+D)=AC2^n+(AD+BC)2^{n\over2}+BD

变形-->
XY=AC2^n+((A-B)(D-C)+AC+BD)2^{n\over2}+BD

该算法的时间复杂度:

由主定理求得
T(n)=O(n^{log3})=O(n^{1.59})

算法的实现:
#include<bits/stdc++.h>
using namespace std;
string big_multiply(string a,string b,int n1,int n2);
int main(){
    string a,b;
    cout<<"输入待乘数:\n";
    while(cin>>a>>b){
        int n1 = a.length();
        int n2 = b.length();
        cout<<big_multiply(a,b,n1,n2)<<endl<<"\n输入待乘数:\n";
    }
    return 0;
}
string one_mul_n(string a,string b){  //1位*n位 
    vector<int> res(b.length()+2);
    string s_res;
    reverse(b.begin(),b.end());
    for(int i=0;i<b.length();i++){
        res[i]+=(a[0]-'0')*(b[i]-'0');
    }
    for(int i=0;i<res.size();i++){
        if(res[i]>9){
            res[i+1]+=res[i]/10;
            res[i]=res[i]%10;
        }
    }
    int index=res.size()-1;
    while(res[index]==0) index--;//找到非0的最高位
    for(int i=index;i>=0;i--){s_res+=res[i]+'0';} 
    return s_res;
}
string add(string a,string b,string c,string d){ //求和
    string s_res;
    int len1 = a.length()-1;
    int len2 = b.length()-1;
    int len3 = c.length()-1;
    int len4 = d.length()-1;
    int num = 0;
    int carry = 0;
    while(len1 >= 0 || len2 >= 0 || len3 >= 0 || len4 >= 0){
        num = carry;
        if(len1 >= 0) num += a[len1]-'0';
        if(len2 >= 0) num += b[len2]-'0';
        if(len3 >= 0) num += c[len3]-'0';
        if(len4 >= 0) num += d[len4]-'0';
        carry = num/10;
        num %= 10;
        s_res += (num+'0');
        --len1;--len2;--len3;--len4;
    }
    while(carry){
        s_res += (carry%10 + '0');
        carry /= 10;
    }
    reverse(s_res.begin(),s_res.end());
    return s_res;
}
string big_multiply(string a,string b,int n1,int n2){
    string s_result="";
    if((a[0]=='-'&&b[0]!='-')||(a[0]!='-'&&b[0]=='-')) s_result="-";//符号判断
    if(a[0]=='-') a.erase(0);
    if(b[0]=='-') b.erase(0);   
    if(n1 == 1){s_result+=one_mul_n(a,b);return s_result;}//1位乘n位 
    if(n2 == 1){s_result+=one_mul_n(b,a);return s_result;}
    int len1 = (n1+1)/2;
    int len2 = (n2+1)/2;
    string A,B,C,D;
    A.append(a,0,len1);//分治 
    B.append(a,len1,n1-len1);
    C.append(b,0,len2);
    D.append(b,len2,n2-len2);
    string num1 = big_multiply(A,C,len1,len2);
    string num2 = big_multiply(A,D,len1,n2-len2);
    string num3 = big_multiply(B,C,n1-len1,len2);
    string num4 = big_multiply(B,D,n1-len1,n2-len2);
    num1.append(n1-len1+n2-len2,'0');
    num2.append(n1-len1,'0');
    num3.append(n2-len2,'0');
    s_result+=add(num1,num2,num3,num4);
    return s_result;
}
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 216,372评论 6 498
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,368评论 3 392
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 162,415评论 0 353
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,157评论 1 292
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,171评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,125评论 1 297
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,028评论 3 417
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,887评论 0 274
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,310评论 1 310
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,533评论 2 332
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,690评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,411评论 5 343
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,004评论 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,659评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,812评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,693评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,577评论 2 353

推荐阅读更多精彩内容