大数相乘 - 整型

 

搞程序的是不是都熬夜的呢?我是的,不过很快就会努力不是了.为了健康,做一个不熬夜的程序员
本来打算了5:00睡的,鬼使神差莫明其妙地放下了PHP写了点C++,于是就有了以下代码的产生,也许是太久没有写C++了.

关于算法我已经写在程序的注释里了.简言之是先相加(不管进位),最后一起进位.

//  $Id: multi.cpp 6 2006-11-01 21:12:48Z JiangMiao $ 
//  JiangMiao's Blog  http://blog.csdn.net/antter
#include  " stdafx.h "
#include  < iostream >
#include  < string >
using   namespace  std;


///  转换字符串为
char *  strToInt( const   string &   in ) 
    {
     char *  rt = new   char [ in .size()];
     for (size_t i = 0 ;i < in .size();i ++ ) 
        {
        rt[i] = in [i] - ' 0 ' ;
        }
     return  rt;
    }

/* *
 * 大数相乘,最大位数 0.04G 即32位int/(9*9)
  */
string  mul( const   string &  sa, const   string &  sb) 
    {
     string  rt;
    size_t la = sa.size();
    size_t lb = sb.size();
    size_t lr = la + lb;  // 最大可能位数
     char *  a = strToInt(sa);
     char *  b = strToInt(sb);
    unsigned   int *  r = new  unsigned  int [lr]; // 分配结果空间
    size_t ia,ib,ir;
    rt.reserve(lr);
    memset(r, 0 ,lr * sizeof ( int ));
     for (ia = 0 ;ia < la;ia ++ )  // 相乘
        {
         for (ib = 0 ;ib < lb;ib ++ )
            {
            ir = ib + ia + 1 ;
            r[ir] += a[ia] * b[ib];
            }
        }
     for (ir = lr - 1 ;ir > 0 ;ir -- )  // 进位
        {
         int  add = r[ir] / 10 ;
        r[ir - 1 ] += add;
        r[ir] %= 10 ;
        }
     for (ir = 0 ;ir < lr;ir ++ )  // 除去前面多余的0
        {
         if (r[ir] != 0 )
             break ;
        }
     for (;ir < lr;ir ++ )  // 生成字串
        {
        rt.append( 1 ,r[ir] + ' 0 ' );
        }
    delete a;
    delete b;
    delete r;
     return  rt;
    }

/*
提升性能的改进:
1. 使用char* 代替string
2. 多数相乘返回非string,最后由intToStr进行字串输出,可极大地节省预处理和生成的时间

实现以上两点性能提升至少45%,乱猜的 :)
*/


///  以下为测试文件
#include  " windows.h "
int  main( int  argc, char **  argv)
    {
     string  rt;
    LARGE_INTEGER fre,begin,end;
    QueryPerformanceFrequency( & fre);
    QueryPerformanceCounter( & begin);

     // 10000阶乘测试
    rt = " 1 " ;
     for ( int  i = 1 ;i <= 10000 ;i ++ )
        {
         char  b[ 6 ];
        _itoa_s(i,b, 10 );
        rt = mul(rt,b);
        }

    QueryPerformanceCounter( & end);
    cout << " Spend  " << (end.QuadPart - begin.QuadPart) * 1000000 / fre.QuadPart << " ns " << endl;
    cout << rt.size() << endl;
     // cout<<rt<<endl;
     return   0 ;
    }


/*  测试10000!的结果
 * C4 2.66MHZ
Spend 16911238ns //17s
35660.
  */

 

欢迎探讨.我的Blog是 http://blog.csdn.net/antter