A245814 -id:A245814 - cp's OEIS frontend

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

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

A245813 Permutation of natural numbers induced when A091205 is restricted to {1} and binary codes for polynomials reducible over GF(2): a(1) = 1, a(n) = A062298(A091205(A091242(n-1))).

Original entry on oeis.org

1, 2, 5, 3, 4, 9, 11, 7, 6, 18, 10, 59, 20, 25, 16, 8, 50, 15, 32, 31, 12, 13, 38, 21, 41, 125, 85, 43, 17, 45, 52, 35, 22, 19, 103, 105, 33, 24, 14, 190, 68, 27, 66, 28, 161, 29, 80, 26, 54, 46, 177, 84, 258, 34, 180, 64, 90, 70, 507, 37, 196, 96, 39, 110, 430, 92, 78, 75, 600, 48, 40, 82, 213, 218, 71, 23, 87, 72, 51, 132, 30

Offset: 1

Antti Karttunen, Aug 16 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: A245814.
Cf. A062298, A091242.
Related permutations: A091205, A245815, A245820.

allocatemem(234567890);
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++; v091226[n] = v091226[n-1]+1, j++; v091242[j] = n; v091226[n] = v091226[n-1]); n++);
A062298(n) = n-primepi(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)}));
A245813(n) = if(n<=1, n, A062298(A091205(A091242(n-1))));
for(n=1, 10001, write("b245813.txt", n, " ", A245813(n)));

(define (A245813 n) (if (<= n 1) n (A062298 (A091205 (A091242 (- n 1))))))

a(1) = 1, and for n > 1, a(n) = A062298(A091205(A091242(n-1))).
As a composition of related permutations:
a(n) = A245815(A245820(n)).

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.

A245819 Permutation of natural numbers induced when A245703 is restricted to nonprime numbers: a(n) = 1+A091245(A245703(A018252(n))).

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 6, 12, 7, 9, 13, 26, 10, 14, 18, 11, 15, 48, 19, 20, 35, 16, 21, 32, 25, 17, 22, 63, 27, 56, 28, 138, 46, 23, 29, 43, 34, 38, 24, 30, 80, 60, 36, 88, 72, 37, 167, 42, 59, 31, 39, 55, 45, 62, 50, 33, 40, 100, 77, 320, 47, 92, 109, 90, 49, 201, 54, 98, 76, 41, 51

Offset: 1

Antti Karttunen, Aug 16 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: A245820.
Related permutations: A245703, A245814, A245815.
Cf. A008578, A018252, A091226, A091245.

\\ The rest of code can be found in A245703.
A245819(n) = if(1==n, 1, 1 + A245703(n-1));

(define (A245819 n) (if (<= n 1) n (+ 1 (A245703 (- n 1)))))

(define (A245819 n) (+ 1 (A091245 (A245703 (A018252 n)))))

(define (A245819 n) (+ 1 (A091226 (A245703 (A008578 n)))))

a(1) = 1, and for n > 1, a(n) = 1 + A245703(n-1).
a(n) = 1+A091245(A245703(A018252(n))). [Induced when A245703 is restricted to nonprime numbers].
a(n) = 1+A091226(A245703(A008578(n))). [Induced also when A245703 is restricted to noncomposite numbers].
As a composition of related permutations:
a(n) = A245814(A245815(n)).

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

A245813 Permutation of natural numbers induced when A091205 is restricted to {1} and binary codes for polynomials reducible over GF(2): a(1) = 1, a(n) = A062298(A091205(A091242(n-1))).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Scheme

Formula

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

A245819 Permutation of natural numbers induced when A245703 is restricted to nonprime numbers: a(n) = 1+A091245(A245703(A018252(n))).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Scheme

Scheme

Scheme

Formula