Прямая проходящая через множество окружностей

program 9;
const n=5;k=3;  {n-Кол-во точек, k-кругов}
type
 TPoint=record
 x,y:real
 end;
 Tcircle=record
 x,y,r:real
 end;
 TPointsSet=array[1..n] of tpoint;
 TcircleSet=array[1..k] of tcircle;
 const
 pointsset:tpointsset=((x:0;y:0),(x:4;y:2),(x:0;y:2),(x:1;y:1),(x:2;y:4));
 circleset:tcircleset=((x:1;y:1;r:4),(x:3;y:2;r:1),(x:0;y:-2;r:3));
 
 
 begin
 end.

type
 TPoint=record
 x,y:real
 end;
_point=record
A:Tpoint;
B:Tpoint;
end;
krug=record
C:Tpoint;
R:integer;
end;
 
const
n=5;
k=3;
kr:array[1..k] of krug=((c:(x:1;y:1);r:4),(c:(x:3;y:2);r:1),(c:(x:0;y:-2);r:3));
ts:array[1..n]of tpoint=((x:0;y:0),(x:4;y:2),(x:0;y:2),(x:1;y:1),(x:2;y:4));
 
 
function Dotr(A1, B1: Tpoint): real;
begin
Dotr:=sqrt(sqr(b1.X-a1.X)+sqr(b1.Y-a1.Y));
end;
 
function Strg(A1, B1, C1: Tpoint): real;
begin
Strg:=abs(((a1.X-c1.X)*(B1.y-c1.Y)-(b1.X-c1.X)*(a1.Y-c1.Y))/2);
end;
 
 
 
procedure rrr(var p:array of _point;m:integer);
var i,j:byte; h:real;q,ch,it:integer;
begin
q:=0;h:=0;it:=0;ch:=0;
for i:=1 to m-1 do
 begin
  for j:=1 to k do
  begin
  h:=(Strg(p[i].A,p[i].B,kr[j].C)*2/Dotr(p[i].a,p[i].B));
 
  if (h<kr[j].R) then ch:=ch+1;
  end;
 
if ch>q then begin
q:=ch;
it:=i;
end;
ch:=0;
end;
 
writeln('number of intersections: ', q);
writeln('  a:    |   b: ');
writeln(p[it].a.x:3:1,',', p[it].a.y:3:1,' | ',p[it].b.x:3:1,',', p[it].b.y:3:1);
end;
 
var
  f,g,m:integer;
  p:array[1..100]of _point;
  begin
  m:=0;
  for f:=1 to n-1 do
  for g:=f+1 to n do
  begin
  m:=m+1;
  p[m].a.x:=ts[f].x;
  p[m].a.y:=ts[f].y;
  p[m].b.x:=ts[g].x;
  p[m].b.y:=ts[g].y;
  end;
 
writeln('Entry Points:');
for f:=1 to n do
writeln(ts[f].x:3:1,',', ts[f].y:3:1);
writeln('Entry Circles:');
writeln('-----------------');
writeln('x    y  |  r');
for f:=1 to k do
writeln(kr[f].c.x:3:1,',', kr[f].c.y:3:1,' | ',kr[f].r:3);
  rrr(p,m);
 
end.