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 21-30 of 35 results. Next

A250247 Permutation of natural numbers: a(1) = 1, a(n) = A083221(a(A055396(n)),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, 42, 43, 44, 33, 46, 47, 48, 49, 50, 75, 52, 53, 54, 65, 56, 99, 58, 59, 60, 61, 62, 57, 64, 95, 66, 67, 68, 111, 70, 71, 72, 103, 74, 51, 76, 77, 78, 79, 80, 45, 82
Offset: 1

Views

Author

Antti Karttunen, Nov 17 2014

Keywords

Comments

The first 7-cycle occurs at: (33 39 63 57 99 81 45), which is mirrored at the cycle (137 167 307 269 523 419 197), consisting of primes (p_33, p_39, p_63, ...).

Examples

			As a(21) = 27, and A000040(21) = 73 and A000040(27) = 103, a(73) = 103.
		

Crossrefs

Inverse: A250248.
Differs from its inverse A250248 for the first time at n = 33, where a(33) = 39, while A250248(33) = 45.
Differs from the "vanilla version" A249817 for the first time at n=73, where a(73) = 103, while A249817(73) = 73.
Differs from "doubly recursed" version A250249 for the first time at n=42, where a(42) = 42, while A250249(42) = 54, thus the first prime where they get different values is p_42 = 181, where a(181) = 181, while A250249(181) = 251 = p_54.

Formula

a(1) = 1, a(n) = A083221(a(A055396(n)),A246277(n)).
Other identities. For all n >= 1:
a(A005843(n)) = A005843(n). [Fixes even numbers].
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"].

A246682 Permutation of natural numbers: a(1) = 0, a(2) = 1, and for n > 1, a(2n) = nthcomposite(a(n)), a(2n-1) = nthprime(a(A064989(2n-1))), where nthprime = A000040, nthcomposite = A002808, and A064989(n) shifts the prime factorization of n one step towards smaller primes.

Original entry on oeis.org

0, 1, 2, 4, 3, 6, 5, 9, 7, 8, 11, 12, 31, 10, 13, 16, 127, 14, 709, 15, 19, 20, 5381, 21, 17, 46, 23, 18, 52711, 22, 648391, 26, 29, 166, 41, 24, 9737333, 858, 71, 25, 174440041, 30, 3657500101, 32, 37, 6186
Offset: 1

Views

Author

Antti Karttunen, Sep 01 2014

Keywords

Comments

Note the indexing: the domain starts from 1, while the range includes also zero.
Has an infinite number of infinite cycles. See comments at A246681.

Crossrefs

Inverse: A246681.
Similar or related permutations: A246376, A246378, A243071, A246368, A064216, A246380.

Programs

Formula

a(1) = 0, a(2) = 1, and for n > 1, a(2n) = nthcomposite(a(n)), a(2n-1) = nthprime(a(A064989(2n-1))), where nthprime = A000040, nthcomposite = A002808, and A064989(n) shifts the prime factorization of n one step towards smaller primes.
As a composition of related permutations:
a(n) = A246378(A243071(n)).
Other identities.
For all n >= 1 the following holds:
a(A000040(n)) = A007097(n-1). [Maps primes to the iterates of primes].
A049076(a(A000040(n))) = n. [Follows from above].
For all n > 1 the following holds:
A010051(a(n)) = A000035(n). [Maps odd numbers larger than one to primes, and even numbers to composites, in some order. Permutations A246378 & A246380 have the same property].

A135044 a(1)=1, then a(c) = p and a(p) = c, where c = T_c(r,k) and p = T_p(r,k), and where T_p contains the primes arranged in rows by the prime index chain and T_c contains the composites arranged in rows by the order of compositeness. See Formula.

Original entry on oeis.org

1, 4, 9, 2, 16, 7, 6, 13, 3, 19, 26, 17, 8, 23, 41, 5, 12, 67, 10, 29, 59, 37, 14, 83, 179, 11, 43, 331, 20, 47, 39, 109, 277, 157, 53, 431, 22, 1063, 31, 191, 15, 2221, 27, 61, 211, 71, 30, 599, 1787, 919, 241, 3001, 35, 73, 8527, 127, 1153, 79, 21, 19577, 44, 89, 283
Offset: 1

Views

Author

Katarzyna Matylla, Feb 11 2008

Keywords

Comments

Exchanges primes with composites, primeth primes with composith composites, etc.
Exchange the k-th prime of order j with the k-th composite of order j and vice versa.
Self-inverse permutation of positive integers.
If n is the composite number A236536(r,k), then a(n) is the corresponding prime A236542(r,k) at the same position (r,k). Vice versa, if n is the prime A236542(r,k), then a(n) is the corresponding composite A236536(r,k) at the same position. - Andrew Weimholt, Jan 28 2014
The original name for this entry did not produce this sequence, but instead A236854, which differs from this permutation for the first time at n=8, where A236854(8)=23, while here a(8)=13. - Antti Karttunen, Feb 01 2014

Examples

			From _Andrew Weimholt_, Jan 29 2014: (Start)
More generally, takes the primes organized in an array according to the sieving process described in the Fernandez paper:
        Row[1](n) = 2, 7, 13, 19, 23, ...
        Row[2](n) = 3, 17, 41, 67, 83, ...
        Row[3](n) = 5, 59, 179, ...
        Row[4](n) = 11, 277, ...
        Lets call this  T_p (n, k)
Also take the composites organized in a similar manner, except we use "composite" numbered positions in our sieve:
        Row[1](n) = 4, 6, 8, 10, 14, 20, 22, ...
        Row[2](n) = 9, 12, 15, 18, 24, ...
        Row[3](n) = 16, 21, 25, ...
        Lets call this T_c (n, k)
If we now take the natural numbers and swap each number (except for 1) with the number which holds the same spot in the other array, then we get the sequence: 1, 4, 9, 2, 16, 7, 6, 13, with for example a(8) = 13 (13 holds the same position in the 'prime' table as 8 does in the 'composite' table). (End)
		

Crossrefs

Programs

  • Maple
    A135044 := proc(n)
        if n = 1 then
            1;
        elif isprime(n) then
            idx := -1 ;
            for r from 1 do
                for c from 1 do
                    if A236542(r,c) = n then
                        idx := [r,c] ;
                    end if;
                    if A236542(r,c) >= n then
                        break;
                    end if;
                end do:
                if type(idx,list)  then
                    break;
                end if;
            end do:
            A236536(r,c) ;
        else
            idx := -1 ;
            for r from 1 do
                for c from 1 do
                    if A236536(r,c) = n then
                        idx := [r,c] ;
                    end if;
                    if A236536(r,c) >= n then
                        break;
                    end if;
                end do:
                if type(idx,list)  then
                    break;
                end if;
            end do:
            A236542(r,c) ;
        end if;
    end proc: # R. J. Mathar, Jan 28 2014
  • Mathematica
    Composite[n_Integer] := Block[{k = n + PrimePi@n + 1}, While[k != n + PrimePi@k + 1, k++ ]; k]; Compositeness[n_] := Block[{c = 1, k = n}, While[ !(PrimeQ@k || k == 1), k = k - 1 - PrimePi@k; c++ ]; c]; Primeness[n_] := Block[{c = 1, k = n}, While[ PrimeQ@k, k = PrimePi@k; c++ ]; c];
    ckj[k_, j_] := Select[ Table[Composite@n, {n, 10000}], Compositeness@# == j &][[k]]; pkj[k_, j_] := Select[ Table[Prime@n, {n, 3000}], Primeness@# == j &][[k]]; f[0]=0; f[1] = 1;
    f[n_] := If[ PrimeQ@ n, pn = Primeness@n; ckj[ Position[ Select[ Table[ Prime@ i, {i, 150}], Primeness@ # == pn &], n][[1, 1]], pn], cn = Compositeness@n; pkj[ Position[ Select[ Table[ Composite@ i, {i, 500}], Compositeness@ # == cn &], n][[1, 1]], cn]]; Array[f, 64] (* Robert G. Wilson v *)

Formula

a(1)=1, a(A236536(r,k))=A236542(r,k), a(A236542(r,k))=A236536(r,k)

Extensions

Edited, corrected and extended by Robert G. Wilson v, Feb 18 2008
Name corrected by Andrew Weimholt, Jan 29 2014

A250248 Permutation of natural numbers: a(1) = 1, a(n) = A246278(a(A055396(n)),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, 42, 43, 44, 81, 46, 47, 48, 49, 50, 75, 52, 53, 54, 125, 56, 63, 58, 59, 60, 61, 62, 39, 64, 55, 66, 67, 68, 135, 70, 71, 72, 103, 74, 51, 76, 77, 78, 79, 80, 99, 82, 83
Offset: 1

Views

Author

Antti Karttunen, Nov 17 2014

Keywords

Crossrefs

Inverse: A250247.
Similar permutations: A250250 for even more recursed variant of A249818.
Differs from the "vanilla version" A249818 for the first time at n=73, where a(73) = 108, while A249818(73) = 73.

Formula

a(1) = 1, a(n) = A246278(a(A055396(n)), A078898(n)).
Other identities. For all n >= 1:
a(A005843(n)) = A005843(n). [Fixes even numbers].
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"].

A333243 Prime numbers with prime indices in A262275.

Original entry on oeis.org

5, 31, 59, 179, 331, 431, 599, 709, 919, 1153, 1297, 1523, 1787, 1847, 2381, 2477, 2749, 3259, 3637, 3943, 4091, 4273, 4549, 5623, 5869, 6113, 6661, 6823, 7607, 7841, 8221, 8527, 8719, 9461, 9739, 9859, 11743, 11953, 12097, 12301, 12547, 13469, 13709, 14177
Offset: 1

Views

Author

Michael P. May, Mar 12 2020

Keywords

Comments

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

Examples

			a(1) = prime(A262275(1)) = prime(3) = 5.
		

Crossrefs

Programs

  • Maple
    b:= proc(n) option remember;
          `if`(isprime(n), 1+b(numtheory[pi](n)), 0)
        end:
    a:= proc(n) option remember; local p;
          p:= `if`(n=1, 1, a(n-1));
          do p:= nextprime(p);
            if (h-> h>1 and h::odd)(b(p)) then break fi
          od; p
        end:
    seq(a(n), n=1..50);  # Alois P. Heinz, Mar 15 2020
  • 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[#>1 && OddQ[#]&[b[p]], Break[]]]; p];
    Array[a, 50] (* Jean-François Alcover, Nov 16 2020, after Alois P. Heinz *)
  • PARI
    b(n)={my(k=0); while(isprime(n), k++; n=primepi(n)); k};
    apply(x->prime(prime(x)), select(n->b(n)%2, [1..500])) \\ Michel Marcus, Nov 18 2022

Formula

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

A333244 Prime numbers with prime indices in A333243.

Original entry on oeis.org

11, 127, 277, 1063, 2221, 3001, 4397, 5381, 7193, 9319, 10631, 12763, 15299, 15823, 21179, 22093, 24859, 30133, 33967, 37217, 38833, 40819, 43651, 55351, 57943, 60647, 66851, 68639, 77431, 80071, 84347, 87803, 90023, 98519, 101701, 103069, 125113, 127643
Offset: 1

Views

Author

Michael P. May, Mar 12 2020

Keywords

Comments

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

Examples

			a(1) = prime(A333243(1)) = prime(5) = 11.
		

Crossrefs

Programs

  • Maple
    b:= proc(n) option remember;
          `if`(isprime(n), 1+b(numtheory[pi](n)), 0)
        end:
    a:= proc(n) option remember; local p;
          p:= `if`(n=1, 1, a(n-1));
          do p:= nextprime(p);
            if (h-> h>2 and h::even)(b(p)) then break fi
          od; p
        end:
    seq(a(n), n=1..42);  # Alois P. Heinz, Mar 15 2020
  • 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[#>2 && EvenQ[#]&[b[p]], Break[]]]; p];
    Array[a, 42] (* Jean-François Alcover, Nov 16 2020, after Alois P. Heinz *)
  • PARI
    b(n)={my(k=0); while(isprime(n), k++; n=primepi(n)); k};
    apply(x->prime(prime(prime(x))), select(n->b(n)%2, [1..500])) \\ Michel Marcus, Nov 18 2022

Formula

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

A322027 Maximum order of primeness among the prime factors of n; a(1) = 0.

Original entry on oeis.org

0, 1, 2, 1, 3, 2, 1, 1, 2, 3, 4, 2, 1, 1, 3, 1, 2, 2, 1, 3, 2, 4, 1, 2, 3, 1, 2, 1, 1, 3, 5, 1, 4, 2, 3, 2, 1, 1, 2, 3, 2, 2, 1, 4, 3, 1, 1, 2, 1, 3, 2, 1, 1, 2, 4, 1, 2, 1, 3, 3, 1, 5, 2, 1, 3, 4, 2, 2, 2, 3, 1, 2, 1, 1, 3, 1, 4, 2, 1, 3, 2, 2, 2, 2, 3, 1, 2
Offset: 1

Views

Author

Gus Wiseman, Nov 24 2018

Keywords

Comments

The order of primeness (A078442) of a prime number p is the number of times one must apply A000720 to obtain a nonprime number.

Examples

			a(105) = 3 because the prime factor of 105 = 3*5*7 with maximum order of primeness is 5, with order 3.
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    p:= proc(n) option remember;
          `if`(isprime(n), 1+p(pi(n)), 0)
        end:
    a:= n-> max(0, map(p, factorset(n))):
    seq(a(n), n=1..120);  # Alois P. Heinz, Nov 24 2018
  • Mathematica
    Table[If[n==1,0,Max@@(Length[NestWhileList[PrimePi,PrimePi[#],PrimeQ]]&/@FactorInteger[n][[All,1]])],{n,100}]

A138947 Square array T[i+1,j] = prime(T[i,j]), T[1,j] = j-th nonprime = A018252(j); read by upward antidiagonals.

Original entry on oeis.org

1, 4, 2, 6, 7, 3, 8, 13, 17, 5, 9, 19, 41, 59, 11, 10, 23, 67, 179, 277, 31, 12, 29, 83, 331, 1063, 1787, 127, 14, 37, 109, 431, 2221, 8527, 15299, 709, 15, 43, 157, 599, 3001, 19577, 87803
Offset: 1

Views

Author

M. F. Hasler, Apr 28 2008

Keywords

Comments

For i>1, T[i,j] = A018252(j)-th number among those not occurring in rows < i.
A permutation of the integers > 0.
Transpose of A114537. See that sequence and the link for more information and references.

Examples

			The first row (1,4,6,8,9,10...) of the array gives the nonprime numbers A018252.
The 2nd row (2,7,13,19,23,29,37,...) of the array gives the primes with nonprime index, A000040(A018252(j)) = A007821(j).
The i-th row is { A000040(k) | A049076(k)=i-1 } = A078442^{-1}(i-1).
Column j is the sequence b(n+1)=prime(b(n)) starting with b(j)=A018252(j): A007097, A057450, A057451, A057452, A057453, A057456, A057457, ...
		

References

  • Alexandrov, Lubomir. "On the nonasymptotic prime number distribution." arXiv preprint math/9811096 (1998). (See Appendix.)

Crossrefs

If the antidiagonals are read in the opposite direction we get A114537.

Programs

  • Mathematica
    t[1, 1] = 1; t[1, 2] = 4; t[1, k_] := (p = t[1, k-1]; If[ PrimeQ[p+1], p+2, p+1]); t[n_ /; n > 1, k_] := Prime[t[n-1, k]]; Flatten[ Table[ t[n, k-n+1], {k, 1, 9}, {n, 1, k}]] (* Jean-François Alcover, Oct 03 2011 *)
  • PARI
    p=c=0; T=matrix( 10,10, i,j, if( i==1, while( isprime(c++),); p=c, p=prime(p))); A138947=concat( vector( vecmin( matsize( T )),i, vector( i,j, T[ j,i+1-j ])))

Formula

T[i,j] = j-th number for which A078442 equals i-1.

A322028 Number of distinct orders of primeness among the prime factors of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 2, 2, 2, 2, 1, 1, 2, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 1, 1, 3, 1, 2, 2, 1, 2, 3, 1, 2, 2, 2, 1, 2, 1, 1, 2, 1, 2, 2, 1, 2, 1, 2, 1, 2, 2, 1, 2
Offset: 1

Views

Author

Gus Wiseman, Nov 24 2018

Keywords

Comments

The order of primeness (A078442) of a prime number p is the number of times one must apply A000720 to obtain a nonprime number.

Examples

			a(105) = 3 because the prime factors of 105 = 3*5*7 have 3 different orders of primeness, namely 2, 3, and 1 respectively.
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    p:= proc(n) option remember;
          `if`(isprime(n), 1+p(pi(n)), 0)
        end:
    a:= n-> nops(map(p, factorset(n))):
    seq(a(n), n=1..120);  # Alois P. Heinz, Nov 24 2018
  • Mathematica
    Table[If[n==1,0,Length[Union[Length[NestWhileList[PrimePi,PrimePi[#],PrimeQ]]&/@FactorInteger[n][[All,1]]]]],{n,100}]

A322030 Numbers whose prime factors all have the same order of primeness.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 14, 16, 17, 19, 23, 25, 26, 27, 28, 29, 31, 32, 37, 38, 41, 43, 46, 47, 49, 51, 52, 53, 56, 58, 59, 61, 64, 67, 71, 73, 74, 76, 79, 81, 83, 86, 89, 91, 92, 94, 97, 98, 101, 103, 104, 106, 107, 109, 112, 113, 116, 121, 122, 123
Offset: 1

Views

Author

Gus Wiseman, Nov 24 2018

Keywords

Comments

The order of primeness (A078442) of a prime number p is the number of times one must apply A000720 to obtain a nonprime number.

Examples

			182 is in the sequence because its prime factors 2, 7, 13 all have order of primeness 1.
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    p:= proc(n) option remember;
          `if`(isprime(n), 1+p(pi(n)), 0)
        end:
    a:= proc(n) option remember; local k; for k from 1+`if`(n=1,
          0, a(n-1)) while nops(map(p, factorset(k)))>1 do od; k
        end:
    seq(a(n), n=1..100);  # Alois P. Heinz, Nov 24 2018
  • Mathematica
    ordpri[n_]:=If[!PrimeQ[n],0,Length[NestWhileList[PrimePi,PrimePi[n],PrimeQ]]];
    Select[Range[200],SameQ@@ordpri/@FactorInteger[#][[All,1]]&]
Previous Showing 21-30 of 35 results. Next