Версия для печати темы
Нажмите сюда для просмотра этой темы в оригинальном формате
Форум программистов > Центр помощи > [Pascal] Определение кратчайшего расстояния


Автор: Gigabyte 15.10.2006, 15:36
с клавиатуры вводяться координаты точки и радиус
надо определить кратчайшее растояние от точки до графика 

Автор: Alexeis 15.10.2006, 16:22
Код

program Project2;

{$N+}

var
  x0, y0, R : Double;

begin
  Write(' x0='); read(x0); writeln;
  Write(' y0='); read(y0); writeln;
  Write(' R=');   read(R);  writeln;
  Write(' L=', Sqrt((x0 * x0 + y0 * y0) * (1 - R / Sqrt(x0 * x0 + y0 * y0))) : 5 : 5);
  Readln;
  Readln;
end.

Автор: Kuvaldis 15.10.2006, 20:15
Gigabyte,
Если решить, что автор задачи не опечатался, когда рисовал окружность без ее части во второй четверти, то вариант решения такой
program Project2;
 
Код

var
  x0, y0, R, L : Double;
begin
  Write(' x0 = '); read(x0); writeln;
  Write(' y0 = '); read(y0); writeln;
  Write(' R = ');  read(R);  writeln;

  // if 2 quater
  if (x0 * y0 < 0) and (x0 < 0) then
  begin
     if (y0 > -x0) then
        L := Sqrt(x0 * x0 + (y0 - R) * (y0 - R))
     else
        L := Sqrt(y0 * y0 + (x0 - R) * (x0 - R))
  end
  else
     L := Abs( Sqrt(x0 * x0 + y0 * y0) - R);

  Writeln(' L=', L : 5 : 5);
  Readln;
  Readln;
end.

Автор: Gigabyte 15.10.2006, 20:41
во 2 четверти график представляет собой 2 отрезка проходящие из точки (0,0)
по осям до окружности

Автор: Kuvaldis 16.10.2006, 00:54
Gigabyte, 
модифицированный вариант
Код

program Project2;
{$N+}
var
  x0, y0, R, L : Double;
  
//*****************************************************************************
function Min(a, b : double) : double;
begin
    if (a <= b) then
       Min := a
    else
       Min := b;
end;
//*****************************************************************************
procedure Process_2_Quater;
begin
    if (y0 > -x0) then
    begin
        if (y0 <= R) then
          L := -x0
        else
          L := Sqrt(x0 * x0 + (y0 - R) * (y0 - R))
    end
    else
    begin
        if (-x0 <= R) then
          L := y0
        else
          L := Sqrt(y0 * y0 + (x0 - R) * (x0 - R))
    end;
end;
//*****************************************************************************
procedure Process_3_Quater;
var
    len : double;
begin
    len := sqrt(x0 * x0 + y0 * y0);

    if (len >= R) then     // if outside of circle
        L := len - R
    else
       L := Min(R - len, -y0);
end;
//*****************************************************************************
procedure Process_1_Quater;
var
    len : double;
begin
    len := sqrt(x0 * x0 + y0 * y0);

    if (len >= R) then     // if outside of circle
        L := len - R
    else
       L := Min(R - len, x0);
end;
//*****************************************************************************
procedure Process_4_Quater;
var
    len : double;
begin
    len := sqrt(x0 * x0 + y0 * y0);

    if (len >= R) then     // if outside of circle
        L := len - R
    else
       L := Min(R - len, len);  
end;
//*****************************************************************************
begin
  Write(' x0 = '); read(x0); writeln;
  Write(' y0 = '); read(y0); writeln;
  Write(' R = ');  read(R);  writeln;
  
  if (x0 * y0 <= 0) and (x0 <= 0) then
    Process_2_Quater
  else if (x0 * y0 >= 0) and (x0 <= 0) then
    Process_3_Quater
  else if (x0 * y0 >= 0) and (x0 > 0) then
    Process_1_Quater
  else
    Process_4_Quater;

  Writeln(' L=', L : 5 : 5);
  Readln;
  Readln;
end.
//*****************************************************************************

Powered by Invision Power Board (http://www.invisionboard.com)
© Invision Power Services (http://www.invisionpower.com)