cp's OEIS Frontend

Square array A(n,k), n >= 1, k >= 1, read by descending antidiagonals.

The group defined by the binary operation A059897(.,.) over the positive integers is commutative with all elements self-inverse, and isomorphic to the additive group of GF(2) polynomial rings such as GF(2)[x,y]. There is a unique isomorphism extending each bijective mapping between respective minimal generating sets. The lexicographically earliest minimal generating set for the A059897 group is A050376, often called the Fermi-Dirac primes. This set has a natural arrangement in a square array, given by A329050(i,j) = prime(i+1)^(2^j), i >= 0, j >= 0. The most meaningful generating set for the additive group of GF(2)[x,y] is {x^i * y^j: i >= 0, j >= 0}, which similarly forms a square array. All this makes A329050(i,j) especially appropriate to be the image (under an isomorphism) of the GF(2) polynomial x^i * y^j.

Using g to denote the intended isomorphism, we specify g(x^i * y^j) = A329050(i,j). This maps minimal generating sets of the additive groups, so the definition of g is completed by specifying g(a+b) = A059897(g(a), g(b)). We then calculate the image under g of polynomial multiplication in GF(2)[x,y], giving us this sequence as the matching multiplicative operator for an isomorphic ring over the positive integers. Using f to denote the inverse of g, A[n,k] = g(f(n) * f(k)).

See the formula section for an alternative definition based on the A329050 array, independent of GF(2)[x,y].

Closely related to A306697 and A297845. If A059897 is replaced in the alternative definition by A059896 (and the definition is supplemented by the derived identity for the absorbing element, shown in the formula section), we get A306697; if A059897 is similarly replaced by A003991 (integer multiplication), we get A297845. This sequence and A306697, considered as multiplicative operators, are carryless arithmetic equivalents of A297845. A306697 uses a method analogous to binary-OR when there would be a multiplicative carry, while this sequence uses a method analogous to binary exclusive-OR. In consequence A(n,k) <> A297845(n,k) exactly when A306697(n,k) <> A297845(n,k). This relationship is not symmetric between the 3 sequences: there are n and k such that A(n,k) = A306697(n,k) <> A297845(n,k). For example A(54,72) = A306697(54,72) = 273375000 <> A297845(54,72) = 22143375000.

From Peter Munn, Apr 02 2019: (Start)

Also column 6 of A059897.

A self-inverse permutation of the positive integers with no fixed points; A073675 composed with A120229.

The permutation swaps pairs of integers whose ratio is 1:6 or 2:3, these ratios corresponding to the factorizations 1*6 = 2*3 = 6. Row 6 is the first row of A059897 to exhibit more than 1 such ratio.

(End)

The integers in the pairs with ratio 1:6 are listed in A036668, the integers in the pairs with ratio 2:3 are listed in A325424. - Peter Munn, Mar 05 2020

Also column 5 of A059897.

A self-inverse permutation of the positive integers with no fixed points that swaps pairs of integers whose ratio is 1:5.

Also column 8 of A059897.

A self-inverse permutation of the positive integers with no fixed points; A073675 composed with A120230.

The permutation swaps pairs of integers whose ratio is 1:8 or 1:2, these ratios corresponding to the factorizations 1*8 = 2*4 = 8.

Also column 24 of A059897.

A self-inverse permutation of the positive integers with no fixed points; a composition of A073675, A120229 and A120230.

The permutation swaps pairs of integers whose ratio is 1:24, 1:6, 3:8 or 2:3, these ratios corresponding to the factorizations 1*24 = 2*12 = 3*8 = 4*6 = 24. Row 24 is the first row of A059897 to exhibit more than 2 such ratios.

When creating A329329, the author realized it was isomorphic to multiplication in the GF(2)[x,y] polynomial ring. However, A329329 was unusual in having A059897(.,.) as additive operator, whereas the equivalent univariate polynomial ring, GF(2)[x], is more commonly mapped (to integers) with bitwise exclusive-or (A003987) representing polynomial addition (and A048720(.,.) representing polynomial multiplication). This sequence shows how multiplication in GF(2)[x] can look when mapped to integers with A059897(.,.) representing polynomial addition.

The group defined by the binary operation A059897(.,.) over the positive integers is commutative with all elements self-inverse, and isomorphic to the additive group of the polynomial ring GF(2)[x]. There is a unique isomorphism extending each bijective mapping between respective minimal generating sets. The lexicographically earliest minimal generating set for the A059897 group is A050376, often called the Fermi-Dirac primes. The most meaningful generating set for the additive group of GF(2)[x] is {x^i: i >= 0}.

Using f to denote the intended isomorphism from GF(2)[x], we specify f(x^i) = A050376(i+1). This maps minimal generating sets of the additive groups, so the definition of f is completed by specifying f(a+b) = A059897(f(a), f(b)). We then calculate the image under f of polynomial multiplication in GF(2)[x], giving us this sequence as the matching multiplication operation for an isomorphic ring over the positive integers. Using g to denote the inverse of f, A(n,k) = f(g(n) * g(k)).

Note that A050376 is closed with respect to A(.,.).

Recall that GF(2)[x] is more usually mapped to integers with A003987(.,.) as addition and A048720(.,.) as multiplication. With this usual mapping, under which A000079(i) is the image of x^i, A052330(.) is the relevant isomorphism from nonnegative integers under A048720(.,.) and A003987(.,.) to positive integers under A(.,.) and A059897(.,.), with A052331(.) its inverse.

A permutation of the squarefree numbers A005117. The missing positive numbers are in A013929. - Alois P. Heinz, Sep 06 2014

From Antti Karttunen, Apr 18 & 19 2017: (Start)

Because a(n) toggles the parity of n there are neither fixed points nor any cycles of odd length.

Conjecture: there are no finite cycles of any length. My grounds for this conjecture: any finite cycle in this sequence, if such cycles exist at all, must have at least one member that occurs somewhere in A285319, the terms that seem already to be quite rare. Moreover, any such a number n should satisfy in addition to A019565(n) < n also that A048675^{k}(n) is squarefree, not just for k=0, 1 but for all k >= 0. As there is on average a probability of only 6/(Pi^2) = 0.6079... that any further term encountered on the trajectory of A048675 is squarefree, the total chance that all of them would be squarefree (which is required from the elements of A019565-cycles) is soon minuscule, especially as A048675 is not very tightly bounded (many trajectories seem to skyrocket, at least initially). I am also assuming that usually there is no significant correlation between the binary expansions of n and A048675(n) (apart from their least significant bits), or, for that matter, between their prime factorizations.

See also the slightly stronger conjecture in A285320, which implies that there would neither be any two-way infinite cycles.

If either of the conjectures is false (there are cycles), then certainly neither sequence A285332 nor its inverse A285331 can be a permutation of natural numbers. (End)

The conjecture made in A087207 (see also A288569) implies the two conjectures mentioned above. A further constraint for cycles is that in any A019565-trajectory which starts from a squarefree number (A005117), every other term is of the form 4k+2, while every other term is of the form 6k+3. - Antti Karttunen, Jun 18 2017

The sequence satisfies the exponential function identity, a(x + y) = a(x) * a(y), whenever x and y do not have a 1-bit in the same position, i.e., when A004198(x,y) = 0. See also A283475. - Antti Karttunen, Oct 31 2019

The above identity becomes unconditional if binary exclusive OR, A003987(.,.), is substituted for addition, and A059897(.,.), a multiplicative equivalent of A003987, is substituted for multiplication. This gives us a(A003987(x,y)) = A059897(a(x), a(y)). - Peter Munn, Nov 18 2019

Also the Heinz number of the binary indices of n, where the Heinz number of a sequence (y_1,...,y_k) is prime(y_1)*...*prime(y_k), and a number's binary indices (A048793) are the positions of 1's in its reversed binary expansion. - Gus Wiseman, Dec 28 2022

Least positive k such that m^k + 1 divides m^n + 1 (with fixed base m). - Vladimir Baltic, Mar 25 2002

To construct the sequence: start with 1, concatenate 1, 1 and double last term gives 1, 2. Concatenate those 2 terms, 1, 2, 1, 2 and double last term 1, 2, 1, 2 -> 1, 2, 1, 4. Concatenate those 4 terms: 1, 2, 1, 4, 1, 2, 1, 4 and double last term -> 1, 2, 1, 4, 1, 2, 1, 8, etc. - Benoit Cloitre, Dec 17 2002

a(n) = gcd(seq(binomial(2*n, 2*m+1)/2, m = 0 .. n - 1)) (odd numbered entries of even numbered rows of Pascal's triangle A007318 divided by 2), where gcd() denotes the greatest common divisor of a set of numbers. Due to the symmetry of the rows it suffices to consider m = 0 .. floor((n-1)/2). - Wolfdieter Lang, Jan 23 2004

Equals the continued fraction expansion of a constant x (cf. A100338) such that the continued fraction expansion of 2*x interleaves this sequence with 2's: contfrac(2*x) = [2; 1, 2, 2, 2, 1, 2, 4, 2, 1, 2, 2, 2, 1, 2, 8, 2, ...].

Simon Plouffe observes that this sequence and A003484 (Radon function) are very similar, the difference being all zeros except for every 16th term (see A101119 for nonzero differences). Dec 02 2004

This sequence arises when calculating the next odd number in a Collatz sequence: Next(x) = (3*x + 1) / A006519, or simply (3*x + 1) / BitAnd(3*x + 1, -3*x - 1). - Jim Caprioli, Feb 04 2005

a(n) = n if and only if n = 2^k. This sequence can be obtained by taking a(2^n) = 2^n in place of a(2^n) = n and using the same sequence building approach as in A001511. - Amarnath Murthy, Jul 08 2005

Also smallest m such that m + n - 1 = m XOR (n - 1); A086799(n) = a(n) + n - 1. - Reinhard Zumkeller, Feb 02 2007

Number of 1's between successive 0's in A159689. - Philippe Deléham, Apr 22 2009

Least number k such that all coefficients of k*E(n, x), the n-th Euler polynomial, are integers (cf. A144845). - Peter Luschny, Nov 13 2009

In the binary expansion of n, delete everything left of the rightmost 1 bit. - Ralf Stephan, Aug 22 2013

The equivalent sequence for partitions is A194446. - Omar E. Pol, Aug 22 2013

Also the 2-adic value of 1/n, n >= 1. See the Mahler reference, definition on p. 7. This is a non-archimedean valuation. See Mahler, p. 10. Sometimes called 2-adic absolute value of 1/n. - Wolfdieter Lang, Jun 28 2014

First 2^(k-1) - 1 terms are also the heights of the successive rectangles and squares of width 2 that are adjacent to any of the four sides of the toothpick structure of A139250 after 2^k stages, with k >= 2. For example: if k = 5 the heights after 32 stages are [1, 2, 1, 4, 1, 2, 1, 8, 1, 2, 1, 4, 1, 2, 1] respectively, the same as the first 15 terms of this sequence. - Omar E. Pol, Dec 29 2020