%I #12 May 08 2021 08:28:55
%S 1,2,3,4,5,7,8,9,11,12,16,17,18,19,20,23,24,25,27,28,31,32,37,40,44,
%T 45,48,49,50,53,54,56,59,61,63,64,67,68,71,72,75,76,80,81,83,88,89,92,
%U 96,97,98,99,103,107,108,112,121,124,125,127,128,131,135,136,144,147,148
%N Numbers whose prime multiplicities are distinct and whose prime indices are term of the sequence.
%C A prime index of n is a number m such that prime(m) divides n.
%C Also Matula-Goebel numbers of rooted trees in which the multiplicities in the multiset of branches directly under any given node are distinct.
%H Andrew Howroyd, <a href="/A319272/b319272.txt">Table of n, a(n) for n = 1..1000</a>
%e 36 is not in the sequence because 36 = 2^2 * 3^2 does not have distinct prime multiplicities.
%e The sequence of terms of the sequence followed by their Matula-Goebel trees begins:
%e 1: o
%e 2: (o)
%e 3: ((o))
%e 4: (oo)
%e 5: (((o)))
%e 7: ((oo))
%e 8: (ooo)
%e 9: ((o)(o))
%e 11: ((((o))))
%e 12: (oo(o))
%e 16: (oooo)
%e 17: (((oo)))
%e 18: (o(o)(o))
%e 19: ((ooo))
%e 20: (oo((o)))
%e 23: (((o)(o)))
%e 24: (ooo(o))
%e 25: (((o))((o)))
%e 27: ((o)(o)(o))
%e 28: (oo(oo))
%e 31: (((((o)))))
%t mgsiQ[n_]:=Or[n==1,And[UnsameQ@@Last/@FactorInteger[n],And@@Cases[FactorInteger[n],{p_,_}:>mgsiQ[PrimePi[p]]]]];
%t Select[Range[100],mgsiQ]
%o (PARI) is(n)={my(f=factor(n)); if(#Set(f[,2])<#f~, 0, for(i=1, #f~, if(!is(primepi(f[i,1])), return(0))); 1)}
%o { select(is, [1..200]) } \\ _Andrew Howroyd_, Mar 01 2020
%Y Cf. A000081, A004111, A007097, A061775, A098859, A130091, A255231, A276625, A316793, A316794, A316795, A316796.
%K nonn
%O 1,2
%A _Gus Wiseman_, Sep 16 2018
%E Terms a(53) and beyond from _Andrew Howroyd_, Mar 01 2020