面试100题系列之13大数的四则运算

关于大数的四则运算，我想网上有很多方法，感觉通用性不强，然后呢，功能不强大，效率不高~~问题多多啊~然后呢，大数这个东西又经常碰到，所以呢，我一咬牙，写了一个比较完整的，花了我一下午的时间啊，不过除法还没有写，有时间再写吧。最近事情太多了，写了之后再补上。
OK，说一下我的思路。首先这里用字符串存数据，这个地球人都知道。第0位是高位。当我们算加减的时候需要从低位开始算起。所以我一开始的思路是将字符串反转一下，然后开始计算，这样也不会担心最后有进位或者借位了，计算好之后再反转回来。OK，这样我为了算一个加法，反转四次啊！！！代价有点大，不能接受。
怎么样能既方便运算又不用反转那么多呢？要是只反转一次就好了。如果把结果的数组指针传进去，算的时候分别从两个字符串尾部开始遍历，将结果从下标0开始一次存入结果数组。最后把结果数组反转一下就可以了。好了，加减搞定了。
对于乘法，str1[i]与str2[j]乘的结果放在结果数组的i+j位，这个结论怎么来的呢？其实也就是我们写的竖式啦，写写看~这里需要一个临时的int数组arr来存放乘的结果，arr[i+j]存放所有满足条件的i和j的元素之和。然后将arr里面的数据拆分到ans数组里面就可以了。
最后说明一下，每个函数加入一个nMaxLen的参数是保证计算的结果不会溢出，也就是不会超过ans数组的最大容量，这里我采用的办法是返回一个NULL，但是这样可能与一些函数（比如说printf）嵌套使用的时候就不是很方便。这里其实可以将返回NULL的操作替换成往结果数组里面打印相应的字符串，来表明错误。
路过的大爷们，转载请注明出处：http://blog.csdn.net/kay_zhyu/article/details/8877032哦~~小女子万分感激~~~
核心代码如下，这里支持负数运算哦，欢迎大家找Bug：

//两个字符串相加
//str1和str2是加数，ans是结果数组，nMaxLen是结果数组的最大空间
char* Add(char* const str1, char* const str2, char *ans, int nMaxLen)
{
	//一些参数检验
	if(!str1 || !str2 || !ans || nMaxLen < 1)
		return NULL;
	//保证是两个整数相加
	if(str1[0] == '-' && str2[0] == '-')//负加负
	{
		ans[0] = '-';
		if(Add(str1 + 1, str2 + 1, ans + 1, nMaxLen - 1))
			return ans;
		else
			return NULL;
	}
	else if(str1[0] == '-')//一正一负
	{
		return Minus(str2, str1 + 1, ans, nMaxLen);
	}
	else if(str2[0] == '-')
	{
		return Minus(str1, str2 + 1, ans, nMaxLen);
	}
	int nLen1 = strlen(str1);
	if(nMaxLen < nLen1)
		return NULL;
	int nLen2 = strlen(str2);
	if(nMaxLen < nLen2)
		return NULL;
	int i,j,k;
	int a,b,carry;
	carry = 0;
	//加公共部分
	for(i = nLen1 - 1, j = nLen2 - 1, k = 0;
		(i >= 0 && j >= 0 && k < nMaxLen); --i,--j,++k)
	{
		a = str1[i] - '0';
		b = str2[j] - '0';
		ans[k] = (a + b + carry) % 10 + '0';
		carry = (a + b + carry) / 10;
	}
	//加str1的多余部分
	for( ; (i >= 0 && k < nMaxLen); --i, ++k)
	{
		a = str1[i] - '0';
		ans[k] = (a + carry) % 10 + '0';
		carry = (a + carry) / 10;
	}
	//加str2的多余部分
	for( ; (j >= 0 && k < nMaxLen); --j, ++k)
	{
		b = str1[j] - '0';
		ans[k] = (b + carry) % 10 + '0';
		carry = (b + carry) / 10;
	}
	if(carry > 0 && k < nMaxLen)
		ans[k++] = carry + '0';
	//结果的越界检查
	if(k >= nMaxLen)
		return NULL;
	//加上结束标记
	ans[k] = '\0';
	Invert(ans, k);
	return ans;
}
//两个大数相减
//str1-str2,ans表示结果，nMaxLen表示ans的最大长度
char* Minus(char* const str1, char * const str2, char *ans, int nMaxLen)
{
	//一些参数检验
	if(!str1 || !str2 || !ans || nMaxLen < 1)
		return NULL;
	//保证是正数减正数
	if(str1[0] != '-' && str2[0] == '-')//正减负
		return Add(str1, str2 + 1, ans, nMaxLen);
	else if(str1[0] == '-' && str2[0] == '-')//负减负
		return Minus(str2 + 1, str1 + 1, ans, nMaxLen);
	else if(str1[0] == '-' && str2[0] != '-')//负减正
	{
		ans[0] = '-';
		if(Add(str1 + 1, str2, ans + 1, nMaxLen - 1))
			return ans;
		else
			return NULL;
	}
	//保证str1大于str2
	int nLen1 = strlen(str1);
	if(nMaxLen < nLen1)
		return NULL;
	int nLen2 = strlen(str2);
	if(nMaxLen < nLen2)
		return NULL;
 //其实这里直接是str1[0]<str2[0]是有bug的，应该是strcmp(str1,str2),确认str1<str2

if(nLen1 < nLen2 || (nLen1 == nLen2 && str1[0] < str2[0])){ans[0] = '-';if(Minus(str2, str1, ans + 1, nMaxLen - 1))return ans;elsereturn NULL;}int i,j,k;int a,b,carry;bool IsPos = true;carry = 0;//减去公共的部分,str1大于str2for(i = nLen1 - 1, j = nLen2 - 1, k = 0;(i
 >= 0 && j >= 0 && k < nMaxLen); --i,--j,++k){a = str1[i] - '0' - carry;b = str2[j] - '0';if(a < b){a += 10;carry = 1;}elsecarry = 0;ans[k] = a - b + '0';}//得到str1剩余的部分for(; (i >= 0 && k < nMaxLen); --i,++k){a = str1[i] - '0' - carry;if(a < carry){a += 10;carry
 = 1;}elsecarry = 0;ans[k] = a + '0';}if(carry > 0 && k < nMaxLen)ans[k++] = '-';if(k >= nMaxLen)return NULL;ans[k] = '\0';Invert(ans, k);return ans;}//大数乘法//str1*str2，ans是结果，nMaxLen是ans的最大容量char *Mutiply(char* const str1, char* const str2, char *ans, int nMaxLen){if(!str1
 || !str2 || !ans || nMaxLen < 1)return NULL;//保证是两个正数相乘if(str1[0] == '-' && str2[0] == '-')return Mutiply(str1 + 1, str2 + 1, ans, nMaxLen);else if(str1[0] == '-'){ans[0] = '-';if(Mutiply(str1 + 1, str2, ans + 1, nMaxLen))return ans;elsereturn NULL;}else if(str2[0]
 == '-'){ans[0] = '-';if(Mutiply(str1, str2 + 1, ans + 1, nMaxLen))return ans;elsereturn NULL;}int nLen1 = strlen(str1);int nLen2 = strlen(str2);if(nMaxLen < nLen1 + nLen2 - 1)//内存不够，返回return NULL;int *arr = (int *)malloc((nLen1 + nLen2 + 1)*sizeof(int));if(!arr)//内存申请失败，返回return
 NULL;memset(arr, 0, (nLen1 + nLen2 + 1)*sizeof(int));int i,j;for(i = 0; i < nLen1; ++i){for(j = 0; j < nLen2; ++j){arr[i+j] += (str1[i] - '0') * (str2[j] - '0');}}int k,carry;carry = 0;//乘积的结果转移到ans中for(k = 0, i = nLen1 + nLen2 - 2; i >= 0 && k < nMaxLen;
 --i,++k){ans[k] = (arr[i] + carry) % 10 + '0';carry = (arr[i] + carry) / 10;}free(arr);//释放辅助内存arr = NULL;while(carry > 0 && k < nMaxLen){ans[k++] = carry % 10 + '0';carry /= 10;}//检查是否溢出if(k > nMaxLen)return NULL;ans[k] = '\0';Invert(ans, k);return ans;}

照旧，给出辅助函数和变量定义，以及main函数的调用，不需要者可以pass~~~~
#include<stdio.h>
#include<string.h>
#include<malloc.h>
char* Add(char* const str1, char* const str2, char *ans, int nMaxLen);
char* Minus(char* const str1, char * const str2, char *ans, int nMaxLen);
//反转字符串
void Invert(char *str, int nLen)
{
	if(!str || nLen < 1)
		return;
	char temp;
	for(int i = 0; i < nLen>>1; ++i)
	{
		temp = str[i];
		str[i] = str[nLen - 1 - i];
		str[nLen - 1 - i] = temp;
	}
}
int main()
{
	const int N = 50;
	char str1[N];
	char str2[N];
	char ans[N];
	char *p;
	while(scanf("%s %s", str1, str2) != EOF)
	{
		p = Add(str1, str2, ans, N);//加
		if(p)
			puts(p);
		p = Minus(str1, str2, ans, N);//减
		if(p)
			puts(p);
		p = Mutiply(str1, str2, ans, N);//乘
		if(p)
			puts(p);	
	}
	return 0;
}