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.

Previous Showing 31-35 of 35 results.

A323161 Lexicographically earliest sequence such that a(i) = a(j) => f(i) = f(j) for all i, j, where f(n<=3) = -n, f(n) = 0 if n-1 is an odd prime, and f(n) = floor((n-1)/2) for all other numbers.

Original entry on oeis.org

1, 2, 3, 4, 5, 4, 6, 4, 7, 7, 8, 4, 9, 4, 10, 10, 11, 4, 12, 4, 13, 13, 14, 4, 15, 15, 16, 16, 17, 4, 18, 4, 19, 19, 20, 20, 21, 4, 22, 22, 23, 4, 24, 4, 25, 25, 26, 4, 27, 27, 28, 28, 29, 4, 30, 30, 31, 31, 32, 4, 33, 4, 34, 34, 35, 35, 36, 4, 37, 37, 38, 4, 39, 4, 40, 40, 41, 41, 42, 4, 43, 43, 44, 4, 45, 45, 46, 46, 47, 4, 48, 48, 49, 49, 50, 50, 51, 4, 52, 52, 53, 4, 54, 4
Offset: 1

Views

Author

Antti Karttunen, Jan 06 2019

Keywords

Comments

For all i, j: A322809(i) = A322809(j) <=> a(i+1) = a(j+1).
For all i, j: a(i) = a(j) => b(i) = b(j), where b can be, but is not limited to, any of the following sequences: A029834, A049084, A062590, A063377, A064891, A078442 (A049076), A175663, A175682, A269668, A292936, A323162, many of which are related to counting primes in certain kinds of chains or iterations.
Why does this work? Consider the function f given in the definition: based on its properties, we can deduce from the value of f(n) the following information about n:
(A) If f(n) = -2, then n is 2, the only even prime,
(B) If f(n) = -3, then n is 3, the first odd prime,
(C) If f(n) is zero, then n is an even composite preceded by a prime, but we don't know which even composite exactly,
(D) If f(n) > 0 and f(1+2*f(n)) = f(2+2*f(n)), then n is either (D1) an odd composite number, or (D2) an even composite number preceded by an odd composite number, and the said composite number in both cases is 1 + 2*f(n),
(E) If f(n) > 0 and f(1+2*f(n)) <> f(2+2*f(n)), then n is an odd prime > 3, specifically, 1 + 2*f(n).
As this sequence is a restricted growth sequence transform of the said function f, we have a(i) = a(j) <=> f(i) = f(j) for all i, j, thus, even without knowing the value of n, but just a(n), we can find the value of f(n) by searching for the minimal k such that a(k) = a(n), then compute f(k) with that k. Furthermore, any function g defined as g(n) = h(f(n)) [where h is any function], clearly satisfies
a(i) = a(j) => g(i) = g(j), for all i, j. [*]
For instances of such functions g, we can consider many sequences like those sequences b(n) listed above, that have g(n) = 0 for all composite numbers, and g(p) > 0 for all primes p. This is usually the pattern, but there are exceptions, like A323162, which is the characteristic function of A005381, composites n such that n-1 is also composite. These are precisely the numbers that occur twice in this sequence, while all other numbers (including primes), occur just once, that is, reside in their own singular equivalence classes. Thus, it is not guaranteed that all sequences g matching to this sequence (i.e. those satisfying the implication *), even if not false positives in strict sense, would necessarily have some consistent relation to primes, instead, they might contain any random values at the positions given by A093515. However, in the current OEIS, such sequences are exceedingly rare.

Crossrefs

Cf. A005381 (numbers that occur twice in this sequence), A093515 (numbers > 1 that occur just once).
Cf. A010051, A029834, A049076, A049084, A062590, A063377, A064891, A078442, A175663, A175682, A269668, A292936, A323162 (some of the matched sequences).

Programs

  • PARI
    up_to = 10000;
    rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
    A323161aux(n) = if(n<=3,-n,if(isprime(n-1),0,((n-1)>>1))); \\ This implements the function f of the definition.
    v323161 = rgs_transform(vector(up_to,n,A323161aux(n)));
    A323161(n) = v323161[n];

Formula

a(1) = 1; for n > 1, a(n) = 1 + A322809(n-1).

A121543 "If k appears then the k-th prime doesn't", with a(1)=1.

Original entry on oeis.org

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

Views

Author

Zak Seidov, Aug 06 2006

Keywords

Comments

Most prime numbers are absent; from first 1000 primes only 138 are terms: 3, 11, 17, 41, 67, 83, 109, 127, 157, 191, 211, 241, 277, 283, 353, ..., all with prime indices. Also, all nonprimes are terms. - Zak Seidov, Sep 12 2015
k is a term if and only if A078442(k) is even. - Pontus von Brömssen, Jul 25 2025

Crossrefs

Programs

  • Mathematica
    s={1};Do[If[ !PrimeQ[n]||(PrimeQ[n] && FreeQ[s,PrimePi[n]]),AppendTo[s,n]],{n,2,100}];s (* Zak Seidov *)

A373338 Characteristic function of A333242: a(n) = 1 if n is a term of A333242.

Original entry on oeis.org

0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0
Offset: 1

Views

Author

Michael P. May, Jun 01 2024

Keywords

Comments

This sequence is the result of applying the N-sieve to generate the prime number subsequence A333242 where 1 indicates a prime number chosen to be included in sequence A333242 and 0 indicates the prime numbers and composites not in A333242.

Crossrefs

Programs

  • Mathematica
    Select[Prime@ Range@ 75, EvenQ@ Length@ NestWhileList[ PrimePi, #, PrimeQ] &]
  • PARI
    A078442(n) = my(k=0); while(isprime(n), k++; n=primepi(n)); k;
    a(n) = A078442(n) % 2; \\ Michel Marcus, Jun 15 2024

Formula

a(n) = A078442(n) mod 2

A358179 Prime numbers with prime indices in A333244.

Original entry on oeis.org

31, 709, 1787, 8527, 19577, 27457, 42043, 52711, 72727, 96797, 112129, 137077, 167449, 173867, 239489, 250751, 285191, 352007, 401519, 443419, 464939, 490643, 527623, 683873, 718807, 755387, 839483, 864013, 985151, 1021271, 1080923, 1128889, 1159901, 1278779, 1323503, 1342907, 1656649, 1693031
Offset: 1

Views

Author

Michael P. May, Nov 11 2022

Keywords

Comments

This sequence can also be generated by the N-sieve.

Examples

			a(1) = prime(A333244(1)) = prime(11) = 31.
		

Crossrefs

Programs

  • Mathematica
    b[n_] := b[n] = If[PrimeQ[n], 1+b[PrimePi[n]], 0];
    a[n_] := a[n] = Module[{p}, p = If[n==1, 1, a[n-1]]; While[True, p = NextPrime[p]; If[#>3 && OddQ[#]&[b[p]], Break[]]]; p];
    Array[a, 50]
  • PARI
    b(n)={my(k=0); while(isprime(n), k++; n=primepi(n)); k};
    apply(x->prime(prime(prime(prime(x)))), select(n->b(n)%2, [1..500])) \\ Michel Marcus, Nov 18 2022

Formula

a(n) = prime(A333244(n)).
a(n) = A049090(A333242(n)).
a(n) = A038580(A262275(n)).
a(n) = A006450(A333243(n)).

A373497 If n is prime, a(n) = 1 if the number of steps in its prime index chain is odd, a(n) = -1 if the number of steps is even, and a(n) = 0 is n is composite or 1.

Original entry on oeis.org

0, 1, -1, 0, 1, 0, 1, 0, 0, 0, -1, 0, 1, 0, 0, 0, -1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, -1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, -1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, -1, 0, 0, 0, 0, 0, 1, 0
Offset: 1

Views

Author

Michael P. May, Jun 06 2024

Keywords

Crossrefs

Programs

  • PARI
    b(n)={my(k=0); while(isprime(n), k++; n=primepi(n)); k}; \\ A078442
    a(n) = if ((n==1) || !isprime(n), return(0)); if (b(n)%2, 1, -1); \\ Michel Marcus, Jun 11 2024

Formula

a(n) = 1 iff n is in A333242.
a(n) = -1 iff n is in A262275.
a(n) = 0 iff n is in A018252.

Extensions

More terms from Michel Marcus, Jun 11 2024
Previous Showing 31-35 of 35 results.