A384991 Order of the permutation of [n] formed by a Josephus elimination variation: take p, skip 1, with p starting at 2 and advancing to the next prime after each skip.
1, 1, 2, 3, 3, 5, 4, 7, 8, 15, 10, 11, 12, 13, 45, 15, 105, 17, 77, 19, 24, 21, 117, 23, 504, 255, 26, 165, 28, 440, 60, 31, 442, 33, 1386, 805, 154, 37, 105, 39, 1020, 216, 208, 43, 40, 45, 2860, 1953, 90, 49, 45, 51, 1092, 120, 184, 55, 56, 150, 58, 6045
Offset: 1
Keywords
Examples
For n=15, the rotations to construct the permutation are 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 \----------------------------------------------/ 1st rotation (p=2) 1, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 2 \----------------------------------------/ 2nd rotation (p=3) 1, 3, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 2, 5 \---------------------------/ 3rd rotation (p=5) 1, 3, 4, 6, 7, 8, 9, 11, 12, 13, 14, 15, 2, 5, 10 \-----/ 4th rotation (p=7) 1, 3, 4, 6, 7, 8, 9, 11, 12, 13, 14, 15, 2, 10, 5 The 4th rotate is an example of an element (5) which was previously rotated to the end, being rotated to the end again. This final permutation has order a(15) = 45 (applying it 45 times reaches the identity permutation again).
Links
- Chuck Seggelin, Table of n, a(n) for n = 1..1000
Crossrefs
Programs
-
Python
from sympy.combinatorics import Permutation from sympy import isprime, prime def apply_transformation(seq): k = 1 p = prime(k) i = p - 1 while i < len(seq): seq.append(seq.pop(i)) k += 1 p = prime(k) i += (p-1) return seq def a(n): seq = list(range(n)) p = apply_transformation(seq.copy()) return Permutation(p).order()
Comments