Sean A. Irvine, <a href="/A049649/b049649_1.txt">Table of n, a(n) for n = 0..500</a> (terms 1..76 from G. C. Greubel)
Sean A. Irvine, <a href="/A049649/b049649_1.txt">Table of n, a(n) for n = 0..500</a> (terms 1..76 from G. C. Greubel)
proposed
approved
editing
proposed
G. CSean A. Greubel, Irvine, <a href="/A049649/b049649_1.txt">Table of n, a(n) for n = 0..75500</a> (terms 1..76 from G. C. Greubel)
proposed
editing
editing
proposed
1, 1, 4, 5, 7, 8, 11, 14, 17, 18, 23, 28, 31, 34, 41, 46, 51, 54, 61, 70, 75, 78, 89, 100, 105, 110, 123, 130, 137, 144, 153, 168, 177, 182, 199, 216, 223, 230, 249, 260, 269, 278, 291, 312, 323, 330, 353, 372, 381, 392, 413, 430, 443, 454, 473, 500, 513, 520, 549, 578, 587, 602
0,23
T[0 | 1, 0 | 1] = 0; T[0 | 1, _] = T[_, 0 | 1] = 1; T[i_, j_] := Module[{slopes, cnt}, slopes = Union@Flatten@Table[(k - j)/(h - i), {h, 0, i - 1}, {k, 0, j - 1}]; cnt[slope_] := Count[Flatten[Table[{h, k}, {h, 0, i - 1}, {k, 0, j - 1}], 1], {h_, k_} /; (k - j)/(h - i) == slope]; Count[cnt /@ slopes, c_ /; c >= 2] + 2]; Table[T[n, n + 3], {n, 1, 0, 50}] (* G. C. Greubel, Dec 05 2017 *)
Missing a(0)=1 inserted by Sean A. Irvine, Aug 06 2021
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed
Cf. A049639.
proposed
editing