Вот моё решение :
var n,i,lik,nod,kol,t,op:longint;
a:array of longint;
b:array of longint; sim:char;
function nodob(a,b:longint):longint;
begin
while (a<>0) and (b<>0) do
if a>b then a:=a mod b else b:=b mod a;
nodob:=a+B;
end;
procedure Swap(Var t,y:longint);
var u:longint;
begin u:=t; t:=y; y:=u; end;
procedure Sort(l,r:longint);
var i,j,x:longint;
begin
x:=a[(l+r) div 2]; i:=l; j:=r;
repeat
while a[i]<x do inc(i);
while a[j]>x do dec(j);
if i<=j then begin Swap(a[i],a[j]); Swap(b[i],b[j]); inc(i); dec(j); end;
until i>j;
if l<j then Sort(l,j);
if i<r then Sort(i,r);
end;
function BIN(lik:longint):longint;
var l,r,sredn:longint;
begin
l:=1; r:=kol; sredn:=(l+r) div 2;
while r-l>1 do
begin
if a[sredn]=lik then begin Bin:=sredn; exit; end
else
if lik<a[sredn] then r:=sredn else l:=sredn;
sredn:=(l+r) div 2;
end;
if a[l]=lik then begin BIN:=l; exit; end;
if a[l+1]=lik then begin BIN:=l+1; exit; end;
BIN:=0;
end;
begin
assign(input,'input.txt');
reset(input);
assign(output,'output.txt');
rewrite(output);
readln(n);
SetLength(a,n+2);
SetLength(b,n+2);
kol:=0; NOD:=0;
for op:= 1 to n do
begin
readln(sim,lik);
case sim of
'+':begin
t:=BIN(lik);
if t=0 then begin
inc(kol);
a[kol]:=lik; b[kol]:=1;
Sort(1,kol);
if kol=1 then nod:=lik else
nod:=nodob(nod,lik);
end else
begin
inc(b[kol]);
end;
end;
'-':begin
t:=BIN(lik);
if b[t]>1 then begin
dec(b[t]);
end else
begin
a[t]:=0; b[t]:=0;
for i:= t to kol-1 do
begin
a[i]:=a[i+1];
b[i]:=b[i+1];
end;
dec(kol);
if kol=0 then begin nod:=1;
end else
begin
nod:=a[1];
for i:= 2 to kol do
begin
nod:=nodob(nod,a[i]);
if nod=1 then break;
end;
end;
end;
end;
end;
writeln(nod);
end;
close(inpuT); close(output);
end.
Почему у меня WA на 4 тесте.

Это как бы уже неправильно. Я читал так:

Проверку в каком смысле? в 5-ом тест получается всё равно количество 0/1 меньше, чем M. Если бы мусор был бы в строке - его можно откинуть, а что делать в таком случае?

Во-первых разных векторов может быть 21^2.(координаты по модулю < 10)
Для каждого такого вектора v, разобьем всех кроликов на классы эквивалентности относительно этого вектора. (отношение эквивалентности: A~B - |A-B| || v.
Отсортируем точки внутри каждого класса эквивалентности по x, потом по y. Теперь мы для запроса мы можем определить: класс эквивалентности в котором лежит точка выстрела за O(1), и сколько точек лежит на этом луче O(log n) - бинарным поиском.

Как по точке (x, y)определить класс эквивалентности относительно вектора (a, b): с = -b * x + a * y;