大數相乘算法

啥也說，直接給代碼

//////////////////////////////////////////////////////////////// //作者：塗遠東 //時間：2010 01 23 //高精度的乘法,乘數的位數沒有限制，小數點後面的位數也沒有限制 //輸入參數：乘數 str1，str2 //輸出參數：str = str1*str2 //返回值類型：string& //////////////////////////////////////////////////////////////// string& large_mult(string str1,string str2,string& str/*輸出參數*/) { str = "0"; int pos1 = str1.find('.'); int pos2 = str2.find('.'); //小數點在字符串的位置 int offset = str1.length() + str2.length() - ((pos1 >0)? (pos1+1):0 ) - ((pos2 >0) ? (pos2+1):0 ); //去除字符串str1和str2中的小數點 if(pos1 != -1) str1.erase(pos1,1); if(pos2 != -1) str2.erase(pos2,1); unsigned char single = 0; //個位數 unsigned char tens = 0; //十位數 unsigned char ret = 0; //乘積 //計算規則，假設str1和str2字符串長度分別爲len1和len2 //乘積str的字符串長度len應該不大於len1+len2 //乘數str1[len1-i]和str2[len2-j]相乘只會乘積中str中str[len-i-j]和str[len-i-j-1]的值 int i = 0; for(string::iterator iter = str1.end()-1; iter >= str1.begin(); iter--,i++) { int index = 0; int offset = 0; int j=0; for(string::iterator iter2=str2.end()-1; iter2 >= str2.begin(); iter2--,j++) { offset = i+j +1; //字符‘0’的ASCII編碼值爲48 ret = (*iter2 -'0') * (*iter -'0'); //乘積 single = ret % 10; //乘積的個位數 tens = ret / 10; //乘積的十位數 //位數不夠左邊補0 index = str.length() - offset; while(index < 1) { str.insert(str.begin(), '0'); index = str.length() - offset; } assert(str.length() > offset); index = str.length() - offset; add(str,index,single); //個位數相加 add(str,index-1,tens); //十位數相加 } } if(str[0] == '0') { str.erase(str.begin()); } if(offset < str.length()) { string::iterator itr = str.end()-offset; str.insert(itr,'.'); } return str; } //////////////////////////////////////////////////////////////// //在大數str的index位處加於num 即 str+= num //////////////////////////////////////////////////////////////// inline bool add(string& str, int index, int num) { assert(index >=0); int p = str[index] + num - '0'; str[index] = p%10 + '0'; if(p >9) add(str,index-1,1); return true; }

測試代碼

#include "stdafx.h" #include <iostream> #include <string> #include <cassert> using namespace std; string& large_mult(string str1,string str2,string& str/*輸出參數*/); bool add(string& str, int index, int num); int main(int argc, char* argv[]) { printf("Please input two large number!/n"); string str1; string str2; std::cout<<"Please input the first large number:/nstr1="; std::cin>> str1; std::cout<<"Please input the second large number::/nstr2="; std::cin>>str2; string str; //str = str1 * str2 large_mult(str1,str2,str); std::cout<<"str=str1*str2/nstr="; std::cout<<str<<std::endl; return 0; }