cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A307223 Irregular table T(n, k) read by rows: n-th row gives number of subsets of the divisors of n which sum to k for 1 <= k <= sigma(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1
Offset: 1

Views

Author

Amiram Eldar, Mar 29 2019

Keywords

Comments

T(n, k) > 0 for all values of k iff n is practical (A005153).

Examples

			Table begins as:
  1
  1,1,1
  1,0,1,1
  1,1,1,1,1,1,1
  1,0,0,0,1,1
  1,1,2,1,1,2,1,1,2,1,1,1
  1,0,0,0,0,0,1,1
  1,1,1,1,1,1,1,1,1,1,1,1,1,1,1
  1,0,1,1,0,0,0,0,1,1,0,1,1
  1,1,1,0,1,1,1,1,0,1,1,1,1,0,1,1,1,1
		

Crossrefs

Cf. A000203 (row lengths), A307224 (row products).

Programs

  • Mathematica
    T[n_,k_] := Module[{d = Divisors[n]}, SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, k}], k]]; Table[T[n, k], {n,1,10}, {k, 1, DivisorSigma[1,n]}] // Flatten

Formula

T(n, n) = A033630(n).
T(n, A030057(n)) = 0 if there is a 0 in the n-th row, i.e. A030057(n) <= sigma(n) or n is not practical.