分治法大数乘法

最新推荐文章于 2021-03-21 22:02:43 发布

oaeo

最新推荐文章于 2021-03-21 22:02:43 发布

阅读量2.5k

点赞数 1

分类专栏：算法文章标签： c struct 算法

本文链接：https://blog.csdn.net/chaoojie/article/details/7944120

版权

算法专栏收录该内容

74 篇文章 0 订阅

订阅专栏

不知道有没有更高效的算法

#include "stdio.h"
 #include "stdlib.h"
 
 #define M 100
 
 typedef struct _Node{
 	int s[M];
 	int l;
 	int c;
 }Node, *pNode;
 
 void cp(pNode src, pNode des, int st, int l){
 	int i, j;
 	for(i=st, j=0; i<st+l; i++, j++){
 		des->s[j] = src->s[i];
 	}
 	des->l = l;
 	des->c = st + src->c;
 }
 
 /*
 分治法 大数乘法
 X = A*10^n + B
 Y = C*10^m + D
 X*Y = A*C*10^(n+m) + A*D*10^n + B*C*10^m + B*D
 */
 
 void add(pNode pa, pNode pb, pNode ans){
 	int i, c;
 	int ta, tb;
 	pNode k;
 
 	if(pa->c < pb->c){
 		k = pa; pa = pb; pb = k;
 	}
 	ans->c = pb->c;
 	c = 0;
 	for(i=0; i<pa->l+pa->c-pb->c; i++){
 		if(i<pa->c-pb->c)
 			ta = 0;
 		else
 			ta = pa->s[i-(pa->c-pb->c)];
 		if(i>=pb->l)
 			tb = 0;
 		else
 			tb = pb->s[i];
 		ans->s[i] = (ta + tb + c)%10;
 		c = (ta + tb + c)/10;
 	}
 	if(c)
 		ans->s[i++] = 1;
 	ans->l = i;
 }
 
 void mul(pNode pa, pNode pb, pNode ans){
 	int i, c, w;
 	int ma = pa->l>>1, mb = pb->l>>1;
 	Node ah, al, bh, bl;
 	Node t1, t2, t3, t4, z;
 	pNode k;
 
 	if(!ma || !mb){
 		if(!ma){
 			k = pa; pa = pb; pb = k;
 		}
 		
 		ans->c = pa->c + pb->c;
 		w = pb->s[0];
 		c = 0;
 		for(i=0; i<pa->l; i++){
 			ans->s[i] = (w*pa->s[i] + c)%10;
 			c = (w*pa->s[i] + c)/10;
 		}
 		if(c) ans->s[i++] = c;
 		ans->l = i;
 		return;
 	}
 
 	cp(pa, &ah, ma, pa->l-ma);
 	cp(pa, &al, 0, ma);
 	cp(pb, &bh, mb, pb->l-mb);
 	cp(pb, &bl, 0, mb);
 	
 	mul(&ah, &bh, &t1);
 	mul(&ah, &bl, &t2);
 	mul(&al, &bh, &t3);
 	mul(&al, &bl, &t4);
 
 	add(&t3, &t4, ans);
 	add(&t2, ans, &z);
 	add(&t1, &z, ans);
 
 }
 
 int main(){
 	Node ans;
 	int i;
 	//54321 * 543
 	Node a = {{1, 2, 3, 4, 5}, 5, 0};
 	Node b = {{3, 4, 5}, 3, 0};
 	mul(&a, &b, &ans);
 	for(i=ans.l-1; i>=0; i--)
 		printf("%d", ans.s[i]);
 	printf("\n");
 	return 0;
 }