厦大C语言上机 1388 高精度乘法

厦大C语言上机 1388 高精度乘法

1388.高精度乘法


时间限制: 1000 MS          内存限制: 65536 K
        
提交数: 564 (0 users)          通过数: 237 (226 users)


问题描述
之前做过了一道关于高精度加法的，想必难不到大家，那么请大家现在尝试一下高精度乘法吧！


输入格式
输入只包含一行，包括两个整数A和B，A、B用空格隔开。(0<=A,B<=10^1000）


输出格式
输出一行，A和B的乘积。


样例输入
111111111111 111111111111


样例输出
12345679012320987654321


来源

xmu

#include <stdio.h>
#include <string.h>void reverse(char *s, int len)
{int i, j;char temp;for (i = 0, j = len-1; i < j; ++i, --j){temp = s[i];s[i] = s[j];s[j] = temp;}
}void multiply(char *a, int len_a, char *b, int len_b, char *product)
{int ptr_a, ptr_b, ptr_int_prod, ptr_prod, len_prod, carry, temp;int int_prod[2005] = { 0 };reverse(a, len_a);reverse(b, len_b);ptr_int_prod = 0;for (ptr_b = 0; ptr_b < len_b; ++ptr_b){ptr_int_prod = ptr_b;for (ptr_a = 0; ptr_a < len_a; ++ptr_a, ++ptr_int_prod)int_prod[ptr_int_prod] += (a[ptr_a] - '0') * (b[ptr_b] - '0');}len_prod = ptr_int_prod;carry = 0;for (ptr_int_prod = 0; ptr_int_prod < len_prod; ++ptr_int_prod){temp = (int_prod[ptr_int_prod] + carry) / 10;int_prod[ptr_int_prod] = (int_prod[ptr_int_prod] + carry) % 10;carry = temp;}while (carry > 0){temp = carry / 10;int_prod[ptr_int_prod] = carry % 10;carry = temp;ptr_int_prod++;len_prod++;}ptr_int_prod = len_prod - 1;while (ptr_int_prod >= 0 && int_prod[ptr_int_prod] == 0){ptr_int_prod--;len_prod--;}if (ptr_int_prod < 0){product[0] = '0';product[1] = '';return;}ptr_prod = 0;while (ptr_int_prod >= 0)product[ptr_prod++] = (char)(int_prod[ptr_int_prod--] + '0');product[ptr_prod] = '';
}int main()
{char a[1005] = { 0 };char b[1005] = { 0 };char product[2005] = { 0 };int len_a, len_b;scanf("%s %s", a, b);len_a = (int)strlen(a);len_b = (int)strlen(b);multiply(a, len_a, b, len_b, product);printf("%sn", product);return 0;
}