|
%I #40 Jun 19 2023 12:40:52
%S 4,7,10,11,13,16,19,21,22,25,28,29,31,34,37,40,41,43,46,49,51,52,55,
%T 56,57,58,61,64,67,70,71,73,76,79,81,82,85,88,91,92,94,97,100,101,103,
%U 106,109,111,112,113,115,118,121,124,127,130,131,133,136,137,139,141,142,145,148,151
%N Nontrivial centered polygonal numbers: numbers of the form A101321(n,k) where n >= 1 and k >= 2.
%C This is a centered polygonal number analog to A090466.
%H Charles R Greathouse IV, <a href="/A275340/b275340.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) ~ c*n with 1.95 < c < 2.29. - _Charles R Greathouse IV_, Jul 28 2016
%p isA275340 := proc(n)
%p local nsearch,ksearch;
%p for nsearch from 1 do
%p if A101321(nsearch,2) > n then
%p return false;
%p end if;
%p for ksearch from 2 do
%p if A101321(nsearch,ksearch) = n then
%p return true;
%p elif A101321(nsearch,ksearch) > n then
%p break;
%p end if;
%p end do:
%p end do:
%p end proc:
%p for n from 1 to 400 do
%p if isA275340(n) then
%p printf("%d,",n) ;
%p end if;
%p end do:
%t maxTerm = 1000;
%t 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 *)
%o (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
%Y Cf. A101321.
%K nonn
%O 1,1
%A _R. J. Mathar_, Jul 28 2016
|
