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.

Showing 1-4 of 4 results.

A342925 a(n) = A003415(sigma(n)), where A003415 is the arithmetic derivative, and sigma is the sum of divisors of n.

Original entry on oeis.org

0, 1, 4, 1, 5, 16, 12, 8, 1, 21, 16, 32, 9, 44, 44, 1, 21, 16, 24, 41, 80, 60, 44, 92, 1, 41, 68, 92, 31, 156, 80, 51, 112, 81, 112, 20, 21, 92, 92, 123, 41, 272, 48, 124, 71, 156, 112, 128, 22, 34, 156, 77, 81, 244, 156, 244, 176, 123, 92, 332, 33, 272, 164, 1, 124, 384, 72, 165, 272, 384, 156, 119, 39, 101, 128, 188
Offset: 1

Views

Author

Antti Karttunen, Apr 07 2021

Keywords

Links

Crossrefs

Cf. A023194 (positions of ones, which is a subsequence of prime powers, A000961).
Cf. A342922, A342923, A342924, A342926, A351568, A351569, A351570, A351571.
Cf. A342021 (fixed points), A343216 [positions k where a(k) < k], A343217 [a(k) >= k], A343218 [a(k) > k].
Cf. A347870 (parity of terms), A347872, A347873, A347877 (positions of odd terms), A347878 (of even terms), A343218, A343220, A344024.

Programs

  • Mathematica
    Array[If[# < 2, 0, # Total[#2/#1 & @@@ FactorInteger[#]]] &@ DivisorSigma[1, #] &, 76] (* Michael De Vlieger, Apr 08 2021 *)
  • PARI
    A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A342925(n) = A003415(sigma(n));

Formula

a(A023194(n)) = 1.
If gcd(m,n) = 1, a(m*n) = sigma(m)*A003415(sigma(n)) + sigma(n)*A003415(sigma(m)) = sigma(m)*a(n) + sigma(n)*a(m).
a(n) = (A351568(n)*A351571(n)) + (A351569(n)*A351570(n)). - Antti Karttunen, Feb 23 2022

A351569 Sum of divisors of the largest unitary divisor of n that is an exponentially odd number.

Original entry on oeis.org

1, 3, 4, 1, 6, 12, 8, 15, 1, 18, 12, 4, 14, 24, 24, 1, 18, 3, 20, 6, 32, 36, 24, 60, 1, 42, 40, 8, 30, 72, 32, 63, 48, 54, 48, 1, 38, 60, 56, 90, 42, 96, 44, 12, 6, 72, 48, 4, 1, 3, 72, 14, 54, 120, 72, 120, 80, 90, 60, 24, 62, 96, 8, 1, 84, 144, 68, 18, 96, 144, 72, 15, 74, 114, 4, 20, 96, 168, 80, 6, 1, 126, 84
Offset: 1

Views

Author

Antti Karttunen, Feb 23 2022

Keywords

Links

Crossrefs

Cf. A000203, A013662, A028982 (positions of odd terms), A268335 (exponentially odd numbers), A350389, A351568, A351571.
Coincides with A001615 on squarefree numbers, A005117.

Programs

  • Mathematica
    f[p_, e_] := If[OddQ[e], (p^(e + 1) - 1)/(p - 1), 1]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Feb 23 2022 *)
  • PARI
    A350389(n) = { my(m=1, f=factor(n)); for(k=1,#f~,if(1==(f[k,2]%2), m *= (f[k,1]^f[k,2]))); (m); };
A351569(n) = sigma(A350389(n));
  • Python
    from math import prod
from sympy import factorint
def A351569(n): return prod((p**(e+1)-1)//(p-1) if e % 2 else 1 for p, e in factorint(n).items()) # Chai Wah Wu, Feb 24 2022

Formula

Multiplicative with a(p^e) = (p^(e+1)-1)/(p-1) if e is odd and 1 otherwise.
a(n) = A000203(A350389(n)).
a(n) = A000203(n) / A351568(n).
Sum_{k=1..n} a(k) ~ c * n^2, where c = zeta(4)/2 = Pi^4/180 = 0.541161... . - Amiram Eldar, Nov 20 2022
Dirichlet g.f.: zeta(2*s) * zeta(2*s-2) * Product_{p prime} (1 + 1/p^(s-1) + 1/p^s - 1/p^(2*s-2)). - Amiram Eldar, Sep 03 2023

A351570 Arithmetic derivative of the sum of the divisors of the largest unitary divisor of n that is a square.

Original entry on oeis.org

0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 20, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 22, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 22, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 22, 1, 38
Offset: 1

Views

Author

Antti Karttunen, Feb 23 2022

Keywords

Comments

Observation: There seems to be no terms in range 2..19 in this sequence.

Links

Crossrefs

Cf. A000203, A003415, A342925, A350388, A351568, A351571, A351572, A351575 (positions of ones).

Programs

  • PARI
    A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A350388(n) = { my(m=1, f=factor(n)); for(k=1,#f~,if(0==(f[k,2]%2), m *= (f[k,1]^f[k,2]))); (m); };
A351570(n) = A003415(sigma(A350388(n)));

Formula

a(n) = A003415(A351568(n)) = A342925(A350388(n)) = A003415(A000203(A350388(n))).

A351573 Arithmetic derivative of the largest unitary divisor of n that is an exponentially odd number.

Original entry on oeis.org

0, 1, 1, 0, 1, 5, 1, 12, 0, 7, 1, 1, 1, 9, 8, 0, 1, 1, 1, 1, 10, 13, 1, 44, 0, 15, 27, 1, 1, 31, 1, 80, 14, 19, 12, 0, 1, 21, 16, 68, 1, 41, 1, 1, 1, 25, 1, 1, 0, 1, 20, 1, 1, 81, 16, 92, 22, 31, 1, 8, 1, 33, 1, 0, 18, 61, 1, 1, 26, 59, 1, 12, 1, 39, 1, 1, 18, 71, 1, 1, 0, 43, 1, 10, 22, 45, 32, 140, 1, 7, 20, 1, 34
Offset: 1

Views

Author

Antti Karttunen, Feb 23 2022

Keywords

Links

Crossrefs

Cf. A003415, A350388, A268335 (exponentially odd numbers), A351571, A351572.

Programs

  • Mathematica
    f1[p_, e_] := If[OddQ[e], p^e, 1]; f2[p_, e_] := If[OddQ[e], e/p, 0]; a[1] = 0; a[n_] := (Times @@ f1 @@@ (f = FactorInteger[n])) * (Plus @@ f2 @@@ f); Array[a, 100] (* Amiram Eldar, Feb 23 2022 *)
  • PARI
    A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A350389(n) = { my(m=1, f=factor(n)); for(k=1,#f~,if(1==(f[k,2]%2), m *= (f[k,1]^f[k,2]))); (m); };
A351573(n) = A003415(A350389(n));

Formula

a(n) = A003415(A350389(n)).
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.