|
|
A124433
|
|
Irregular array {a(n,m)} read by rows where (sum{n>=1} sum{m=1 to A001222(n)+1} a(n,m)*y^m/n^x) = 1/(zeta(x)-1+1/y) for all x and y where the double sum converges.
|
|
4
|
|
|
1, 0, -1, 0, -1, 0, -1, 1, 0, -1, 0, -1, 2, 0, -1, 0, -1, 2, -1, 0, -1, 1, 0, -1, 2, 0, -1, 0, -1, 4, -3, 0, -1, 0, -1, 2, 0, -1, 2, 0, -1, 3, -3, 1, 0, -1, 0, -1, 4, -3, 0, -1, 0, -1, 4, -3, 0, -1, 2, 0, -1, 2, 0, -1, 0, -1, 6, -9, 4, 0, -1, 1, 0, -1, 2, 0, -1, 2, -1, 0, -1, 4, -3, 0, -1, 0, -1, 6, -6, 0, -1, 0, -1, 4, -6, 4, -1, 0, -1, 2, 0, -1, 2, 0, -1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,13
|
|
COMMENTS
|
Row n has A001222(n)+1 terms. The polynomial P_n(y) = (sum{m=1 to A001222(n)+1} a(n,m)*y^m) is a generalization of the Mobius (Moebius) function, where P_n(1) = A008683(n).
Up to sign, also the number of strict length-k chains of divisors from n to 1, 1 <= k <= 1 + A001222(n). For example, row n = 36 counts the following chains (empty column indicated by dot):
. 36/1 36/2/1 36/4/2/1 36/12/4/2/1
36/3/1 36/6/2/1 36/12/6/2/1
36/4/1 36/6/3/1 36/12/6/3/1
36/6/1 36/9/3/1 36/18/6/2/1
36/9/1 36/12/2/1 36/18/6/3/1
36/12/1 36/12/3/1 36/18/9/3/1
36/18/1 36/12/4/1
36/12/6/1
36/18/2/1
36/18/3/1
36/18/6/1
36/18/9/1
(End)
|
|
LINKS
|
Mohammad K. Azarian, A Double Sum, Problem 440, College Mathematics Journal, Vol. 21, No. 5, Nov. 1990, p. 424. Solution published in Vol. 22. No. 5, Nov. 1991, pp. 448-449.
|
|
FORMULA
|
a(1,1)=1. a(n,1) = 0 for n>=2. a(n,m+1) = -sum{k|n,k < n} a(k,m), where, for the purpose of this sum, a(k,m) = 0 if m > A001222(k)+1.
|
|
EXAMPLE
|
1/(zeta(x) - 1 + 1/y) = y - y^2/2^x - y^2/3^x + ( - y^2 + y^3)/4^x - y^2/5^x + ( - y^2 + 2y^3)/6^x - y^2/7^x + ...
The sequence of rows begins:
1: 1 16: 0 -1 3 -3 1 31: 0 -1
2: 0 -1 17: 0 -1 32: 0 -1 4 -6 4 -1
3: 0 -1 18: 0 -1 4 -3 33: 0 -1 2
4: 0 -1 1 19: 0 -1 34: 0 -1 2
5: 0 -1 20: 0 -1 4 -3 35: 0 -1 2
6: 0 -1 2 21: 0 -1 2 36: 0 -1 7 -12 6
7: 0 -1 22: 0 -1 2 37: 0 -1
8: 0 -1 2 -1 23: 0 -1 38: 0 -1 2
9: 0 -1 1 24: 0 -1 6 -9 4 39: 0 -1 2
10: 0 -1 2 25: 0 -1 1 40: 0 -1 6 -9 4
11: 0 -1 26: 0 -1 2 41: 0 -1
12: 0 -1 4 -3 27: 0 -1 2 -1 42: 0 -1 6 -6
13: 0 -1 28: 0 -1 4 -3 43: 0 -1
14: 0 -1 2 29: 0 -1 44: 0 -1 4 -3
15: 0 -1 2 30: 0 -1 6 -6 45: 0 -1 4 -3
(End)
|
|
MATHEMATICA
|
f[l_List] := Block[{n = Length[l] + 1, c}, c = Plus @@ Last /@ FactorInteger[n]; Append[l, Prepend[ -Plus @@ Pick[PadRight[ #, c] & /@ l, Mod[n, Range[n - 1]], 0], 0]]]; Nest[f, {{1}}, 34] // Flatten(* Ray Chandler, Feb 13 2007 *)
chnsc[n_]:=If[n==1, {{}}, Prepend[Join@@Table[Prepend[#, n]&/@chnsc[d], {d, DeleteCases[Divisors[n], 1|n]}], {n}]];
Table[(-1)^k*Length[Select[chnsc[n], Length[#]==k&]], {n, 30}, {k, 0, PrimeOmega[n]}] (* Gus Wiseman, Aug 24 2020 *)
|
|
CROSSREFS
|
A008480 gives rows ends (up to sign).
A008683 gives row sums (the Moebius function).
A097805 is the restriction to powers of 2 (up to sign).
A251683 is the unsigned version with zeros removed.
A334996 is the unsigned version (except with a(1) = 0).
A334997 is an unsigned non-strict version.
A337107 is the restriction to factorial numbers.
A001222 counts prime factors with multiplicity.
A067824 counts strict chains of divisors starting with n.
A074206 counts strict chains of divisors from n to 1.
A122651 counts strict chains of divisors summing to n.
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts strict chains of divisors.
A337105 counts strict chains of divisors from n! to 1.
|
|
KEYWORD
|
sign,tabf
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|