OFFSET
1,1
COMMENTS
A semiprime (A001358) is a product of any two prime numbers.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
The sequence of terms together with their prime indices begins:
16: {1,1,1,1}
32: {1,1,1,1,1}
64: {1,1,1,1,1,1}
81: {2,2,2,2}
96: {1,1,1,1,1,2}
128: {1,1,1,1,1,1,1}
160: {1,1,1,1,1,3}
192: {1,1,1,1,1,1,2}
224: {1,1,1,1,1,4}
243: {2,2,2,2,2}
256: {1,1,1,1,1,1,1,1}
288: {1,1,1,1,1,2,2}
320: {1,1,1,1,1,1,3}
352: {1,1,1,1,1,5}
384: {1,1,1,1,1,1,1,2}
416: {1,1,1,1,1,6}
448: {1,1,1,1,1,1,4}
486: {1,2,2,2,2,2}
For example, a complete list of all factorizations of 192 into primes or semiprimes is:
(2*2*2*2*2*2*3)
(2*2*2*2*2*6)
(2*2*2*2*3*4)
(2*2*2*4*6)
(2*2*3*4*4)
(2*4*4*6)
(3*4*4*4)
Since none of these is strict, 192 is in the sequence.
MAPLE
filter:= proc(n)
g(map(t -> t[2], ifactors(n)[2]))
end proc;
g:= proc(L) option remember; local x, i, j, t, s, Cons, R;
if nops(L) = 1 then return L[1] > 3
elif nops(L) = 2 then return max(L) > 4
fi;
Cons:= {seq(x[i] + x[i, i] + add(x[j, i], j=1..i-1)
+ add(x[i, j], j=i+1..nops(L)) = L[i], i=1..nops(L))};
R:= traperror(Optimization:-LPSolve(0, Cons, assume=binary));
type(R, string)
end proc:
select(filter, [$2..2000]); # Robert Israel, Dec 28 2020
MATHEMATICA
facs[n_]:=If[n<=1, {{}}, Join@@Table[Map[Prepend[#, d]&, Select[facs[n/d], Min@@#>=d&]], {d, Rest[Divisors[n]]}]];
Select[Range[1000], Select[facs[#], UnsameQ@@#&&SubsetQ[{1, 2}, PrimeOmega/@#]&]=={}&]
CROSSREFS
Allowing only primes gives A013929.
Removing all squares of primes gives A339740.
These are the positions of zeros in A339839.
The complement is A339889.
A002100 counts partitions into squarefree semiprimes.
A293511 are a product of distinct squarefree numbers in exactly one way.
A320663 counts non-isomorphic multiset partitions into singletons or pairs.
A339841 have exactly one factorization into primes or semiprimes.
The following count factorizations:
- A001055 into all positive integers > 1.
- A320655 into semiprimes.
- A320656 into squarefree semiprimes.
- A320732 into primes or semiprimes.
- A322353 into distinct semiprimes.
- A339661 into distinct squarefree semiprimes.
- A339742 into distinct primes or squarefree semiprimes.
- A339839 into distinct primes or semiprimes.
The following count vertex-degree partitions and give their Heinz numbers:
- A321728 is conjectured to count non-half-loop-graphical partitions of n.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Dec 20 2020
STATUS
approved