OFFSET
1,16
COMMENTS
LINKS
Michel Marcus, Table of n, a(n) for n = 1..5050
EXAMPLE
The first 6 terms of A017666 are 1, 2, 3, 4, 5, 1 where 1 appears twice, 2 to 5 appear once and 6 is absent; giving the 6th row: 2, 1, 1, 1, 1, 0.
Triangle starts
1;
1, 1;
1, 1, 1;
1, 1, 1, 1;
1, 1, 1, 1, 1;
2, 1, 1, 1, 1, 0;
2, 1, 1, 1, 1, 0, 1;
2, 1, 1, 1, 1, 0, 1, 1;
2, 1, 1, 1, 1, 0, 1, 1, 1;
2, 1, 1, 1, 2, 0, 1, 1, 1, 0;
...
MATHEMATICA
Table[Length@ Select[Range@ n, Denominator[DivisorSigma[-1, #]] == k &], {n, 13}, {k, n}] // Flatten (* Michael De Vlieger, Sep 22 2015 *)
PROG
(PARI) tabl(nn) = {vds = vector(nn, n, denominator(sigma(n, -1))); for (n=1, nn, vin = vector(n, k, vds[k]); rown = vector(n, k, #select(x->x==k, vin)); for(k=1, n, print1(rown[k], ", ")); print(); ); }
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Michel Marcus, Sep 22 2015
STATUS
approved