login
Octagonal numbers of the form C*(3C - 2) with 3C - 2 = k^2 and C a composite number.
1

%I #12 Aug 19 2020 09:19:45

%S 96,225,1408,3400,9633,12936,22016,43681,53600,78408,93633,130625,

%T 152776,205408,308481,350208,446216,696008,854400,943041,1250656,

%U 1493896,1628033,1923201,2085000,2439008,3280256,3774408,4041441

%N Octagonal numbers of the form C*(3C - 2) with 3C - 2 = k^2 and C a composite number.

%H Robert Israel, <a href="/A125511/b125511.txt">Table of n, a(n) for n = 1..10000</a>

%p From _R. J. Mathar_, Apr 22 2010: (Start)

%p A000567 := proc(n) n*(3*n-2) ; end proc:

%p for n from 4 to 1600 do if not isprime(n) then if issqr(3*n-2) then printf("%d,", A000567(n) ); end if; end if; end do: (End)

%t Reap[Do[If[CompositeQ[n], If[IntegerQ[Sqrt[3n-2]], Sow[n(3n-2)]]], {n, 4, 10000}]][[2, 1]] (* _Jean-François Alcover_, Aug 19 2020 *)

%Y Cf. A124106.

%K nonn

%O 1,1

%A _Giovanni Teofilatto_, Jan 19 2007

%E Extended beyond 22016 by _R. J. Mathar_, Apr 22 2010