|
%I #17 Jan 25 2024 02:40:44
%S 20801,22781,37411,47531,55651,75031,80011,100111,120901,133661,
%T 161101,177661,191101,199001,201001,230051,236531,240901,245311,
%U 263351,279661,289201,313751,323851,326401,368561,376751,436601,439561,445511,472781,475861,488281,507211,539561
%N Centered 10-gonal numbers which are sphenic numbers.
%H Robert Israel, <a href="/A368426/b368426.txt">Table of n, a(n) for n = 1..5080</a>
%e A062786(65) = 20801 = 5 * 65 * (65-1) + 1 = 11 * 31 * 61.
%e A062786(68) = 22781 = 5 * 68 * (68-1) + 1 = 11 * 19 * 109.
%p isc10:= proc(n) issqr(20*n+5) end proc:
%p P:= select(isprime, [seq(seq(10*i+j,j=[-1,1]),i=1..1000)]): nP:= nops(P):
%p M:= P[1]*P[2]*P[-1]:
%p A:= NULL:
%p for i from 1 to nP-2 while P[i]^3 < M do
%p for j from i+1 to nP-1 while P[i]*P[j]^2 < M do
%p for k from j+1 to nP do
%p q:= P[i]*P[j]*P[k];
%p if q > M then break fi;
%p if isc10(q) then A:= A,q fi
%p od od od:
%p sort([A]); # _Robert Israel_, Dec 24 2023
%t Select[Table[5*n*(n+1) + 1, {n, 1, 330}], FactorInteger[#][[;; , 2]] == {1, 1, 1} &] (* _Amiram Eldar_, Dec 24 2023 *)
%Y Intersection of A007304 and A062786.
%K nonn
%O 1,1
%A _Massimo Kofler_, Dec 24 2023
|
