OFFSET
1,7
EXAMPLE
This triangle is defined by:
* T(n,k) = A100995(n/k) if k|n, 0 otherwise.
Sequence A100995 is defined by:
* A100995(n) = m if n = p^m for some prime p, 0 otherwise.
* A054525(n,k) = MoebiusMu(n/k) if k|n, 0 otherwise.
Triangle A051731 is defined by:
* A051731(n,k) = 1 if k|n, 0 otherwise.
The matrix log of triangle A051731 begins:
0;
1, 0;
1, 0, 0;
1/2, 1, 0, 0;
1, 0, 0, 0, 0;
0, 1, 1, 0, 0, 0;
1, 0, 0, 0, 0, 0, 0;
1/3, 1/2, 0, 1, 0, 0, 0, 0;
1/2, 0, 1, 0, 0, 0, 0, 0, 0;
0, 1, 0, 0, 1, 0, 0, 0, 0, 0;
1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
0, 0, 1/2, 1, 0, 1, 0, 0, 0, 0, 0, 0;
1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0;
0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
1/4, 1/3, 0, 1/2, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0; ...
where all nonzero elements are positive unit fractions.
PROG
(PARI) T(n, k)=if(n%k==0, if(#(factor(n/k)~)==1, factor(n/k)[1, 2], 0), 0)
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Paul D. Hanna, Jan 13 2006
STATUS
approved