发段C程序，练练手：用64位保存2个无符号整数的乘法结果

xdsupermanli

闲来无事，写了段小C代码，贴上来玩玩，哈哈哈～
/*****************************************************************
* 实现两个unsigned int 类型的整数相乘，结果保存在两个unsigned int变量
* 中。其基本思想是把unsigned int分拆为两个unsigned short分别相乘，就和
* 做2位十进制数的乘法一样一样的，只不过一个是逢10进位，一个是逢0x10000进位
*****************************************************************/

typedef unsigned long u32;

void u32_multiply(u32 a, u32 b, u32 *result_hi, u32 *result_lw)
{
        u32 result[4];
        u32 g; //进位
        u32 t; //余数
        int len;

        //b的低16位与a相乘，结果保存在result数组中
        memset(result, 0, sizeof(result));
        result[0] = (a & 0xffff) * (b & 0xffff);
        g = result[0] >> 16;
        t = result[0] & 0xffff;
        result[0] = t;       
        len = 1;

        while(g)
        {
                result[len++] += g & 0xffff;
                g >>= 16;
        }

        result[1] += (a >> 16) * (b & 0xffff);
        g = result[1] >> 16;
        t = result[1] & 0xffff;
        result[1] = t;
        len = 2;

        while(g)
        {
                result[len++] += g & 0xffff;
                g >>= 16;
        }
       
        //b的高16位与a相乘，结果累加在result数组中
        result[1] += (a & 0xffff) * (b >> 16);
        g = result[1] >> 16;
        t = result[1] & 0xffff;
        result[1] = t;
       
        len = 2;
        while(g)
        {
                result[len++] += g & 0xffff;
                g >>= 16;
        }

        result[2] += (a >> 16) * (b >> 16);
        g = result[2] >> 16;
        t = result[2] & 0xffff;
        result[2] = t;

        len = 3;
        while(g)
        {
                result[len++] += g & 0xffff;
                g >>= 16;
        }

        *result_lw = ((result[1] & 0xffff) << 16) | (result[0] & 0xffff);
        *result_hi = ((result[3] & 0xffff) << 16) | (result[2] & 0xffff);

}