Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #24 Jun 07 2019 11:05:00
%S 1,2,2,3,6,6,3,4,12,6,6,12,4,5,20,12,15,10,10,15,12,20,5,6,30,20,12,
%T 15,10,10,15,12,20,30,6,7,42,30,20,28,21,15,35,14,14,35,15,21,28,20,
%U 30,42,7,8,56,42,30,20,28,21,24,40,35,14,14,35,40,24,21,28,20,30,42,56,8
%N Irregular triangle read by rows: T(n,k) = Farey(n,k+1) - Farey(n,k) where Farey(n,k) = A006842(n,k)/A006843(n,k).
%C This is also the product of the denominators of pairs of consecutive terms in the Farey sequence.
%C Each term of this sequence is an integer: (Proof by induction)
%C Assume that the reciprocal of Farey differences of order n are the product of the consecutive denominators, i.e., if x/y and c/d are adjacent, then |x/y - c/d| = 1/dy. Let a/b and p/q be adjacent in Farey sequence up to n, such that n+1 = b+q (so only their mediant is in the middle).
%C As |a/b - p/q| = 1/bq, |aq - bp| = 1, so |aq - bp + ab - ab| = 1, so |a/b - (a+p)/(b+q)| = 1. The base case is trivial. QED
%F T(n,k) = Farey(n,k+1) - Farey(n,k) with Farey(n,k) = A006842(n,k)/A006843(n,k).
%F T(n,k) = A006843(n,k)*A006843(n,k+1).
%e T(1,1) = 1/(1 - 0);
%e T(2,1) = 1/(1/2 - 0);
%e T(2,2) = 1/(1 - 1/2);
%e T(3,1) = 1/(1/3 - 0);
%e T(3,2) = 1/(1/2 - 1/3);
%e T(3,3) = 1/(2/3 - 1/2);
%e T(3,4) = 1/(1 - 2/3);
%e ...
%e If written as an array:
%e 1;
%e 2, 2;
%e 3, 6, 6, 3;
%e 4, 12, 6, 6, 12, 4;
%e 5, 20, 12, 15, 10, 10, 15, 12, 20, 5;
%e ...
%o (PARI) rowf(n) = {my(vf = [0]); for (k=1, n, for (m=1, k, vf = concat(vf, m/k); ); ); vecsort(Set(vf));} \\ A006842/A006843
%o row(n) = my(vf = rowf(n)); vector(#vf-1, k, 1/(vf[k+1] - vf[k])); \\ _Michel Marcus_, Jun 07 2019
%Y Cf. A006842, A006843.
%K nonn,tabf
%O 1,2
%A _Isaac Kaufmann_, May 30 2019
%E More terms from _Michel Marcus_, Jun 07 2019