%I #94 Aug 16 2024 23:09:36
%S 0,0,14,78,252,620,1290,2394,4088,6552,9990,14630,20724,28548,38402,
%T 50610,65520,83504,104958,130302,159980,194460,234234,279818,331752,
%U 390600,456950,531414,614628,707252,809970,923490,1048544,1185888,1336302,1500590,1679580
%N a(n) = n^4 - n.
%C a(n) is the number of ways to assign 4 different students to n different dorm rooms, each of which can hold at most 3 students. In other words, a(n) is the number of functions f:[4]->[n] with the size of the pre-image set of each element of the codomain at most 3. - _Dennis P. Walsh_, Mar 21 2013
%C a(n) are the values of m that yield integer solutions to this family of equations: x = sqrt(m + sqrt(x)), which may also be viewed as an infinitely recursive radical. The real solutions for x at each m = a(n) is n^2, except at n = 1 (m = 0) where x = 0 or 1 is a solution. - _Richard R. Forberg_, Oct 15 2014
%H Harry J. Smith, <a href="/A058895/b058895.txt">Table of n, a(n) for n = 0..2000</a>
%H Milan Janjic, <a href="http://www.pmfbl.org/janjic/">Enumerative Formulas for Some Functions on Finite Sets</a>
%H <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (5,-10,10,-5,1).
%F a(n) = n*(n-1)*(n^2+n+1) = A000583(n) - n = A002061(n+1) * A002378(n-1) = (n-1) * A027444(n) = -n * A024001(n).
%F a(n) = 2*A027482(n). - _Zerinvary Lajos_, Jan 28 2008
%F a(n) = floor(n^7/(n^3+1)). - _Gary Detlefs_, Feb 11 2010
%F a(n)^3 = (a(n)/n)^4 + (a(n)/n)^3. - _Vincenzo Librandi_, Feb 23 2012
%F a(n)^3 + A068601(n)^3 + A033562(n)^3 = A185065(n)^3, for n > 0. - _Vincenzo Librandi_, Mar 13 2012
%F G.f.: 2*x^2*(7+4*x+x^2)/(1-x)^5. - _Colin Barker_, Apr 23 2012
%F a(n) = 14*C(n,2) + 36*C(n,3) + 24*C(n,4). - _Dennis P. Walsh_, Mar 21 2013
%F Sum_{n>=2} (-1)^n/a(n) = (Pi/3)*sech(Pi*sqrt(3)/2) + 4*log(2)/3 - 1 = 0.06147271494... . - _Amiram Eldar_, Jul 04 2020
%F Sum_{n>=2} 1/a(n) = A339605. - _R. J. Mathar_, Jan 08 2021
%F E.g.f.: exp(x)*x^2*(7 + 6*x + x^2). - _Stefano Spezia_, Jul 09 2021
%F a(n) = 12*A000332(n+2) + 2*A000537(n-1). - _Yasser Arath Chavez Reyes_, Apr 05 2024
%p seq(n*(n^3-1),n=0..25) ; # _R. J. Mathar_, Dec 10 2015
%t Table[n^4 - n, {n, 0, 40}] (* _Vladimir Joseph Stephan Orlovsky_, Feb 20 2012 *)
%o (PARI) { for (n = 0, 2000, write("b058895.txt", n, " ", n^4-n); ) } \\ _Harry J. Smith_, Jun 23 2009
%o (Magma) [n^4-n: n in [0..40]]; // _Vincenzo Librandi_, Feb 23 2012
%o (Sage) (2*x^2*(7+4*x+x^2)/(1-x)^5).series(x, 37).coefficients(x, sparse=False) # _Stefano Spezia_, Jul 09 2021
%Y Cf. A000583, A002061, A002378, A024001, A027444, A027482, A068601, A033562, A185065, A339605.
%K easy,nonn
%O 0,3
%A _Henry Bottomley_, Jan 08 2001