OFFSET
1,2
COMMENTS
Old name: 'Positions of zeros in A268387: numbers n such that when the exponents e_1 .. e_k in their prime factorization n = p_1^e_1 * ... * p_k^e_k are bitwise-xored together, the result is zero.
From Peter Munn, Sep 14 2019 and Dec 01 2019: (Start)
When trailing zeros are removed from the terms written in base p, for any prime p, every positive integer not divisible by p appears exactly once. This is the lexicographically earliest sequence with this property.
The closure of A238748 with respect to the commutative binary operation A059897(.,.). As integers are self-inverse under A059897(.,.), the sequence thereby forms a subgroup, denoted H, of the positive integers under A059897(.,.). H is a subgroup of A000379.
(The symbol ^ can take on a meaning in relation to a group operation. However, in this comment ^ denotes the power operator for standard integer multiplication.) For any prime p, the subgroup {p^k : k >= 0} and H are each a (left and right) transversal of the other. For k >= 0 and primes p_1 and p_2, the cosets (p_1^k)H and (p_2^k)H are the same.
(End)
From Peter Munn, Dec 01 2021: (Start)
If we take the square root of the square terms we reproduce the sequence itself. The set of all products of a square term and a squarefree term is the sequence as a set.
The terms are the elements of the ideal generated by {6} in the ring defined in A329329. Similarly, the ideal generated by {8} gives A262675. 6 and 8 are images of each other under A225546(.), which is an automorphism of the ring. So this sequence and A262675, as sets, are images of each other under A225546(.). The elements of the ideal generated by {6,8} form the notable set A000379.
(End)
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
OEIS Wiki, Ideal.
Eric Weisstein's World of Mathematics, Closure, Group, Left Transversal, Right Transversal, Square part, Squarefree part.
FORMULA
From Peter Munn, Oct 30 2019: (Start)
For k >= 0, prime p_1, prime p_2, {m : m = A059897(p_1^k, a(n)), n >= 1} = {m : m = A059897(p_2^k, a(n)), n >= 1}.
(End)
From Peter Munn, Nov 24 2021: (Start)
{a(n) : n >= 1} = {A000188(a(n)) : n >= 1}.
(End)
EXAMPLE
1 has an empty factorization, and as XOR of an empty set is zero, 1 is included.
6 = 2^1 * 3^1 and as XOR(1,1) = 0, 6 is included.
30 = 2^1 * 3^1 * 5^1 is NOT included, as XOR(1,1,1) = 1.
360 = 2^3 * 3^2 * 5^1 is included, as the bitwise-XOR of exponents 3, 2 and 1 ("11", "10" and "01" in binary) results zero.
MATHEMATICA
Select[Range[200], # == 1 || BitXor @@ Last /@ FactorInteger[#] == 0 &] (* Amiram Eldar, Nov 27 2020 *)
PROG
(Scheme, with Antti Karttunen's IntSeq-library)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Antti Karttunen, Feb 05 2016
EXTENSIONS
New name from Peter Munn, Jul 15 2024
STATUS
approved