login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
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).
From Gus Wiseman, Aug 24 2020: (Start)
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 + ...
From Gus Wiseman, Aug 24 2020: (Start)
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).
A073093 gives row lengths.
A074206 gives unsigned row sums.
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.
A001055 counts factorizations.
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.
Sequence in context: A352427 A257024 A334996 * A287104 A190483 A090239
KEYWORD
sign,tabf
AUTHOR
Leroy Quet, Dec 15 2006
EXTENSIONS
Extended by Ray Chandler, Feb 13 2007
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 18 22:56 EDT 2024. Contains 370952 sequences. (Running on oeis4.)