login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A323712 Number of the following described shuffles required to return a deck of n cards to its original state. Create two piles by alternating a card from the top of the deck left-right-left-right until the deck is exhausted. Then, placing the left pile on top of the right pile constitutes one shuffle. 1
1, 1, 3, 4, 4, 6, 6, 3, 9, 5, 5, 12, 12, 4, 12, 8, 8, 9, 9, 6, 6, 22, 22, 20, 20, 9, 27, 28, 28, 10, 10, 5, 15, 12, 12, 36, 36, 12, 12, 20, 20, 7, 7, 12, 36, 46, 46, 42, 42, 8, 24, 52, 52, 20, 20, 9, 9, 29, 29, 60, 60, 6, 18, 12, 12, 33, 33, 22, 66, 70, 70, 18 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,3

COMMENTS

The card shuffling procedure is the same for even n and odd n.

Here are a few conjectures.

a(n) <= n for all n.

a(p)=a(p-1) and a(p)|p-1 when p is a prime >= 5.

a(n)=a(n-1) and a(n)|n-1 for nonprimes 341=31*11 and 22369621=8191*2731 and probably other pseudoprimes of the form p*((p+2)/3) where p is a Mersenne prime and (p+2)/3 is prime.

n cards are returned to their original state after n shuffles when n=1, 3, 4, 6, 9, 12, 22, 27, 28, 36, 46, 52, 60, 70, 78, 81, ... . These values of n are either of the form p-1 where p is an odd prime number or 3^i, i >= 0.

a(c) is relatively small (compared with nearby values) when c is a Catalan number.

a(2n+1)=3*a(2n) or a(2n+1)=a(2n) for all n.

LINKS

Alois P. Heinz, Table of n, a(n) for n = 1..1024

FORMULA

a(2^m) = m if m is odd, a(2^m) = 2m if m is even. - Alois P. Heinz, Feb 15 2019

EXAMPLE

For n=4, {a1,a2,a3,a4}-->{a3,a1,a4,a2}-->{a4,a3,a2,a1}-->{a2,a4,a1,a3}-->{a1,a2,a3,a4}, so a(4)=4.

For n=5, {a1,a2,a3,a4,a5}-->{a5,a3,a1,a4,a2}-->{a2,a1,a5,a4,a3}-->{a3,a5,a2,a4,a1}-->{a1,a2,a3,a4,a5}, so a(5)=4.

PROG

(PARI) perm(n, vn) = {my(va = List(), vb = List()); for (k=1, n, if (k % 2, listput(va, vn[k]), listput(vb, vn[k])); ); Vec(concat(Vecrev(va), Vecrev(vb))); }

a(n) = {my(vn = vector(n, k, k), vs = perm(n, vn), nb = 1); while (vs != vn, vs = perm(n, vs); nb++); nb; } \\ Michel Marcus, Feb 06 2019

CROSSREFS

Cf. A024222, A022998, A163776.

Sequence in context: A309530 A061117 A255171 * A215250 A229022 A014683

Adjacent sequences:  A323709 A323710 A323711 * A323713 A323714 A323715

KEYWORD

nonn

AUTHOR

David Lovler, Jan 24 2019

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 12 17:17 EDT 2020. Contains 336439 sequences. (Running on oeis4.)