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 11-20 of 77 results. Next

A091203 Factorization-preserving isomorphism from binary codes of GF(2) polynomials to integers.

Original entry on oeis.org

0, 1, 2, 3, 4, 9, 6, 5, 8, 15, 18, 7, 12, 11, 10, 27, 16, 81, 30, 13, 36, 25, 14, 33, 24, 17, 22, 45, 20, 21, 54, 19, 32, 57, 162, 55, 60, 23, 26, 63, 72, 29, 50, 51, 28, 135, 66, 31, 48, 35, 34, 243, 44, 39, 90, 37, 40, 99, 42, 41, 108, 43, 38, 75, 64, 225, 114, 47, 324
Offset: 0

Views

Author

Antti Karttunen, Jan 03 2004

Keywords

Comments

E.g. we have the following identities: A000040(n) = a(A014580(n)), A091219(n) = A008683(a(n)), A091220(n) = A000005(a(n)), A091221(n) = A001221(a(n)), A091222(n) = A001222(a(n)), A091225(n) = A010051(a(n)), A091227(n) = A049084(a(n)), A091247(n) = A066247(a(n)).

Links

Crossrefs

Inverse: A091202.
Several variants exist: A235042, A091205, A106443, A106445, A106447.
Cf. also A000005, A000040, A001221, A001222, A004247, A008683, A010051, A014580, A048720, A049084, A066247, A091219, A091220, A091221, A091222, A091225, A091227, A091247, A234741, A234742, A245703, A245704, A278233.
Cf. also A302024, A302026, A305418, A305428 for other similar permutations.

Programs

  • PARI
    A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
A091225(n) = polisirreducible(Pol(binary(n))*Mod(1, 2));
A305419(n) = if(n<3,1, my(k=n-1); while(k>1 && !A091225(k),k--); (k));
A305422(n) = { my(f = subst(lift(factor(Pol(binary(n))*Mod(1, 2))),x,2)); for(i=1,#f~,f[i,1] = Pol(binary(A305419(f[i,1])))); fromdigits(Vec(factorback(f))%2,2); };
A091203(n) = if(n<=1,n,if(!(n%2),2*A091203(n/2),A003961(A091203(A305422(n))))); \\ Antti Karttunen, Jun 10 2018

Formula

a(0)=0, a(1)=1. For n's coding an irreducible polynomial ir_i, that is if n=A014580(i), we have a(n) = A000040(i) and for composite polynomials a(ir_i X ir_j X ...) = p_i * p_j * ..., where p_i = A000040(i) and X stands for carryless multiplication of GF(2)[X] polynomials (A048720) and * for the ordinary multiplication of integers (A004247).
Other identities. For all n >= 1, the following holds:
A010051(a(n)) = A091225(n). [After a(1)=1, maps binary representations of irreducible GF(2) polynomials, A014580, to primes and the binary representations of corresponding reducible polynomials, A091242, to composite numbers. The permutations A091205, A106443, A106445, A106447, A235042 and A245704 have the same property.]
From Antti Karttunen, Jun 10 2018: (Start)
For n <= 1, a(n) = n, for n > 1, a(n) = 2*a(n/2) if n is even, and if n is odd, then a(n) = A003961(a(A305422(n))).
a(n) = A005940(1+A305418(n)) = A163511(A305428(n)).
A046523(a(n)) = A278233(n).
(End)

A091214 Composite numbers whose binary representation encodes a polynomial irreducible over GF(2).

Original entry on oeis.org

25, 55, 87, 91, 115, 117, 143, 145, 171, 185, 203, 213, 247, 253, 285, 299, 301, 319, 333, 351, 355, 357, 361, 369, 375, 391, 395, 415, 425, 445, 451, 471, 477, 501, 505, 515, 529, 535, 539, 545, 623, 637, 665, 675, 687, 695, 721, 731, 789, 799, 803, 817
Offset: 1

Views

Author

Antti Karttunen, Jan 03 2004

Keywords

Comments

"Encoded in binary representation" means that a polynomial a(n)*X^n+...+a(0)*X^0 over GF(2) is represented by the binary number a(n)*2^n+...+a(0)*2^0 in Z (where each coefficient a(k) = 0 or 1).

Links

Crossrefs

Intersection of A002808 and A014580.
Subsequence of A235033, A236834 and A236838.
Left inverse: A235044.
Cf. A091206 (Primes whose binary expansion encodes a polynomial irreducible over GF(2)), A091209 (Primes that encode a polynomial reducible over GF(2)), A091212 (Composite, and reducible over GF(2)).
Cf. A091215, A235027, A235046, A236841, A236845, A236850, A236851, A236861.
Cf. also A235041-A235042.

Programs

  • Mathematica
    fQ[n_] := Block[{ply = Plus @@ (Reverse@ IntegerDigits[n, 2] x^Range[0, Floor@ Log2@ n])}, ply == Factor[ply, Modulus -> 2] && n != 2^Floor@ Log2@ n && ! PrimeQ@ n]; Select[ Range@ 840, fQ] (* Robert G. Wilson v, Aug 12 2011 *)
  • PARI
    isA014580(n)=polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from Charles R Greathouse IV
isA091214(n) = (!isprime(n) && isA014580(n));
n = 0; i = 0; while(n < 2^20, n++; if(isA091214(n), i++; write("b091214.txt", i, " ", n)));
\\ The b-file was computed with this program. Antti Karttunen, May 17 2015

Formula

Other identities. For all n >= 1:
A235044(a(n)) = n. [A235044 works as a left inverse of this sequence.]
a(n) = A014580(A091215(n)). - Antti Karttunen, May 17 2015

Extensions

Entry revised and name corrected by Antti Karttunen, May 17 2015

A235042 Factorization-preserving bijection from GF(2)[X]-polynomials to nonnegative integers, version which fixes the elements that are irreducible in both semirings.

Original entry on oeis.org

0, 1, 2, 3, 4, 9, 6, 7, 8, 21, 18, 11, 12, 13, 14, 27, 16, 81, 42, 19, 36, 49, 22, 39, 24, 5, 26, 63, 28, 33, 54, 31, 32, 93, 162, 91, 84, 37, 38, 99, 72, 41, 98, 15, 44, 189, 78, 47, 48, 77, 10, 243, 52, 57, 126, 17, 56, 117, 66, 59, 108, 61, 62, 147, 64, 441
Offset: 0

Views

Author

Antti Karttunen, Jan 02 2014

Keywords

Comments

Like A091203 this is a factorization-preserving isomorphism from GF(2)[X]-polynomials to integers. The former are encoded in the binary representation of n like this: n=11, '1011' in binary, stands for polynomial x^3+x+1, n=25, '11001' in binary, stands for polynomial x^4+x^3+1. However, this version does not map the irreducible GF(2)[X] polynomials (A014580) straight to the primes (A000040), but instead fixes the intersection of those two sets (A091206), and maps the elements in their set-wise difference A014580 \ A000040 (= A091214) in numerical order to the set-wise difference A000040 \ A014580 (= A091209).
The composite values are defined by the multiplicativity. E.g., we have a(A048724(n)) = 3*a(n) and a(A001317(n)) = A000244(n) = 3^n for all n.
This map satisfies many of the same identities as A091203, e.g., we have A091219(n) = A008683(a(n)), A091220(n) = A000005(a(n)), A091221(n) = A001221(a(n)), A091222(n) = A001222(a(n)), A091225(n) = A010051(a(n)) and A091247(n) = A066247(a(n)) for all n >= 1.

Examples

			Here (t X u) = A048720(t,u):
a(2)=2, a(3)=3 and a(7)=7, as 2, 3 and 7 are all in A091206.
a(4) = a(2 X 2) = a(2)*a(2) = 2*2 = 4.
a(5) = a(3 X 3) = a(3)*a(3) = 3*3 = 9.
a(9) = a(3 X 7) = a(3)*a(7) = 3*7 = 21.
a(10) = a(2 X 3 X 3) = a(2)*a(3)*a(3) = 2*3*3 = 18.
a(15) = a(3 X 3 X 3) = a(3)^3 = 3^3 = 27.
a(17) = a(3 X 3 X 3 X 3) = a(3)^4 = 3^4 = 81.
a(21) = a(7 X 7) = a(7)*a(7) = 7*7 = 49.
a(25) = 5, as 25 is the first term of A091214 and 5 is the first term of A091209.
a(50) = a(2 X 25) = a(2)*a(25) = 2*5 = 10.

Links

Crossrefs

Inverse: A235041. Fixed points: A235045.
Cf. A010051, A004247, A091225, A091209, A235044.
Similar cross-multiplicative permutations: A091203, A091205, A106443, A106445, A106447.

Formula

a(0)=0, a(1)=1, a(p) = p for those irreducible GF(2)[X]-polynomials whose binary encoding is a prime (i.e., p is in A091206), and for the rest of irreducible GF(2)[X]-polynomials (those which are encoded by a composite natural number, i.e., q is in A091214), a(q) = A091209(A235044(q)), and for reducible polynomials, a(i X j X k X ...) = a(i) * a(j) * a(k) * ..., where each i, j, k, ... is in A014580, X stands for carryless multiplication of GF(2)[X] polynomials (A048720) and * for the ordinary multiplication of integers (A004247).

A091245 Number of reducible GF(2)[X]-polynomials in range [0,n].

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Jan 03 2004

Keywords

Comments

Analogous to A065855.

Examples

			In range [0,8] there are the following four reducible polynomials: 4,5,6,8 thus a(8) = 4.

Links

Crossrefs

Partial sums of A091247. Cf. A091242.

Programs

  • PARI
    first(n)=my(s); concat([0,0], vector(n-1,k, s += !polisirreducible(Pol(binary(k+1))*Mod(1,2)))) \\ Charles R Greathouse IV, Sep 02 2015

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

Views

Author

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

Keywords

Comments

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.

Links

Crossrefs

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.

Programs

  • PARI
    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

Formula

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

A106490 Total number of bases and exponents in Quetian Superfactorization of n, excluding the unity-exponents at the tips of branches.

Original entry on oeis.org

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

Views

Author

Antti Karttunen, May 09 2005 based on Leroy Quet's message ('Super-Factoring' An Integer) posted to SeqFan-mailing list on Dec 06 2003

Keywords

Comments

Quetian Superfactorization proceeds by factoring a natural number to its unique prime-exponent factorization (p1^e1 * p2^e2 * ... pj^ej) and then factoring recursively each of the (nonzero) exponents in similar manner, until unity-exponents are finally encountered.

Examples

			a(64) = 3, as 64 = 2^6 = 2^(2^1*3^1) and there are three non-1 nodes in that superfactorization. Similarly, for 360 = 2^(3^1) * 3^(2^1) * 5^1 we get a(360) = 5. a(65536) = a(2^(2^(2^(2^1)))) = 4.

Links

Crossrefs

Cf. A028234, A064372, A067029, A106444, A106491, A106492, A106493.
Cf. A276230 (gives first k such that a(k) = n, i.e., this sequence is a left inverse of A276230).
After n=1 differs from A038548 for the first time at n=24, where A038548(24)=4, while a(24)=3.

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=1, 0,
      add(1+a(i[2]), i=ifactors(n)[2]))
    end:
seq(a(n), n=1..100);  # Alois P. Heinz, Nov 07 2014
  • Mathematica
    a[n_] := a[n] = If[n == 1, 0, Sum[1 + a[i[[2]]], {i,FactorInteger[n]}]]; Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Nov 11 2015, after Alois P. Heinz *)
  • PARI
    A067029(n) = if(n<2, 0, factor(n)[1,2]);
A028234(n) = my(f = factor(n)); if (#f~, f[1, 1] = 1); factorback(f); /* after Michel Marcus */
a(n) = if(n<2, 0, 1 + a(A067029(n)) + a(A028234(n)));
for(n=1, 150, print1(a(n),", ")) \\ Indranil Ghosh, Mar 23 2017, after formula by Antti Karttunen

Formula

Additive with a(p^e) = 1 + a(e).
a(1) = 0; for n > 1, a(n) = 1 + a(A067029(n)) + a(A028234(n)). - Antti Karttunen, Mar 23 2017
Other identities. For all n >= 1:
a(A276230(n)) = n.
a(n) = A106493(A106444(n)).
a(n) = A106491(n) - A064372(n).

A091221 Number of distinct irreducible polynomials dividing n-th GF(2)[X]-polynomial.

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Jan 03 2004

Keywords

Links

Crossrefs

a(n) = A001221(A091203(n)) = A001221(A091205(n)). A000374(n) = a(A000051(n)).

Programs

  • Maple
    f:= proc(n) local L,P,R,i;
  L:= convert(n,base,2);
  P:= add(L[i]*X^(i-1),i=1..nops(L));
  R:= Factors(P) mod 2;
  nops(R[2]);
end proc:
map(f, [$1.200]); # Robert Israel, Oct 11 2024

A091250 List of primitive irreducible polynomials over GF(2), interpreted as binary vectors, then written in base 10.

Original entry on oeis.org

3, 7, 11, 13, 19, 25, 37, 41, 47, 55, 59, 61, 67, 91, 97, 103, 109, 115, 131, 137, 143, 145, 157, 167, 171, 185, 191, 193, 203, 211, 213, 229, 239, 241, 247, 253, 285, 299, 301, 333, 351, 355, 357, 361, 369, 391, 397, 425, 451, 463, 487, 501, 529, 539, 545
Offset: 1

Views

Author

Antti Karttunen, Jan 03 2004

Keywords

Comments

Or, primitive irreducible polynomials over GF(2) listed in lexicographic order and evaluated at X=2.

Links

Crossrefs

a(n) = A014580(A091249(n)). Cf. A011260, A091252. A007088(a(n)) = A058947(n) (same sequence in binary).

A091212 Composite numbers whose binary representation encodes a polynomial reducible over GF(2).

Original entry on oeis.org

4, 6, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 22, 24, 26, 27, 28, 30, 32, 33, 34, 35, 36, 38, 39, 40, 42, 44, 45, 46, 48, 49, 50, 51, 52, 54, 56, 57, 58, 60, 62, 63, 64, 65, 66, 68, 69, 70, 72, 74, 75, 76, 77, 78, 80, 81, 82, 84, 85, 86, 88, 90, 92, 93, 94, 95, 96, 98, 99
Offset: 1

Views

Author

Antti Karttunen, Jan 03 2004

Keywords

Comments

"Encoded in binary representation" means that a polynomial a(n)*X^n+...+a(0)*X^0 over GF(2) is represented by the binary number a(n)*2^n+...+a(0)*2^0 in Z (where each coefficient a(k) = 0 or 1).
From Reinhard Zumkeller, Jul 05-12 2011, values for maximum n corrected by Antti Karttunen, May 18 2015: (Start)
a(n) = A192506(n) for n <= 36.
a(n) = A175526(n) for n <= 36.
(End)

Links

Crossrefs

Intersection of A002808 and A091242.
Cf. A257688 (complement, either 1, irreducible in GF(2)[X] or prime), A091206 (prime and irreducible), A091209 (prime and reducible), A091214 (nonprime and irreducible).
Cf. A091213, A236861, A235036 (a subsequence, apart from 1).
Differs from both A175526 and A192506 for the first time at n=37, where a(37) = 56, while A175526(37) = A192506(37) = 55, a term missing from here (as 55 encodes a polynomial which is irreducible in GF(2)[X]).
Differs from its subsequence A205783(n+1) for the first time at n=47, where a(47) = 69, while 69 is missing from A205783.

Programs

  • PARI
    isA014580(n)=polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from Charles R Greathouse IV
isA091212(n) = ((n > 1) && !isprime(n) && !isA014580(n));
n = 0; i = 0; while(n < 2^16, n++; if(isA091212(n), i++; write("b091212.txt", i, " ", n)));

Formula

a(n) = A091242(A091213(n)).

A091220 Number of divisors of the n-th GF(2)[X]-polynomial.

Original entry on oeis.org

1, 2, 2, 3, 3, 4, 2, 4, 4, 6, 2, 6, 2, 4, 4, 5, 5, 8, 2, 9, 3, 4, 4, 8, 2, 4, 6, 6, 4, 8, 2, 6, 4, 10, 4, 12, 2, 4, 6, 12, 2, 6, 4, 6, 8, 8, 2, 10, 4, 4, 6, 6, 4, 12, 2, 8, 6, 8, 2, 12, 2, 4, 6, 7, 9, 8, 2, 15, 3, 8, 4, 16, 2, 4, 8, 6, 4, 12, 4, 15, 3, 4, 8, 9, 7, 8, 2, 8, 4, 16, 2, 12, 4, 4, 6, 12, 2
Offset: 1

Views

Author

Antti Karttunen, Jan 03 2004

Keywords

Links

Crossrefs

a(n) = A000005(A091203(n)) = A000005(A091205(n)). Cf. A091257.

Programs

  • PARI
    a(n)=local(p,fm,k);while(n>0,p+=Mod(n,2)*x^k;n\=2;k++);fm=factor(p);prod(k=1,matsize(fm)[1],fm[k,2]+1) \\ Franklin T. Adams-Watters, Jun 22 2010
Previous Showing 11-20 of 77 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.