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.

A141151 L.g.f.: A(x) = log( Sum_{n>=0} n^n*x^n ) = Sum_{n>=1} a(n)*x^n/n.

This page as a plain text file.
%I A141151 #14 Sep 08 2022 08:45:35
%S A141151 1,7,70,899,14001,255532,5342541,125876003,3300437302,95338188007,
%T A141151 3009043615073,103043811158864,3805827820399125,150819894172935183,
%U A141151 6383815674758486310,287459477551898694403,13721584934214631377921
%N A141151 L.g.f.: A(x) = log( Sum_{n>=0} n^n*x^n ) = Sum_{n>=1} a(n)*x^n/n.
%H A141151 Seiichi Manyama, <a href="/A141151/b141151.txt">Table of n, a(n) for n = 1..385</a>
%F A141151 a(n) ~ n^(n+1). - _Vaclav Kotesovec_, May 30 2019
%e A141151 L.g.f.: A(x) = x + 7*x^2/2 + 70*x^3/3 + 899*x^4/4 + 14001*x^5/5 + ...
%e A141151 exp(A(x)) = 1 + x + 4*x^2 + 27*x^3 + 256*x^4 + 3125*x^5 + 46656*x^6 + ...
%t A141151 With[{m=20}, Drop[CoefficientList[Series[Log[Sum[If[n==0, 1, n^n*x^n], {n, 0, m+2}]], {x, 0, m}], x], 1]*Range[1, m]] (* _G. C. Greubel_, May 30 2019 *)
%o A141151 (PARI) {a(n)=polcoeff(x*deriv(log(Ser(concat(1,vector(n+1,k,k^k))))),n)}
%o A141151 (Magma) m:=20; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Log( (&+[n^n*x^n: n in [0..m+2]]) ) )); [n*b[n]: n in [1..m-1]]; // _G. C. Greubel_, May 30 2019
%o A141151 (Sage) m = 20; T = taylor(log( sum(n^n*x^n for n in (0..m+2)) ), x, 0, m); [n*T.coefficient(x, n) for n in (1..m)] # _G. C. Greubel_, May 30 2019
%Y A141151 Cf. A141152, A306628.
%K A141151 nonn
%O A141151 1,2
%A A141151 _Paul D. Hanna_, Jun 11 2008