|
%I #5 Mar 30 2012 18:57:40
%S 1,3,2,7,4,5,12,8,9,6,19,13,14,10,11,27,20,21,15,16,17,36,28,29,22,23,
%T 24,18,47,37,38,30,31,32,25,26,59,48,49,39,40,41,33,34,35,73,60,61,50,
%U 51,52,42,43,44,45,88,74,75,62,63,64,53,54,55,56,46,104,89,90
%N Natural interspersion of A194102; a rectangular array, by antidiagonals.
%C See A194029 for definitions of natural fractal sequence and natural interspersion. Every positive integer occurs exactly once (and every pair of rows intersperse), so that as a sequence, A194100 is a permutation of the positive integers; its inverse is A194101.
%e Northwest corner:
%e 1...3...7...12...19
%e 2...4...8...13...20
%e 5...9...14..21...29
%e 6...10..15..22...30
%e 11..16..23..31...40
%t z = 40; g = Sqrt[2];
%t c[k_] := Sum[Floor[j*g], {j, 1, k}];
%t c = Table[c[k], {k, 1, z}] (* A194102 *)
%t f[n_] := If[MemberQ[c, n], 1, 1 + f[n - 1]]
%t f = Table[f[n], {n, 1, 800}] (* A194103 new *)
%t r[n_] := Flatten[Position[f, n]]
%t t[n_, k_] := r[n][[k]]
%t TableForm[Table[t[n, k], {n, 1, 8}, {k, 1, 7}]]
%t p = Flatten[Table[t[k, n - k + 1], {n, 1, 16}, {k, 1, n}]] (* A194104 *)
%t q[n_] := Position[p, n]; Flatten[Table[q[n], {n, 1, 80}]] (* A194105 *)
%Y Cf. A194029, A194102, A194103, A194105.
%K nonn,tabl
%O 1,2
%A _Clark Kimberling_, Aug 15 2011
|
