OFFSET
1,2
LINKS
Colin Barker, Table of n, a(n) for n = 1..1000
FORMULA
Conjectures:
a(n) = a(n-1)+2*a(n-3)-2*a(n-4)-a(n-6)+a(n-7) for n>7.
G.f. x*(1+3*x+5*x^2-x^3-x^5+x^6) / ((1-x)^3*(1+x+x^2)^2).
PROG
(PARI)
pg(m, n) = (n^2*(m-2)-n*(m-4))/2 \\ The n-th m-gonal number
pgr(m, r) = n=1; L=List(); while((t=pg(m, n))<r, listput(L, t); n++); Vec(L)
pgpgs(p, q, r) = setunion(pgr(p, r), pgr(q, r))
pgpgs(4, 10, 2000)
CROSSREFS
KEYWORD
nonn
AUTHOR
Colin Barker, Apr 24 2016
STATUS
approved