Проблема с алгоритмом возведения в степень

expt ?

expt ?

Не понял.

(print (expt 2 5))
>>> 32

upd. не заметил, что её и нужно заменить)

(expt 2 3) 8

Если самому хочется - то в случае с if-ами будет удобнее использовать cond. Ещё вариант - рекурсивно, без setf.

Я не фанат лиспянки, но - фуфуфу писать рекурсивные факториалы и возведение в степень вне университетского пруда для перваков.

Тут вот госпадин Дибвиг поделился кодом

(defun power (x n)
"Power raises x to the nth power. Ν must be an integer >= 0.
This executes in log n time, because of the check for even n.
(cond ((= n 0) 1)
((evenp n) (expt (power x (/ n 2 )) 2))
(t (* x (power x (- n 1))))))

хотелось бы увидеть (trace (power 2 3))

но не получается, подскажешь как?

(defun pow(x n)
    (if(< n 0)
       1
       (if(= n 1)
          x
          (* x (pow x (- n 1))))
       )
    )

clisp переваривает

(defun pow(x n) (if(< n 0) 1 (if(= n 1) x (* x (pow x (- n 1)))) ) )

Может подскажешь почему

(trace-function ’pow)

выдаёт такое

Lisp nesting exceeds ‘max-lisp-eval-depth’

но код работает (pow 2 3) => 8

есть ограничение на глубину вложенности(рекурсии), можно изменить ограничение так

(setq max-lisp-eval-depth 1000)

Но я бы рекомендовал переписать по старинке пошаговым умножением или оптимизированным(при чётном n возводим временный результат) в квадрат, а n делим на 2, т.к. так можно будет отвязаться от размера стека

(setq max-lisp-eval-depth 1000)

На этом я уже горел, у меня стек потек, пришлось всё переписывать в итеративном стиле.

Но я бы рекомендовал переписать по старинке пошаговым умножением или оптимизированным

У меня степени в основном float, и отрицательные числа возводятся в степень. Т.е. там где expt x y

If x is a finite negative number and y is a finite non-integer, expt returns a NaN.

И вот это надо как то обрулить.

Нужно поменять expt x y в коде из-за её

А конкретно-то тебе что нужно? Почему просто не кастануть x y к float?

проблема в том что отрицательные нецелые степени порождают НЧ

Почему просто не кастануть x y к float?

Не понял, это как?

А конкретно-то тебе что нужно?

У меня NaN вылетала.

проблема в том что отрицательные нецелые степени порождают НЧ

Класс, что такое НЧ?

If x is a finite negative number and y is a finite non-integer, expt returns a NaN.

Так это насколько я понимаю для случаев когда результат - мнимое число, под это попадает -1 в степени 1/2, т.е. квадратный корень из -1, а он не определён и если хочется считать именно такое, то придётся в комплексных числах считать

Не Число)

порождают НЧ

Ну так правильно, результат-то комплексный.

ну так прикол в том, что тогда придётся либо забить на такое, либо навелосипедить комплексные числа

Не понял, это как?

(expt (float 2) (float 3))
8.0
(expt (float -2) (float 3))
-8.0

Я может, что то упустил, попробуй код на лиспе, то же самое или число выдаст

[code] (defun pow (x n) (exp (* n (log x))))

(pow -2 0.3) => 0.0e+NaN

[/code]

clisp

(print (expt -2 0.3)) ; #C(0.7236485 0.9960168) 

придётся либо забить на такое

Спасибо, товарищ капитан.

так вы логарифм от отрицательного берёте и тут при всём желании вы не придумаете действительной степени возведя в которую положительной число получится отрицательное

рядовой - смирно!

(print (expt -2 0.3)) ; #C(0.7236485 0.9960168)

Как тогда на clisp такое пролезло?

может, что то упустил

Ты упустил, что результат это не целое и даже не вещественное число. Если хочешь упороться в математику на emacs смотри сюда

ну так полагаю elisp он просто для другого и они выбросили поддержку комплексных, т.к. не нужно.

Если хочешь упороться в математику на emacs смотри сюда

Ты намекаешь, что GNU Emacs Calculator можно использовать прям в elisp code.

на clisp

То другой лисп, коммунальный. Там чего только нет, дажё чёрт в ступе. elisp сильно проще.

Программируется, класс, спасибо.

можно использовать прям в elisp code.

А что мешает? Я хз что ты там делаешь, конечно.

В моём emacs нету evenp. И ещё вместо trace - trace-function.