A245822 -id:A245822 - cp's OEIS frontend

A245823 Fixed points of A245821 and A245822.

1, 2, 3, 4, 5, 7, 11, 17, 24, 31, 59, 89, 127, 184, 277, 461, 669, 709, 1097, 1787, 1995, 3259, 4999, 5381, 8807, 15299, 17351, 30133, 48593, 52711, 60810, 91081, 167449, 192263

Offset: 1

Antti Karttunen, Aug 02 2014

nonn

If n is a member, then the n-th prime, p_n (= A000040(n)) is also a member. If p_n is a member, then its index n is also a member. Thus the sequence is completely determined by its nonprime terms: 1, 4, 24, 184, 669, 1995, 60810, ..., and is obtained as a union of prime recurrences that start with those values: A007097 U A057450 U ..., etc.

A007097 and A057450 are subsequences.

Cf. A000040, A245821, A245822.

default(primelimit, 2^30);for(n=1, 2^18, if(A245821(n) == n, print1(n, ", "))); \\ Other code as in A245821.

;; With Antti Karttunen's IntSeq-library.
(define A245823 (FIXED-POINTS 1 1 A245821))

A245816 Permutation of natural numbers induced when A245822 is restricted to nonprime numbers: a(n) = A062298(A245822(A018252(n))).

Original entry on oeis.org

1, 2, 4, 5, 3, 10, 6, 22, 7, 16, 9, 23, 27, 51, 15, 17, 35, 13, 37, 11, 39, 56, 69, 38, 14, 18, 48, 78, 33, 120, 20, 19, 46, 67, 24, 62, 42, 34, 28, 73, 25, 103, 31, 206, 40, 55, 68, 92, 300, 26, 76, 50, 99, 65, 157, 281, 165, 184, 8, 121, 134, 277, 423, 30, 47, 36, 223, 70, 514, 75, 101, 116, 236, 139, 74

Offset: 1

Antti Karttunen, Aug 02 2014

nonn

This permutation is induced when A245822 is restricted to nonprimes, A018252, the first column of A114537, but equally, when it is restricted to column 2 (A007821), column 3 (A049078), etc. of that square array, or alternatively, to the successive rows of A236542.

The sequence of fixed points f(n) begins as 1, 2, 15, 142, 548, 1694, 54681. A018252(f(n)) gives the nonprime terms of A245823.

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

Inverse: A245815.
Related permutations: A245814, A245820, A245822.
Cf. A007821, A018252, A049078, A062298, A114537, A236542, A245823.

allocatemem(123456789);
default(primelimit, 2^22)
A014580 = vector(2^18);
A091226 = vector(2^22);
A002808(n)={my(k); for(k=0, primepi(n), isprime(n++)&&k--); n}; \\ This function from M. F. Hasler, Oct 31 2008
A062298(n) = n-primepi(n);
A018252(n) = if(1==n,1,A002808(n-1));
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++; A014580[i] = n; A091226[n] = A091226[n-1]+1, A091226[n] = A091226[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));
A245816(n) = A062298(A245822(A018252(n)));
for(n=1, 10001, write("b245816.txt", n, " ", A245816(n)));

(define (A245816 n) (A062298 (A245822 (A018252 n))))

a(n) = A062298(A245822(A018252(n))).
As a composition of related permutations:
a(n) = A245820(A245814(n)).
Also following holds for all n >= 1:
a(n) = A062298(A049084(A245822(A007821(n)))).
a(n) = A062298(A049084(A049084(A245822(A049078(n))))).
etc.

A245818 Difference in size between rooted trees which are encoded with Matula-Goebel numbers A245822(n) and n: a(n) = A061775(A245822(n)) - A061775(n).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 2, -1, 2, 0, 0, 1, 0, 1, 2, -1, 0, -1, 2, -1, 1, 0, -1, 0, 2, 1, 2, 1, -1, 0, 1, 1, 2, 0, 0, 2, 2, -1, -1, -1, 0, 1, -1, 1, 0, 2, -1, -1, 1, 0, 3, 0, 0, 0, 3, -1, 0, 0, 0, 1, 3, -2, 0, 0, 0, 1, 2, -1, 2, 3, 0, 2, 2, -4, 1, -1, 1, 3, 4, -1, 0, 0, -2, 1, 0, 1, 1, 0, 0, -1, 1, 1, 0

Offset: 1

Antti Karttunen, Aug 16 2014

sign

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

Cf. A000040, A061775, A245817, A245822, A245823.

\\ Execute first the code given in A061775 and A245822.
A245818(n) = A061775(A245822(n)) - A061775(n);
for(n=1, 10001, write("b245818.txt", n, " ", A245818(n)));

(define (A245818 n) (- (A061775 (A245822 n)) (A061775 n)))

a(n) = A061775(A245822(n)) - A061775(n).
Other identities. For all n >= 1, the following holds:
a(A000040(n)) = a(n). [The result for the n-th prime is same as for n itself].
a(A245823(n)) = 0. [A245823 gives a (proper) subsequence of the positions of the zeros].

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

A245703 Permutation of natural numbers: a(1) = 1, a(p_n) = A014580(a(n)), a(c_n) = A091242(a(n)), where p_n = n-th prime, c_n = n-th composite number and A014580(n) and A091242(n) are binary codes for n-th irreducible and n-th reducible polynomials over GF(2), respectively.

Original entry on oeis.org

1, 2, 3, 4, 7, 5, 11, 6, 8, 12, 25, 9, 13, 17, 10, 14, 47, 18, 19, 34, 15, 20, 31, 24, 16, 21, 62, 26, 55, 27, 137, 45, 22, 28, 42, 33, 37, 23, 29, 79, 59, 35, 87, 71, 36, 166, 41, 58, 30, 38, 54, 44, 61, 49, 32, 39, 99, 76, 319, 46, 91, 108, 89, 48, 200, 53, 97, 75, 40, 50, 203, 70, 67, 57, 78, 64, 43, 51

Offset: 1

Antti Karttunen, Aug 02 2014

All the permutations A091202, A091204, A106442, A106444, A106446, A235041 share the same property that primes (A000040) are mapped bijectively to the binary representations of irreducible GF(2) polynomials (A014580) but while they determine the mapping of composites (A002808) to the corresponding binary codes of reducible polynomials (A091242) by a simple multiplicative rule, this permutation employs index-recursion also in that case.

Inverse: A245704.
Cf. A000040, A002808, A000720, A007097, A010051, A014580, A065855, A091225, A091230, A091242.
Similar or related permutations: A091202, A091204, A106442, A106444, A106446, A235041, A135141, A245701, A245702, A245821, A245822, A244987, A245450.

allocatemem(123456789);
a014580 = vector(2^18);
a091242 = 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++; a014580[i] = n, j++; a091242[j] = n); n++)
A245703(n) = if(1==n, 1, if(isprime(n), a014580[A245703(primepi(n))], a091242[A245703(n-primepi(n)-1)]));
for(n=1, 10001, write("b245703.txt", n, " ", A245703(n)));

;; With memoization-macro definec.
(definec (A245703 n) (cond ((= 1 n) n) ((= 1 (A010051 n)) (A014580 (A245703 (A000720 n)))) (else (A091242 (A245703 (A065855 n))))))

a(1) = 1, a(p_n) = A014580(a(n)) and a(c_n) = A091242(a(n)), where p_n is the n-th prime, A000040(n) and c_n is the n-th composite, A002808(n).
a(1) = 1, after which, if A010051(n) is 1 [i.e. n is prime], then a(n) = A014580(a(A000720(n))), otherwise a(n) = A091242(a(A065855(n))).
As a composition of related permutations:
a(n) = A245702(A135141(n)).
a(n) = A091204(A245821(n)).
Other identities. For all n >= 1, the following holds:
a(A007097(n)) = A091230(n). [Maps iterates of primes to the iterates of A014580. Permutation A091204 has the same property]
A091225(a(n)) = A010051(n). [Maps primes to binary representations of irreducible GF(2) polynomials, A014580, and nonprimes to union of {1} and the binary representations of corresponding reducible polynomials, A091242. The permutations A091202, A091204, A106442, A106444, A106446 and A235041 have the same property.]

A245704 Permutation of natural numbers: a(1) = 1, a(A014580(n)) = A000040(a(n)), a(A091242(n)) = A002808(a(n)), where A000040(n) = n-th prime, A002808(n) = n-th composite number, and A014580(n) and A091242(n) are binary codes for n-th irreducible and n-th reducible polynomial over GF(2), respectively.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 5, 9, 12, 15, 7, 10, 13, 16, 21, 25, 14, 18, 19, 22, 26, 33, 38, 24, 11, 28, 30, 34, 39, 49, 23, 55, 36, 20, 42, 45, 37, 50, 56, 69, 47, 35, 77, 52, 32, 60, 17, 64, 54, 70, 78, 94, 66, 51, 29, 105, 74, 48, 41, 84, 53, 27, 88, 76, 95, 106, 73, 125, 91, 72, 44, 140, 97, 100, 68, 58, 115, 75, 40

Offset: 1

Antti Karttunen, Aug 02 2014

nonn

All the permutations A091203, A091205, A106443, A106445, A106447, A235042 share the same property that the binary representations of irreducible GF(2) polynomials (A014580) are mapped bijectively to the primes (A000040) but while they determine the mapping of corresponding reducible polynomials (A091242) to the composite numbers (A002808) by a simple multiplicative rule, this permutation employs index-recursion also in that case.

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

Inverse: A245703.
Cf. A000040, A002808, A007097, A010051, A014580, A091225, A091226, A091230, A091242, A091245.
Similar or related permutations: A245701, A245822, A227413, A091203, A091205, A106443, A106445, A106447, A235042, A244987, A245450.

allocatemem(123456789);
default(primelimit, 2^22)
A091226 = vector(2^22);
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
j=0; n=2; while((n < 2^22), if(isA014580(n), A091226[n] = A091226[n-1]+1, A091226[n] = A091226[n-1]); n++)
A091245(n) = ((n-A091226[n])-1);
A245704(n) = if(1==n, 1, if(isA014580(n), prime(A245704(A091226[n])), A002808(A245704(A091245(n)))));
for(n=1, 10001, write("b245704.txt", n, " ", A245704(n)));

;; With memoization-macro definec.
(definec (A245704 n) (cond ((= 1 n) n) ((= 1 (A091225 n)) (A000040 (A245704 (A091226 n)))) (else (A002808 (A245704 (A091245 n))))))

a(1) = 1, after which, if A091225(n) is 1 [i.e. n is in A014580], then a(n) = A000040(a(A091226(n))), otherwise a(n) = A002808(a(A091245(n))).
As a composition of related permutations:
a(n) = A227413(A245701(n)).
a(n) = A245822(A091205(n)).
Other identities. For all n >= 1, the following holds:
a(A091230(n)) = A007097(n). [Maps iterates of A014580 to the iterates of primes. Permutation A091205 has the same property].
A010051(a(n)) = A091225(n). [After a(1)=1, maps binary representations of irreducible GF(2) polynomials (= A014580) to primes and the corresponding representations of reducible polynomials to composites].

A091204 Factorization and index-recursion preserving isomorphism from nonnegative integers to polynomials over GF(2).

Original entry on oeis.org

0, 1, 2, 3, 4, 7, 6, 11, 8, 5, 14, 25, 12, 19, 22, 9, 16, 47, 10, 31, 28, 29, 50, 13, 24, 21, 38, 15, 44, 61, 18, 137, 32, 43, 94, 49, 20, 55, 62, 53, 56, 97, 58, 115, 100, 27, 26, 37, 48, 69, 42, 113, 76, 73, 30, 79, 88, 33, 122, 319, 36, 41, 274, 39, 64, 121, 86, 185

Offset: 0

Antti Karttunen, Jan 03 2004. Name changed Aug 16 2014

nonn

This "deeply multiplicative" isomorphism is one of the deep variants of A091202 which satisfies most of the same identities as the latter, but it additionally preserves also the structures where we recurse on prime's index. E.g. we have: A091230(n) = a(A007097(n)) and A061775(n) = A091238(a(n)). This is because the permutation induces itself when it is restricted to the primes: a(n) = A091227(a(A000040(n))).

On the other hand, when this permutation is restricted to the nonprime numbers (A018252), permutation A245814 is induced.

Antti Karttunen, Table of n, a(n) for n = 0..8192
A. Karttunen, Scheme-program for computing this sequence.
Index entries for sequences operating on GF(2)[X]-polynomials
Index entries for sequences that are permutations of the natural numbers

Inverse: A091205.
Similar or related permutations: A091202, A106442, A106444, A106446, A235041, A245703, A245814, A245822.
Cf. A000040, A007097, A010051, A014580, A018252, A048720, A061775, A091225, A091230, A091238, A091242.

v014580 = vector(2^18); A014580(n) = v014580[n];
isA014580(n)=polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from Charles R Greathouse IV
i=0; n=2; while((n < 2^22), if(isA014580(n), i++; v014580[i] = n); 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))}));
for(n=0, 8192, write("b091204.txt", n, " ", A091204(n)));
\\ Antti Karttunen, Aug 16 2014

a(0)=0, a(1)=1, a(p_i) = A014580(a(i)) for primes with index i and for composites a(p_i * p_j * ...) = a(p_i) X a(p_j) X ..., where X stands for carryless multiplication of GF(2)[X] polynomials (A048720).
As a composition of related permutations:
a(n) = A245703(A245822(n)).
Other identities.
For all n >= 0, the following holds:
a(A007097(n)) = A091230(n). [Maps iterates of primes to the iterates of A014580. Permutation A245703 has the same property]
For all n >= 1, the following holds:
A091225(a(n)) = A010051(n). [Maps primes bijectively to binary representations of irreducible GF(2) polynomials, A014580, and nonprimes to union of {1} and the binary representations of corresponding reducible polynomials, A091242, in some order. The permutations A091202, A106442, A106444, A106446, A235041 and A245703 have the same property.]

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].

cp's OEIS Frontend

A245823 Fixed points of A245821 and A245822.

Views

Author

Keywords

Comments

Crossrefs

Programs

PARI

Scheme

A245816 Permutation of natural numbers induced when A245822 is restricted to nonprime numbers: a(n) = A062298(A245822(A018252(n))).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Scheme

Formula

A245818 Difference in size between rooted trees which are encoded with Matula-Goebel numbers A245822(n) and n: a(n) = A061775(A245822(n)) - A061775(n).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Scheme

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

A245703 Permutation of natural numbers: a(1) = 1, a(p_n) = A014580(a(n)), a(c_n) = A091242(a(n)), where p_n = n-th prime, c_n = n-th composite number and A014580(n) and A091242(n) are binary codes for n-th irreducible and n-th reducible polynomials over GF(2), respectively.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Scheme

Formula

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Scheme

Formula

A091204 Factorization and index-recursion preserving isomorphism from nonnegative integers to polynomials over GF(2).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

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

Views

Author

Keywords