如何用Pascal找2-1000中的亲密数对

var
i,a,he,b:longint;
begin
for a:=2 to 1000 do
begin
he:=1;
for i:=2 to trunc(sqrt(a)) do
if a mod i=0 then
he:=he+i+a div i;{he 是a的因数和}
b:=1;
for i:=2 to trunc(sqrt(he)) do
if he mod i=0 then
b:=b+i+he div i;{b 是he的因数和}
if (a=b)and(a<he){防止重复} then
writeln(a:6,he:6);
end;
end.
答案是220和284