OFFSET
1,1
COMMENTS
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
By Hall's marriage theorem, k is a term if and only if there is a sub-multiset S of the prime indices of k such that fewer than |S| numbers are divisors of a member of S. Equivalently, k is divisible by a member of A370348. - Robert Israel, Feb 15 2024
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Wikipedia, Cartesian product.
EXAMPLE
The terms together with their prime indices begin:
4: {1,1}
8: {1,1,1}
12: {1,1,2}
16: {1,1,1,1}
18: {1,2,2}
20: {1,1,3}
24: {1,1,1,2}
27: {2,2,2}
28: {1,1,4}
32: {1,1,1,1,1}
36: {1,1,2,2}
40: {1,1,1,3}
44: {1,1,5}
48: {1,1,1,1,2}
For example, the choices of a divisor of each prime index of 90 are: (1,1,1,1), (1,1,1,3), (1,1,2,1), (1,1,2,3), (1,2,1,1), (1,2,1,3), (1,2,2,1), (1,2,2,3). But none of these has all distinct elements, so 90 is in the sequence.
MAPLE
filter:= proc(n) uses numtheory, GraphTheory; local B, S, F, D, E, G, t, d;
F:= ifactors(n)[2];
F:= map(t -> [pi(t[1]), t[2]], F);
D:= `union`(seq(divisors(t[1]), t = F));
F:= map(proc(t) local i; seq([t[1], i], i=1..t[2]) end proc, F);
if nops(D) < nops(F) then return false fi;
E:= {seq(seq({t, d}, d=divisors(t[1])), t = F)};
S:= map(t -> convert(t, name), [op(F), op(D)]);
E:= map(e -> map(convert, e, name), E);
G:= Graph(S, E);
B:= BipartiteMatching(G);
B[1] = nops(F);
end proc:
remove(filter, [$1..200]); # Robert Israel, Feb 15 2024
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Select[Range[100], Select[Tuples[Divisors/@primeMS[#]], UnsameQ@@#&]=={}&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jul 22 2022
STATUS
approved