cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A152009 (L)-sieve transform of {1,4,7,10,...,3n-2,...} (A016777).

Original entry on oeis.org

1, 3, 6, 10, 16, 25, 39, 60, 91, 138, 208, 313, 471, 708, 1063, 1596, 2395, 3594, 5392, 8089, 12135, 18204, 27307, 40962, 61444, 92167
Offset: 1

Views

Author

John W. Layman, Nov 19 2008

Keywords

Comments

The (L)-sieve transform of the sequence {a(n)} of positive integers is defined as follows: Denote the sequence of natural numbers by N. Remove the first term of N, which we denote by s(1) and then from the resulting sequence delete all terms whose index is a term of {a(n)}, to obtain the sequence N'.
Then remove the first term of N', denoted by s(2) and then from the resulting sequence delete all terms whose index is a term of {a(n)}, to obtain N''. Repeat this process indefinitely to obtain the transform LST({a(n)}) = {s(1), s(2),...}, the sequence of initial terms removed at each stage.
The (L)-sieve transform is quite different from the transform introduced by N. J. A. Sloane in A099361 and used by T. D. Noe in A100424 - A100426 and seems to lead to more interesting results and relationships among sequences. An interesting property of the (L)-sieve transform is that the (L)-sieve transform of the sequence {1,3,6,10,...,n(n+1)/2,...} of triangular numbers is again the triangular numbers. Another (conjectured) connection with the triangular numbers is given in the following.
Conjecture. Let x(0) be a random sequence of positive integers and, for n>0, let x(n)=S[x(n-1)], where S is the (L)-sieve transform. Then the limit of {x(n)} as n goes to infinity is the sequence of triangular numbers {1,3,6,10,...,n(n+1)/2,...}.
Illustration of the conjecture:
x(0)={3,8,12,14,18,22,25,31,34,39,42,45,...} (A random initial sequence.)
x(1)={1,2,3,5,7,10,14,20,28,38,51,69,...}
x(2)={1,5,12,20,30,41,53,65,78,91,105,119,...}
x(3)={1,3,5,8,11,15,19,24,29,35,41,48,...}
x(4)={1,3,7,13,21,31,43,56,71,88,107,127,...}
x(5)={1,3,6,10,15,20,26,33,40,48,56,65,...}
x(6)={1,3,6,10,15,22,30,39,50,62,75,90,...}
x(7)={1,3,6,10,15,21,28,36,45,55,66,78,...} ...
t={1,3,6,10,15,21,28,36,45,55,66,78,...} (Triangular numbers)

Programs

Formula

It appears that {a(n)} is given by a(n)=floor[(3*a(n-1)+3)/2], with a(1)=1.

A100424 A sieve transform applied three times to the positive integers.

Original entry on oeis.org

2, 3, 13, 37, 107, 139, 223, 251, 359, 397, 503, 647, 683, 857, 887, 1033, 1151, 1249, 1291, 1429, 1493, 1601, 1667, 1783, 1831, 2003, 2053, 2267, 2377, 2459, 2593, 2677, 2753, 2801, 2903, 3119, 3209, 3347, 3461, 3557, 3607, 3727, 3851, 4079, 4139, 4243
Offset: 1

Views

Author

T. D. Noe, Nov 19 2004

Keywords

Comments

The process is described in A099361. The first application of the sieve transform produces the prime numbers A000040. The second application yields sequence A099361. In principle, the sieve transform can be applied to any sequence of positive integers. For instance, the sieve transform of the positive even numbers is 2, 4, 8, 16,.... Also note that the transform can produce a finite sequence. See A100425 and A100426 for more examples.

Crossrefs

Programs

  • Mathematica
    SieveTransform[b_List] := Module[{d, nn=Length[b], a=b}, Do[d=a[[i]]; If[d>1, Do[a[[j]]=-1, {j, i+d, nn, d}]], {i, nn}]; DeleteCases[a, -1]]; SieveTransform[SieveTransform[SieveTransform[Range[2, 5000]]]]

A100425 A sieve transform applied to the composite numbers.

Original entry on oeis.org

4, 6, 8, 9, 12, 14, 18, 21, 25, 26, 35, 36, 50, 52, 56, 66, 68, 74, 77, 95, 96, 98, 108, 118, 141, 143, 146, 162, 166, 176, 185, 186, 187, 203, 207, 228, 235, 245, 250, 254, 273, 275, 285, 290, 295, 302, 322, 340, 350, 361, 375, 380, 402, 404, 416, 417, 436, 445
Offset: 1

Views

Author

T. D. Noe, Nov 19 2004

Keywords

Crossrefs

Programs

  • Mathematica
    (* see A100424 for the definition of SieveTransform *) SieveTransform[Complement[Range[2, 1000], Prime[Range[PrimePi[1000]]]]]
Showing 1-3 of 3 results.