Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A143210
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A143210 Gray code applied to Eulerian numbers: p(x,n)=(1 - x)^(n + 1)*Sum[k^n*x^k, {k, 0, Infinity}]/x; t0*n,m)=Coefficients(t0(n,m)); t(n,m)=GrayCode(t0(n,m)). +0
1
1, 1, 1, 1, 6, 1, 1, 14, 14, 1, 1, 23, 99, 23, 1, 1, 37, 441, 441, 37, 1, 1, 68, 244, 456, 244, 68, 1, 1, 140, 167, 386, 386, 167, 140, 1, 1, 269, 408, 255, 785, 255, 408, 269, 1, 1, 527, 912, 788, 987, 987, 788, 912, 527, 1 (list; graph; listen)
OFFSET

1,5

COMMENT

Row sums are:

{1, 2, 8, 30, 147, 958, 1082, 1388, 2651, 6430}.

REFERENCES

http://mathworld.wolfram.com/notebooks/Combinatorics/GrayCode.nb; Eric's MathWorld entry http://mathworld.wolfram.com/GrayCode.html.

FORMULA

p(x,n)=(1 - x)^(n + 1)*Sum[k^n*x^k, {k, 0, Infinity}]/x; t0*n,m)=Coefficients(t0(n,m)); t(n,m)=GrayCode(t0(n,m)).

EXAMPLE

{1},

{1, 1},

{1, 6, 1},

{1, 14, 14, 1},

{1, 23, 99, 23, 1},

{1, 37, 441, 441, 37, 1},

{1, 68, 244, 456, 244, 68, 1},

{1, 140, 167, 386, 386, 167, 140, 1},

{1, 269, 408, 255, 785, 255, 408, 269, 1},

{1, 527, 912, 788, 987, 987, 788, 912, 527, 1}

MATHEMATICA

Code GrayCodeList[k_] := Module[{b = IntegerDigits[k, 2], i}, Do[ If[b[[i - 1]] == 1, b[[i]] = 1 - b[[i]]], {i, Length[b], 2, -1} ]; b ]; FromGrayCodeList[d_] := Module[{b = d, i, j}, Do[ If[Mod[Sum[b[[j]], {j, i - 1}], 2] == 1, b[[i]] = 1 - b[[i]]], {i, n = Length[d], 2, -1} ]; FromDigits[b, 2] ]; GrayCode[i_, n_] := FromDigits[BitXor @@@ Partition[Prepend[ IntegerDigits[i, 2, n], 0], 2, 1], 2] FromGrayCode[i_, n_] := FromDigits[BitXor[IntegerDigits[i, 2, n], FoldList[ BitXor, 0, Most[IntegerDigits[i, 2, n]]]], 2]; Clear[f, a, n, m, x]; (*A123125*) f[x_, n_] := f[x, n] = (1 - x)^(n + 1)*Sum[k^n*x^k, {k, 0, Infinity}]/x; Table[FullSimplify[ExpandAll[f[x, n]]], {n, 0, 10}]; a = Table[CoefficientList[FullSimplify[ExpandAll[f[x, n]]], x], {n, 1, 10}]; Table[Flatten[Table[GrayCode[a[[n]][[m]], 10], {m, 1, n}]], {n, 1, Length[ a]}]; Flatten[%]

CROSSREFS

Cf. A123125.

Sequence in context: A081579 A146772 A082105 this_sequence A152238 A086645 A154980

Adjacent sequences: A143207 A143208 A143209 this_sequence A143211 A143212 A143213

KEYWORD

nonn,uned,probation

AUTHOR

Roger L. Bagula and Gary W. Adamson (rlbagulatftn(AT)yahoo.com), Oct 20 2008

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