A127170 Triangle read by rows: T(n,k) is the number of divisors of n that are divisible by k, with 1 <= k <= n.

1, 2, 1, 2, 0, 1, 3, 2, 0, 1, 2, 0, 0, 0, 1, 4, 2, 2, 0, 0, 1, 2, 0, 0, 0, 0, 0, 1, 4, 3, 0, 2, 0, 0, 0, 1, 3, 0, 2, 0, 0, 0, 0, 0, 1, 4, 2, 0, 0, 2, 0, 0, 0, 0, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 6, 4, 3, 2, 0, 2, 0, 0, 0, 0, 0, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 4, 2, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 1

Offset

1,2

Comments

Column k lists the terms of A000005 interleaved with k - 1 zeros.

Eigensequence of the triangle = A007557; i.e., sequence A007557 shifts to the left upon multiplication by A127170. - Gary W. Adamson, Apr 27 2009

Links

Michael De Vlieger, Table of n, a(n) for n = 1..11325 (rows n = 1..150, flattened)

Formula

A007429(n) = Sum_{i=1..n} i*a(i).

T(n,k) = A000005(n/k), if k divides n, otherwise 0, with n >= 1 and 1 <= k <= n. - Omar E. Pol, Apr 01 2015

Example

First 10 rows of the triangle:
  1;
  2, 1;
  2, 0, 1;
  3, 2, 0, 1;
  2, 0, 0, 0, 1;
  4, 2, 2, 0, 0, 1;
  2, 0, 0, 0, 0, 0, 1;
  4, 3, 0, 2, 0, 0, 0, 1;
  3, 0, 2, 0, 0, 0, 0, 0, 1;
  4, 2, 0, 0, 2, 0, 0, 0, 0, 1;
  ...

Maple

T:= (n, k)-> `if`(irem(n, k)=0, numtheory[tau](n/k), 0):
seq(seq(T(n, k), k=1..n), n=1..14);  # Alois P. Heinz, Feb 16 2022

Mathematica

Table[Function[D, Table[Count[D, _?(Mod[#, k] == 0 &)], {k, n}]]@ Divisors[n], {n, 12}] // Flatten (* Michael De Vlieger, Feb 16 2022 *)

Prog

(PARI) tabl(nn) = {m = matrix(nn, nn, n, k, if ((n % k) == 0, 1, 0)); m = m^2; for (n=1, nn, for (k=1, n, print1(m[n, k], ", "); ); print(); ); } \\ Michel Marcus, Apr 01 2015

Crossrefs

Row sums give A007425.

Cf. A000005, A007557, A051731, A126988, A007429, A127170, A195050.

Sequence in context: A190491 A364835 A143352 * A143335 A099505 A156837

Adjacent sequences: A127167 A127168 A127169 * A127171 A127172 A127173

Keyword

nonn,tabl

Author

Gary W. Adamson, Jan 06 2007

Extensions

8 terms taken from Example section and then corrected in Data section by Omar E. Pol, Mar 30 2015

Extended beyond a(21) by Omar E. Pol, Apr 01 2015

New name (which was a comment dated Mar 30 2015) from Omar E. Pol, Feb 16 2022

Status

approved