OFFSET
1,2
COMMENTS
The positive integers are partitioned between this sequence, A332821 and A332822, which list the integers in respective cosets of the subgroup.
As the sequence lists the integers in a multiplicative subgroup of the positive rationals, the sequence is closed under multiplication and, provided the result is an integer, under division.
It follows that for any n in this sequence, all powers n^k are present (k >= 0), as are all cubes.
If we take each odd term of this sequence and replace each prime in its factorization by the next smaller prime, the resulting numbers are a permutation of the full sequence; and if we take the square root of each square term we get the full sequence.
There are no primes in the sequence, therefore if k is present and p is a prime, k*p and k/p are absent (noting that k/p might not be an integer). This property extends from primes to all terms of A050376 (often called Fermi-Dirac primes), therefore to squares of primes, 4th powers of primes etc.
The terms are the even numbers in A332821 halved. The terms are also the numbers m such that 5m is in A332821, and so on for alternate primes: 11, 17, 23 etc. Likewise, the terms are the numbers m such that 3m is in A332822, and so on for alternate primes: 7, 13, 19 etc.
The numbers that are half of the even terms of this sequence are in A332822, which consists exactly of those numbers. The numbers that are one third of the terms that are multiples of 3 are in A332821, which consists exactly of those numbers. These properties extend in a pattern of alternating primes as described in the previous paragraph.
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.
If m and n are in this sequence then so is m*n (the definition of "multiplicative semigroup"), while if n is in this sequence, and x is in the complement A359830, then n*x is in A359830. This essentially follows from the fact that A048675 is totally additive sequence. Compare to A329609. - Antti Karttunen, Jan 17 2023
FORMULA
From Peter Munn, Mar 17 2021: (Start)
{a(n) : n >= 1} = {k : k >= 1, 3|A048675(k)}.
{a(n) : n >= 1} = {k : k >= 1, 3|A195017(k)}.
(End)
MATHEMATICA
Select[Range@ 225, Or[Mod[Total@ #, 3] == 0 &@ Map[#[[-1]]*2^(PrimePi@ #[[1]] - 1) &, FactorInteger[#]], # == 1] &] (* Michael De Vlieger, Mar 15 2020 *)
PROG
(PARI) isA332820(n) = { my(f = factor(n)); !((sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2)%3); };
CROSSREFS
Cf. A048675, A195017, A332821, A332822, A353350 (characteristic function), A353348 (its Dirichlet inverse), A359830 (complement).
Comparable 2 or 3-way classifications: A000379/A000028, A001969/A000069, A003159/A036554, A005843/A005408, A028260/A026424, A191257/A067368/A213258, A325431/A325432, A329609/A329604/A332812.
Subsequences: A000578\{0}, A006094, A090090, A099788, A245630 (A191002 in ascending order), A244726\{0}, A325698, A338471, A338556, A338907.
Subsequence of {1} U A268388.
KEYWORD
nonn
AUTHOR
Antti Karttunen and Peter Munn, Feb 25 2020
EXTENSIONS
New name from Peter Munn, Mar 08 2021
STATUS
approved