cp's OEIS Frontend

A206706 Triangle read by rows, T(n,k) n>=0, 0<=k<=n; T(0,0) = -1 and for n > 0 T(n,k) = moebius(n,k+1) - moebius(n,k) where moebius(n,k) = mu(floor(n/k)) if k<>0 and k divides n, 0 otherwise; mu=A008683.

%I A206706 #23 Mar 07 2020 09:04:02
%S A206706 -1,1,-1,-1,2,-1,-1,1,1,-1,0,-1,1,1,-1,-1,1,0,0,1,-1,1,-2,0,1,0,1,-1,
%T A206706 -1,1,0,0,0,0,1,-1,0,0,0,-1,1,0,0,1,-1,0,0,-1,1,0,0,0,0,1,-1,1,-2,1,0,
%U A206706 -1,1,0,0,0,1,-1,-1,1,0,0,0,0,0,0,0,0,1,-1,0
%N A206706 Triangle read by rows, T(n,k) n>=0, 0<=k<=n; T(0,0) = -1 and for n > 0 T(n,k) = moebius(n,k+1) - moebius(n,k) where moebius(n,k) = mu(floor(n/k)) if k<>0 and k divides n, 0 otherwise; mu=A008683.
%C A206706 This is a variant of Paul D. Hanna's A123706 which uses a definition given by Mats Granvik. It adds the column T(n,0) = mu(n) at the left hand side of the triangle.
%C A206706 The value T(0,0) was set to -1 to make the triangle invertible as a matrix with uniform signs of the entries of the inverse.
%e A206706 [ 0]                    -1,
%e A206706 [ 1]                   1, -1,
%e A206706 [ 2]                 -1, 2, -1,
%e A206706 [ 3]                -1, 1, 1, -1,
%e A206706 [ 4]              0, -1, 1, 1, -1,
%e A206706 [ 5]             -1, 1, 0, 0, 1, -1,
%e A206706 [ 6]           1, -2, 0, 1, 0, 1, -1,
%e A206706 [ 7]          -1, 1, 0, 0, 0, 0, 1, -1,
%e A206706 [ 8]        0, 0, 0, -1, 1, 0, 0, 1, -1,
%e A206706 [ 9]       0, 0, -1, 1, 0, 0, 0, 0, 1, -1,
%e A206706 The inverse of this triangle as a matrix begins
%e A206706 [-1,  0,  0,  0,  0,  0,  0]
%e A206706 [-1, -1,  0,  0,  0,  0,  0]
%e A206706 [-1, -2, -1,  0,  0,  0,  0]
%e A206706 [-1, -3, -1, -1,  0,  0,  0]
%e A206706 [-1, -4, -2, -1, -1,  0,  0]
%e A206706 [-1, -5, -2, -1, -1, -1,  0]
%e A206706 [-1, -6, -3, -2, -1, -1, -1]
%p A206706 with(numtheory): A206706 := proc(n,k) local moebius;
%p A206706 moebius := (n, k) -> `if`(k<>0 and irem(n,k) = 0, mobius(iquo(n,k)), 0);
%p A206706 moebius(n, k+1) - moebius(n, k) end:
%t A206706 mu[n_, k_] := If[k != 0 && Divisible[n, k], MoebiusMu[n/k], 0];
%t A206706 T[0, 0] = -1; T[n_, k_] /; 0 <= k <= n := mu[n, k+1] - mu[n, k];
%t A206706 Table[T[n, k], {n, 0, 11}, {k, 0, n}] // Flatten (* _Jean-François Alcover_, Jul 12 2019 *)
%o A206706 (Sage)
%o A206706 def mur(n,k): return moebius(n//k) if k != 0 and n%k == 0 else 0
%o A206706 def A206706(n,k) : return -1 if n==0 and k==0 else mur(n,k+1) - mur(n,k)
%Y A206706 Cf. A123706, A008683.
%K A206706 sign,tabl
%O A206706 0,5
%A A206706 _Peter Luschny_, Feb 11 2012