|
|
A014899
|
|
a(n) = (16^(n+1) - 15*n - 16)/225.
|
|
3
|
|
|
0, 1, 18, 291, 4660, 74565, 1193046, 19088743, 305419896, 4886718345, 78187493530, 1250999896491, 20015998343868, 320255973501901, 5124095576030430, 81985529216486895, 1311768467463790336, 20988295479420645393, 335812727670730326306
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
FORMULA
|
a(n) = 16*a(n-1) + n = 18*a(n-1) - 33*a(n-2) + 16*a(n-3).
|
|
MAPLE
|
a:=n->sum((16^(n-j)-1)/15, j=0..n): seq(a(n), n=1..16); # Zerinvary Lajos, Jan 05 2007
n0:=20: tabl:=array(1..n0-1): for n from 0 to n0 do: tabl[n+1]:=(4^(2*n+2) - 15*n - 16)/225:od:print( tabl): # Michel Lagneau, Apr 26 2010
|
|
MATHEMATICA
|
Table[(16^(n+1)-15*n-16)/225, {n, 0, 20}] (* Harvey P. Dale, Dec 20 2010 *)
|
|
PROG
|
(Magma) I:=[0, 1, 18]; [n le 3 select I[n] else 18*Self(n-1) - 33*Self(n-2) + 16*Self(n-3): n in [1..20]]; // Vincenzo Librandi, Oct 20 2012
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|