Версия для печати темы

Нажмите сюда для просмотра этой темы в обычном формате

Форум «Всё о Паскале» _ Задачи _ Проблемы в Pascal (для всех компиляторов!)

Автор: ONi++ 10.11.2009 21:11

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

Автор: andriano 10.11.2009 21:22

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

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

Автор: ONi++ 10.11.2009 21:46

знаешь формулу для нахождения площади фигуры по координатам? вот она: 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 10.11.2009 22:08

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

Автор: ONi++ 11.11.2009 17:21


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 11.11.2009 17:34

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

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

Автор: ONi++ 11.11.2009 17:44

-2 -2
3 1
0 1
0 1

Автор: volvo 11.11.2009 17:52

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

Автор: ONi++ 11.11.2009 17:54

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


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.