|
|
A005917
|
|
Rhombic dodecahedral numbers: a(n) = n^4 - (n - 1)^4.
(Formerly M4968)
|
|
59
|
|
|
1, 15, 65, 175, 369, 671, 1105, 1695, 2465, 3439, 4641, 6095, 7825, 9855, 12209, 14911, 17985, 21455, 25345, 29679, 34481, 39775, 45585, 51935, 58849, 66351, 74465, 83215, 92625, 102719, 113521, 125055, 137345, 150415, 164289, 178991
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Final digits of a(n), i.e., a(n) mod 10, are repeated periodically with period of length 5 {1,5,5,5,9}. There is a symmetry in this list since the sum of two numbers equally distant from the ends is equal to 10 = 1 + 9 = 5 + 5 = 2*5. Last two digits of a(n), i.e., a(n) mod 100, are repeated periodically with period of length 50. - Alexander Adamchuk, Aug 11 2006
If Y is a 3-subset of a 2n-set X then, for n >= 2, a(n-2) is the number of 4-subsets of X intersecting Y. - Milan Janjic, Nov 18 2007
The numbers are the constant number found in magic squares of order n, where n is an odd number, see the comment in A006003. A Magic Square of side 1 is 1; 3 is 15; 5 is 65 and so on. - David Quentin Dauthier, Nov 07 2008
Two times the area of the triangle with vertices at (0,0), ((n - 1)^2, n^2), and (n^2, (n - 1)^2). - J. M. Bergot, Jun 25 2013
Construct an array M with M(0,n) = 2*n^2 + 4*n + 1 = A056220(n+1), M(n,0) = 2*n^2 + 1 = A058331(n) and M(n,n) = 2*n*(n+1) + 1 = A001844(n). Row(n) begins with all the increasing odd numbers from A058331(n) to A001844(n) and column(n) begins with all the decreasing odd numbers from A056220(n+1) to A001844(n). The sum of the terms in row(n) plus those in column(n) minus M(n,n) equals a(n+1). The first five rows of array M are [1, 7, 17, 31, 49, ...]; [3, 5, 15, 29, 47, ...]; [9, 11, 13, 27, 45, ...]; [19, 21, 23, 25, 43, ...]; [33, 35, 37, 39, 41, ...]. - J. M. Bergot, Jul 16 2013 [This contribution was moved here from A047926 by Petros Hadjicostas, Mar 08 2021.]
For n>=2, these are the primitive sides s of squares of type 2 described in A344332. - Bernard Schott, Jun 04 2021
|
|
REFERENCES
|
J. H. Conway and R. K. Guy, The Book of Numbers, p. 53.
E. Deza and M. M. Deza, Figurate Numbers, World Scientific Publishing, 2012, pp. 123-124.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
|
|
FORMULA
|
a(n) = (2*n - 1)*(2*n^2 - 2*n + 1).
G.f.: x*(1+x)*(1+10*x+x^2)/(1-x)^4. - Simon Plouffe in his 1992 dissertation
More generally, g.f. for n^m - (n - 1)^m is Euler(m, x)/(1 - x)^m, where Euler(m, x) is Eulerian polynomial of degree m (cf. A008292). E.g.f.: x*(exp(y/(1 - x)) - exp(x*y/(1 - x)))/(exp(x*y/(1 - x))-x*exp(y/(1 - x))). - Vladeta Jovovic, May 08 2002
a(n) = sum of the next (2*n - 1) odd numbers; i.e., group the odd numbers so that the n-th group contains (2*n - 1) elements like this: (1), (3, 5, 7), (9, 11, 13, 15, 17), (19, 21, 23, 25, 27, 29, 31), ... E.g., a(3) = 65 because 9 + 11 + 13 + 15 + 17 = 65. - Xavier Acloque, Oct 11 2003
a(n) = 2*n - 1 + 12*Sum_{i = 1..n} (i - 1)^2. - Xavier Acloque, Oct 16 2003
a(n) = (4*binomial(n,2) + 1)*sqrt(8*binomial(n,2) + 1). - Paul Barry, Mar 14 2004
Binomial transform of [1, 14, 36, 24, 0, 0, 0, ...], if the offset is 0. - Gary W. Adamson, Dec 20 2007
Sum_{i=1..n-1}(a(i) + a(i+1)) = 8*Sum_{i=1..n}(i^3 + i) = 16*A002817(n-1) for n > 1. - Bruno Berselli, Mar 04 2011
|
|
MATHEMATICA
|
#[[2]]-#[[1]]&/@Partition[Range[0, 40]^4, 2, 1] (* More efficient than the above Mathematica program because it only has to calculate each 4th power once *) (* Harvey P. Dale, Feb 07 2015 *)
|
|
PROG
|
(Haskell)
a005917 n = a005917_list !! (n-1)
a005917_list = map sum $ f 1 [1, 3 ..] where
f x ws = us : f (x + 2) vs where (us, vs) = splitAt x ws
(Python)
A005917_list, m = [], [24, -12, 2, 1]
for _ in range(10**2):
for i in range(3):
|
|
CROSSREFS
|
(1/12)*t*(2*n^3 - 3*n^2 + n) + 2*n - 1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A063493, A063494, A063495, A063496.
|
|
KEYWORD
|
nonn,easy,nice
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|