cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-5 of 5 results.

A267969 Conjectured list of positive numbers k such that the Collatz-like map T: if x > 1 and x odd, x -> 5*x+1 and, if x even, x -> x/2, when started at k, eventually reaches a cycle that does not contain 1 (cf. A232711).

Original entry on oeis.org

5, 10, 13, 17, 20, 26, 27, 33, 34, 40, 43, 52, 54, 66, 68, 80, 83, 86, 104, 105, 108, 132, 136, 160, 166, 172, 181, 185, 208, 210, 211, 215, 216, 245, 263, 264, 269, 272, 275, 320, 329, 332, 344, 362, 370, 416, 420, 422, 430, 432, 435, 453, 457, 463, 490, 526
Offset: 1

Views

Author

Michel Lagneau, Jan 19 2016

Keywords

Comments

A companion to A232711 (which conjecturally lists the numbers whose trajectory reaches 1), and A267970 (which conjecturally lists the numbers whose trajectory diverges).
This is conjectural in that there is no proof that the list is complete as far as it goes. Some of the terms in A267970 could belong to A232711 or this sequence if the trajectory is extended far enough. - N. J. A. Sloane, Jan 23 2016
It appears that the trajectories of all terms in this sequence reach one of two length 10 loops, one containing 13 and the other 17. This has been checked for terms up to 10^4 assuming trajectories with more than 10000 odd terms are infinite and also up to 10^6 assuming trajectories with more than 1000 odd terms are infinite. - Gary Detlefs, Jan 25 2022

Crossrefs

Programs

  • Maple
    f:= proc(m,b,n) if n mod 2 = 1 then return m*n+1 else return n/2 fi end proc
    F:= proc(m,b,n,i) option remember; if i=1 then return f(m,b,n) else return f(m,b,F(m,b,n,i-1)) fi end proc
    for x from 1 to 1000 do for y from 1 to 1000 do if F(5,1,x,y)= 86 or F(5,1,x,y)=26 then print(x): x=x+1; y:=y+1 fi od od
    # use print(x,y) to give the number of iterations needed to reach the cycle point
    # Gary Detlefs, Jan 25 2022

Extensions

Entry revised by N. J. A. Sloane, Jan 23 2016
a(16)-a(55) added by Gary Detlefs, Jan 25 2022

A267970 Conjectured list of positive numbers n such that the Collatz-like map T: if x > 1 and x odd, x -> 5*x+1 and, if x even, x -> x/2, when started at n, never repeats.

Original entry on oeis.org

7, 9, 11, 14, 18, 21, 22, 23, 25, 28, 29, 31, 35, 36, 37, 39, 41, 42, 44, 45, 46, 47, 49, 50, 53, 55, 56, 57, 58, 59, 61, 62, 63, 67, 69, 70, 71, 72, 73, 74, 75, 77, 78, 79, 81, 82, 84, 85, 87, 88, 89, 90, 91, 92, 93, 94, 95, 98, 99, 100
Offset: 1

Views

Author

Michel Lagneau, Jan 19 2016

Keywords

Comments

A companion to A232711 (which conjecturally lists the numbers whose trajectory reaches 1), and A267969 (which conjecturally lists the numbers whose trajectory cycles but does not reach 1). The present sequences conjecturally lists the numbers whose trajectory diverges.
The numbers listed have trajectories that have not cycled in the first 10000 steps. (But this does not mean much. There are quite small numbers that only cycle after 13000 iterations. - N. J. A. Sloane, Jan 23 2016)

Crossrefs

Extensions

Entry revised by N. J. A. Sloane, Jan 23 2016
a(18)-a(60) from Dmitry Kamenetsky, Jun 26 2024

A267703 Conjectured list of numbers whose trajectory under the '7x+1' map eventually reaches 1.

Original entry on oeis.org

1, 2, 4, 5, 8, 9, 10, 16, 18, 20, 32, 36, 40, 41, 64, 72, 73, 80, 82, 128, 144, 146, 160, 164, 167, 256, 288, 292, 320, 328, 329, 334, 512, 576, 584, 585, 640, 656, 658, 668, 1024, 1152, 1168, 1170, 1280, 1312, 1316, 1336, 1337, 1965, 2048, 2304, 2336, 2340, 2560
Offset: 1

Views

Author

Michel Lagneau, Jan 19 2016

Keywords

Comments

This is conjectural in that there is no known proof that the missing numbers 3, 6, 7, ... are really missing. It may be that after a very large number of iterations they will cycle. - N. J. A. Sloane, Jan 23 2016
Note that the computer program does not actually calculate a complete list of "numbers k such that the Collatz-like map T: if x odd, x -> 7*x+1 and if x even, x -> x/2, when started at k, eventually reaches 1".

Examples

			5 is in the sequence because the trajectory of 5 is 5 -> 36 -> 18 -> 9 -> 64 -> 32 -> 16 -> 8 -> 4 -> 2 -> 1.
		

Crossrefs

Programs

  • Maple
    nn:=10000:
    for n from 1 to 2340 do:
      m:=n:cyc:={n}:
        for i from 1 to nn do:
         if irem(m,2)=0
          then
           m:=m/2:
          else
          m:=7*m+1:
         fi:
        cyc:=cyc union {m}:
        od:
        n0:=nops(cyc):
        if n0N. J. A. Sloane, Jan 23 2016)

Extensions

Entry revised by N. J. A. Sloane, Jan 23 2016
a(19)-a(55) from Dmitry Kamenetsky, Jun 24 2024

A268338 Numbers that cycle under the following transformation: if m is even, divide by 2, if m is congruent to 1 mod 4, multiply by 3 and add 1; if m is congruent to 3 mod 4, multiply by 7 and add 1.

Original entry on oeis.org

1, 2, 4, 19, 23, 31, 38, 41, 46
Offset: 1

Views

Author

David Consiglio, Jr., Feb 01 2016

Keywords

Comments

Some numbers appear to grow indefinitely under these rules, but it is possible that they may eventually cycle at some point. All numbers up to 50 either cycle or transform to another number that cycles (typically 1). 51 is the first open case: it may eventually cycle or may continue to grow indefinitely.

Examples

			23 is a member of this sequence. 23 is congruent to 3 mod 4.  As a result, 23 transforms to 23*7+1 = 162.  From there 162 -> 81 -> 244 -> 122 -> 61 -> 184 -> 92 -> 46 -> 23.  23 is the least member of this cycle.
49 is not a member of this sequence because it eventually reduces to 19, which cycles.
		

Crossrefs

Programs

  • Python
    a = 1
    b = 1
    prev = []
    keep = []
    count = 0
    while b < 51:
        keep.append(a)
        flag1 = False
        flag2 = False
        if a % 2 == 0:
            a /= 2
        elif a % 4 == 1:
            a = a*3+1
        else:
            a = a*7+1
        if count > 50:
            b += 1
            a = b
            count = 0
            keep = []
        if keep.count(a) == 2 and a not in prev and a <= 50:
            prev.append(a)
            count = 0
            keep = []
            b += 1
            a = b
        count += 1
    print(sorted(prev))
    # David Consiglio, Jr., Feb 01 2016

Extensions

Corrected and edited by David Consiglio, Jr., Apr 20 2016

A270968 Reduced 5x+1 function R applied to the odd integers: a(n) = R(2n-1), where R(k) = (5k+1)/2^r, with r as large as possible.

Original entry on oeis.org

3, 1, 13, 9, 23, 7, 33, 19, 43, 3, 53, 29, 63, 17, 73, 39, 83, 11, 93, 49, 103, 27, 113, 59, 123, 1, 133, 69, 143, 37, 153, 79, 163, 21, 173, 89, 183, 47, 193, 99, 203, 13, 213, 109, 223, 57, 233, 119, 243, 31, 253, 129, 263, 67, 273, 139, 283, 9, 293, 149, 303
Offset: 1

Views

Author

Michel Lagneau, Mar 27 2016

Keywords

Comments

The odd-indexed terms a(2i+1) = 10i+3 = A017305(i), i>=0;
a(4i+4) = 10i+9 = A017377(i), i>=0;
a(8i+6) = 10i+7 = A017353(i), i>=0;
a(16i+2) = 10i+1 = A017281(i), i>=0.
Note that a(n) = a(16n-6) = a(6n-2)/3. No multiple of 5 is in this sequence.
a(n) = R(2n-1) < 2n-1 for n = 2, 6, 10, ..., 2+4i,...

Examples

			a(4)=9 because (2*4-1) = 7  -> (5*7+1)/2^2 = 9.
		

Crossrefs

Programs

  • Mathematica
    nextOddK[n_] := Module[{m=5n+1}, While[EvenQ[m], m=m/2]; m]; (* assumes odd n *) Table[nextOddK[n], {n, 1, 200, 2}]
  • PARI
    a(n) = my(m = 2*n-1, c = 5*m+1); c/2^valuation(c, 2); \\ Michel Marcus, Mar 27 2016

Formula

a(n) = A000265(A017341(n-1)). - Michel Marcus, Mar 27 2016
Showing 1-5 of 5 results.