login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A257571
Triangular array read by rows: d(h,k) = distance between h and k in the tree at A232558, for h >=0, k = 0..h.
3
0, 1, 0, 2, 1, 0, 3, 2, 1, 0, 3, 2, 1, 2, 0, 4, 3, 2, 3, 1, 0, 4, 3, 2, 1, 3, 4, 0, 5, 4, 3, 2, 4, 5, 1, 0, 4, 3, 2, 3, 1, 2, 4, 5, 0, 5, 4, 3, 4, 2, 3, 5, 6, 1, 0, 5, 4, 3, 4, 2, 1, 5, 6, 3, 4, 0, 6, 5, 4, 5, 3, 2, 6, 7, 4, 5, 1, 0, 5, 4, 3, 2, 4, 5, 1, 2
OFFSET
1,4
COMMENTS
The distance between h and k is the length of the path from h to k in the tree defined from the root 0 by edges from x to x+1 and x to 2x if x is even, and an edge from x to x+1 if x is odd. This is the tree defined at A232558; it is a subtree of the tree defined at A257569.
LINKS
EXAMPLE
First ten rows:
0
1 0
2 1 0
3 2 1 0
3 2 1 2 0
4 3 2 3 1 0
4 3 2 1 3 4 0
5 4 3 2 4 5 1 0
4 3 2 3 1 2 4 5 0
5 4 3 5 2 3 6 7 1 0
d(6,4) = 3 counts the edges in the path 6,3,2,4;
d(46,21) = 6 counts the edges in the path 46,23,22,11,10,20,21.
MATHEMATICA
f[{x_, y_}] := f[x, y] = If[EvenQ[x], {y, x/2}, {x - 1, y}];
g[{x_, y_}] := g[x, y] = Drop[FixedPointList[f, {x, y}], -1];
s[n_] := s[n] = Reverse[Select[Sort[Flatten[Select[g[{n, 0}], #[[2]] == 0 &]]], # > 0 &]];
m[h_, k_] := m[h, k] = Max[Intersection[s[h], s[k]]];
j[h_, k_] := j[h, k] = Join[Select[s[h], # >= m[h, k] &], Reverse[Select[s[k], # > m[h, k] &]]];
d[h_, k_] := d[h, k] = If[k*h == 0, Length[j[h, k]], -1 + Length[j[h, k]]];
TableForm[Table[d[h, k], {h, 0, 59}, {k, 0, 59}]]; (* A257570 array *)
Flatten[Table[d[h - k, k], {h, 0, 59}, {k, 0, h}]]; (* A257570 sequence *)
CROSSREFS
Sequence in context: A074984 A112658 A190693 * A219649 A292160 A025581
KEYWORD
nonn,tabl,easy
AUTHOR
Clark Kimberling, May 01 2015
STATUS
approved