A122081 Duplicate of A007140.
1, 1, 3, 14, 115, 2086, 101791, 14835870, 6852422567, 10338780759514, 51804974736769271
Offset: 0
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.
Triangles A106498 and A123547 begin: n = 0 k = 0 : 1, 1 Total = 1, 1 n = 1 k = 0 : 1, 0 k = 1 : 1, 1 Total = 2, 1 n = 2 k = 0 : 1, 0 k = 1 : 1, 0 k = 2 : 2, 1 k = 3 : 1, 1 k = 4 : 1, 1 Totals = 6, 3 n = 3 k = 0 : 1, 0 k = 1 : 1, 0 k = 2 : 2, 0 k = 3 : 4, 1 k = 4 : 5, 2 k = 5 : 5, 4 k = 6 : 4, 3 k = 7 : 2, 2 k = 8 : 1, 1 k = 9 : 1, 1 Totals = 26, 14
permcount[v_] := Module[{m = 1, s = 0, k = 0, t}, For[i = 1, i <= Length[v], i++, t = v[[i]]; k = If[i > 1 && t == v[[i - 1]], k + 1, 1]; m *= t*k; s += t]; s!/m]; edges[v_] := Sum[GCD[v[[i]], v[[j]]], {i, 2, Length[v]}, {j, 1, i - 1}] + Total @ Quotient[v + 1, 2]; A122082[n_] := Module[{s = 0}, Do[s += permcount[p]*2^edges[p], {p, IntegerPartitions[n]}]; s/n!]; a[n_] := A122082[n] - A122082[n-1]; a /@ Range[0, 17] (* Jean-François Alcover, Sep 05 2019, after Andrew Howroyd in A122082 *)
Comments