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-2 of 2 results.

A366308 The number of infinitary divisors of n that are terms of A366242.

Original entry on oeis.org

1, 2, 2, 1, 2, 4, 2, 2, 1, 4, 2, 2, 2, 4, 4, 2, 2, 2, 2, 2, 4, 4, 2, 4, 1, 4, 2, 2, 2, 8, 2, 4, 4, 4, 4, 1, 2, 4, 4, 4, 2, 8, 2, 2, 2, 4, 2, 4, 1, 2, 4, 2, 2, 4, 4, 4, 4, 4, 2, 4, 2, 4, 2, 2, 4, 8, 2, 2, 4, 8, 2, 2, 2, 4, 2, 2, 4, 8, 2, 4, 2, 4, 2, 4, 4, 4, 4
Offset: 1

Views

Author

Amiram Eldar, Oct 06 2023

Keywords

Crossrefs

Programs

  • Mathematica
    s[0] = 0; s[n_] := s[n] = s[Floor[n/4]] + If[OddQ[Mod[n, 4]], 1, 0]; f[p_, e_] := 2^s[e]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    s(e) = if(e > 3, s(e\4)) + e%2 \\ after Charles R Greathouse IV at A139351
    a(n) = vecprod(apply(x -> 2^s(x), factor(n)[, 2]));

Formula

Multiplicative with a(p^e) = 2^A139351(e).
a(n) = 2^A366246(n).
a(n) = A037445(n)/A366309(n).
a(n) = A037445(A366244(n)).
a(n) >= 1, with equality if and only if n is in A366243.
a(n) <= A037445(n), with equality if and only if n is in A366242.

A370650 Numbers whose number of infinitary divisors that are terms of A366242 is equal to the number of infinitary divisors that are terms of A366243.

Original entry on oeis.org

1, 8, 12, 18, 20, 27, 28, 44, 45, 50, 52, 63, 64, 68, 75, 76, 92, 98, 99, 116, 117, 124, 125, 144, 147, 148, 153, 164, 171, 172, 175, 188, 207, 212, 216, 236, 242, 244, 245, 261, 268, 275, 279, 284, 292, 316, 324, 325, 332, 333, 338, 343, 356, 360, 363, 369, 387, 388, 400
Offset: 1

Views

Author

Amiram Eldar, Feb 24 2024

Keywords

Comments

Numbers k such that A366308(k) = A366309(k).
Numbers k such that A366246(k) = A366247(k) = A064547(k)/2.
If k is a term, then all the numbers with the same prime signature as k are terms. The least terms with each prime signature are listed in A370651.
p^A039004(k) is a term for all primes p and all k >= 1.

Crossrefs

Programs

  • Mathematica
    s1[0] = 0; s1[n_] := s1[n] = s1[Floor[n/4]] + If[OddQ[Mod[n, 4]], 1, 0]; f1[p_, e_] := s1[e]; a1[1] = 0; a1[n_] := Plus  @@ f1 @@@ FactorInteger[n];
    s2[0] = 0; s2[n_] := s2[n] = s2[Floor[n/4]] + If[Mod[n, 4] > 1, 1, 0]; f2[p_, e_] := s2[e]; a2[1] = 0; a2[n_] := Plus @@ f2 @@@ FactorInteger[n];
    q[n_] := a1[n] == a2[n]; Select[Range[400], q]
Showing 1-2 of 2 results.