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.

A385379 The maximum possible number of distinct composite prime powers (A246547) in the factorization of n into prime powers.

Original entry on oeis.org

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

Views

Author

Amiram Eldar, Jun 27 2025

Keywords

Comments

Differs from A376679 at n = 1, 48, 72, 80, ... .
The factorization includes primes if n is not a powerful number (A001694) that is larger than 1.
a(n) depends only on the prime signature of n (A118914).

Examples

			         n | a(n) | factorization
  ---------+------+----------------------------------------
         4 |  1   | 2^2
        32 |  2   | 2^2 * 2^3
       288 |  3   | 2^2 * 2^3 * 3^2
      4608 |  4   | 2^2 * 2^3 * 3^2 * 2^4
    115200 |  5   | 2^2 * 2^3 * 3^2 * 2^4 * 5^2
   3110400 |  6   | 2^2 * 2^3 * 3^2 * 2^4 * 5^2 * 3^3
  99532800 |  7   | 2^2 * 2^3 * 3^2 * 2^4 * 5^2 * 3^3 * 2^5
		

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := Floor[(Sqrt[8*e + 9] - 3)/2]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = vecsum(apply(x -> (sqrtint(8*x+9)-1)\2 , factor(n)[, 2]));

Formula

Additive with a(p^e) = A052146(e+1).
a(n) = 0 if and only if n is squarefree (A005117).
a(A385380(n)) = n-1.
Sum_{k=1..n} a(k) ~ n * (log(log(n)) + B + C), where B is Mertens's constant (A077761), C = Sum_{k>=1} P(k*(k+3)/2) = 0.49006911093767425812..., and P is the prime zeta function.

A376657 Number of integer factorizations of n into nonsquarefree factors > 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 07 2024

Keywords

Examples

			The a(n) factorizations for n = 16, 64, 72, 144, 192, 256, 288:
  (16)   (64)     (72)    (144)    (192)     (256)      (288)
  (4*4)  (8*8)    (8*9)   (4*36)   (4*48)    (4*64)     (4*72)
         (4*16)   (4*18)  (8*18)   (8*24)    (8*32)     (8*36)
         (4*4*4)          (9*16)   (12*16)   (16*16)    (9*32)
                          (12*12)  (4*4*12)  (4*8*8)    (12*24)
                          (4*4*9)            (4*4*16)   (16*18)
                                             (4*4*4*4)  (4*8*9)
                                                        (4*4*18)
		

Crossrefs

For prime-powers we have A000688.
Positions of zeros are A005117 (squarefree numbers), complement A013929.
For squarefree instead of nonsquarefree we have A050320, strict A050326.
For nonprime numbers we have A050370.
The version for partitions is A114374.
For perfect-powers we have A294068.
For non-perfect-powers we have A303707.
For non-prime-powers we have A322452.
The strict case is A376679.
Nonsquarefree numbers:
- A078147 (first differences)
- A376593 (second differences)
- A376594 (inflections and undulations)
- A376595 (nonzero curvature)
A000040 lists the prime numbers, differences A001223.
A001055 counts integer factorizations, strict A045778.
A005117 lists squarefree numbers, differences A076259.
A317829 counts factorizations of superprimorials, strict A337069.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],NoneTrue[SquareFreeQ]]],{n,100}]
Showing 1-2 of 2 results.