A192504 -id:A192504 - cp's OEIS frontend

A003309 Ludic numbers: apply the same sieve as Eratosthenes, but cross off every k-th remaining number.

1, 2, 3, 5, 7, 11, 13, 17, 23, 25, 29, 37, 41, 43, 47, 53, 61, 67, 71, 77, 83, 89, 91, 97, 107, 115, 119, 121, 127, 131, 143, 149, 157, 161, 173, 175, 179, 181, 193, 209, 211, 221, 223, 227, 233, 235, 239, 247, 257, 265, 277, 283, 287, 301, 307, 313

Offset: 1

N. J. A. Sloane

The definition can obviously only be applied from k = a(2) = 2 on: for k = 1, all remaining numbers would be deleted. - M. F. Hasler, Nov 02 2024

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Donovan Johnson, Table of n, a(n) for n = 1..100000
David Applegate, C program for A003309.
Donovan Johnson, Ludic numbers computed up to A003309(1236290) = 23000711.
OEIS Wiki, Ludic numbers.
Popular Computing (Calabasas, CA), Sieves: Problem 43, Vol. 2 (No. 13, Apr 1974), pp. 6-7. This is Sieve #1. [Annotated and scanned copy]
Rosettacode Wiki, Ludic numbers.
Index entries for sequences generated by sieves

Without the initial 1 occurs as the leftmost column in arrays A255127 and A260717.
Cf. A003310, A003311, A100464, A100585, A100586 (variants).
Cf. A192503 (primes in sequence), A192504 (nonprimes), A192512 (number of terms <= n).
Cf. A192490 (characteristic function).
Cf. A192607 (complement).
Cf. A260723 (first differences).
Cf. A255420 (iterates of f(n) = A003309(n+1) starting from n=1).
Subsequence of A302036.
Cf. A237056, A237126, A237427, A235491, A255407, A255408, A255421, A255422, A260435, A260436, A260741, A260742 (permutations constructed from Ludic numbers).
Cf. also A000959, A008578, A255324, A254100, A272565 (Ludic factor of n), A297158, A302032, A302038.
Cf. A376237 (ludic factorial: cumulative product), A376236 (ludic Fortunate numbers).

a003309 n = a003309_list !! (n - 1)
a003309_list = 1 : f [2..] :: [Int]
   where f (x:xs) = x : f (map snd [(u, v) | (u, v) <- zip [1..] xs,
                                             mod u x > 0])
-- Reinhard Zumkeller, Feb 10 2014, Jul 03 2011

ludic:= proc(N) local i, k,S,R;
  S:= {$2..N};
  R:= 1;
  while nops(S) > 0 do
    k:= S[1];
    R:= R,k;
    S:= subsop(seq(1+k*j=NULL, j=0..floor((nops(S)-1)/k)),S);
  od:
[R];
end proc:
ludic(1000); # Robert Israel, Feb 23 2015

t = Range[2, 400]; r = {1}; While[Length[t] > 0, k = First[t]; AppendTo[r, k]; t = Drop[t, {1, -1, k}];]; r (* Ray Chandler, Dec 02 2004 *)

t=vector(399,x,x+1); r=[1]; while(length(t)>0, k=t[1];r=concat(r,[k]);t=vector((length(t)*(k-1))\k,x,t[(x*k+k-2)\(k-1)])); r \\ Phil Carmody, Feb 07 2007

A3309=[1]; next_A003309(n)=nn && break); n+!if(n=setsearch(A3309,n+1,1),return(A3309[n])) \\ Should be made more efficient if n >> max(A3309). - M. F. Hasler, Nov 02 2024
{A003309(n) = while(n>#A3309, next_A003309(A3309[#A3309])); A3309[n]} \\ Should be made more efficient in case n >> #A3309. - M. F. Hasler, Nov 03 2024

upto(nn)= my(r=List([1..nn]), p=1); while(p++<#r, my(k=r[p], i=p); while((i+=k)<=#r, listpop(~r, i); i--)); Vec(r); \\ Ruud H.G. van Tol, Dec 13 2024

remainders = [0]
ludics = [2]
N_MAX = 313
for i in range(3, N_MAX) :
    ludic_index = 0
    while ludic_index < len(ludics) :
        ludic = ludics[ludic_index]
        remainder = remainders[ludic_index]
        remainders[ludic_index] = (remainder + 1) % ludic
        if remainders[ludic_index] == 0 :
            break
        ludic_index += 1
    if ludic_index == len(ludics) :
        remainders.append(0)
        ludics.append(i)
ludics = [1] + ludics
print(ludics)
# Alexandre Herrera, Aug 10 2023

def A003309(): # generator of the infinite list of ludic numbers
    L = [2, 3]; yield 1; yield 2; yield 3
    while k := len(L)//2: # could take min{k | k >= L[-1-k]-1}
        for j in L[-1-k::-1]: k += 1 + k//(j-1)
        L.append(k+2); yield k+2
A003309_upto = lambda N=99: [t for t,_ in zip(A003309(),range(N))]
# M. F. Hasler, Nov 02 2024

(define (A003309 n) (if (= 1 n) n (A255127bi (- n 1) 1))) ;; Code for A255127bi given in A255127.
;; Antti Karttunen, Feb 23 2015

Complement of A192607; A192490(a(n)) = 1. - Reinhard Zumkeller, Jul 05 2011
From Antti Karttunen, Feb 23 2015: (Start)
a(n) = A255407(A008578(n)).
a(n) = A008578(n) + A255324(n).
(End)

More terms from David Applegate and N. J. A. Sloane, Nov 23 2004

A192506 Numbers that are neither ludic nor prime.

Original entry on oeis.org

4, 6, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 22, 24, 26, 27, 28, 30, 32, 33, 34, 35, 36, 38, 39, 40, 42, 44, 45, 46, 48, 49, 50, 51, 52, 54, 55, 56, 57, 58, 60, 62, 63, 64, 65, 66, 68, 69, 70, 72, 74, 75, 76, 78, 80, 81, 82, 84, 85, 86, 87, 88, 90, 92, 93, 94

Offset: 1

Reinhard Zumkeller, Jul 05 2011

nonn

Intersection of A002808 and A192607; (1-A010051(a(n)))*(1-A192490(a(n)))=1;
a(n) = A091212(n) for n <= 60.
a(n) = A175526(n) for n <= 53. - Reinhard Zumkeller, Jul 12 2011
In other words, composite numbers that are nonludic. - Antti Karttunen, May 11 2015

Reinhard Zumkeller (first 1000 terms) & Antti Karttunen, Table of n, a(n) for n = 1..10000

Cf. A257689 (complement, either ludic or prime), A192503 (ludic and prime), A192504 (ludic and nonprime), A192505 (nonludic and prime).
Cf. A010051, A192490, A002808, A192607.
a(n) differs from A091212(n) and A205783(n+1) for the first time at n=37, where a(37) = 55, while 55 is missing from both A091212 and A205783.
Differs from A175526 for the first time at n=54, where a(54) = 78, while A175526(54) = 77, a term which is missing from here.

a192506 n = a192506_list !! (n-1)
a192506_list = filter ((== 0) . a010051) a192607_list
(Scheme, with Antti Karttunen's IntSeq-library)
(define A192506 (MATCHING-POS 1 1 (lambda (n) (and (zero? (A192490 n)) (zero? (A010051 n))))))
;; Antti Karttunen, May 07 2015

a3309[nmax_] := a3309[nmax] = Module[{t = Range[2, nmax], k, r = {1}}, While[Length[t] > 0, k = First[t]; AppendTo[r, k]; t = Drop[t, {1, -1, k}]]; r];
ludicQ[n_, nmax_] /; 1 <= n <= nmax := MemberQ[a3309[nmax], n];
terms = 1000;
f[nmax_] := f[nmax] = Select[Range[nmax], !ludicQ[#, nmax] && !PrimeQ[#]&] // PadRight[#, terms]&;
f[nmax = terms];
f[nmax = 2 nmax];
While[f[nmax] != f[nmax/2], nmax = 2 nmax];
seq = f[nmax] (* Jean-François Alcover, Dec 10 2021, after Ray Chandler in A003309 *)

A257689 Numbers that are either ludic or prime.

Original entry on oeis.org

1, 2, 3, 5, 7, 11, 13, 17, 19, 23, 25, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 77, 79, 83, 89, 91, 97, 101, 103, 107, 109, 113, 115, 119, 121, 127, 131, 137, 139, 143, 149, 151, 157, 161, 163, 167, 173, 175, 179, 181, 191, 193, 197, 199, 209, 211, 221, 223, 227, 229, 233, 235, 239, 241, 247, 251, 257, 263, 265

Offset: 1

Antti Karttunen, May 07 2015

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10181

Union of primes (A000040) and ludic numbers (A003309).
Cf. A192506 (complement, neither ludic nor prime), A192503 (ludic and prime), A192504 (ludic and nonprime), A192505 (nonludic and prime).
Differs from A206074(n-1), A186891(n) and A257688(n) for the first time at n=19, where a(19) = 59, while A206074(18) = A186891(19) = A257688(19) = 55, a term missing from here.
Differs from A257691 for the first time at n=24, where a(24) = 77, while A257691(24) = 79.

a3309[nmax_] := a3309[nmax] = Module[{t = Range[2, nmax], k, r = {1}}, While[Length[t] > 0, k = First[t]; AppendTo[r, k]; t = Drop[t, {1, -1, k}]]; r];
ludicQ[n_, nmax_] /; 1 <= n <= nmax := MemberQ[a3309[nmax], n];
terms = 1000;
f[nmax_] := f[nmax] = Select[Range[nmax], ludicQ[#, nmax] || PrimeQ[#]&] // PadRight[#, terms]&;
f[nmax = terms];
f[nmax = 2 nmax];
While[f[nmax] != f[nmax/2], nmax = 2 nmax];
seq = f[nmax] (* Jean-François Alcover, Dec 10 2021, after Ray Chandler in A003309 *)

A276448 Numbers n for which A272565(n) > A020639(n).

Original entry on oeis.org

25, 55, 77, 85, 91, 115, 119, 121, 133, 143, 145, 161, 175, 187, 203, 205, 209, 221, 235, 247, 253, 265, 287, 295, 301, 325, 329, 341, 343, 355, 361, 371, 377, 385, 403, 407, 413, 415, 437, 445, 473, 475, 481, 493, 497, 505, 511, 517, 527, 535, 539, 551, 553, 565, 581, 595, 623, 625, 655, 667, 671, 685, 697, 703, 707, 713, 715, 721

Offset: 1

Antti Karttunen, Sep 11 2016

nonn

			85 is the fourth number of the fourth row of A255127, which starts with 7: 7, 31, 59, 85, ..., thus A272565(85)=7, while on A083221 it occurs as the sixth term on the third row that starts with 5, thus A020639(85)=5, and 7 > 5, thus 85 is included in this sequence.

Antti Karttunen, Table of n, a(n) for n = 1..2000

Cf. A276447 (complement in A276437), A276347.

Cf. A020639, A272565, A192504, A083221, A255127.

cp's OEIS Frontend

A003309 Ludic numbers: apply the same sieve as Eratosthenes, but cross off every k-th remaining number.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

Python

Python

Scheme

Formula

Extensions

A192506 Numbers that are neither ludic nor prime.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

A257689 Numbers that are either ludic or prime.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

A276448 Numbers n for which A272565(n) > A020639(n).

Views

Author

Keywords

Examples

Links

Crossrefs