%I #33 Jan 03 2022 19:25:25
%S 2523,3112819,5656019,10132171,12167825,16639567,25302173,31995475,
%T 35158921,37334419,43890719,44816821,47715269,53548223,55534523,
%U 90526075,90533525,127558319,142929025,143167073,144989575,147182225
%N Numbers n such that n, n+2, n+4, n+6 are of the form p^2*q where p and q are distinct primes.
%C All terms are odd. Proof: if n is even then out of the 4 numbers n, n+2, n+4, n+6, 2 of them must be either both of the form 2*p^2, 2*q^2, or both of the form 4*p, 4*q. In either case, for p != q and p, q prime, the difference between these 2 numbers are more than 6, reaching a contradiction. - _Chai Wah Wu_, Jun 24 2019
%H Giovanni Resta, <a href="/A308736/b308736.txt">Table of n, a(n) for n = 1..10000</a> (first 123 terms from Ray Chandler)
%e 2523 = 3*29*29, 2525 = 5*5*101, 2527 = 7*19*19, 2529 = 3*3*281.
%t psx = Table[{0}, {7}]; nmax = 150000000; n = 1; lst = {};
%t While[n < nmax, n++;
%t psx = RotateRight[psx];
%t psx[[1]] = Sort[Last /@ FactorInteger[n]];
%t If[Union[{psx[[1]], psx[[3]], psx[[5]], psx[[7]]}] == {{1, 2}}, AppendTo[lst, n - 6]];];
%t lst
%o (Python)
%o from sympy import factorint
%o A308736_list, n, mlist = [], 3, [False]*4
%o while len(A308736_list) < 100:
%o if mlist[0] and mlist[1] and mlist[2] and mlist[3]:
%o A308736_list.append(n)
%o n += 2
%o f = factorint(n+6)
%o mlist = mlist[1:] + [(len(f),sum(f.values())) == (2,3)] # _Chai Wah Wu_, Jun 24 2019, Jan 03 2022.
%Y Cf. A074173, A308735.
%K nonn
%O 1,1
%A _Ray Chandler_, Jun 24 2019