|
Search: id:A060452
|
|
|
| A060452 |
|
Let v = (1,4,9,...,n^2), x = (0,1,2,4,6,...) [first n terms of A002620]; a(n) = v.v * x.x - (v.x)^2. |
|
+0
3
|
|
| 0, 1, 6, 38, 107, 350, 728, 1752, 3090, 6215, 9878, 17654, 26117, 42924, 60256, 93024, 125460, 184509, 241110, 341110, 434511, 595562, 742808, 991640, 1215110, 1586403, 1914822, 2452646, 2922185, 3681560, 4337024, 5385600, 6281704, 7701561, 8904294, 10793862, 12381939, 14858038, 16924440, 20124440, 22778042, 26862143, 30229430, 35383062, 39609933, 46046276, 51299936, 59262560, 65733500, 75499125
(list; graph; listen)
|
|
|
OFFSET
|
1,3
|
|
|
REFERENCES
|
N. J. A. Sloane and V. Vaishampayan, in preparation, 2001.
|
|
MAPLE
|
fv := n->1/30*n*(1+n)*(2*n+1)*(3*n^2+3*n-1); # this is A000538
f1 := n->1/160*(n-1)*(1+n)*(2*n^3+5*n^2+2*n-5);
f2 := n->1/160*n*(n+2)*(2*n^3+n^2-2*n+4);
f7 := n->if n mod 2 = 0 then f2(n) else f1(n) end if; # this is A059859
f3 := n->1/20*n^5+1/8*n^4+1/24*n^3-11/120*n-1/8*n^2;
f4 := n->1/20*n^5+1/8*n^4+1/24*n^3+1/30*n;
f5:-n-> if `mod`(n, 2) = 0 then f4(n) else f3(n) end if; # this is A060453
A060452 := n->f7(n)*fv(n)-f5(n)^2;
|
|
CROSSREFS
|
Cf. A002620, A000538, A059859. Agrees with A060453 for first 37 terms.
Sequence in context: A039293 A055713 A060454 this_sequence A045949 A026296 A037499
Adjacent sequences: A060449 A060450 A060451 this_sequence A060453 A060454 A060455
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
njas and Vinay Vaishampayan, Apr 09 2001
|
|
|
Search completed in 0.002 seconds
|
