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

A275340
Nontrivial centered polygonal numbers: numbers of the form A101321(n,k) where n >= 1 and k >= 2.
4
4, 7, 10, 11, 13, 16, 19, 21, 22, 25, 28, 29, 31, 34, 37, 40, 41, 43, 46, 49, 51, 52, 55, 56, 57, 58, 61, 64, 67, 70, 71, 73, 76, 79, 81, 82, 85, 88, 91, 92, 94, 97, 100, 101, 103, 106, 109, 111, 112, 113, 115, 118, 121, 124, 127, 130, 131, 133, 136, 137, 139, 141, 142, 145, 148, 151
OFFSET
1,1
COMMENTS
This is a centered polygonal number analog to A090466.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
FORMULA
a(n) ~ c*n with 1.95 < c < 2.29. - Charles R Greathouse IV, Jul 28 2016
MAPLE
isA275340 := proc(n)
local nsearch, ksearch;
for nsearch from 1 do
if A101321(nsearch, 2) > n then
return false;
end if;
for ksearch from 2 do
if A101321(nsearch, ksearch) = n then
return true;
elif A101321(nsearch, ksearch) > n then
break;
end if;
end do:
end do:
end proc:
for n from 1 to 400 do
if isA275340(n) then
printf("%d, ", n) ;
end if;
end do:
MATHEMATICA
maxTerm = 1000;
Table[1+n*k*(k+1)/2, {n, 1, maxTerm-1}, {k, 2, Sqrt[2maxTerm] // Ceiling}] // Flatten // Union // Select[#, # <= maxTerm&]& (* Jean-François Alcover, Jun 17 2023 *)
PROG
(PARI) list(lim)=my(v=List(), t); lim\=1; for(k=2, sqrt(8*lim-7)/2, t=k*(k+1)/2; forstep(a=t+1, lim, t, listput(v, a))); Set(v) \\ Charles R Greathouse IV, Jul 28 2016
CROSSREFS
Cf. A101321.
Sequence in context: A072125 A370759 A223024 * A082206 A115566 A364840
KEYWORD
nonn
AUTHOR
R. J. Mathar, Jul 28 2016
STATUS
approved