cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A354341 a(n) = n! * Sum_{k=1..n} ( Sum_{d|k} 1/(d * ((k/d)!)^d) )/(n-k)!.

Original entry on oeis.org

1, 4, 12, 38, 130, 557, 2877, 18314, 136458, 1180457, 11389081, 122833207, 1446973931, 18594740348, 257507754524, 3835059283282, 60937544854850, 1030871972064485, 18469079943443229, 349656695460113159, 6969526853682012755, 145958486484692023936
Offset: 1

Views

Author

Seiichi Manyama, Aug 15 2022

Keywords

Crossrefs

Programs

  • PARI
    a182926(n) = n!*sumdiv(n, d, 1/(d*(n/d)!^d));
    a(n) = sum(k=1, n, a182926(k)*binomial(n, k));
    
  • PARI
    my(N=30, x='x+O('x^N)); Vec(serlaplace(-exp(x)*sum(k=1, N, log(1-x^k/k!))))

Formula

a(n) = Sum_{k=1..n} A182926(k) * binomial(n,k).
E.g.f.: -exp(x) * Sum_{k>0} log(1-x^k/k!).

A356597 Expansion of e.g.f. ( Product_{k>0} 1/(1 - x^k/k) )^exp(x).

Original entry on oeis.org

1, 1, 5, 26, 172, 1354, 12403, 127945, 1471006, 18589503, 255951308, 3808299648, 60871219649, 1039240205691, 18868377309780, 362838034712928, 7364831540699076, 157305165900364641, 3526069495916583260, 82744901973286823822, 2028396974232995349291
Offset: 0

Views

Author

Seiichi Manyama, Aug 15 2022

Keywords

Crossrefs

Programs

  • PARI
    my(N=30, x='x+O('x^N)); Vec(serlaplace(1/prod(k=1, N, 1-x^k/k)^exp(x)))
    
  • PARI
    a354339(n) = n!*sum(k=1, n, sumdiv(k, d, 1/(d*(k/d)^d))/(n-k)!);
    a_vector(n) = my(v=vector(n+1)); v[1]=1; for(i=1, n, v[i+1]=sum(j=1, i, a354339(j)*binomial(i-1, j-1)*v[i-j+1])); v;

Formula

a(0) = 1; a(n) = Sum_{k=1..n} A354339(k) * binomial(n-1,k-1) * a(n-k).
Showing 1-2 of 2 results.