zju 1041 雷达覆盖
Description
以雷达心为圆心的半圆形雷达覆盖范围有多个点 雷达可旋转,求最多覆盖数(含在边界的)
Input
Output
题解
先把雷达不可能覆盖的点排除,每个点与中点都组成一个直线方程完后线性规划看其他的点在这条直线左边的inc(l)在右边的inc(r)在直线的inc(m),完后比较m+l与max的大小和r+m与max的大小。这题就完了。
代码
var
xn,yn,n,num:longint;
rn:real;
x,y:array [0..1001] of longint;
procedure init;
var
i,xx,yy,nn:longint;
begin
n:=0;
readln(xn,yn,rn);
if rn<=0 then halt;
readln(nn);
for i:=1 to nn do
begin
readln(xx,yy);
if sqrt(sqr(xx-xn)+sqr(yy-yn))<=rn then
begin
inc(n);
x[n]:=xx; y[n]:=yy;
end;
end;
end;
function max(o,p:longint):longint;
begin
if o>p then exit(o);
exit(p);
end;
procedure main;
var
i,j,r,l,m:longint;
begin
num:=0;
for i:=1 to n do
begin
r:=0; l:=0; m:=0;
for j:=1 to n do
begin
if (x[i]-xn)*(y[j]-yn)-(x[j]-xn)*(y[i]-yn)=0 then inc(m);
if (x[i]-xn)*(y[j]-yn)-(x[j]-xn)*(y[i]-yn)<0 then inc(r);
if (x[i]-xn)*(y[j]-yn)-(x[j]-xn)*(y[i]-yn)>0 then inc(l);
end;
num:=max(max(num,m+r),m+l);
end;
end;
begin
while 1=1 do
begin
init;
main;
writeln(num);
end;
end.