Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A158202
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A158202 Coefficients of polynomials from matrix characteristic polynomials: m(n,m,d)=If[ m <= n, Mod[Binomial[n, m], 2], 0]; M(n)=m(n,m,d).Transpose[m(n,m,d)].Transpose[m(n,m,d)].m(n,m,d) +0
1
1, 1, -1, 1, -2, 1, 1, -11, 11, -1, 1, -12, 22, -12, 1, 1, -21, 106, -106, 21, -1, 1, -30, 255, -708, 255, -30, 1, 1, -91, 1065, -3963, 3963, -1065, 91, -1, 1, -92, 1156, -5028, 7926, -5028, 1156, -92, 1, 1, -101, 1880, -12688, 34482, -34482, 12688, -1880, 101 (list; table; graph; listen)
OFFSET

0,5

COMMENT

Row sums are:

{1, 0, 0, 0, 0, 0, -256, 0, 0, 0, 0,...}.

Absolute value row sums are:

{1, 2, 4, 24, 48, 256, 1280, 10240, 20480, 98304, 405504,...}.

Matrix example is:

M(3)={{3, 2, 2},

{2, 3, 2},

{6, 6, 5}}.

I've only managed to get three levels of these:

(1,5,5,1),(1,7,7,1),(1,11,11,1)

but they seem stable that you get them with combinations of the binomial

matrix, the Transpose[] and Reverse[].

Plotting them shows that they aren't Sierpinski-Pascal at modulo two.

FORMULA

m(n,m,d)=If[ m <= n, Mod[Binomial[n, m], 2], 0];

M(n)=m(n,m,d).Transpose[m(n,m,d)].Transpose[m(n,m,d)].m(n,m,d);

Out_(n,m)=coefficients(characteristicpolynomial(M(n),x),x)

EXAMPLE

{1},

{1, -1},

{1, -2, 1},

{1, -11, 11, -1},

{1, -12, 22, -12, 1},

{1, -21, 106, -106, 21, -1},

{ 1, -30, 255, -708, 255, -30, 1},

{1, -91, 1065, -3963, 3963, -1065, 91, -1},

{1, -92, 1156, -5028, 7926, -5028, 1156, -92, 1},

{1, -101, 1880, -12688, 34482, -34482, 12688, -1880, 101, -1},

{1, -110, 2669, -25128, 98706, -152276, 98706, -25128, 2669, -110, 1}

MATHEMATICA

Clear[M, T, d, a, x, a0];

T[n_, m_, d_] := If[ m <= n, Mod[Binomial[n, m], 2], 0];

M[d_] := Table[T[n, m, d], { n, 1, d}, {m, 1, d}].Transpose[Table[T[ n, m, d], {n, 1, d}, { m, 1, d}]].Transpose[Table[T[n, m, d], {n, 1, d}, {m, 1, d}]].Table[ T[n, m, d], {n, 1, d}, {m, 1, d}];

a0 = Table[M[d], {d, 1, 10}];

Table[Det[M[d]], {d, 1, 10}];

Table[CharacteristicPolynomial[M[d], x], {d, 1, 10}];

a = Join[{{1}}, Table[CoefficientList[ Expand[CharacteristicPolynomial[M[n], x]], x], {n, 1, 10}]];

Flatten[a]

Join[{1}, Table[Apply[Plus, CoefficientList[ Expand[CharacteristicPolynomial[M[n], x]], x]], {n, 1, 10}]];

CROSSREFS

Sequence in context: A064307 A165883 A110905 this_sequence A066094 A010246 A156885

Adjacent sequences: A158199 A158200 A158201 this_sequence A158203 A158204 A158205

KEYWORD

sign,tabl,uned

AUTHOR

Roger L. Bagula (rlbagulatftn(AT)yahoo.com), Mar 14 2009

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified November 25 20:09 EST 2009. Contains 167514 sequences.


AT&T Labs Research