A049076 -id:A049076 - cp's OEIS frontend

A007821 Primes p such that pi(p) is not prime.

2, 7, 13, 19, 23, 29, 37, 43, 47, 53, 61, 71, 73, 79, 89, 97, 101, 103, 107, 113, 131, 137, 139, 149, 151, 163, 167, 173, 181, 193, 197, 199, 223, 227, 229, 233, 239, 251, 257, 263, 269, 271, 281, 293, 307, 311, 313, 317, 337, 347, 349, 359, 373

Offset: 1

Monte J. Zerger (mzerger(AT)cc4.adams.edu), Clark Kimberling

nonn

Primes prime(k) such that A049076(k)=1, sorted along increasing k. - R. J. Mathar, Jan 28 2014

The complement of A006450 (primes with prime index) within the primes A000040.

C. Kimberling, Fractal sequences and interspersions, Ars Combinatoria, vol. 45 p 157 1997.

R. Zumkeller, Table of n, a(n) for n = 1..1000
Lubomir Alexandrov, On the nonasymptotic prime number distribution, arXiv:math/9811096 [math.NT], 1998.
N. Fernandez, An order of primeness, F(p)
N. Fernandez, An order of primeness [cached copy, included with permission of the author]

Cf. A049076, A049078, A049079, A049080, A049081, A058322, A058324, A058325, A058326, A058327, A058328, A093046, A006450.

Let A = primes A000040, B = nonprimes A018252. The 2-level compounds are AA = A006450, AB = A007821, BA = A078782, BB = A102615. The 3-level compounds AAA, AAB, ..., BBB are A038580, A049078, A270792, A102617, A270794, A270795, A270796, A102616.

a007821 = a000040 . a018252
a007821_list = map a000040 a018252_list
-- Reinhard Zumkeller, Jan 12 2013

A007821 := proc(n) ithprime(A018252(n)) ; end proc: # R. J. Mathar, Jul 07 2012

Prime[ Select[ Range[75], !PrimeQ[ # ] &]] (* Robert G. Wilson v, Mar 15 2004 *)
With[{nn=100},Pick[Prime[Range[nn]],Table[If[PrimeQ[n],0,1],{n,nn}],1]] (* Harvey P. Dale, Aug 14 2020 *)

forprime(p=2, 1e3, if(!isprime(primepi(p)), print1(p, ", "))) \\ Felix Fröhlich, Aug 16 2014

from sympy import primepi
def A007821(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-(p:=primepi(x))+primepi(p)
    return bisection(f,n,n) # Chai Wah Wu, Oct 19 2024

A137588(a(n)) = n; a(n) = A000040(A018252(n)). - Reinhard Zumkeller, Jan 28 2008
A000040 = A007821 U A006450. - Juri-Stepan Gerasimov, Sep 24 2009
A175247 U { a(n); n > 1 } = A000040. { a(n) } = { 2 } U { primes (A000040) with composite index (A002808) }. - Jaroslav Krizek, Mar 13 2010
G.f. over nonprime powers: Sum_{k >= 1} prime(k)*x^k-prime(prime(k))*x^prime(k). - Benedict W. J. Irwin, Jun 11 2016

Edited by M. F. Hasler, Jul 31 2015

A038580 Primes with indices that are primes with prime indices.

Original entry on oeis.org

5, 11, 31, 59, 127, 179, 277, 331, 431, 599, 709, 919, 1063, 1153, 1297, 1523, 1787, 1847, 2221, 2381, 2477, 2749, 3001, 3259, 3637, 3943, 4091, 4273, 4397, 4549, 5381, 5623, 5869, 6113, 6661, 6823, 7193, 7607, 7841, 8221, 8527, 8719, 9319, 9461, 9739

Offset: 1

Brian Galebach

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from T. D. Noe)
Robert G. Batchko, A prime fractal and global quasi-self-similar structure in the distribution of prime-indexed primes, arXiv preprint arXiv:1405.2900 [math.GM], 2014.
Robert E. Dressler and S. Thomas Parker, Primes with a prime subscript, J. ACM, Volume 22 Issue 3, July 1975, 380-381.
Neil Fernandez, An order of primeness, F(p).
Neil Fernandez, An order of primeness. [cached copy, included with permission of the author]
Neil Fernandez, More terms of this and other sequences related to A049076.
Michael P. May, Properties of Higher-Order Prime Number Sequences, Missouri J. Math. Sci. (2020) Vol. 32, No. 2, 158-170; and arXiv version, arXiv:2108.04662 [math.NT], 2021.

Primes p for which A049076(p) > 3.
Cf. A049090, A049202, A049203, A057847, A057849, A057850, A057851, A058332, A086749, A093047.
Second differences give A245175.
Let A = primes A000040, B = nonprimes A018252. The 2-level compounds are AA = A006450, AB = A007821, BA = A078782, BB = A102615. The 3-level compounds AAA, AAB, ..., BBB are A038580, A049078, A270792, A102617, A270794, A270795, A270796, A102616.

[NthPrime(NthPrime(NthPrime(n))): n in [1..50]]; // Vincenzo Librandi, Jul 17 2016

a:= ithprime@@3;
seq(a(n), n=1..50);  # Alois P. Heinz, Jun 14 2015
# For Maple code for the prime/nonprime compound sequences (listed in cross-references) see A003622. - N. J. A. Sloane, Mar 30 2016

Table[ Prime[ Prime[ Prime[ n ] ] ], {n, 1, 60} ]
Nest[Prime, Range[45], 3] (* Robert G. Wilson v, Mar 15 2004 *)

a(n) = prime(prime(prime(n))) \\ Charles R Greathouse IV, Apr 28 2015

list(lim)=my(v=List(),q,r); forprime(p=2,lim, if(isprime(q++) && isprime(r++), listput(v,p))); Set(v) \\ Charles R Greathouse IV, Feb 14 2017

a(n) = prime(prime(prime(n))).

a(n) ~ n*log(n)^3. - Ilya Gutkovskiy, Jul 17 2016

A078442 a(p) = a(n) + 1 if p is the n-th prime, prime(n); a(n)=0 if n is not prime.

Original entry on oeis.org

0, 1, 2, 0, 3, 0, 1, 0, 0, 0, 4, 0, 1, 0, 0, 0, 2, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 5, 0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 3, 0, 1, 0, 0, 0, 0, 0, 2, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0

Offset: 1

Henry Bottomley, Dec 31 2002

nonn

Fernandez calls this the order of primeness of n.
a(A007097(n))=n, for any n >= 0. - Paul Tek, Nov 12 2013
When a nonoriented rooted tree is encoded as a Matula-Goebel number n, a(n) tells how many edges needs to be climbed up from the root of the tree until the first branching vertex (or the top of the tree, if n is one of the terms of A007097) is encountered. Please see illustrations at A061773. - Antti Karttunen, Jan 27 2014
Zero-based column index of n in the Kimberling-style dispersion table of the primes (see A114537). - Allan C. Wechsler, Jan 09 2024

			a(1) = 0 since 1 is not prime;
a(2) = a(prime(1)) = a(1) + 1 = 1 + 0 = 1;
a(3) = a(prime(2)) = a(2) + 1 = 1 + 1 = 2;
a(4) = 0 since 4 is not prime;
a(5) = a(prime(3)) = a(3) + 1 = 2 + 1 = 3;
a(6) = 0 since 6 is not prime;
a(7) = a(prime(4)) = a(4) + 1 = 0 + 1 = 1.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
N. Fernandez, An order of primeness, F(p)
N. Fernandez, An order of primeness [cached copy, included with permission of the author]
Index entries for sequences related to Matula-Goebel numbers

A left inverse of A007097.
One less than A049076.
a(A000040(n)) = A049076(n).
Cf. A373338 (mod 2), A018252 (positions of zeros).
Cf. A000720, A049084, A061773.
Cf. permutations A235489, A250247/A250248, A250249/A250250, A245821/A245822 that all preserve a(n).
Cf. also array A114537 (A138947) and permutations A135141/A227413, A246681.

a078442 n = fst $ until ((== 0) . snd)
                        (\(i, p) -> (i + 1, a049084 p)) (-2, a000040 n)
-- Reinhard Zumkeller, Jul 14 2013

A078442 := proc(n)
    if not isprime(n) then
        0 ;
    else
        1+procname(numtheory[pi](n)) ;
    end if;
end proc: # R. J. Mathar, Jul 07 2012

a[n_] := a[n] = If[!PrimeQ[n], 0, 1+a[PrimePi[n]]]; Array[a, 105] (* Jean-François Alcover, Jan 26 2018 *)

A078442(n)=for(i=0,n, isprime(n) || return(i); n=primepi(n)) \\ M. F. Hasler, Mar 09 2010

a(n) = A049076(n)-1.
a(n) = if A049084(n) = 0 then 0 else a(A049084(n)) + 1. - Reinhard Zumkeller, Jul 14 2013
For all n, a(n) = A007814(A135141(n)) and a(A227413(n)) = A007814(n). Also a(A235489(n)) = a(n). - Antti Karttunen, Jan 27 2014

A317713 Number of distinct terminal subtrees of the rooted tree with Matula-Goebel number n.

Original entry on oeis.org

1, 2, 3, 2, 4, 3, 3, 2, 3, 4, 5, 3, 4, 3, 4, 2, 4, 3, 3, 4, 4, 5, 4, 3, 4, 4, 3, 3, 5, 4, 6, 2, 5, 4, 5, 3, 4, 3, 4, 4, 5, 4, 4, 5, 4, 4, 5, 3, 3, 4, 5, 4, 3, 3, 5, 3, 4, 5, 5, 4, 4, 6, 4, 2, 5, 5, 4, 4, 4, 5, 5, 3, 5, 4, 4, 3, 6, 4, 6, 4, 3, 5, 5, 4, 6, 4, 5, 5, 4, 4, 5, 4, 6, 5, 5, 3, 5, 3, 5, 4, 5, 5, 4, 4, 5, 3, 4, 3

Offset: 1

Gus Wiseman, Aug 05 2018

nonn

			20 is the Matula-Goebel number of the tree (oo((o))), which has 4 distinct terminal subtrees: {(oo((o))), ((o)), (o), o}. So a(20) = 4.
See also illustrations in A061773.

Antti Karttunen, Table of n, a(n) for n = 1..16384
Index entries for sequences related to Matula-Göbel numbers

One more than A324923.

Cf. A000081, A007097, A049076, A061773, A061775, A076146, A109082, A109129, A206491, A303431, A316476.

ids[n_]:=Union@@FixedPointList[Union@@(Cases[If[#==1,{},FactorInteger[#]],{p_,_}:>PrimePi[p]]&/@#)&,{n}];
Table[Length[ids[n]],{n,100}]

A006530(n) = if(1==n, n, my(f=factor(n)); f[#f~, 1]);
A324923(n) = { my(lista = List([]), gpf, i); while(n > 1, gpf=A006530(n); i = primepi(gpf); n /= gpf; n *= i; listput(lista,i)); #Set(lista); }; \\ Antti Karttunen, Oct 23 2023
A317713(n) = (1+A324923(n)); \\ Antti Karttunen, Oct 23 2023

a(n) = 1+A324923(n). - Antti Karttunen, Oct 23 2023

Data section extended up to a(108) by Antti Karttunen, Oct 23 2023

A250249 Permutation of natural numbers: a(1) = 1, a(n) = A083221(a(A055396(n)), a(A246277(n))).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 27, 22, 23, 24, 25, 26, 21, 28, 29, 30, 31, 32, 39, 34, 35, 36, 37, 38, 63, 40, 41, 54, 43, 44, 33, 46, 47, 48, 49, 50, 75, 52, 53, 42, 65, 56, 99, 58, 59, 60, 61, 62, 57, 64, 95, 78, 67, 68, 111, 70, 71, 72, 103, 74, 51, 76, 77, 126, 79, 80, 45, 82

Offset: 1

Antti Karttunen, Nov 17 2014

nonn

This is a "doubly-recursed" version of A249817.
For primes p_n, a(p_n) = p_{a(n)}.
The first 7-cycle occurs at: (33 39 63 57 99 81 45), which is mirrored by the cycle (66 78 126 114 198 162 90) with terms double the size and also by the cycle (137 167 307 269 523 419 197), consisting of primes (p_33, p_39, p_63, ...).

			For n = 42 = 2*3*7, we see that it occurs as the 21st term on the top row of A246278 (A055396(42) = 1 and A246277(42) = 21), recursing on both yields a(1) = 1, a(21) = 27, thus we find A083221(1,27), the 27th term on A083221's topmost row (also A005843) which is 54, thus a(42) = 54.
Examples for cases where n is a prime:
a(3709) = a(p_518) = p_{a(518)} = A000040(1162) = 9397.
a(3719) = a(p_519) = p_{a(519)} = A000040(1839) = 15767.

Antti Karttunen, Table of n, a(n) for n = 1..3718
Index entries for sequences that are permutations of the natural numbers

Inverse: A250250.
Fixed points: A250251, their complement: A249729.
Cf. A000035, A000040, A055396, A078442, A083221, A049076.
Differs from A250250 for the first time at n=33, where a(33) = 39, while A250250(33) = 45.
Differs from the "vanilla version" A249817 for the first time at n=42, where a(42) = 54, while A249817(42) = 42.

a(1) = 1, a(n) = A083221(a(A055396(n)), a(A246277(n))).
Other identities. For all n >= 1:
a(2n) = 2*a(n), or equally, a(n) = a(2n)/2. [The even bisection halved gives the sequence back].
a(p_n) = p_{a(n)}, or equally, a(n) = A049084(a(A000040(n))). [Restriction to primes induces the same sequence].
A078442(a(n)) = A078442(n), A049076(a(n)) = A049076(n). [Preserves the "order of primeness of n"].
A000035(n) = A000035(a(n)). [Preserves the parity].

A245821 Permutation of natural numbers: a(n) = A091205(A245703(n)).

Original entry on oeis.org

1, 2, 3, 4, 5, 9, 7, 6, 8, 12, 11, 15, 23, 81, 18, 10, 17, 30, 13, 162, 27, 36, 19, 24, 16, 25, 38, 46, 37, 45, 31, 135, 14, 20, 50, 57, 47, 69, 21, 55, 83, 115, 419, 87, 60, 210, 61, 42, 54, 26, 90, 28, 29, 35, 32, 63, 171, 52, 59, 138, 113, 180, 111, 48, 88, 39, 41, 621, 72, 22, 953, 230, 103, 207, 126, 64, 33, 243

Offset: 1

Antti Karttunen, Aug 02 2014

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10001
Index entries for sequences that are permutations of the natural numbers

Inverse: A245822.
Other related permutations:  A091205, A245703, A245815.
Fixed points: A245823.
Cf. A000040, A049084, A078442, A049076, A018252, A062298.

allocatemem(234567890);
v014580 = vector(2^18);
v091226 = vector(2^22);
v091242 = vector(2^22);
isA014580(n)=polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from Charles R Greathouse IV
i=0; j=0; n=2; while((n < 2^22), if(isA014580(n), i++; v014580[i] = n; v091226[n] = v091226[n-1]+1, j++; v091242[j] = n; v091226[n] = v091226[n-1]); n++);
A014580(n) = v014580[n];
A091226(n) = v091226[n];
A091242(n) = v091242[n];
A091205(n) = if(n<=1, n, if(isA014580(n), prime(A091205(A091226(n))), {my(irfs, t); irfs=subst(lift(factor(Mod(1, 2)*Pol(binary(n)))), x, 2); irfs[,1]=apply(t->A091205(t), irfs[,1]); factorback(irfs)}));
A245703(n) = if(1==n, 1, if(isprime(n), A014580(A245703(primepi(n))), A091242(A245703(n-primepi(n)-1))));
A245821(n) = A091205(A245703(n));
for(n=1, 10001, write("b245821.txt", n, " ", A245821(n)));

(define (A245821 n) (A091205 (A245703 n)))

a(n) = A091205(A245703(n)).
Other identities. For all n >= 1, the following holds:
A078442(a(n)) = A078442(n), A049076(a(n)) = A049076(n). [Preserves "the order of primeness of n"].
a(p_n) = p_{a(n)} where p_n is the n-th prime, A000040(n).
a(n) = A049084(a(A000040(n))). [Thus the same permutation is induced also when it is restricted to primes].
A245815(n) = A062298(a(A018252(n))). [While restriction to nonprimes induces another permutation].

A250250 Permutation of natural numbers: a(1) = 1, a(n) = A246278(a(A055396(n)),a(A078898(n))).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 27, 22, 23, 24, 25, 26, 21, 28, 29, 30, 31, 32, 45, 34, 35, 36, 37, 38, 33, 40, 41, 54, 43, 44, 81, 46, 47, 48, 49, 50, 75, 52, 53, 42, 125, 56, 63, 58, 59, 60, 61, 62, 39, 64, 55, 90, 67, 68, 135, 70, 71, 72, 103, 74, 51, 76, 77, 66, 79, 80, 99, 82, 83

Offset: 1

Antti Karttunen, Nov 17 2014

nonn

This is a "doubly-recursed" version of A249818.

Antti Karttunen, Table of n, a(n) for n = 1..3432
Index entries for sequences that are permutations of the natural numbers

Inverse: A250249.
Fixed points: A250251, their complement: A249729.
Cf. A000035, A000040, A055396, A078442, A078898, A049076, A246278.
See also other (somewhat) similar permutations: A245821, A057505.
Differs from the "vanilla version" A249818 for the first time at n=42, where a(42) = 54, while A249818(42) = 42.

a(1) = 1, a(n) = A246278(a(A055396(n)), a(A078898(n))).
Other identities. For all n >= 1:
a(2n) = 2*a(n), or equally, a(n) = a(2n)/2. [The even bisection halved gives the sequence back].
a(p_n) = p_{a(n)}, or equally, a(n) = A049084(a(A000040(n))). [Restriction to primes induces the same sequence].
A078442(a(n)) = A078442(n), A049076(a(n)) = A049076(n). [Preserves the "order of primeness of n"].
A000035(n) = A000035(a(n)). [Preserves the parity].

A236542 Array T(n,k) read along descending antidiagonals: row n contains the primes with n steps in the prime index chain.

Original entry on oeis.org

2, 7, 3, 13, 17, 5, 19, 41, 59, 11, 23, 67, 179, 277, 31, 29, 83, 331, 1063, 1787, 127, 37, 109, 431, 2221, 8527, 15299, 709, 43, 157, 599, 3001, 19577, 87803, 167449, 5381, 47, 191, 919, 4397, 27457, 219613, 1128889, 2269733, 52711

Offset: 1

R. J. Mathar, Jan 28 2014

Row n contains the primes A000040(j) for which A049076(j) = n.

			The array starts:
    2,    7,   13,   19,   23,   29,   37,   43,   47,   53,...
    3,   17,   41,   67,   83,  109,  157,  191,  211,  241,...
    5,   59,  179,  331,  431,  599,  919, 1153, 1297, 1523,...
   11,  277, 1063, 2221, 3001, 4397, 7193, 9319,10631,12763,...
   31, 1787, 8527,19577,27457,42043,72727,96797,112129,137077,...

N. Fernandez, An order of primeness, F(p).

Cf. A007821 (row 1), A049078 (row 2), A049079 (row 3), A007097 (column 1), A058010 (diagonal), A057456 - A057457 (columns), A135044, A236536.

A236542 := proc(n,k)
    option remember ;
    if n = 1 then
        A007821(k) ;
    else
        ithprime(procname(n-1,k)) ;
    end if:
end proc:
for d from 2 to 10 do
    for k from d-1 to 1 by -1 do
            printf("%d,",A236542(d-k,k)) ;
    end do:
end do:

A007821 = Prime[Select[Range[15], !PrimeQ[#]&]];
T[n_, k_] := T[n, k] = If[n == 1, If[k <= Length[A007821], A007821[[k]], Print["A007821 must be extended"]; Abort[]], Prime[T[n-1, k]]];
Table[T[n-k+1, k], {n, 1, 9}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Apr 16 2020 *)

T(1,k) = A007821(k).

T(n,k) = prime( T(n-1,k) ), n>1 .

A245822 Permutation of natural numbers: a(n) = A245704(A091204(n)).

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 7, 9, 6, 16, 11, 10, 19, 33, 12, 25, 17, 15, 23, 34, 39, 70, 13, 24, 26, 50, 21, 52, 53, 18, 31, 55, 77, 93, 54, 22, 29, 27, 66, 105, 67, 48, 137, 156, 30, 28, 37, 64, 91, 35, 85, 58, 97, 49, 40, 98, 36, 135, 59, 45, 47, 261, 56, 76, 92, 122, 83, 374, 38, 102, 139, 69, 167, 130, 88, 203, 351, 212, 349, 235, 14

Offset: 1

Antti Karttunen, Aug 02 2014

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10001
Index entries for sequences that are permutations of the natural numbers

Inverse: A245821.
Other related permutations: A091204, A245704, A245816.
Fixed points: A245823.
Cf. A000040, A049084, A078442, A049076, A018252, A062298.

allocatemem(123456789);
default(primelimit, 2^22)
v014580 = vector(2^18); A014580(n) = v014580[n];
v091226 = vector(2^22); A091226(n) = v091226[n];
A002808(n)={ my(k=-1); while( -n + n += -k + k=primepi(n), ); n}; \\ This function from M. F. Hasler
isA014580(n)=polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from Charles R Greathouse IV
i=0; j=0; n=2; while((n < 2^22), if(isA014580(n), i++; v014580[i] = n; v091226[n] = v091226[n-1]+1, v091226[n] = v091226[n-1]); n++)
A091204(n) = if(n<=1, n, if(isprime(n), A014580(A091204(primepi(n))), {my(pfs, t, bits, i); pfs=factor(n); pfs[,1]=apply(t->Pol(binary(A091204(t))), pfs[,1]); sum(i=1, #bits=Vec(factorback(pfs))%2, bits[i]<<(#bits-i))}));
A091245(n) = ((n-A091226(n))-1);
A245704(n) = if(1==n, 1, if(isA014580(n), prime(A245704(A091226(n))), A002808(A245704(A091245(n)))));
A245822(n) = A245704(A091204(n));
for(n=1, 10001, write("b245822.txt", n, " ", A245822(n)));

(define (A245822 n) (A245704 (A091204 n)))

a(n) = A245704(A091204(n)).
Other identities. For all n >= 1, the following holds:
A078442(a(n)) = A078442(n), A049076(a(n)) = A049076(n). [Preserves "the order of primeness of n"].
a(p_n) = p_{a(n)} where p_n is the n-th prime, A000040(n).
a(n) = A049084(a(A000040(n))). [Thus the same permutation is induced also when it is restricted to primes].
A245816(n) = A062298(a(A018252(n))). [While restriction to nonprimes induces another permutation].

A050439 Fifth-order composites.

Original entry on oeis.org

39, 49, 55, 56, 60, 69, 74, 77, 78, 84, 93, 94, 95, 100, 105, 106, 110, 115, 119, 124, 125, 126, 130, 133, 140, 141, 145, 152, 155, 156, 159, 162, 164, 165, 170, 174, 180, 183, 184, 188, 189, 198, 201, 202, 203, 206, 207, 209, 212, 213, 218, 222, 225, 231

Offset: 1

Michael Lugo (mlugo(AT)thelabelguy.com), Dec 22 1999

			C(C(C(C(C(8))))) = C(C(C(C(15)))) = C(C(C(25))) = C(C(38)) = C(55) = 77. So 77 is in the sequence.

N. Fernandez, An order of primeness, F(p)
N. Fernandez, An order of primeness [cached copy, included with permission of the author]

Cf. A049076-A049081, A006450, A050435, A050436, A050438, A050440.

C := remove(isprime,[$4..1000]): seq(C[C[C[C[C[n]]]]],n=1..100);

Let C(n) be the n-th composite number, with C(1)=4. Then these are numbers C(C(C(C(C(n))))).

More terms from Asher Auel Dec 15 2000

cp's OEIS Frontend

A007821 Primes p such that pi(p) is not prime.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

Extensions

A038580 Primes with indices that are primes with prime indices.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Formula

A078442 a(p) = a(n) + 1 if p is the n-th prime, prime(n); a(n)=0 if n is not prime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A317713 Number of distinct terminal subtrees of the rooted tree with Matula-Goebel number n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A250249 Permutation of natural numbers: a(1) = 1, a(n) = A083221(a(A055396(n)), a(A246277(n))).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A245821 Permutation of natural numbers: a(n) = A091205(A245703(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Scheme

Formula

A250250 Permutation of natural numbers: a(1) = 1, a(n) = A246278(a(A055396(n)),a(A078898(n))).

Views

Author

Keywords

Comments