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

A135769
Pentagonal numbers > 0 which are not the difference of two other pentagonal numbers > 0.
5
1, 5, 12, 51, 92, 117, 176, 330, 477, 852, 1080, 4347, 9801, 13776, 24512, 27270, 39285, 69876, 88452, 124272, 137562, 220992, 268182, 315792, 354051, 403782, 523626, 612162, 629532, 699392, 796797, 919242, 1119312, 1158522, 1190376
OFFSET
1,2
COMMENTS
A subsequence of A136113, obtained by omitting A136113(A135771(k)), k=1,2,3,... ; i.e. those which are not the difference of two larger pentagonal numbers, but the difference of a larger and a smaller pentagonal number.
The definition ("...two other...") excludes the case P(n) = P(m)-P(n), cf. comment by R. J. Mathar in A000326.
LINKS
FORMULA
P(n)=n*(3*n-1)/2 <=> n*(n-1/3) = (2/3)*P(n), thus m = P(n) <=> m = P([sqrt(2m/3)]+1)
and m = P(n) <=> 24m+1 = (6n-1)^2, useful for investigating the possibility of writing P(n)=P(n')+P(n"): this is possible whenever (6n-1)^2=(6n'-1)^2+(6n"-1)^2-1.
EXAMPLE
See A135768 for a list of P(n) which are in A136113 but not in A135769.
PROG
(PARI) P(n)=n*(3*n-1)/2
isPent(t)=P(sqrtint((t*2)\3)+1)==t
for( i=1, 999, for( j=1, (P(i)-1)\3, isPent(P(i)+P(j))&j!=i&next(2)); print1(P(i)", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
R. J. Mathar and M. F. Hasler, Feb 07 2008
STATUS
approved