|
|
A008133
|
|
a(n) = floor(n/3)*floor((n+1)/3).
|
|
4
|
|
|
0, 0, 0, 1, 1, 2, 4, 4, 6, 9, 9, 12, 16, 16, 20, 25, 25, 30, 36, 36, 42, 49, 49, 56, 64, 64, 72, 81, 81, 90, 100, 100, 110, 121, 121, 132, 144, 144, 156, 169, 169, 182, 196, 196, 210, 225, 225, 240, 256, 256, 272, 289, 289, 306, 324, 324, 342, 361, 361, 380, 400, 400, 420, 441, 441, 462, 484, 484, 506
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,6
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
Partial sums of A087509. a(n+1)=sum{j=0..n, sum{k=0..j, if (mod(jk, 3)=2, 1, 0) }}. - Paul Barry, Sep 14 2003
Empirical g.f.: -x^3*(x^2+1) / ((x-1)^3*(x^2+x+1)^2). - Colin Barker, Mar 31 2013
|
|
PROG
|
(Haskell)
a008133 n = a008133_list !! n
a008133_list = zipWith (*) (tail ts) ts where ts = map (`div` 3) [0..]
(PARI) a(n) = floor(n/3)*floor((n+1)/3); /* Joerg Arndt, Mar 31 2013 */
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|