高精度加减乘除

一。高精度乘法（敲了好久。。头秃）

//str1,str2是存数字的字符串，返回最后的位数，c存储结果int highmult( char str1[], char str2[], int c[] ){    int len1 = strlen(str1), len2 = strlen(str2);    int a[len1 + 1] = {0}, b[len2 + 1] = {0};    int i, j;    memset(a, 0, sizeof(a));    memset(b, 0, sizeof(b));    memset(c, 0, sizeof(c));
/*
*用这个函数写1000阶乘时这句话没用了，
*不知道为啥，可把我无奈坏了，找了半天错，
*好像memset对int数组不好用，有时候出问题
*/
    //把字符变成数字，方便运算    for( i = len1 - 1, j = 0; i >= 0; i--)        a[j++] = str1[i] - '0';    for( i = len2 - 1, j = 0; i >= 0; i--)        b[j++] = str2[i] - '0';    for( i = 0; i < len1; i++)//i, j是权重，i+j就是c[]里面的权重        for( j = 0; j < len2; j++)        {            c[i + j] += a[i] * b[j];            printf("%d\n", c[i + j]);            if( c[i + j] >= 10 )//因为9 + 9 * 9 < 100，所以只进一位            {                int k = c[i + j];                c[i + j + 1] += k / 10;                c[i + j] = k % 10;            }        }    
    //找位数。。。这个大概可以在前面的代码里优化，，懒，，    for( i = len1 + len2 + 2; (c[i] == 0) && (i >= 0); i-- )        ;    return i + 1;}
二。高精度加法
这个是我还不知道高精度是啥东西时百度的算法，记不得找的哪儿的了，一直存着,稍微改了点，，，原创看见了误喷一定加上链接。
恩！原创棒棒的！
//加起来存进c
int highadd( char num1[], char num2[], int c[] ){    int lena = strlen(num1), lenb = strlen(num2), lenc, i, x;//x 是进位用的    int a[100], b[100];    memset(a, 0, sizeof(a));    memset(b, 0, sizeof(b));    memset(c, 0, sizeof(c));
    //存数字    for( i = 0; i < lena; i++)        a[lena - i - 1] = num1[i] - '0';    for( i = 0; i < lenb; i++)        b[lenb - i - 1] = num2[i] - '0';    lenc = 0;    x = 0;    while( lenc <= lena || lenc <= lenb )    {        c[lenc] = a[lenc] + b[lenc] + x;        x = c[lenc] / 10;        c[lenc] %= 10;        lenc++;    }    c[lenc] = x;    if( c[lenc] == 0 )        lenc--;    return lenc;}
三。高精度减法
//贼他妈难写，
int highsub( char a[], char b[], int c[] )
{    for(int i = 0; i < MAXLEN; ++i)        c[i] = 0;    int la = strlen(a), lb = strlen(b);    int sign(0), m(0);    int com = compare(a, b);
    
    //记录符号    if( com > 0 )    {        sign = 1;        m = la;    }    else if( com < 0)    {        sign = -1;        m = lb;    }    else        return 1;//一位以便于输出0    if(sign == 1)    {        int sub = la - lb;        for(int i = 0; i < sub; ++i)//先把前面的权重高的解决了            c[i] = a[i] - '0';        for(int i = sub; i < m; ++i)//和我写的乘法一样，算一次就处理一次进位        {            if( a[i] > b[i - sub])                c[i] = a[i] - b[i - sub];            else if( a[i] == b[i - sub] )                ;            else//小于时向上一位要一个，后面再集中处理如果上一位也是零的情况            {                --c[i - 1];                c[i] = (a[i] - '0') + 10 - (b[i - sub] - '0');            }        }    }    else//一样的    {        int sub = lb - la;        for(int i = 0; i < sub; ++i)            c[i] = b[i] - '0';        for(int i = sub; i < m; ++i)        {            if( b[i] > a[i - sub])                c[i] = b[i] - a[i - sub];            else if( a[i - sub] == b[i])                ;            else            {                --c[i - 1];                c[i] = (b[i] - '0') + 10 - (a[i - sub] - '0');            }        }    }    for(int i = m - 1; i >= 0; --i)//处理退位时产生的负数问题        if( c[i] < 0 )        {            --c[i - 1];            c[i] += 10;        }    if(c[0] == 0)//处理a， b中前几位相同以至于c前几位是0的问题        for(int i = 1; i < MAXLEN; ++i)            if( c[i] != 0 )            {                for(int j = i; j < m; ++j)                    c[j - i] = c[j];                m -= i;                break;            }    c[0] *= sign;    return m;}
四。高精度除法
这个代码对于大数不能算位数相差超过六位的，太慢了，因为是用减法实现的，循环减，减了几次答案就是几
以后灵光一闪优化了（偷偷百度一下）再修改吧
（最难写的高精度。。）
//要存小数点就用字符数组c了，n是确定小数点后几位,默认c里面都是0，懒得加memset
int highdiv( char a[], char b[], char c[], int n ){    int d[MAXLEN];    int j(0), k(0);    int zs;//记录小数点位置，便于小数部分运算
    
    //代码核心，把除降解为减    while( compare( a, b ) >= 0 )    {        int len = highsub(a, b, d);        memset(a, 0, sizeof(0));        for(int i = 0; i < len; ++i)            a[i] = d[i] + '0';        a[len] = '\0';        ++c[j];        if(c[j] >= 10)//烦人的进位问题        {            for(int i = j; i > 0; --i)                if( c[i] >= 10 )                {                    ++c[i - 1];                    c[i] %= 10;                }            if(c[0] >= 10)//进到头了就把c整体往后挪一位            {                for(int i = j; i >= 0; --i)                    c[i + 1] = c[i];                c[0] = 1;                c[1] = 0;                ++j;            }        }    }    zs = j + 2;//实际上是指向了小数点后一位
    //把减的比b小的a膨胀开来，添加n + 1个0，模拟手算的加0过程    for(int i = strlen(a); i < MAXLEN && k <= n; ++i, ++k)    {        a[i] = '0';        if( k == n )            a[i + 1] = '\0';    }    j += 2;//j指向第一个小数位，再来一波和上面差不多的操作    while( compare( a, b ) >= 0 )    {        int len = highsub(a, b, d);        memset(a, 0, sizeof(0));        for(int i = 0; i < len; ++i)            a[i] = d[i] + '0';        a[len] = '\0';        ++c[j];        if(c[j] >= 10)        {            for(int i = j; i > zs; --i)                if( c[i] >= 10 )                {                    ++c[i - 1];                    c[i] %= 10;                }            if(c[zs] >= 10)            {                for(int i = j; i >= zs; --i)                    c[i + 1] = c[i];                c[zs] = 1;                c[zs + 1] = 0;                ++j;            }        }    }    k = j;    if( c[k] >= 5 )//四舍五入，循环是怕一入就出现了10        do        {            ++c[k - 1];            c[k] = 0;        }        while( c[--k] >= 10 );    j = zs + n;    for(int i = 0; i <= j; ++i)//数字变字符，便于小数点打印        c[i] += '0';    c[zs - 1] = '.';    c[j] = '\0';    return j + 1;}
最后把compare也贴出来吧，不复杂
int compare(char a[], char b[]){    int la = strlen(a), lb = strlen(b);    if( la > lb )        return 1;    else if( la < lb )        return -1;    else        for(int i = 0; i < la; ++i)            if( a[i] > b[i] )                return 1;            else if( a[i] < b[i] )                return -1;    return 0;}

还是初学者，自己敲可能有bug，也应该还可以优化，比如处理进位写成函数什么的。
要是以后变大佬了就回来优化一波算法。