pascal高精度除法程序

<code class="pascal">program HighPrecisionDivision;

{$APPTYPE CONSOLE}

var
  dividend, divisor, quotient, remainder: array[0..1000] of integer;
  dividend_length, divisor_length, quotient_length, remainder_length: integer;
  i, j, k: integer;

begin
  // 初始化变量
  dividend_length := 0;
  divisor_length := 0;
  quotient_length := 0;
  remainder_length := 0;
  for i := 0 to 1000 do
  begin
    dividend[i] := 0;
    divisor[i] := 0;
    quotient[i] := 0;
    remainder[i] := 0;
  end;

  // 输入被除数
  write('Enter dividend: ');
  readln(dividend_length);
  for i := 1 to dividend_length do
  begin
    write('Dividend digit ', i, ': ');
    readln(dividend[i]);
  end;

  // 输入除数
  write('Enter divisor: ');
  readln(divisor_length);
  for i := 1 to divisor_length do
  begin
    write('Divisor digit ', i, ': ');
    readln(divisor[i]);
  end;

  // 执行除法
  for i := dividend_length downto 1 do
  begin
    remainder[i] := remainder[i] * 10 + dividend[i];
    quotient[i] := remainder[i] div divisor[divisor_length];
    remainder[i] := remainder[i] mod divisor[divisor_length];
    for j := divisor_length downto 1 do
    begin
      dividend[i-j+divisor_length] := dividend[i-j+divisor_length] - divisor[j] * quotient[i];
    end;
  end;

  // 输出商
  quotient_length := dividend_length - divisor_length + 1;
  write('Quotient: ');
  for i := quotient_length downto 1 do
    write(quotient[i]);
  writeln;

  // 输出余数
  remainder_length := divisor_length - 1;
  write('Remainder: ');
  for i := remainder_length downto 1 do
    write(remainder[i]);
  writeln;
end.</code>