模板：高精度浮点数相乘

 1 #include <string>

 2 #include <vector>

 3 #include <algorithm>

 4 

 5 string mul(string str1,string str2)

 6 {

 7     vector<int> v_res(str1.size()+str2.size(),0);

 8     string::size_type i,j;

 9     vector<int>::size_type k,p;

10 

11     reverse(str1.begin(),str1.end());

12     reverse(str2.begin(),str2.end());

13     for(i = 0; i != str1.size(); ++i)

14     {

15         for(j = 0; j != str2.size(); ++j)

16         {

17             v_res[i+j] += (str1[i]-'0') * (str2[j] - '0');

18         }

19     }

20     for(k = 0; k != v_res.size() - 1; ++k)

21     {

22         v_res[k+1] += v_res[k] / 10;

23         v_res[k] = v_res[k] % 10;

24     }

25 

26     for(p = v_res.size() - 1; p != -1; --p)

27     {

28         if(v_res[p] != 0) break;

29     }

30     if(p == -1) p = 0;

31 

32     string s_res(p+1,'0');

33     for(k = p; k != -1; --k) s_res[p-k] = char(v_res[k] + '0');

34     

35 

36     return s_res;

37 

38 }

39 

40 string real_mul(string str1,string str2)

41 {

42     string::size_type idx_str1_point = str1.find(".");

43     if(idx_str1_point == string::npos)

44     {

45         str1 += ".0";

46         idx_str1_point = str1.find(".");

47     }

48     str1.erase(idx_str1_point,1);

49 

50     string::size_type idx_str2_point = str2.find(".");

51     if(idx_str2_point == string::npos)

52     {

53         str2 += ".0";

54         idx_str2_point = str2.find(".");

55     }

56     str2.erase(idx_str2_point,1);

57 

58     string::size_type dec_res_len = (str1.size() - idx_str1_point) + (str2.size() - idx_str2_point);

59 

60     string res = mul(str1,str2);

61 

62     if(res.size() < dec_res_len + 1)

63     {

64         res = string(dec_res_len + 1 - res.size(),'0') + res;

65     }

66 

67     res.insert(res.size() - dec_res_len,".");

68 

69     string::size_type idx_res_tail = res.find_last_not_of("0");

70     res = res.substr(0,idx_res_tail+1);

71 

72     if(res[res.size() - 1] == '.') res.erase(res.size()-1);

73 

74     return res;

75 }