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”).

A339010
a(n) is the number of ways to write n as the difference of two centered k-gonal numbers for k >= 3.
1
0, 0, 1, 1, 1, 2, 1, 2, 3, 2, 1, 5, 1, 2, 5, 3, 1, 6, 1, 5, 5, 2, 1, 8, 3, 2, 6, 5, 1, 10, 1, 4, 5, 2, 5, 12, 1, 2, 5, 8, 1, 10, 1, 5, 12, 2, 1, 11, 3, 6, 5, 5, 1, 12, 5, 8, 5, 2, 1, 19, 1, 2, 12, 5, 5, 10, 1, 5, 5, 10, 1, 18, 1, 2, 12, 5, 5, 10, 1, 11, 10, 2
OFFSET
1,6
COMMENTS
Records occur at indices n = 1, 3, 6, 9, 12, 18, 24, 30, 36, 60, 90, 120, 180, 270, 360, 420, 540, 630, 840, 1080, ...
LINKS
Code Golf Stack Exchange, Uncentered Polygons
Eric Weisstein's World of Mathematics, Centered Polygonal Number
FORMULA
a(n) = Sum_{d|n, 3*d <= n} A001227(d).
EXAMPLE
For n = 35, the a(35) = 5 differences are:
A101321( 5,4) - A101321( 5,2) = 51 - 16 = 35,
A101321( 5,7) - A101321( 5,6) = 141 - 106 = 35,
A101321( 7,3) - A101321( 7,1) = 43 - 8 = 35,
A101321( 7,5) - A101321( 7,4) = 106 - 71 = 35, and
A101321(36,1) - A101321(36,0) = 36 - 1 = 35.
PROG
(PARI) a(n) = sumdiv(n, d, if (3*d <= n, numdiv(d>>valuation(d, 2)))); \\ Michel Marcus, Nov 19 2020
CROSSREFS
Cf. A333822 (polygonal numbers), A333836 (positive polygonal numbers), A333868 (binomial coefficients), A333880 (perfect powers).
Sequence in context: A332509 A336157 A190167 * A332842 A352696 A171565
KEYWORD
nonn
AUTHOR
Peter Kagey, Nov 18 2020
STATUS
approved