A260426 a(1) = 1, a(A206074(n)) = A014580(a(n)), a(A205783(1+n)) = A091242(a(n)), where A014580 [respectively A091242] give binary codes for irreducible [resp. reducible] polynomials over GF(2), while A206074 and A205783 give similar codes for polynomials with coefficients 0 or 1 that are irreducible [resp. reducible] over Q.
1, 2, 3, 4, 7, 5, 11, 6, 8, 12, 25, 9, 13, 17, 10, 14, 47, 18, 19, 34, 15, 20, 31, 24, 55, 16, 21, 62, 137, 26, 37, 27, 45, 22, 28, 42, 59, 33, 71, 23, 87, 29, 41, 79, 166, 35, 61, 49, 36, 58, 30, 38, 319, 54, 91, 76, 44, 89, 97, 32, 203, 108, 39, 53, 99, 200, 67, 46, 103, 78, 185, 64, 131, 48, 75, 40, 379, 50, 73, 373, 109, 70, 433, 113, 95, 57, 1123, 111, 143, 121
Offset: 1
Links
Crossrefs
Programs
-
PARI
allocatemem(234567890); vecsize = (2^24)-4; uplim = 2^25; v014580 = vector(vecsize); A014580 = n -> v014580[n]; v091242 = vector(vecsize); A091242 = n -> v091242[n]; v255574 = vector(vecsize); A255574 = n -> v255574[n]; A255572 = n -> (n - A255574(n) - 1); isA014580(n) = polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from Charles R Greathouse IV isA206074(n) = polisirreducible(Pol(binary(n))); A257000 = n -> isA206074(n); v255574[1] = 0; i=0; j=0; n=2; while((n < uplim), if(!(n%65536),print1(n,", ")); v255574[n] = v255574[n-1]+A257000(n); if(isA014580(n), i++; v014580[i] = n, j++; v091242[j] = n); n++); print(n); A260426(n) = if(1==n, 1, if(isA206074(n), A014580(A260426(A255574(n))), A091242(A260426(A255572(n))))); for(n=1, 7968, write("b260426.txt", n, " ", A260426(n)));
Comments