%I #23 Jan 18 2021 02:41:03
%S 3,4,7,10,13,18,19,22,24,25,29,32,34,37,42,43,45,46,53,55,56,60,61,62,
%T 71,78,79,80,81,82,85,89,94,98,99,101,104,105,107,108,113,114,115,118,
%U 121,131,132,134,139,140,144,146,150,151,152,153,155,163,166,173,174,176,181,182,187,189,192,193,194,195,199,200,204
%N One part of a 3-way classification of the positive integers. Numbers n for which A048675(n) == 2 (mod 3).
%C The positive integers are partitioned between A332820, A332821 and this sequence.
%C For each prime p, the terms include exactly one of p and p^2. The primes alternate between this sequence and A332821. This sequence has the primes with even indexes, those in A031215.
%C The terms are the even numbers in A332820 halved. The terms are also the numbers m such that 5m is in A332820, and so on for alternate primes: 11, 17, 23 etc. Likewise, the terms are the numbers m such that 3m is in A332821, and so on for alternate primes: 7, 13, 19, 29 etc.
%C If we take each odd term of this sequence and replace each prime in its factorization by the next smaller prime, we get the same set of numbers as we get from halving the even terms of this sequence, and A332821 consists exactly of those numbers. The numbers that are one third of the terms that are multiples of 3 are in A332820, which consists exactly of those numbers. The numbers that are one fifth of the terms that are multiples of 5 constitute A332821, and for larger primes, an alternating pattern applies as described in the previous paragraph.
%C The product of any 2 terms of this sequence is in A332821, the product of any 3 terms is in A332820, and the product of a term of A332820 and a term of this sequence is in this sequence. So if a number k is present, k^2 is in A332821, k^3 is in A332820, and k^4 is in this sequence.
%C If k is an even number, exactly one of {k/2, k, 2k} is in the sequence (cf. A191257 / A067368 / A213258); and generally if k is a multiple of a prime p, exactly one of {k/p, k, k*p} is in the sequence.
%F {a(n) : n >= 1} = {2 * A332821(k) : k >= 1} U {A003961(A332821(k)) : k >= 1}.
%F {a(n) : n >= 1} = {A332821(k)^2 : k >= 1} U {A331590(2, A332821(k)) : k >= 1}.
%t Select[Range@ 204, Mod[Total@ #, 3] == 2 &@ Map[#[[-1]]*2^(PrimePi@ #[[1]] - 1) &, FactorInteger[#]] &] (* _Michael De Vlieger_, Mar 15 2020 *)
%o (PARI) isA332822(n) = { my(f = factor(n)); (2==((sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2)%3)); };
%Y Positions of terms valued -1 in A332823; equivalently, numbers in row 3k-1 of A277905 for some k >= 1.
%Y Cf. A048675, A332820, A332821.
%Y Comparable 2 or 3-way classifications: A000379/A000028, A001969/A000069, A003159/A036554, A005843/A005408, A028260/A026424, A191257/A067368/A213258, A325431/A325432, A329609/A329604/A332812.
%Y Subsequences: intersection of A026478 and A066207, A031215 (prime terms), A033430\{0}, A117642\{0}, A169604, A244727\{0}, A244729\{0}, A338910 (semiprime terms).
%K nonn
%O 1,1
%A _Antti Karttunen_ and _Peter Munn_, Feb 25 2020