A372351 -id:A372351 - cp's OEIS frontend

A371094 a(n) = m*(2^e) + ((4^e)-1)/3, where m = 3n+1, and e is the 2-adic valuation of m.

1, 21, 7, 21, 13, 341, 19, 45, 25, 117, 31, 69, 37, 341, 43, 93, 49, 213, 55, 117, 61, 5461, 67, 141, 73, 309, 79, 165, 85, 725, 91, 189, 97, 405, 103, 213, 109, 1877, 115, 237, 121, 501, 127, 261, 133, 1109, 139, 285, 145, 597, 151, 309, 157, 5461, 163, 333, 169, 693, 175, 357, 181, 1493, 187, 381, 193, 789, 199

Offset: 0

Antti Karttunen (proposed by Ali Sada), Apr 19 2024

Construction: take the binary expansion of 3n+1 (A016777(n)), and substitute "01" for all trailing 0-bits that follow after its odd part (= A067745(1+n)), of which there are A371093(n) in total. See the examples.

			For n=1, 3*n+1 = 4, "100" in binary, when we substitute 01's for the two trailing 0's, we obtain 21, "10101" in binary, therefore a(1) = 21.
For n=6, 3*6+1 = 19, "10011" in binary, and there are no trailing 0's, and no changes, therefore a(6) = 19.
For n=7, 3*7+1 = 22, "10110" in binary, with one trailing 0, which when replaced with 01 gives us 45, "101101" in binary, therefore a(7) = 45.
For n=229, there are e=4 trailing bit expansions 0 -> 01,
  3n+1 = binary  101011  0 0 0 0
  a(n) = binary  101011 01010101

Antti Karttunen, Table of n, a(n) for n = 0..16384
Michael De Vlieger, Plot binary expansion of a(n) arranged with smallest bit at bottom and largest at top, n increasing from left to right for n = 0..3071. Black indicates 1, white indicates 0.
Index entries for sequences related to 3x+1 (or Collatz) problem
Index entries for sequences related to binary expansion of n

Cf. A016777, A067745, A075677, A086893, A096773, A372289, A371093.
Cf. A016921, A372351 (even and odd bisection), A372290 (numbers occurring in the latter).
Cf. arrays A372282, A372283 and A371096, A371100, A371102.
Cf. also A302338.

Array[#2*(2^#3) + ((4^#3) - 1)/3 & @@ {#1, #2, IntegerExponent[#2, 2]} & @@ {#, 3 #1 + 1} &, 67, 0] (* Michael De Vlieger, Apr 19 2024 *)

A371094(n) = { my(m=1+3*n, e=valuation(m,2)); ((m*(2^e)) + (((4^e)-1)/3)); };

def A371094(n): return ((m:=3*n+1)<<(e:=(~m & m-1).bit_length()))+((1<<(e<<1))-1)//3 # Chai Wah Wu, Apr 28 2024

a(n) = A372289(A016777(n)).

a(2n) = A016777(2n) = A016921(n).

A372282 Array read by upward antidiagonals: A(n, k) = A371094(A(n-1, k)) for n > 1, k >= 1; A(1, k) = 2*k-1.

Original entry on oeis.org

1, 21, 3, 5461, 21, 5, 357913941, 5461, 341, 7, 1537228672809129301, 357913941, 1398101, 45, 9, 28356863910078205288614550619314017621, 1537228672809129301, 23456248059221, 1109, 117, 11, 9649340769776349618630915417390658987772498722136713669954798667326094136661, 28356863910078205288614550619314017621, 6602346876188694799461995861, 873813, 11605, 69, 13

Offset: 1

Antti Karttunen, Apr 28 2024

			Array begins:
n\k|    1     2        3     4      5     6        7     8      9     10
---+----------------------------------------------------------------------
1  |    1,    3,       5,    7,     9,   11,      13,   15,    17,    19,
2  |   21,   21,     341,   45,   117,   69,     341,   93,   213,   117,
3  | 5461, 5461, 1398101, 1109, 11605, 3413, 1398101, 2261, 87381, 11605,

Cf. A005408 (row 1), A372351 (row 2, bisection of A371094), A372444 (column 14).
Arrays derived from this one:
  A372283,
  A372285 the number of terms of A086893 in the interval [A(n, k), A(1+n, k)],
  A372287 the column index of A(n, k) in array A257852,
  A372288 the sum of digits of A(n, k) in "Jacobsthal greedy base",
  A372353 differences between A(n,k) and the largest term of A086893 <= A(n,k),
  A372354 floor(log_2(.)) of terms, A372356 (and their columnwise first differences),
  A372359 terms xored with binary words of the same length, either of the form 10101...0101 or 110101...0101, depending on whether the binary length is odd or even.
Cf. also arrays A371096, A371102 that give subsets of columns of this array, and array A371100 that gives the terms of the row 2 in different order.

up_to = 28;
A371094(n) = { my(m=1+3*n, e=valuation(m,2)); ((m*(2^e)) + (((4^e)-1)/3)); };
A372282sq(n,k) = if(1==n,2*k-1,A371094(A372282sq(n-1,k)));
A372282list(up_to) = { my(v = vector(up_to), i=0); for(a=1,oo, for(col=1,a, i++; if(i > up_to, return(v)); v[i] = A372282sq((a-(col-1)),col))); (v); };
v372282 = A372282list(up_to);
A372282(n) = v372282[n];

A371100 Array A read by upward antidiagonals in which the entry A(n,k) in row n and column k is defined by A(n, k) = 4^n(6k - 3 - 2*(-1)^n) + (4^n - 1)/3, n,k >= 1.

Original entry on oeis.org

21, 21, 45, 341, 117, 69, 341, 725, 213, 93, 5461, 1877, 1109, 309, 117, 5461, 11605, 3413, 1493, 405, 141, 87381, 30037, 17749, 4949, 1877, 501, 165, 87381, 185685, 54613, 23893, 6485, 2261, 597, 189, 1398101, 480597, 283989, 79189, 30037, 8021, 2645, 693, 213, 1398101, 2970965, 873813, 382293, 103765, 36181, 9557, 3029, 789, 237

Offset: 1

Antti Karttunen and Ali Sada, Apr 18 2024

			The top left corner of the array:
n\k|      1       2       3        4        5        6        7        8
---+--------------------------------------------------------------------------
1  |     21,     45,     69,      93,     117,     141,     165,     189, ...
2  |     21,    117,    213,     309,     405,     501,     597,     693, ...
3  |    341,    725,   1109,    1493,    1877,    2261,    2645,    3029, ...
4  |    341,   1877,   3413,    4949,    6485,    8021,    9557,   11093, ...
5  |   5461,  11605,  17749,   23893,   30037,   36181,   42325,   48469, ...
6  |   5461,  30037,  54613,   79189,  103765,  128341,  152917,  177493, ...
7  |  87381, 185685, 283989,  382293,  480597,  578901,  677205,  775509, ...
8  |  87381, 480597, 873813, 1267029, 1660245, 2053461, 2446677, 2839893, ...
...

Paolo Xausa, Table of n, a(n) for n = 1..11325 (first 150 antidiagonals, flattened).

Cf. A007283, A079319, A257852, A371094, A371101.
Cf. A372351 (same terms, in different order), A372290 (sorted into ascending order, without duplicates), A372293 (odd numbers that do not occur here).
Leftmost column is A144864 duplicated, without its initial 1.
Row 1: A102603.

A371100[n_, k_] := 4^n*(6*k - 3 - 2*(-1)^n) + (4^n - 1)/3;
Table[A371100[n - k + 1, k], {n, 10}, {k, n}] (* Paolo Xausa, Apr 21 2024 *)

up_to = 55;
A371100sq(n,k) = 4^n*(6*k - 3 - 2*(-1)^n) + (4^n - 1)/3;
A371100list(up_to) = { my(v = vector(up_to), i=0); for(a=1,oo, for(col=1,a, i++; if(i > up_to, return(v)); v[i] = A371100sq((a-(col-1)),col))); (v); };
v371100 = A371100list(up_to);
A371100(n) = v371100[n];

A(n, k) = A007283(n)*A257852(n,k) + A079319(n).
A(n, k) = A371094(A257852(n,k)).
A(n+2, k) = 5 + 16*A(n,k).

A372288 Array read by upward antidiagonals: A(n, k) = A265745(A372282(n, k)), n,k >= 1, where A265745(n) is the sum of digits in "Jacobsthal greedy base".

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 3, 3, 1, 1, 1, 3, 3, 1, 1, 1, 1, 3, 3, 3, 3, 1, 1, 1, 1, 3, 3, 1, 3, 1, 1, 1, 1, 3, 1, 1, 3, 3, 1, 1, 1, 1, 1, 1, 1, 5, 3, 3, 1, 1, 1, 1, 1, 1, 1, 3, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 5, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 3, 5, 3

Offset: 1

Antti Karttunen, Apr 28 2024

Collatz conjecture is equal to the claim that each column will eventually settle to constant 1's, somewhere under its topmost row. This works as only the bisection A002450 of Jacobsthal numbers (A001045) contains numbers of the form 4k+1, while the other bisection contains only numbers of the form 4k+3, which do not occur among the range of A372351. See also the comments in A371094.

			Array begins:
n\k| 1  2  3  4  5  6  7  8  9 10 11 12 13     14 15    16 17 18 19 20    21 22
---+----------------------------------------------------------------------------
1  | 1, 1, 1, 3, 3, 1, 3, 3, 3, 3, 1, 3, 3,     3, 3,    3, 3, 3, 3, 5,    5, 1,
2  | 1, 1, 1, 3, 3, 3, 1, 3, 3, 3, 1, 5, 5,     5, 3,    5, 3, 3, 3, 5,    5, 3,
3  | 1, 1, 1, 3, 3, 3, 1, 5, 1, 3, 1, 3, 3,     5, 3,    5, 5, 1, 3, 3,    5, 3,
4  | 1, 1, 1, 3, 3, 1, 1, 3, 1, 3, 1, 1, 3,     5, 3,    3, 3, 1, 3, 5,    5, 3,
5  | 1, 1, 1, 1, 3, 1, 1, 1, 1, 3, 1, 1, 3,     5, 1,    5, 3, 1, 3, 3,    3, 3,
6  | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3,     3, 1,    5, 3, 1, 1, 5,    5, 3,
7  | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     5, 1,    3, 3, 1, 1, 3,    5, 3,
8  | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     5, 1,    5, 1, 1, 1, 3,    3, 3,
9  | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     3, 1,    5, 1, 1, 1, 3,    5, 1,
10 | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     5, 1,    5, 1, 1, 1, 3,    5, 1,
11 | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     5, 1, 2155, 1, 1, 1, 1,    5, 1,
12 | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     5, 1,    5, 1, 1, 1, 1, 6251, 1,
13 | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 10347, 1,    5, 1, 1, 1, 1,    5, 1,
14 | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     5, 1,    5, 1, 1, 1, 1,    5, 1,
15 | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     5, 1,    7, 1, 1, 1, 1,    5, 1,
16 | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     5, 1,    5, 1, 1, 1, 1,    7, 1,

Index entries for sequences related to 3x+1 (or Collatz) problem

Cf. A001045, A002450, A265745, A265747, A371094, A372351, A372282, A372283, A372287.

Cf. also array A372561 (formed by columns whose indices in this array are given by A372443).

up_to = 105;
A130249(n) = (#binary(3*n+1)-1);
A001045(n) = (2^n - (-1)^n) / 3;
A265745(n) = { my(s=0); while(n,s++; n -= A001045(A130249(n))); (s); };
A371094(n) = { my(m=1+3*n, e=valuation(m,2)); ((m*(2^e)) + (((4^e)-1)/3)); };
A372282sq(n,k) = if(1==n,2*k-1,A371094(A372282sq(n-1,k)));
A372288sq(n,k) = A265745(A372282sq(n,k));
A372288list(up_to) = { my(v = vector(up_to), i=0); for(a=1,oo, for(col=1,a, i++; if(i > up_to, return(v)); v[i] = A372288sq((a-(col-1)),col))); (v); };
v372288 = A372288list(up_to);
A372288(n) = v372288[n];

cp's OEIS Frontend

A371094 a(n) = m*(2^e) + ((4^e)-1)/3, where m = 3n+1, and e is the 2-adic valuation of m.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A372282 Array read by upward antidiagonals: A(n, k) = A371094(A(n-1, k)) for n > 1, k >= 1; A(1, k) = 2*k-1.

Views

Author

Keywords

Examples

Crossrefs

Programs

PARI

A371100 Array A read by upward antidiagonals in which the entry A(n,k) in row n and column k is defined by A(n, k) = 4^n*(6*k - 3 - 2*(-1)^n) + (4^n - 1)/3, n,k >= 1.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A372288 Array read by upward antidiagonals: A(n, k) = A265745(A372282(n, k)), n,k >= 1, where A265745(n) is the sum of digits in "Jacobsthal greedy base".

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

A371100 Array A read by upward antidiagonals in which the entry A(n,k) in row n and column k is defined by A(n, k) = 4^n(6k - 3 - 2*(-1)^n) + (4^n - 1)/3, n,k >= 1.