OFFSET
1,2
COMMENTS
The number 1 is in the sequence by convention.
All primes p are trivially in the sequence.
All semiprimes greater than 4 are in the sequence because they have only two unordered factorizations pq = p*q whose sums are distinct. They are distinct because the only solution to p*q = p+q is p=q=2.
If a number m is not in the sequence, then all multiples of m are not in the sequence. For example, multiples of 4 are not in the sequence because there always exist at least two factorizations 4*k = 2*2*k whose factors sum to the same value 4+k = 2+2+k.
The complement is in A337080.
LINKS
Eric Weisstein's World of Mathematics, Unordered Factorization.
EXAMPLE
All unordered factorization of 30 are 30 = 2*15 = 3*10 = 5*6 = 2*3*5. Corresponding sums of factors are distinct: 30, 17 = 15+2, 13 = 10+3, 11 = 6+5, 10 = 2+3+5. Therefore 30 is in the sequence.
All unordered factorization of 90 are 90 = 45*2 = 30*3 = 18*5 = 15*6 = 15*3*2 = 10*9 = 9*5*2 = 10*3*3 = 6*5*3 = 5*3*3*2. Corresponding sums of factors are not all distinct: 90, 57, 33, 23, 21, 20, 19, 16, 16, 14, 13 because the sum 16 = 10+3+3 = 9+5+2 appears twice. Therefore 90 is not in the sequence.
PROG
(PARI) factz(n, minn) = {my(v=[]); fordiv(n, d, if ((d>=minn) && (d<=sqrtint(n)), w = factz(n/d, d); for (i=1, #w, w[i] = concat([d], w[i]); ); v = concat(v, w); ); ); concat(v, [[n]]); }
factorz(n) = factz(n, 2);
isok(n) = my(vs = apply(x->vecsum(x), factorz(n))); #vs == #Set(vs); \\ Michel Marcus, Aug 13 2020
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Matej Veselovac, Aug 12 2020
STATUS
approved