用c语言编写高精度除法

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

void high_precision_divide(char *dividend, char *divisor, char **quotient, char **remainder)
{
    // 1. 初始化
    int dividend_length = strlen(dividend);
    int divisor_length = strlen(divisor);
    int max_length = dividend_length > divisor_length ? dividend_length : divisor_length;
    *quotient = (char *)malloc(max_length + 1);
    *remainder = (char *)malloc(max_length + 1);
    int quotient_index = 0;
    int remainder_index = 0;
    int i;

    // 2. 归一化
    int padding = max_length - dividend_length;
    for (i = 0; i < padding; i++) {
        dividend[i] = '0';
    }
    padding = max_length - divisor_length;
    for (i = 0; i < padding; i++) {
        divisor[i] = '0';
    }

    // 3-4. 循环除法
    for (i = max_length - 1; i >= 0; i--) {
        // 将余数左移一位，并将下一个被除数数字添加到余数中
        int remainder_int = remainder[remainder_index] - '0';
        remainder_int *= 10;
        remainder_int += dividend[i] - '0';
        remainder[remainder_index] = remainder_int + '0';

        // 计算余数除以除数的数字
        int quotient_digit = remainder_int / (divisor[i] - '0');

        // 将商添加到商中
        (*quotient)[quotient_index] = quotient_digit + '0';
        quotient_index++;

        // 将除数乘以商的数字，从被除数中减去该乘积
        int product = quotient_digit * (divisor[i] - '0');
        int j;
        for (j = 0; j <= i; j++) {
            int dividend_int = dividend[j] - '0';
            dividend_int -= product % 10;
            dividend[j] = dividend_int + '0';
            product /= 10;
        }

        // 更新余数索引
        remainder_index++;
    }

    // 5. 删除前导零
    while ((*quotient)[0] == '0' && quotient_index > 0) {
        (*quotient)++;
        quotient_index--;
    }
    while ((*remainder)[0] == '0' && remainder_index > 0) {
        (*remainder)++;
        remainder_index--;
    }
}