常量系数多项式除法(Delphi版)

最近研究二维码时发现有的地方需要进行多项式除法，于是自己写了一个算法，以活动一下脑子。

//初始化数组，赋值0.

procedure InitArray(var bb : array of double);
var
a : integer;
begin
for a := 0 to high(bb) do
begin
bb[a] := 0;
end;
end;

//多项式计算.

procedure CountMultiExpress(var dXiShu1,dXiShu2,dXiShu3,dXiShu4 : array of double;AMaxMiCi:integer);
var
//多项式表达规则：认为输入的都是标准的表达式f(x)=a0*x^0+a1*x^1+a2*x^2+...+an*x^n
//系数从0-20存a0-a20,幂为对应的数组索引
//dXiShu1为分母,dXishu2为分子, 3为商 4为余数。
//这里的计算仅仅支持常数系数的。
dXiShuTmp : array of double;
dTmp1,dTmp2,dMaxMiCi,dMaxXiShu : double;
i,j,iTmpMiCi,iMaxMiCi:integer;

procedure InitArray(var bb : array of double);
var
a : integer;
begin
for a := 0 to high(bb) do
begin
bb[a] := 0;
end;
end;
begin
Setlength(dXiShuTmp,AMaxMiCi);
//获取多项式分子最大幂次和系数。
for i:= 20 downto 0 do
begin
if dXiShu2[i]<>0.0 then
begin
dMaxMiCi := i;
dMaxXiShu := dXiShu2[i];
break;
end;
end;
//计算，从高到低
iMaxMiCi := trunc(dMaxMiCi);
for i:= AMaxMiCi downto iMaxMiCi do
begin
if dXiShu1[i]<> 0.0 then
begin
dTmp1 := dXiShu1[i] / dMaxXiShu;
iTmpMiCi := i - iMaxMiCi;
dXiShu3[iTmpMiCi] := dTmp1;//商.
InitArray(dXiShuTmp);
//分子乘dTmp1
for j := iMaxMiCi downto 0 do
begin
dXishuTmp[j+iTmpMiCi] := dXiShu2[j] * dTmp1;
end;
//做减法分母-dXishuTmp
for j := i downto 0 do
begin
dXishu1[j] := dXishu1[j]-dXishuTmp[j];
end;
end;
end;
//余数
for j := 0 to AMaxMiCi do
begin
dXiShu4[j] := dXiShu1[j];
end;
Setlength(dXiShuTmp,0);
end;

使用：

procedure TForm1.Button5Click(Sender: TObject);
var
dXiShu1,dXiShu2,dXiShu3,dXiShu4 : array[0..20] of double;
strOut1,strOut2:string;
j : integer;
begin
//初始化
InitArray(dXiShu1);
InitArray(dXiShu2);
InitArray(dXiShu3);
InitArray(dXiShu4);
//分母，分子赋值
dXiShu1[0]:= 7;
dXiShu1[1]:= 4;
dXiShu1[2]:= 6;
dXiShu1[3]:= 0;
dXiShu1[4]:= 0;
dXiShu1[5]:= 0;

dXiShu2[0]:= 3;
dXiShu2[1]:= 2;
dXiShu2[2]:= 0;
dXiShu2[3]:= 0;
dXiShu2[4]:= 0;
dXiShu2[5]:= 0;
CountMultiExpress(dXiShu1,dXiShu2,dXiShu3,dXiShu4,20);
//输出结果
strOut1 := '';
strOut2 := '';
for j := 0 to 20 do
begin
//商
if dXiShu3[j] > 0.0 then
begin
strOut2 := ' +'+floattostr(dXiShu3[j])+'x^'+inttostr(j)+' ' + strOut2 ;
end;
if dXiShu3[j] < 0.0 then
begin
strOut2 := ' '+floattostr(dXiShu3[j])+'x^'+inttostr(j)+' ' + strOut2 ;
end;
//余数
if dXiShu4[j] > 0.0 then
begin
strOut1 := ' +'+floattostr(dXiShu4[j])+'x^'+inttostr(j)+' ' + strOut1 ;
end;
if dXiShu4[j] < 0.0 then
begin
strOut1 := ' '+floattostr(dXiShu4[j])+'x^'+inttostr(j)+' ' + strOut1 ;
end;
end;
strOut1 := trim(strOut1);
strOut2 := trim(strOut2);
if strOut1[1]='+' then
strOut1 := copy(strOut1,2,length(strOut1)-1);
if strOut2[1]='+' then
strOut2 := copy(strOut2,2,length(strOut2)-1);
memo3.Lines.Clear();
memo3.Lines.Append('商数：'+ strOut2);
memo3.Lines.Append('余数：'+ strOut1);

end;