login
This site is supported by donations to The OEIS Foundation.
Logo

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A157926 Coefficients of the first modulo two factor of polynomials of the type: p(x,n)=(x^Prime[n] + 1)/(x + 1). 0
1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1 (list; graph; refs; listen; history; internal format)
OFFSET

2,1

COMMENTS

Row sums are:

{3, 5, 7, 3, 15, 11, 15, 15, 3, 23, 5, 17, 27, 21, 15,...}.

This procedure gives more than the types with two factors:

Mod[(x^Prime[n] + 1)/(x + 1),2]=f1(x)*f2(x);

In each case the program picks just the first factor.

The classic Golay factorization:

Factor[PolynomialMod[(x^23 + 1)/((x + 1)), 2], Modulus -> 2]

(1 + x + x^5 + x^6 + x^7 + x^9 + x^11)

(1 + x^2 + x^4 + x^5 + x^6 + x^10 + x^11)

The other one mentioned by Sloane and Conway in "Sphere Packings":

Factor[PolynomialMod[(x^47 + 1)/((x + 1)), 2], Modulus -> 2]

(1 + x + x^2 + x^3 + x^5 + x^6 + x^7 + x^9 + x^10 + x^12 + x^13 + x^14 + x^18 + x^19 + x^23)

(1 + x^4 + x^5 + x^9 +x^10 + x^11 + x^13 + x^14 + x^16 + x^17 + x^18 + x^20 +x^21 + x^22 + x^23)

REFERENCES

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, pp. 231.

FORMULA

Mod[(x^Prime[n] + 1)/(x + 1),2]=Product[fi(x),{i,0,n}];

Out_(n,m)=coefficients(f1(x)).

EXAMPLE

{1, 1, 0, 1},

{1, 0, 0, 1, 1, 1, 0, 0, 1},

{1, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1},

{1, 0, 1, 0, 0, 1},

{1, 0, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 0, 1},

{1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1},

{1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1},

{1, 1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 0, 1},

{1, 1, 0, 0, 0, 0, 0, 0, 0, 1},

{1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1},

{1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1},

{1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1},

{1, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 1},

{1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1},

{1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1}

MATHEMATICA

a = Flatten[Table[If[ ( Factor[PolynomialMod[(x^Prime[n] + 1)/((x + 1)), 2], Modulus -> 2] == PolynomialMod[(x^Prime[n] + 1)/((x + 1)), 2]) /. x -> 2, {}, FactorList[PolynomialMod[( x^Prime[n] + 1)/((x + 1)), 2], Modulus -> 2][[2]][[1]]], {n, 2, 30}]];

Table[CoefficientList[a[[n]], x], {n, 1, Length[a]}];

Flatten[%] Table[Apply[Plus, CoefficientList[a[[n]], x]], {n, 1, Length[a]}];

CROSSREFS

Sequence in context: A079559 A175480 A014577 * A131377 A077049 A124895

Adjacent sequences:  A157923 A157924 A157925 * A157927 A157928 A157929

KEYWORD

nonn,tabf,uned,obsc

AUTHOR

Roger L. Bagula and Gary W. Adamson (rlbagulatftn(AT)yahoo.com), Mar 09 2009

EXTENSIONS

I started to edit this, changing the definition to read: "Irregular triangle read by rows: factorize (x^prime(n)+1)/(x+1) mod 2, then write down coefficients of lowest degree factor, with exponents in increasing order." This would make a perfectly reasonable sequence, which would begin {1,1}, {1,1,1}, {1,1,1,1,1}, {1,1,0,1}, ... Unfortunately this does not match the terms of the present sequence. Possibly the exponents should be 2^(prime(n)-1) rather than prime(n). I did not investiagte this further. - N. J. A. Sloane, Dec 16 2010

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

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

Last modified February 15 03:33 EST 2012. Contains 205694 sequences.