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.

A385006 The sum of the biquadratefree divisors of n.

Original entry on oeis.org

1, 3, 4, 7, 6, 12, 8, 15, 13, 18, 12, 28, 14, 24, 24, 15, 18, 39, 20, 42, 32, 36, 24, 60, 31, 42, 40, 56, 30, 72, 32, 15, 48, 54, 48, 91, 38, 60, 56, 90, 42, 96, 44, 84, 78, 72, 48, 60, 57, 93, 72, 98, 54, 120, 72, 120, 80, 90, 60, 168, 62, 96, 104, 15, 84, 144
Offset: 1

Views

Author

Amiram Eldar, Jun 15 2025

Keywords

Comments

First differs from A365682 and A366992 at n = 32.
The number of these divisors is A252505(n), and the largest of them is A058035(n).

Links

Crossrefs

Cf. A046100, A048250, A058035, A073185, A252505, A365682, A366992.
The sum of divisors d of n such that d is: A000593 (odd), A033634 (exponentially odd), A035316 (square), A038712 (power of 2), A048250 (squarefree), A072079 (3-smooth), A073185 (cubefree), A113061 (cube), A162296 (nonsquarefree), A183097 (powerful), A186099 (5-rough), A353900 (exponentially 2^n), A385005 (cubefull), this sequence (biquadratefree).

Programs

  • Mathematica
    f[p_, e_] := (p^Min[e+1, 4] - 1)/(p - 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n), p, e); prod(i = 1, #f~, p = f[i,1]; e = f[i,2]; (p^min(e+1, 4) - 1)/(p - 1));}

Formula

Multiplicative with a(p^e) = (p^min(e+1, 4) - 1)/(p - 1).
In general, the sum of the k-free (numbers that are not divisible by a k-th power larger than 1) divisors of n is multiplicative with a(p^e) = (p^min(e+1, k) - 1)/(p - 1).
Dirichlet g.f.: zeta(s) * zeta(s-1) /zeta(4*s-4).
In general, the sum of the k-free divisors of n has Dirichlet g.f.: zeta(s)*zeta(s-1)/zeta(k*s-k).
Sum_{k=1..n} a(k) ~ (15/(2*Pi^2)) * n^2.
In general, the sum of the k-free divisors of n has an average order (Pi^2/(12*zeta(k))) * n^2.

A366989 The number of prime powers p^q dividing n, where p is prime and q is either 1 or prime (A334393 without the first term 1).

Original entry on oeis.org

0, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 3, 1, 2, 2, 3, 1, 3, 1, 3, 2, 2, 1, 4, 2, 2, 3, 3, 1, 3, 1, 4, 2, 2, 2, 4, 1, 2, 2, 4, 1, 3, 1, 3, 3, 2, 1, 4, 2, 3, 2, 3, 1, 4, 2, 4, 2, 2, 1, 4, 1, 2, 3, 4, 2, 3, 1, 3, 2, 3, 1, 5, 1, 2, 3, 3, 2, 3, 1, 4, 3, 2, 1, 4, 2, 2, 2
Offset: 1

Views

Author

Amiram Eldar, Oct 31 2023

Keywords

Comments

First differs from A122810 at n = 48, and from A318322 at n = 64.

Links

Crossrefs

Cf. A000720, A001221, A005117, A077761, A334393, A366988, A366991, A366992, A366994.
Cf. A122810, A318322.

Programs

  • Mathematica
    f[p_, e_] := PrimePi[e] + 1; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n)); sum(i = 1, #f~, 1 + primepi(f[i, 2]));}

Formula

Additive with a(p^e) = A000720(e) + 1.
a(n) = 1 is and only if n is squarefree (A005117) > 1.
a(n) = A366988(n) + A001221(n).
Sum_{k=1..n} a(k) ~ n * (log(log(n)) + B + C), where B is Mertens's constant (A077761), C = Sum_{p prime} P(p) = 0.67167522222173297323..., and P(s) is the prime zeta function.

A366991 The number of divisors of n that are not terms of A322448.

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 4, 3, 4, 2, 6, 2, 4, 4, 4, 2, 6, 2, 6, 4, 4, 2, 8, 3, 4, 4, 6, 2, 8, 2, 5, 4, 4, 4, 9, 2, 4, 4, 8, 2, 8, 2, 6, 6, 4, 2, 8, 3, 6, 4, 6, 2, 8, 4, 8, 4, 4, 2, 12, 2, 4, 6, 5, 4, 8, 2, 6, 4, 8, 2, 12, 2, 4, 6, 6, 4, 8, 2, 8, 4, 4, 2, 12, 4, 4
Offset: 1

Views

Author

Amiram Eldar, Oct 31 2023

Keywords

Comments

First differs from A365680 at n = 64.
The number of divisors of n whose prime factorization has exponents that are all either 1 or primes.
The sum of these divisors is A366992(n) and the largest of them is A366994(n).

Links

Crossrefs

Cf. A000005, A000720, A046100, A365680, A366989, A366992, A366994.

Programs

  • Mathematica
    f[p_, e_] := PrimePi[e] + 2; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, primepi(f[i, 2]) + 2);}

Formula

Multiplicative with a(p^e) = A000720(e) + 2.
a(n) <= A000005(n), with equality if and only if n is a biquadratefree number (A046100).

A366994 The largest divisor of n that is not a term of A322448.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 8, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 24, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 32, 65, 66, 67, 68
Offset: 1

Views

Author

Amiram Eldar, Oct 31 2023

Keywords

Comments

First differs from A365683 at n = 64.
The largest divisor of n whose prime factorization has exponents that are all either 1 or primes.
The number of these divisors is A366991(n) and their sum is A366992(n).

Links

Crossrefs

Cf. A007917, A322448, A365683, A366989, A366991, A366992.

Programs

  • Mathematica
    f[p_, e_] := p^If[e == 1, 1, NextPrime[e+1, -1]]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i, 2] == 1, f[i, 1], f[i, 1]^precprime(f[i, 2])));}

Formula

Multiplicative with a(p) = p and a(p^e) = p^A007917(e) for e >= 2.
a(n) <= n, with equality if and only if n is not in A322448.
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{p prime} f(1/p) = 0.48535795387619596052..., where f(x) = (1 - x) * (1 + Sum_{k>=1} x^(2*k-s(k))), s(k) = A007917(k) for k >= 2, and s(1) = 1.
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.