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.

A249770 Irregular triangle read by rows: T(n,k) is the number of Abelian groups of order n with k invariant factors (2 <= n, 1 <= k).

This page as a plain text file.
%I A249770 #38 Dec 29 2018 13:02:17
%S A249770 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,1,1,1,1,
%T A249770 1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,1,1,1,1,1,1,3,1,1,1,1,1,1,1,1,1,1,1,1,
%U A249770 1,1,1
%N A249770 Irregular triangle read by rows: T(n,k) is the number of Abelian groups of order n with k invariant factors (2 <= n, 1 <= k).
%C A249770 The length of n-th row is A051903(n) and its last element is A249773(A101296(n)).
%C A249770 T(n,k) depends only on k and the prime signature of n.
%H A249770 Álvar Ibeas, <a href="/A249770/b249770.txt">Rows n=2..58654, flattened</a>
%F A249770 T(n,k) = A249771(A101296(n),k).
%F A249770 T(n,1) = 1. If k > 1 and n = Product(p_i^e_i), T(n,k) = Sum(Product(A008284(e_i,k), i in I) * Product(A026820(e_i,k-1), i not in I)), where the sum is taken over nonempty subsets I of {1,...,omega(n)}.
%F A249770 If p is prime and gcd(p,n) = 1, T(pn,k) = T(n,k).
%F A249770 Dirichlet g.f. of column sums: zeta(s)zeta(2s)···zeta(ms) = 1 + Sum_{n >= 2} (Sum_{k=1..m} T(n,k)) / n^s.
%F A249770 T(n,1) + T(n,2) = A046951(n)
%e A249770 First rows:
%e A249770 1;
%e A249770 1;
%e A249770 1,1;
%e A249770 1;
%e A249770 1;
%e A249770 1;
%e A249770 1,1,1;
%e A249770 1,1;
%e A249770 1;
%e A249770 1;
%e A249770 1,1;
%e A249770 1;
%e A249770 1;
%e A249770 1;
%e A249770 1,2,1,1;
%e A249770 1;
%e A249770 ...
%t A249770 f[{x_, y_}] := x^IntegerPartitions[y];
%t A249770 g[n_] := FactorInteger[n][[1, 1]];
%t A249770 h[list_] := Apply[Times,Map[PadRight[#, Max[Map[Length, SplitBy[list, g]]], 1] &,SplitBy[list, g]]]; t[list_] := Tally[Map[Length, list]][[All, 2]];
%t A249770 Map[t, Table[Map[h, Join @@@ Tuples[Map[f, FactorInteger[n]]]], {n, 2, 50}]] // Grid (* _Geoffrey Critzer_, Nov 26 2015 *)
%Y A249770 Refinement of A000688.
%Y A249770 Cf. A008284, A051903, A026820, A101296, A249771, A249773, A046951, A264809.
%K A249770 nonn,tabf
%O A249770 2,21
%A A249770 _Álvar Ibeas_, Nov 06 2014