login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A292072 Main diagonal of A292068. 2
1, -1, -3, -20, 66, 4439, 454420, 4873175, -3803048954, -7320203267692, -1403057989033446, 6669491545211096686, 78492109668913945526447, 69591502229308312804788424, -6243846072108996200105800383026, -604234376454072219680822138902122079 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
LINKS
FORMULA
a(n) = [x^n] Product_{k=1..n} 1/(1 + k^n*x^k).
MAPLE
b:= proc(n, i, k) option remember; (m->
`if`(m<n, 0, `if`(n=m, i!^k, b(n, i-1, k)+
`if`(i>n, 0, i^k*b(n-i, i-1, k)))))(i*(i+1)/2)
end:
g:= proc(n, k) option remember; `if`(n=0, 1,
-add(b(n-i$2, k)*g(i, k), i=0..n-1))
end:
a:= n-> g(n$2):
seq(a(n), n=0..15); # Alois P. Heinz, Sep 12 2017
MATHEMATICA
b[n_, i_, k_] := b[n, i, k] = Function[m, If[m < n, 0, If[n == m, i!^k, b[n, i - 1, k] + If[i > n, 0, i^k*b[n - i, i - 1, k]]]]][i*(i + 1)/2];
g[n_, k_] := g[n, k] = If[n == 0, 1, -Sum[b[n-i, n-i, k]*g[i, k], {i, 0, n-1}]];
a[n_] := g[n, n];
Table[a[n], {n, 0, 15}] (* Jean-François Alcover, Jun 03 2018, after Alois P. Heinz *)
PROG
(PARI) {a(n) = polcoeff(1/prod(k=1, n, 1+k^n*x^k+x*O(x^n)), n)}
(Python)
from sympy.core.cache import cacheit
from sympy import factorial as f
@cacheit
def b(n, i, k):
m=i*(i + 1)/2
return 0 if m<n else f(i)**k if n==m else b(n, i - 1, k) + (0 if i>n else i**k*b(n - i, i - 1, k))
@cacheit
def g(n, k): return 1 if n==0 else -sum([b(n - i, n - i, k)*g(i, k) for i in range(n)])
def a(n): return g(n, n)
print([a(n) for n in range(16)]) # Indranil Ghosh, Sep 14 2017, after Maple program
CROSSREFS
Sequence in context: A099721 A024402 A345689 * A183377 A067600 A348208
KEYWORD
sign
AUTHOR
Seiichi Manyama, Sep 12 2017
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 05:48 EDT 2024. Contains 371265 sequences. (Running on oeis4.)