login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A161806
A trisection of A161804: a(n) = A161804(3n+1) for n>=0.
3
3, 30, 141, 513, 1815, 5727, 15882, 42417, 108165, 255831, 585258, 1302966, 2762349, 5705829, 11577633, 22708053, 43675938, 83011398, 153929484, 281210994, 509494515, 905832642, 1591395774, 2778237765, 4776943011
OFFSET
0,1
COMMENTS
G.f. of A161804 is exp( Sum_{n>=1} A002129(n) * 3*A038500(n) * q^n/n ),
where A002129 forms the l.g.f. of log[ Sum_{n>=0} x^(n(n+1)/2) ], and
A038500(n) is the highest power of 3 dividing n.
EXAMPLE
G.f.: T_1(q) = 3 + 30*q + 141*q^2 + 513*q^3 + 1815*q^4 + 5727*q^5 +...
Terms are divisible by 3:
A/3=[1,10,47,171,605,1909,5294,14139,36055,85277,195086,434322,...].
PROG
(PARI) {a(n)=local(L=sum(m=1, 3*n+1, 3*3^valuation(m, 3)*sumdiv(m, d, -(-1)^d*d)*x^m/m)+x*O(x^(3*n+1))); polcoeff(exp(L), 3*n+1)}
CROSSREFS
Cf. A161804, other trisections: A161805 (T_0), A161807 (T_2).
Sequence in context: A221516 A174774 A020874 * A003689 A127868 A002463
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 20 2009
STATUS
approved