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.

A167008 a(n) = Sum_{k=0..n} C(n,k)^k.

This page as a plain text file.
%I A167008 #37 Apr 14 2023 07:14:45
%S A167008 1,2,4,14,106,1732,66634,5745700,1058905642,461715853196,
%T A167008 461918527950694,989913403174541980,5009399946447021173140,
%U A167008 60070720443204091719085184,1548154498059133199618813305334,92346622775540905956057053976278584
%N A167008 a(n) = Sum_{k=0..n} C(n,k)^k.
%C A167008 Row sums of A219206.
%H A167008 Reinhard Zumkeller, <a href="/A167008/b167008.txt">Table of n, a(n) for n = 0..75</a>
%H A167008 Vaclav Kotesovec, <a href="https://oeis.org/wiki/User:Vaclav_Kotesovec">Interesting asymptotic formulas for binomial sums</a>, Jun 09 2013.
%F A167008 Limit_{n->oo} a(n)^(1/n^2) = (1-r)^(-r/2) = 1.533628065110458582053143..., where r = A220359 = 0.70350607643066243... is the root of the equation (1-r)^(2*r-1) = r^(2*r). - _Vaclav Kotesovec_, Dec 12 2012
%t A167008 Flatten[{1,Table[Sum[Binomial[n, k]^k, {k,0,n}], {n,20}]}]
%t A167008 (* Program for numerical value of the limit a(n)^(1/n^2) *) (1-r)^(-r/2)/.FindRoot[(1-r)^(2*r-1)==r^(2*r),{r,1/2},WorkingPrecision->100] (* _Vaclav Kotesovec_, Dec 12 2012 *)
%t A167008 Total/@Table[Binomial[n,k]^k,{n,0,20},{k,0,n}] (* _Harvey P. Dale_, Oct 19 2021 *)
%o A167008 (PARI) a(n)=sum(k=0,n,binomial(n,k)^k)
%o A167008 (Haskell) a167008 = sum . a219206_row  -- _Reinhard Zumkeller_, Feb 27 2015
%o A167008 (Magma) [(&+[Binomial(n,j)^j: j in [0..n]]): n in [0..20]]; // _G. C. Greubel_, Aug 26 2022
%o A167008 (SageMath) [sum(binomial(n,j)^j for j in (0..n)) for n in (0..20)] # _G. C. Greubel_, Aug 26 2022
%Y A167008 Cf. A000169, A014062, A167009, A167010, A184731, A219206, A220359, A362288.
%K A167008 nonn,nice
%O A167008 0,2
%A A167008 _Paul D. Hanna_, Nov 17 2009