login
a(n,m) = (binomial(n,m) mod 2)*Gray_Code(n,m).
0

%I #12 Jan 30 2023 07:40:12

%S 1,1,1,1,0,1,1,0,0,1,1,0,0,0,1,1,1,0,0,1,1,1,0,0,0,0,0,1,1,0,0,0,0,0,

%T 0,1,1,0,0,0,0,0,0,0,1,1,1,0,0,0,0,0,0,1,1,1,0,1,0,0,0,0,0,1,0,1,1,0,

%U 0,1,0,0,0,0,1,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,1,1,1,0,0,0,0,0,0,0,0,0,0,1,1

%N a(n,m) = (binomial(n,m) mod 2)*Gray_Code(n,m).

%F a(n,m) = (binomial(n,m) mod 2)*Gray_Code(n,m).

%e {1},

%e {1, 1},

%e {1, 0, 1},

%e {1, 0, 0, 1},

%e {1, 0, 0, 0, 1},

%e {1, 1, 0, 0, 1, 1},

%e {1, 0, 0, 0, 0, 0, 1},

%e {1, 0, 0, 0, 0, 0, 0, 1},

%e {1, 0, 0, 0, 0, 0, 0, 0, 1},

%e {1, 1, 0, 0, 0, 0, 0, 0, 1, 1},

%e {1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1},

%e {1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1},

%e {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1},

%e {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},

%e {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1},

%e {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1}

%t c[i_, k_] := Floor[Mod[i/2^k, 2]];

%t b[i_, k_] = If[c[i, k] == 0 && c[i, k + 1] == 0, 0, If[c[i, k] == 1 && c[i, k + 1] == 1, 0, 1]];

%t n = 15

%t a0 = Table[If[Sum[b[i, k]*b[j, k], {k, 0, n}] == 0, 1, 0], {j, 0, n}, {i, 0, n}];

%t ListDensityPlot[a0, Mesh -> False];

%t c = Delete[Table[Reverse[Table[a0[[n, l - n]], {n, 1, l - 1}]], {l, 1, Dimensions[a0][[1]] + 1}], 1];

%t Flatten[c];

%t Dimensions[c];

%t d = Table[Table[Mod[Binomial[n0, m], 2], {m, 0, n0}], {n0, 0, n}]

%t e = Table[Table[c[[n0, m]]*d[[n0, m]], {m, 1, n0}], {n0, 1, n + 1}]

%t Flatten[e]

%Y Cf. A047999.

%K nonn,tabl,uned

%O 1,1

%A _Roger L. Bagula_, Sep 27 2007