POJ 1328

 

题目大意：求在x轴上放一些半径为d的圆，问至少多少个能将n个点覆盖，多组询问。

解：原来想着：线段树+heap，可是看着这一个个0ms，不敢切，以后有时间切切看虐自己。正解是先把可行圆心区间求出来，然后求最多的不互斥区间（因为区间相斥的话，说明有共点，于是不用分配多一个雷达了），证明是首先如果当前l>rrr,inc(s), 若有l<rrr and r<rrr则新预设雷达于r，rrr:=r;

好像有种隐约敢，说不出来，好难受，还是以后想透彻了再补

本题的trick点也很多，例如，d<0, y<0，注意y=0竟然是合法的！！wa了两次，囧

//RadarInstallation
const
        maxn=1111;
        inf='1.txt';
        def=1e-10;
type
        data=record
          l, r: double;
        end;
var
        line: array[0..maxn]of data;
        test, n, d, ans: longint;
procedure qsort(b, e: longint);
var
        i, j: longint;
        x: double;
        k: data;
begin
  i := b; j := e; x := line[(i+j)>>1].l;
  repeat
    while line[i].l<x do inc(i);
    while line[j].l>x do dec(j);
    if i<=j then begin
      k := line[i]; line[i] := line[j]; line[j] := k;
      inc(i); dec(j);
    end;
  until i>j;
  if j>b then qsort(b, j);
  if i<e then qsort(i, e);
end;

procedure init;
var
        i, x, y: longint;
        tmp: double;
begin
  ans := 0;
  for i := 1 to n do begin
    readln(x, y);
    with line[i] do begin
      if (y<0)or(y>d) then begin ans := -1; continue; end;
      if y=d then tmp := 0
      else tmp := sqrt(sqr(d)-sqr(y));
      l := x - tmp; r := x + tmp;
    end;
  end;
  if ans=-1 then exit;
  qsort(1, n);
end;

procedure main;
var
        rrr: double;
        i: longint;
begin
  rrr := -maxlongint;
  for i := 1 to n do with line[i] do begin
    if l{+def}>rrr then begin
      inc(ans);
      rrr := r;
    end
    else begin
      if r<rrr then rrr := r;
    end;
  end;
end;

procedure print;
begin
  writeln('Case ',test,': ', ans);
end;

begin
  assign(input,inf); reset(input);
  readln(n, d);
  test := 0;
  while (n<>0)or(d<>0) do begin
    inc(test);
    init;
    if ans=0 then main;
    print;
    readln;
    readln(n, d);
  end;
end.