Реализации алгоритмов/Числа Эйлера первого рода: различия между версиями
Содержимое удалено Содержимое добавлено
Maxal (обсуждение | вклад) перенесено из w:Числа Эйлера I рода |
Maxal (обсуждение | вклад) мНет описания правки |
||
Строка 10: | Строка 10: | ||
{ E(n, k) = sum(j=0, k, (-1)^j * (k-j)^n * binomial(n+1,j) ) } |
{ E(n, k) = sum(j=0, k, (-1)^j * (k-j)^n * binomial(n+1,j) ) } |
||
== [[w:Python]] == |
== [[w:Python|Python]] == |
||
<source lang="python"> |
<source lang="python"> |
||
from math import factorial |
from math import factorial |
Версия от 17:54, 7 апреля 2011
PARI/GP
\\ рекуррентная формула { E(n, k) = if(k<1 || k>n, 0, if(n==1, 1, k*E(n-1,k) + (n-k+1)*E(n-1,k-1) ); ) } \\ явная формула { E(n, k) = sum(j=0, k, (-1)^j * (k-j)^n * binomial(n+1,j) ) }
Python
from math import factorial
# вспомогательные функции
def binomial(n, k):
if n >= k:
return factorial(n) / (factorial(k) * factorial(n-k))
else:
return 0
def table(f, n):
return [[f(i, j) for j in range(n)] for i in range(n)]
def format(table):
return "\n".join([" ".join(["%6s" % col for col in row]) for row in table])
# рекуррентная формула
def euler_rec(n, k):
if k<1:
return 0
elif k>n:
return 0
elif n==1:
return 1
else:
return k * euler_rec(n-1, k) + (n-k+1) * euler_rec(n-1, k-1)
# явная формула
def euler_iter(n, k):
if k<1:
return 0
elif k>n:
return 0
elif n==1:
return 1
else:
return sum([binomial(n+1, j) * (-1)**j * (k-j)**n for j in range(k)])
# демонстрация работы
print(format(table(euler_rec, 10)))
print(format(table(euler_iter, 10)))
"""
0 0 0 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0 0
0 1 1 0 0 0 0 0 0 0
0 1 4 1 0 0 0 0 0 0
0 1 11 11 1 0 0 0 0 0
0 1 26 66 26 1 0 0 0 0
0 1 57 302 302 57 1 0 0 0
0 1 120 1191 2416 1191 120 1 0 0
0 1 247 4293 15619 15619 4293 247 1 0
0 1 502 14608 88234 156190 88234 14608 502 1
"""