login
Lexicographically earliest infinite sequence of distinct positive terms such that, for any prime p, any run of consecutive multiples of p has length exactly 2.
55

%I #216 Jul 28 2024 13:10:24

%S 1,2,4,3,6,8,5,10,12,9,7,14,16,11,22,18,15,20,24,21,28,26,13,17,34,30,

%T 45,19,38,32,23,46,36,27,25,35,42,48,29,58,40,55,33,39,52,44,77,49,31,

%U 62,50,65,78,54,37,74,56,63,51,68,60,75,41,82,64,43,86

%N Lexicographically earliest infinite sequence of distinct positive terms such that, for any prime p, any run of consecutive multiples of p has length exactly 2.

%C In other words, each multiple of a prime p has exactly one neighbor that is also a multiple of p.

%C This sequence is similar to A280866; the first difference occurs at n=42: a(42)=55 whereas A280866(42)=50.

%C Conjectured to be a permutation of the positive integers.

%C Sometimes referred to as the "cup of coffee" sequence, since it feels as if just one more cup of coffee is all it would take to prove that this is indeed a permutation of the positive integers. - _N. J. A. Sloane_, Nov 04 2020

%C There are several short cycles, and apparently at least two infinite cycles. For a list see the attached file "Properties of A280864". - _N. J. A. Sloane_, Feb 03 2017

%C Properties (For proofs, see the attached file "Properties of A280864")

%C Theorem 1: This sequence contains every prime and every even number. (Added by _N. J. A. Sloane_, Jan 15 2017)

%C Theorem 2: The sequence contains infinitely many odd composite numbers. (Added by _N. J. A. Sloane_, Feb 14 2017)

%C Theorem 3: If p is an odd prime, the sequence contains infinitely many odd multiples of p. (Added by _N. J. A. Sloane_, Mar 12 2017, with corrected proof Apr 03 2017)

%C There are two types of primes in this sequence: Type I, the first time a term a(n) is divisible by p is when a(n)=p for some n; Type II, the first time a term a(n) is divisible by p is when a(n)=k*p for some n and some k>1 (the Type II primes are listed in A280745).

%C Conjecture 4: If a prime p divides a(n) then p <= n. - _N. J. A. Sloane_, Apr 07 2017 and Apr 16 2017

%C Theorem 5: The sequence is a permutation of the natural numbers iff it contains every square. - _N. J. A. Sloane_, Apr 14 2017

%C From _Bob Selcoe_, Apr 03 2017: (Start)

%C Define the "radical class" C_R to be the set of numbers which have the same radical R (or the same largest squarefree divisor - i.e., the same product of their prime factors). These are the columns in A284311. So for example C_10 is the set of numbers with radical 10 or prime factors {2,5}: {10, 20, 40, 50 80, 100, 160, ...).

%C If the sequence contains any members of C_R, then those members must appear in order; so for example, if 160 has appeared, {10, 20, 40, 50, 80} will have already appeared, in that order. Naturally, this holds for prime powers; for example, C_5: if 3125 has appeared, {5, 25, 125, 625} will have appeared earlier, in that order.

%C After seeing a(n), let S be smallest missing number (A280740) and let prime(G) be largest prime already appearing in the sequence. Conjecture: Prime(G) < S <= prime(G+1), and a(35) = 25 = S is the only nonprime S term (following a(31) = 23, preceding a(39) = 29). (End)

%H N. J. A. Sloane, <a href="/A280864/b280864.txt">Table of n, a(n) for n = 1..100000</a> (First 10000 terms from Rémy Sigrist)

%H Dana G. Korssjoen, Biyao Li, Stefan Steinerberger, Raghavendra Tripathi, and Ruimin Zhang, <a href="https://arxiv.org/abs/2012.04625">Finding structure in sequences of real numbers via graph theory: a problem list</a>, arXiv:2012.04625, Dec 08, 2020

%H Rémy Sigrist, <a href="/A280864/a280864.gp.txt">PARI program for A280864</a>

%H N. J. A. Sloane, <a href="/A280864/a280864_5.txt">Properties of A280864</a> [Revised, Apr 25 2017]

%H N. J. A. Sloane, <a href="/A280864/a280864.txt">Table of n, a(n) for n = 1..1000000</a>, computed using Sigrist's PARI program.

%H N. J. A. Sloane, <a href="/A195264/a195264.pdf">Confessions of a Sequence Addict (AofA2017)</a>, slides of invited talk given at AofA 2017, Jun 19 2017, Princeton. Mentions this sequence.

%e The first terms, alongside their required and forbidden prime factors are:

%e n a(n) Required Forbidden

%e -- ---- -------- ---------

%e 1 1 none none

%e 2 2 none none

%e 3 4 2 none

%e 4 3 none 2

%e 5 6 3 none

%e 6 8 2 3

%e 7 5 none 2

%e 8 10 5 none

%e 9 12 2 5

%e 10 9 3 2

%e 11 7 none 3

%e 12 14 7 none

%e 13 16 2 7

%e 14 11 none 2

%e 15 22 11 none

%e 16 18 2 11

%e 17 15 3 2

%e 18 20 5 3

%e 19 24 2 5

%e 20 21 3 2

%e 21 28 7 3

%e 22 26 2 7

%e 23 13 13 2

%e 24 17 none 13

%e 25 34 17 none

%e 26 30 2 17

%e 27 45 3, 5 2

%e 28 19 none 3, 5

%e 29 38 19 none

%e 30 32 2 19

%e 31 23 none 2

%e 32 46 23 none

%e 33 36 2 23

%e 34 27 3 2

%e 35 25 none 3

%e 36 35 5 none

%e 37 42 7 5

%e 38 48 2, 3 7

%e 39 29 none 2, 3

%e 40 58 29 none

%e 41 40 2 29

%e 42 55 5 2

%p N:= 1000: # to get all terms until the first term > N

%p A[1]:= 1:

%p A[2]:= 2:

%p G:= {}:

%p Avail:= [$3..N]:

%p found:= true:

%p lastn:= 2:

%p for n from 3 while found and nops(Avail)>0 do

%p found:= false;

%p H:= G;

%p G:= numtheory:-factorset(A[n-1]);

%p r:= convert(G minus H,`*`);

%p s:= convert(G intersect H, `*`);

%p for j from 1 to nops(Avail) do

%p if Avail[j] mod r = 0 and igcd(Avail[j],s) = 1 then

%p found:= true;

%p A[n]:= Avail[j];

%p Avail:= subsop(j=NULL,Avail);

%p lastn:= n;

%p break

%p fi

%p od;

%p od:

%p seq(A[i],i=1..lastn); # _Robert Israel_, Mar 22 2017

%t terms = 100;

%t rad[n_] := Times @@ FactorInteger[n][[All, 1]];

%t A280864 = Reap[present = 0; p = 1; pp = 1; Do[forbidden = GCD[p, pp]; mandatory = p/forbidden; a = mandatory; While[BitGet[present, a] > 0 || GCD[forbidden, a] > 1, a += mandatory]; Sow[a]; present += 2^a; pp = p; p = rad[a], terms]][[2, 1]] (* _Jean-François Alcover_, Nov 23 2017, translated from _Rémy Sigrist_'s PARI program *)

%Y See A280738, A280740, A280741 (inverse), A280742, A280743, A280744, A280745, A280746, A280755, A280770, A280771, A280773, A280774, A283832, A284724, A284725, A284726, A284785, A285181 for various subsidiary sequences.

%Y A280754 gives fixed points.

%Y Cf. A280866.

%Y In the same spirit as A064413 and A098550.

%Y A338338, A338444, and A375029 are variants.

%Y A373797 is a finite version.

%K nonn,nice

%O 1,2

%A _Rémy Sigrist_, Jan 09 2017

%E Added "infinite" to definition. - _N. J. A. Sloane_, Sep 28 2019