login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

The Judd Trump's "infinite plant" sequence: prime numbers become nonprime numbers by striking the cue ball 2 with a cue stick to the right (see the Comments section).
1

%I #10 Dec 14 2020 05:19:02

%S 2,11,3,23,29,5,13,31,7,41,43,37,47,53,59,17,61,67,71,83,89,19,97,73,

%T 101,103,107,79,109,127,113,149,131,151,157,137,139,163,167,173,181,

%U 179,211,191,193,197,223,199,227,229,233,239,241,251,257,263,269,271,281,283,277,293,307,311,331,337,313,347

%N The Judd Trump's "infinite plant" sequence: prime numbers become nonprime numbers by striking the cue ball 2 with a cue stick to the right (see the Comments section).

%C There is a non-snooker description of this sequence: first erase all spaces between terms; then move every comma 1 position to the left; the new sequence is now made by nonprimes only (with duplicates, sometimes); the starting sequence (this one) is the lexicographically earliest with this property that has no duplicates and no nonprimes.

%e Striking 2 to the right pushes 2 against 11;

%e the last digit of 11 is then pushed against 3 (leaving 21 behind - a nonprime);

%e the last digit of 3 is then pushed against 23 (leaving 1 behind - a nonprime);

%e the last digit of 23 is then pushed against 29 (leaving 32 behind - a nonprime);

%e the last digit of 29 is then pushed against 5 (leaving 32 behind - a nonprime);

%e the last digit of 5 is then pushed against 13 (leaving 9 behind - a nonprime);

%e etc.

%e This is the lexicographically earliest sequence of distinct positive terms with this property.

%o (Python)

%o from sympy import isprime

%o def aupto(n):

%o alst, used, strakm1 = [0, 2], {2}, "2"

%o for k in range(2, n+1):

%o ball = (str(alst[k-1]))[-1]

%o ak = 1

%o ball_left = ball + (str(ak))[:-1]

%o while isprime(int(ball_left)) or ak in used or not isprime(ak):

%o ak += 2 # continue to only test odds

%o ball_left = ball + (str(ak))[:-1]

%o alst.append(ak)

%o used.add(ak)

%o return alst[1:] # use alst[n] for a(n) function

%o print(aupto(70)) # _Michael S. Branicky_, Dec 11 2020

%Y Cf. A339467 (the Ronnie O'Sullivan sequence).

%K base,nonn

%O 1,1

%A _Eric Angelini_ and _Carole Dubois_, Dec 10 2020