Проблемы в Pascal (для всех компиляторов!) Опции

[ONi++]

Программу одну составил по олимпиадным задачам. Суть задачи состоит в том, что нужно узнать по данным координатам, находится ли точка в треугольнике. Я сравниваю площадь треугольника и площадь четырехугольника (получившегося из соединения линиями треугольника с точкой) с помощью функции round(). Но, когда она [функция] округляет 4.5 два раза, то сначала получается 5, а потом 4. Что делать? Не хотелось бы менять весь алгоритм из-за какой-то ошибки в компиляторах (кстати, я проверял и на PascalABC, и на TurboPascal, и на FreePascal).

[andriano]

Это распространенная болезнь начинающих - в собственных ошибках винить компилятор.
Вероятнее всего, что во втором случае у тебя вместо 4.5 что-то вроде 4.49999999999.

Кстати, подозреваю, что и реализованный алгоритм также не дает правильного ответа на поставленный вопрос.
Зачем, кстати, ты округляешь?

[ONi++]

знаешь формулу для нахождения площади фигуры по координатам? вот она: s=s+sqrt(p*(p-a)*(p-b)*(p-c), где
p=(a+b+c)/2, a=sqrt((x[1]-x[i])^2+(y[1]-y[i])^2),b=sqrt((x[i]-x[i+1])^2+(y[i]-y[i+1])^2),c=sqrt((x[i+1]-x[1])^2+(y[i+1]-y[1])^2), 2<=i<=n-1. Так мы находим площадь треугольника, а чтобы найти площадь фигуры "треугольник+точка", просто меняем x[1] и y[1] на xx и yy (координаты точки), x[n+1]=x[1],y[n+1]=y[1],1<=i<=n. Получается два вещественных числа, которые надо округлить, потому что без округления числа не равны в любом случае (кстати, почему???), округляем. На первых тестах все проходит, но следущих тестах одинаковые площади двух фигур (допустим, 4.5 и 4.5) округляются по разному - 5 и 4. А без округления они не равны. В чем же глюк??

[volvo]

Ты код приведи, который показывает ошибку (вот именно так: чтобы программа при сравнении одинаковых - по твоему - площадей выдавала сообщение, что площади не равны), а не разливай тут теоретическую воду, договорились? А то я тоже много чего могу сказать...

[ONi++]

var
input,output:text;
x,y:array[1..4]of integer;
i,n,xx,yy:integer;
a,b,c,p,s1,s2:real;
begin
  assign(input,'tria-pt.in');
  assign(output,'tria-pt.out');
  reset(input);
  rewrite(output);
  n:=3;
  i:=1;
  while i<=n do
	 begin
		readln(input,x[i],y[i]);
		i:=i+1;
	 end;
  readln(input,xx,yy);
  a:=sqrt(sqr(x[1]-x[2])+sqr(y[1]-y[2]));
  b:=sqrt(sqr(x[2]-x[3])+sqr(y[2]-y[3]));
  c:=sqrt(sqr(x[3]-x[1])+sqr(y[3]-y[1]));
  p:=(a+b+c)/2;
  s1:=sqrt(p*(p-a)*(p-b)*(p-c));
  x[n+1]:=x[1];
  y[n+1]:=y[1];
  i:=1;
  while i<=n do
	 begin
		a:=sqrt(sqr(xx-x[i])+sqr(yy-y[i]));
		b:=sqrt(sqr(x[i]-x[i+1])+sqr(y[i]-y[i+1]));
		c:=sqrt(sqr(x[i+1]-xx)+sqr(y[i+1]-yy));
		p:=(a+b+c)/2;
		s2:=s2+sqrt(p*(p-a)*(p-b)*(p-c));
		i:=i+1;
	 end;
  if round(s1)=round(s2) then
	 writeln(output,'In')
  else
	 writeln(output,'Out');
  close(input);
  close(output);
end.

[volvo]

Ну, и что? Запустил программу под FPC, треугольник: (1, 1), (4, 8), (6, 4) и точка (5, 6). S1 = S2, выдало ответ "In", что я делаю не так? Значит, проблема не в компиляторе? С какими значениями у тебя неверный ответ?

P.S. Прочитай уже, что тебе написали выше: выброси Round и сравнивай s1 и s2 как и положено сравнивать вещественные числа, с какой-то погрешностью: Abs(s1 - s2) < eps

[ONi++]

-2 -2
3 1
0 1
0 1

[volvo]

-2 -2
3 1
0 1
0 1

И с этими значениями проходит равенство, и выдает In (и FPC 2.2.4, и Турбо-Паскаль 7 со включенным сопроцессором. Без сопроцессора Round работает чуть по-другому и результат иной)...

[ONi++]

Спасибо, действительно, разница помогла. Вот код:

program triapt;
var
input,output:text;
x,y:array[1..4]of integer;
i,n,xx,yy:integer;
a,b,c,p,s1,s2:real;
begin
  assign(input,'tria-pt.in');
  assign(output,'tria-pt.out');
  reset(input);
  rewrite(output);
  n:=3;
  i:=1;
  while i<=n do
	 begin
		readln(input,x[i],y[i]);
		i:=i+1;
	 end;
  readln(input,xx,yy);
  a:=sqrt(sqr(x[1]-x[2])+sqr(y[1]-y[2]));
  b:=sqrt(sqr(x[2]-x[3])+sqr(y[2]-y[3]));
  c:=sqrt(sqr(x[3]-x[1])+sqr(y[3]-y[1]));
  p:=(a+b+c)/2;
  s1:=sqrt(p*(p-a)*(p-b)*(p-c));
  x[n+1]:=x[1];
  y[n+1]:=y[1];
  i:=1;
  while i<=n do
	 begin
		a:=sqrt(sqr(xx-x[i])+sqr(yy-y[i]));
		b:=sqrt(sqr(x[i]-x[i+1])+sqr(y[i]-y[i+1]));
		c:=sqrt(sqr(x[i+1]-xx)+sqr(y[i+1]-yy));
		p:=(a+b+c)/2;
		s2:=s2+sqrt(p*(p-a)*(p-b)*(p-c));
		i:=i+1;
	 end;
  if round(abs(s1-s2))=0 then //вот здесь я изменил - и все работает!
	 writeln(output,'In')
  else
	 writeln(output,'Out');
  close(input);
  close(output);
end.

Сообщение отредактировано: ONi++ - 11.11.2009 17:54