Алгебра у программистов. Преподаватель пишет на доске: sin X = 1.
- Ну, кто может найти X?
Выбегает программист и так радостно:
- Вот, вот X, - показывает на X.
Это может звучать тривиально, но как мне возвести число в степень? Например,
2^12 = 4095.
На самом деле вопрос далеко не тривиальный. Проблема в том, что сам алгоритм
функции далеко не прост. Функцией Power(X, N) (т.е. X^N) должны четко
отслеживаться несколько возможных ситуаций:
X любое число, N = 0
X = 1, N любое число
X = 0 и N > 0
X = 0 и N < 0
X > 0
X < 0 и N нечетное целое
X < 0 и N целое
X < 0 и N нецелое
Посмотрите на следующую, абсолютно правильно
работающую функцию (тем не менее она может быть и не самой эффективной!):
interfacetype
EPowerException = class(Exception)
end;
implementationfunction Power(X, N: real): extended;
var
t: longint;
r: real;
isInteger: boolean;
beginif N = 0 thenbegin
result := 1.0;
exit;
end;
if X = 1.0 thenbegin
result := 1.0;
exit;
end;
if X = 0.0 thenbeginif N > 0.0 thenbegin
result := 0.0;
exit;
endelseraise EPowerException.Create('Результат - бесконечность');
end;
if (X > 0) thentry
result := exp(N * ln(X));
exit;
exceptraise
EPowerException.Create('Результат - переполнение или потеря значимости');
end;
{ X - отрицательный, но мы все еще можем вычислить результат, если n целое. }{ пытаемся получить целую часть n с использованием типа longint, вычисление }{ четности n не займет много времени }try
t := trunc(n);
if (n - t) = 0 then
isInteger := trueelse
isInteger := False;
except{ Лишний бит может вызвать переполнение или потерю значимости }
r := int(n);
if (n - r) = 0 thenbegin
isInteger := true;
if frac(r / 2) = 0.5 then
t := 1
else
t := 2;
endelse
isInteger := False;
end;
if isInteger thenbegin{n целое}if odd(t) then{n нечетное}try
result := -exp(N * ln(-X));
exit;
exceptraise
EPowerException.Create('Результат - переполнение или потеря значимости');
endelse{n четное}try
result := exp(N * ln(-X));
exit;
exceptraise
EPowerException.Create('Результат - переполнение или потеря значимости');
end;
endelseraise EPowerException.Create('Результат невычисляем');
end;
Статья Возведение числа в степень раздела Синтаксис Математика может быть полезна для разработчиков на Delphi и FreePascal.
Комментарии и вопросы
Материалы статей собраны из открытых источников, владелец сайта не претендует на авторство. Там где авторство установить не удалось, материал подаётся без имени автора. В случае если Вы считаете, что Ваши права нарушены, пожалуйста, свяжитесь с владельцем сайта.