Modern Computing for Mathematicians (Math 587)
A generalization of the Syracuse iterator

February 24, 2009

Recall that the Syracuse function s is defined for integers n by sn=1ifn13n+1ifn>1is oddn2ifn>1is even

Many questions about the iterative behavior of s – in particular, the question of whether the value of some iterate from any starting point is 1 – are unaffected if s is replaced with the function s1 defined as follows: s1n=1ifn13n+1ifn>1is oddn2kifn=2kmwheremis odd andk1 or, as one might more informally write, for n1, s1n=3n+1made coprime to2.

Generalizing this, for given pairwise coprime integers a,b,m>0 with m2, one defines for n1 fabmn=an+bmade coprime tom. Here, for an integer x, the phrase “x made coprime to m” means that for any common prime divisor p of x and m the highest power of p dividing x is removed as a factor. Note that the meaning of fabm is not changed when m is replaced by the product of the distinct primes dividing m; that is, without loss of generality one may restrict to the case where m is square-free.

Example: s1=f312.

Exercises:

  1. Write code for gp to investigate the iterates of fabm from a given integer. In particular, the code should be able to determine whether from a given starting integer n a cycle is formed within the first N iterates.

  2. Determine what cycles, if any, are formed and whether there seems to be a pattern of unbounded growth in the iterates for n,N10000 in the following cases:

    1. f325.

    2. f513.

    3. f17130030.