|
|
A341526
|
|
Numerator of ratio n*sigma(A003961(n)) / sigma(n)*A003961(n), where sigma is the sum of divisors of n, and A003961 shifts the prime factorization of n one step towards larger primes.
|
|
11
|
|
|
1, 8, 9, 52, 20, 4, 21, 64, 279, 160, 77, 26, 117, 28, 6, 1936, 170, 248, 114, 1040, 189, 308, 115, 32, 1425, 104, 1053, 26, 464, 16, 589, 1664, 231, 1360, 10, 124, 777, 304, 1053, 1280, 902, 42, 516, 22, 372, 230, 423, 968, 343, 3800, 17, 676, 530, 468, 110, 224, 513, 3712, 1829, 104, 2074, 589, 5859, 69952, 780, 154
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
|
|
MATHEMATICA
|
f[p_, e_] := NextPrime[p]^e; g[1] = 1; g[n_] := Times @@ f @@@ FactorInteger[n]; a[n_] := Numerator[n*DivisorSigma[1, (gn = g[n])]/(DivisorSigma[1, n] * gn)]; Array[a, 100] (* Amiram Eldar, Feb 17 2021 *)
|
|
PROG
|
(PARI)
A003961(n) = { my(f=factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); }; \\ From A003961
A341526(n) = { my(s=A003961(n)); numerator((sigma(s)*n)/(sigma(n)*s)); };
|
|
CROSSREFS
|
Cf. A341626 (same sequence as applied onto prime shift array A246278).
|
|
KEYWORD
|
nonn,frac
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|