%I #38 Jun 18 2021 02:37:36
%S 2,4,6,8,10,12,14,15,16,18,20,22,24,26,28,30,32,33,34,35,36,38,40,42,
%T 44,45,46,48,50,51,52,54,55,56,58,60,62,64,66,68,69,70,72,74,75,76,77,
%U 78,80,82,84,85,86,88,90,92,93,94,95,96,98,99,100,102,104
%N Numbers k such that the gcd of the indices j for which the j-th prime prime(j) divides k is 1.
%C Any integer k in the sequence encodes (by 'Heinz encoding' cf. A056239) a multiset of integers whose gcd is 1, namely the multiset containing r_j copies of j if k factors as Product_j prime(j)^{r_j} with gcd_j j = 1.
%C Clearly the sequence contains all even numbers and no odd primes or odd prime powers. It also clearly contains all numbers that are divisible by consecutive primes.
%C The sequence is the list of those k such that A289508(k) = 1.
%C It is also the list of those k such that A289506(k) = A289507(k).
%C Heinz numbers of integer partitions with relatively prime parts, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - _Gus Wiseman_, Apr 13 2018
%H Alois P. Heinz, <a href="/A289509/b289509.txt">Table of n, a(n) for n = 1..20000</a>
%e 6 is a term because 6 = p_1*p_2 and gcd(1,2) = 1.
%e From _Gus Wiseman_, Apr 13 2018: (Start)
%e Sequence of integer partitions with relatively prime parts begins:
%e 02 : (1)
%e 04 : (11)
%e 06 : (21)
%e 08 : (111)
%e 10 : (31)
%e 12 : (211)
%e 14 : (41)
%e 15 : (32)
%e 16 : (1111)
%e 18 : (221)
%e 20 : (311)
%e 22 : (51)
%e 24 : (2111)
%e 26 : (61)
%e 28 : (411)
%e 30 : (321)
%e 32 : (11111)
%e 33 : (52)
%e 34 : (71)
%e 35 : (43)
%e 36 : (2211)
%e 38 : (81)
%e 40 : (3111)
%e (End)
%p p:=1:for ind to 10000 do p:=nextprime(p);primeindex[p]:=ind;od:
%p out:=[]:for n from 2 to 100 do m:=[];f:=ifactors(n)[2];g:=0;
%p for k to nops(f) do mk:=primeindex[f[k][1]];m:=[op(m),mk];
%p g:=gcd(g,mk);od; if g=1 then out:=[op(out),n];fi;od:out;
%t Select[Range[200],GCD@@PrimePi/@FactorInteger[#][[All,1]]===1&] (* _Gus Wiseman_, Apr 13 2018 *)
%o (PARI) isok(n) = my(f=factor(n)); gcd(apply(x->primepi(x), f[,1])) == 1; \\ _Michel Marcus_, Jul 19 2017
%o (Python)
%o from sympy import gcd, primepi, primefactors
%o def ok(n): return gcd([primepi(p) for p in primefactors(n)]) == 1
%o print([n for n in range(1, 151) if ok(n)]) # _Indranil Ghosh_, Aug 06 2017
%Y Cf. A001222, A007359, A051424, A056239, A289506, A289507, A289508, A296150, A302696, A302697, A302698, A302796.
%K nonn
%O 1,1
%A _Christopher J. Smyth_, Jul 11 2017