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

A384914 The number of unordered factorizations of n into numbers of the form p^(k^2) where p is prime and k >= 0 (A323520).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Amiram Eldar, Jun 12 2025

Keywords

Comments

First differs from A203640, A295658 and A365333 at n = 64, from A043289 and A053164 at n = 81, and from A063775 at n = 512.

Examples

			a(16) = 2 since 4 has 2 factorizations: 2^1 * 2^1 * 2^1 * 2^1 and 2^4, with exponents 1 and 4 that are squares.
		

Crossrefs

Programs

  • Mathematica
    s[n_] := s[n] = If[n == 0, 1, Sum[Sum[d * Boole[IntegerQ[Sqrt[d]]], {d, Divisors[j]}] * s[n-j], {j, 1, n}] / n];
    f[p_, e_] := s[e]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    s(n) = if(n < 1, 1, sum(j = 1, n, sumdiv(j, d, d*issquare(d)) * s(n-j))/n);
    a(n) = vecprod(apply(s, factor(n)[, 2]));

Formula

Multiplicative with a(p^e) = A001156(e).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} f(1/p) = 1.08451356983124311685..., where f(x) = (1-x) / Product_{k>=1} (1-x^(k^2)).

A295657 Multiplicative with a(p^e) = p^floor((e-1)/2).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 2, 1
Offset: 1

Views

Author

Antti Karttunen, Nov 28 2017

Keywords

Crossrefs

Programs

  • Mathematica
    Array[Apply[Times, FactorInteger[#] /. {p_, e_} /; p > 0 :> p^Floor[(e - 1)/2]] &, 105] (* Michael De Vlieger, Nov 28 2017 *)
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, f[i,1]^floor((f[i,2]-1)/2));} \\ Amiram Eldar, Nov 30 2022

Formula

a(1) = 1; for n > 1, a(n) = A020639(n)^A004526(A067029(n)-1) * a(A028234(n)).
a(n) = A000188(A003557(n)).
a(n) = 1 iff A212793(n) = 1.
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 15/Pi^2 = 1.519817... (A082020). - Amiram Eldar, Nov 30 2022

A365333 The number of exponentially odd coreful divisors of the largest square dividing n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Amiram Eldar, Sep 01 2023

Keywords

Comments

First differs from A043289, A053164, A063775, A203640 and A295658 at n = 64.
The number of squares dividing the largest exponentially odd divisor of n is A325837(n).
The sum of the exponentially odd divisors of the largest square dividing n is A365334(n). [corrected, Sep 08 2023]
The number of exponentially odd divisors of the largest square dividing n is the same as the number of squares dividing n, A046951(n). - Amiram Eldar, Sep 08 2023

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := Max[1, Floor[e/2]]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = vecprod(apply(x -> max(1, x\2), factor(n)[, 2]));

Formula

a(n) = A325837(A008833(n)).
a(n) = 1 if and only if n is a biquadratefree number (A046100).
Multiplicative with a(p^e) = max(1, floor(e/2)).
Dirichlet g.f.: zeta(s) * zeta(4*s) * zeta(6*s) / zeta(12*s).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 15015/(1382*Pi^2) = 1.100823... .

Extensions

Name corrected by Amiram Eldar, Sep 08 2023
Showing 1-3 of 3 results.