大数相乘解决办法

#include "stdafx.h"
#include <iostream>
using namespace std;
// 大数相乘解决办法
int main(int argc, char* argv[])
{
    unsigned char s1[] = "987654321987654321";
    unsigned char s2[] = "987654321987654321";
    int is1 = strlen((const char *)s1);
    int is2 = strlen((const char *)s2);
    unsigned char *s = new unsigned char[is1 + is2 + 1];

    memset(s, 0, is1 + is2 + 1);
    for (int i = 0; i < is1; i++)
    {
        if ((s1[i] > '9') || (s1[i] < '0'))
        {
            return 0;
        }

        s1[i] -= '0';
    }

    for (i = 0; i < is2; i++)
    {
        if ((s2[i] > '9') || (s2[i] < '0'))
        {
            return 0;
        }

        s2[i] -= '0';
    }

    for (int k = is2 - 1; k >= 0; k--)
    {
        for (int j = is1 - 1; j >= 0; j--)
        {
            int t = s2[k] * s1[j];
            s[k + j + 1] += t % 10;

            if (k + j >= 0)
            {
                s[k + j] += t / 10;
                s[k + j] += s[k + j + 1] / 10;
                s[k + j + 1] = s[k + j + 1] % 10;
            }
        }
    }

    for (i = 0; i < (is1 + is2); i++)
    {
        s[i] += '0';
    }

    cout << "s= " << s << endl;
    delete [] s;
    double d1 = 987654321987654321.0;
    double d2 = 987654321987654321.0;
    cout << d1 * d2 << endl;
    return 0;
}