A356065 -id:A356065 - cp's OEIS frontend

A050361 Number of factorizations into distinct prime powers greater than 1.

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, 2, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 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, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1

Offset: 1

Christian G. Bower, Oct 15 1999

a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24 = 2^3*3 and 375 = 3*5^3 both have prime signature (3,1).

The number of unordered factorizations of n into 1 and exponentially odd prime powers, i.e., p^e where p is a prime and e is odd (A246551). - Amiram Eldar, Jun 12 2025

			From _Gus Wiseman_, Jul 30 2022: (Start)
The A000688(216) = 9 factorizations of 216 into prime powers are:
  (2*2*2*3*3*3)
  (2*2*2*3*9)
  (2*2*2*27)
  (2*3*3*3*4)
  (2*3*4*9)
  (2*4*27)
  (3*3*3*8)
  (3*8*9)
  (8*27)
Of these, the a(216) = 4 strict cases are:
  (2*3*4*9)
  (2*4*27)
  (3*8*9)
  (8*27)
(End)

Antti Karttunen, Table of n, a(n) for n = 1..100000 (first 10000 terms from Reinhard Zumkeller)
Index entries for sequences computed from exponents in factorization of n.

Cf. A000009, A050360, A050362, A050363, A050364, A101296, A246551.
Cf. A124010.
This is the strict case of A000688.
Positions of 1's are A004709, complement A046099.
The case of primes (instead of prime-powers) is A008966, non-strict A000012.
The non-strict additive version allowing 1's A023893, ranked by A302492.
The non-strict additive version is A023894, ranked by A355743.
The additive version (partitions) is A054685, ranked by A356065.
The additive version allowing 1's is A106244, ranked by A302496.
A001222 counts prime-power divisors.
A005117 lists all squarefree numbers.
A034699 gives maximal prime-power divisor.
A246655 lists all prime-powers (A000961 includes 1), towers A164336.
A296131 counts twice-factorizations of type PQR, non-strict A295935.
Cf. A001970, A002110, A025487, A055887, A063834, A076610, A085970, A279786, A302590, A302601, A354911, A355742.

a050361 = product . map a000009 . a124010_row
-- Reinhard Zumkeller, Aug 28 2014

A050361 := proc(n)
    local a,f;
    if n = 1 then
        1;
    else
        a := 1 ;
        for f in ifactors(n)[2] do
            a := a*A000009(op(2,f)) ;
        end do:
    end if;
end proc: # R. J. Mathar, May 25 2017

Table[Times @@ PartitionsQ[Last /@ FactorInteger[n]], {n, 99}] (* Arkadiusz Wesolowski, Feb 27 2017 *)

A000009(n,k=(n-!(n%2))) = if(!n,1,my(s=0); while(k >= 1, if(k<=n, s += A000009(n-k,k)); k -= 2); (s));
A050361(n) = factorback(apply(A000009,factor(n)[,2])); \\ Antti Karttunen, Nov 17 2019

Dirichlet g.f.: Product_{n is a prime power >1}(1 + 1/n^s).
Multiplicative with a(p^e) = A000009(e).
a(A002110(k))=1.
a(n) = A050362(A101296(n)). - R. J. Mathar, May 26 2017
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} f(1/p) = 1.26020571070524171076..., where f(x) = (1-x) * Product_{k>=1} (1 + x^k). - Amiram Eldar, Oct 03 2023

A023893 Number of partitions of n into prime power parts (1 included); number of nonisomorphic Abelian subgroups of symmetric group S_n.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 10, 14, 20, 27, 36, 48, 63, 82, 105, 134, 171, 215, 269, 335, 415, 511, 626, 764, 929, 1125, 1356, 1631, 1953, 2333, 2776, 3296, 3903, 4608, 5427, 6377, 7476, 8744, 10205, 11886, 13818, 16032, 18565, 21463, 24768, 28536

Offset: 0

Olivier Gérard

nonn

			From _Gus Wiseman_, Jul 28 2022: (Start)
The a(0) = 1 through a(6) = 10 partitions:
  ()  (1)  (2)   (3)    (4)     (5)      (33)
           (11)  (21)   (22)    (32)     (42)
                 (111)  (31)    (41)     (51)
                        (211)   (221)    (222)
                        (1111)  (311)    (321)
                                (2111)   (411)
                                (11111)  (2211)
                                         (3111)
                                         (21111)
                                         (111111)
(End)

Seiichi Manyama, Table of n, a(n) for n = 0..10000
Index entries for sequences related to groups

Cf. A009490, A023894 (first differences), A062297 (number of Abelian subgroups).
Cf. A018819, A062051, A131995.
The multiplicative version (factorizations) is A000688.
Not allowing 1's gives A023894, strict A054685, ranked by A355743.
The version for just primes (not prime-powers) is A034891, strict A036497.
The strict version is A106244.
These partitions are ranked by A302492.
A000041 counts partitions, strict A000009.
A001222 counts prime-power divisors.
A072233 counts partitions by sum and length.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
Cf. A001970, A055887, A063834, A085970, A279784, A295935, A356065.

Table[Length[Select[IntegerPartitions[n],Count[Map[Length,FactorInteger[#]], 1] == Length[#] &]], {n, 0, 35}] (* Geoffrey Critzer, Oct 25 2015 *)
nmax = 50; Clear[P]; P[m_] := P[m] = Product[Product[1/(1-x^(p^k)), {k, 1, m}], {p, Prime[Range[PrimePi[nmax]]]}]/(1-x)+O[x]^nmax // CoefficientList[ #, x]&; P[1]; P[m=2]; While[P[m] != P[m-1], m++]; P[m] (* Jean-François Alcover, Aug 31 2016 *)

lista(m) = {x = t + t*O(t^m); gf = prod(k=1, m, if (isprimepower(k), 1/(1-x^k), 1))/(1-x); for (n=0, m, print1(polcoeff(gf, n, t), ", "));} \\ Michel Marcus, Mar 09 2013

from functools import lru_cache
from sympy import factorint
@lru_cache(maxsize=None)
def A023893(n):
    @lru_cache(maxsize=None)
    def c(n): return sum((p**(e+1)-p)//(p-1) for p,e in factorint(n).items())+1
    return (c(n)+sum(c(k)*A023893(n-k) for k in range(1,n)))//n if n else 1 # Chai Wah Wu, Jul 15 2024

G.f.: (Product_{p prime} Product_{k>=1} 1/(1-x^(p^k))) / (1-x).

A023894 Number of partitions of n into prime power parts (1 excluded).

Original entry on oeis.org

1, 0, 1, 1, 2, 2, 3, 4, 6, 7, 9, 12, 15, 19, 23, 29, 37, 44, 54, 66, 80, 96, 115, 138, 165, 196, 231, 275, 322, 380, 443, 520, 607, 705, 819, 950, 1099, 1268, 1461, 1681, 1932, 2214, 2533, 2898, 3305, 3768, 4285, 4872, 5530, 6267, 7094, 8022, 9060

Offset: 0

Olivier Gérard

nonn

			From _Gus Wiseman_, Jul 28 2022: (Start)
The a(0) = 1 through a(9) = 7 partitions:
  ()  .  (2)  (3)  (4)   (5)   (33)   (7)    (8)     (9)
                   (22)  (32)  (42)   (43)   (44)    (54)
                               (222)  (52)   (53)    (72)
                                      (322)  (332)   (333)
                                             (422)   (432)
                                             (2222)  (522)
                                                     (3222)
(End)

Seiichi Manyama, Table of n, a(n) for n = 0..10000
E. Grosswald, Partitions into prime powers

The multiplicative version (factorizations) is A000688, coprime A354911.
Allowing 1's gives A023893, strict A106244, ranked by A302492.
The strict version is A054685.
The version for just primes is ranked by A076610, squarefree A356065.
Twice-partitions of this type are counted by A279784, factorizations A295935.
These partitions are ranked by A355743.
A000041 counts partitions, strict A000009.
A001222 counts prime-power divisors.
A072233 counts partitions by sum and length.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
Cf. A001970, A055887, A063834, A085970.

Table[Length[Select[IntegerPartitions[n],And@@PrimePowerQ/@#&]],{n,0,30}] (* Gus Wiseman, Jul 28 2022 *)

is_primepower(n)= {ispower(n, , &n); isprime(n)}
lista(m) = {x = t + t*O(t^m); gf = prod(k=1, m, if (is_primepower(k), 1/(1-x^k), 1)); for (n=0, m, print1(polcoeff(gf, n, t), ", "));}
\\ Michel Marcus, Mar 09 2013

from functools import lru_cache
from sympy import factorint
@lru_cache(maxsize=None)
def A023894(n):
    @lru_cache(maxsize=None)
    def c(n): return sum((p**(e+1)-p)//(p-1) for p,e in factorint(n).items())
    return (c(n)+sum(c(k)*A023894(n-k) for k in range(1,n)))//n if n else 1 # Chai Wah Wu, Jul 15 2024

G.f.: Prod(p prime, Prod(k >= 1, 1/(1-x^(p^k))))

A355743 Numbers whose prime indices are all prime-powers.

Original entry on oeis.org

1, 3, 5, 7, 9, 11, 15, 17, 19, 21, 23, 25, 27, 31, 33, 35, 41, 45, 49, 51, 53, 55, 57, 59, 63, 67, 69, 75, 77, 81, 83, 85, 93, 95, 97, 99, 103, 105, 109, 115, 119, 121, 123, 125, 127, 131, 133, 135, 147, 153, 155, 157, 159, 161, 165, 171, 175, 177, 179, 187

Offset: 1

Gus Wiseman, Jul 24 2022

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Also MM-numbers of multiset partitions into constant multisets, where the multiset of multisets with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset of multisets with MM-number 78 is {{},{1},{1,2}}.

			The terms together with their prime indices begin:
   1: {}
   3: {2}
   5: {3}
   7: {4}
   9: {2,2}
  11: {5}
  15: {2,3}
  17: {7}
  19: {8}
  21: {2,4}
  23: {9}
  25: {3,3}
  27: {2,2,2}
  31: {11}
  33: {2,5}
  35: {3,4}
  41: {13}
  45: {2,2,3}

Robert Price, Table of n, a(n) for n = 1..1410

The multiplicative version is A000688, strict A050361, coprime A354911.
The case of only primes (not all prime-powers) is A076610, strict A302590.
Allowing prime index 1 gives A302492.
These are the products of elements of A302493.
Requiring n to be a prime-power gives A302601.
These are the positions of 1's in A355741.
The squarefree case is A356065.
The complement is A356066.
A001222 counts prime-power divisors.
A023894 counts ptns into prime-powers, strict A054685, with 1's A023893.
A034699 gives maximal prime-power divisor.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
A355742 chooses a prime-power divisor of each prime index.
Cf. A085970, A106244, A279784, A295935, A355731, A356064.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],And@@PrimePowerQ/@primeMS[#]&]

A085970 Number of integers ranging from 2 to n that are not prime-powers.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 3, 3, 4, 5, 5, 5, 6, 6, 7, 8, 9, 9, 10, 10, 11, 11, 12, 12, 13, 13, 13, 14, 15, 16, 17, 17, 18, 19, 20, 20, 21, 21, 22, 23, 24, 24, 25, 25, 26, 27, 28, 28, 29, 30, 31, 32, 33, 33, 34, 34, 35, 36, 36, 37, 38, 38, 39, 40, 41, 41, 42, 42, 43

Offset: 1

Reinhard Zumkeller, Jul 06 2003

nonn

For n > 2, a(n) gives the number of duplicate eliminations performed by the Sieve of Eratosthenes when sieving the interval [2, n]. - Felix Fröhlich, Dec 10 2016
Number of terms of A024619 <= n. - Felix Fröhlich, Dec 10 2016
First differs from A082997 at n = 30. - Gus Wiseman, Jul 28 2022

			The a(30) = 13 numbers: 6, 10, 12, 14, 15, 18, 20, 21, 22, 24, 26, 28, 30. - _Gus Wiseman_, Jul 28 2022

Felix Fröhlich, Table of n, a(n) for n = 1..10000

Cf. A000720, A024619, A085972.
The complement is counted by A065515, without 1's A025528.
For primes instead of prime-powers we have A065855, with 1's A062298.
Partial sums of A143731.
The version not treating 1 as a prime-power is A356068.
A000688 counts factorizations into prime-powers.
A001222 counts prime-power divisors.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
Cf. A023893, A023894, A036234, A069513, A207375, A354911, A355743, A356064, A356065.

With[{nn = 75}, Table[n - Count[#, k_ /; k < n] - 1, {n, nn}] &@ Join[{1}, Select[Range@ nn, PrimePowerQ]]] (* Michael De Vlieger, Dec 11 2016 *)

a(n) = my(i=0); forcomposite(c=4, n, if(!isprimepower(c), i++)); i \\ Felix Fröhlich, Dec 10 2016

from sympy import primepi, integer_nthroot
def A085970(n): return n-1-sum(primepi(integer_nthroot(n,k)[0]) for k in range(1,n.bit_length())) # Chai Wah Wu, Aug 20 2024

a(n) = Max{A024619(k)<=n} k;

a(n) = n - A065515(n) = A085972(n) - A000720(n).

Name modified by Gus Wiseman, Jul 28 2022. Normally 1 is not considered a prime-power, cf. A000961, A246655.

A356068 Number of integers ranging from 1 to n that are not prime-powers (1 is not a prime-power).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 4, 4, 5, 6, 6, 6, 7, 7, 8, 9, 10, 10, 11, 11, 12, 12, 13, 13, 14, 14, 14, 15, 16, 17, 18, 18, 19, 20, 21, 21, 22, 22, 23, 24, 25, 25, 26, 26, 27, 28, 29, 29, 30, 31, 32, 33, 34, 34, 35, 35, 36, 37, 37, 38, 39, 39, 40, 41, 42

Offset: 1

Gus Wiseman, Jul 31 2022

nonn

			The a(30) = 14 numbers: 1, 6, 10, 12, 14, 15, 18, 20, 21, 22, 24, 26, 28, 30.

The complement is counted by A025528, with 1's A065515.
For primes instead of prime-powers we have A062298, with 1's A065855.
The version treating 1 as a prime-power is A085970.
One more than the partial sums of A143731.
A000688 counts factorizations into prime-powers.
A001222 counts prime-power divisors.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
Cf. A000720, A023894, A036234, A050361, A054685, A069513, A355743, A356064, A356065, A356066.

Table[Length[Select[Range[n],!PrimePowerQ[#]&]],{n,100}]

a(n) = A085970(n) + 1.

A354911 Number of factorizations of n into relatively prime prime-powers.

Original entry on oeis.org

1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 2, 0, 1, 1, 0, 0, 2, 0, 2, 1, 1, 0, 3, 0, 1, 0, 2, 0, 1, 0, 0, 1, 1, 1, 4, 0, 1, 1, 3, 0, 1, 0, 2, 2, 1, 0, 5, 0, 2, 1, 2, 0, 3, 1, 3, 1, 1, 0, 2, 0, 1, 2, 0, 1, 1, 0, 2, 1, 1, 0, 6, 0, 1, 2, 2, 1, 1, 0, 5, 0, 1, 0, 2, 1, 1, 1

Offset: 1

Gus Wiseman, Jul 25 2022

nonn

			The a(n) factorizations for n = 6, 12, 24, 36, 48, 72, 96:
  2*3  3*4    3*8      4*9      3*16       8*9        3*32
       2*2*3  2*3*4    2*2*9    2*3*8      2*4*9      3*4*8
              2*2*2*3  3*3*4    3*4*4      3*3*8      2*3*16
                       2*2*3*3  2*2*3*4    2*2*2*9    2*2*3*8
                                2*2*2*2*3  2*3*3*4    2*3*4*4
                                           2*2*2*3*3  2*2*2*3*4
                                                      2*2*2*2*2*3

Wikipedia, Coprime integers.

This is the relatively prime case of A000688, partitions A023894.
Positions of 0's are A246655 (A000961 includes 1).
For strict instead of relatively prime we have A050361, partitions A054685.
Positions of 1's are A000469 (A120944 excludes 1).
For pairwise coprime instead of relatively prime we have A143731.
The version for partitions instead of factorizations is A356067.
A000005 counts divisors.
A001055 counts factorizations.
A001221 counts distinct prime divisors, with sum A001414.
A001222 counts prime-power divisors.
A289509 lists numbers whose prime indices are relatively prime.
A295935 counts twice-factorizations with constant blocks (type PPR).
A355743 lists numbers with prime-power prime indices, squarefree A356065.
Cf. A000837, A023893, A076610, A085970, A106244, A279784, A318721, A355737, A355742, A356064, A356066.

ufacs[s_,n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&, Select[ufacs[Select[s,Divisible[n/d,#]&],n/d],Min@@#>=d&]],{d,Select[s,Divisible[n,#]&]}]];
Table[Length[Select[ufacs[Select[Divisors[n],PrimePowerQ[#]&],n],GCD@@#<=1&]],{n,100}]

a(n) = A000688(n) if n is nonprime, otherwise a(n) = 0.

A356064 Numbers with a prime index other than 1 that is not a prime-power. Complement of A302492.

Original entry on oeis.org

13, 26, 29, 37, 39, 43, 47, 52, 58, 61, 65, 71, 73, 74, 78, 79, 86, 87, 89, 91, 94, 101, 104, 107, 111, 113, 116, 117, 122, 129, 130, 137, 139, 141, 142, 143, 145, 146, 148, 149, 151, 156, 158, 163, 167, 169, 172, 173, 174, 178, 181, 182, 183, 185, 188, 193

Offset: 1

Gus Wiseman, Jul 25 2022

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

These are numbers divisible by a prime number not of the form prime(q^k) where q is a prime number and k >= 1.

			The terms together with their prime indices begin:
   13: {6}
   26: {1,6}
   29: {10}
   37: {12}
   39: {2,6}
   43: {14}
   47: {15}
   52: {1,1,6}
   58: {1,10}
   61: {18}
   65: {3,6}
   71: {20}
   73: {21}
   74: {1,12}
   78: {1,2,6}
   79: {22}
   86: {1,14}
   87: {2,10}

Heinz numbers of the partitions counted by A023893.
Allowing prime index 1 gives A356066.
A000688 counts factorizations into prime-powers, strict A050361.
A001222 counts prime-power divisors.
A023894 counts partitions into prime-powers, strict A054685.
A034699 gives the maximal prime-power divisor.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
A355742 chooses a prime-power divisor of each prime index.
A355743 = numbers whose prime indices are prime-powers, squarefree A356065.
Cf. A076610, A085970, A106244, A302492, A302493, A302601, A330946, A354911.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],!And@@PrimePowerQ/@DeleteCases[primeMS[#],1]&]

A356066 Numbers with a prime index that is not a prime-power. Complement of A355743.

Original entry on oeis.org

2, 4, 6, 8, 10, 12, 13, 14, 16, 18, 20, 22, 24, 26, 28, 29, 30, 32, 34, 36, 37, 38, 39, 40, 42, 43, 44, 46, 47, 48, 50, 52, 54, 56, 58, 60, 61, 62, 64, 65, 66, 68, 70, 71, 72, 73, 74, 76, 78, 79, 80, 82, 84, 86, 87, 88, 89, 90, 91, 92, 94, 96, 98, 100, 101

Offset: 1

Gus Wiseman, Jul 31 2022

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

			The terms together with their prime indices begin:
    2: {1}
    4: {1,1}
    6: {1,2}
    8: {1,1,1}
   10: {1,3}
   12: {1,1,2}
   13: {6}
   14: {1,4}
   16: {1,1,1,1}
   18: {1,2,2}
   20: {1,1,3}
   22: {1,5}
   24: {1,1,1,2}

The complement is A355743, counted by A023894.
The squarefree complement is A356065, counted by A054685.
Allowing prime index 1 gives A356064, complement A302492.
A000688 counts factorizations into prime-powers, strict A050361.
A001222 counts prime-power divisors.
A034699 gives the maximal prime-power divisor.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
A355742 chooses a prime-power divisor of each prime index.
Cf. A076610, A085970, A106244, A302493, A302601, A330946, A354911.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],!And@@PrimePowerQ/@primeMS[#]&]

Union of A299174 and A356064.

A382524 Number of ways to choose a different constant partition of each part of a constant partition of n.

Original entry on oeis.org

1, 1, 2, 2, 5, 2, 6, 2, 10, 3, 6, 2, 24, 2, 6, 4, 17, 2, 36, 2, 18, 4, 6, 2, 86, 3, 6, 10, 18, 2, 44, 2, 50, 4, 6, 4, 159, 2, 6, 4, 62, 2, 44, 2, 18, 30, 6, 2, 486, 3, 12, 4, 18, 2, 140, 4, 62, 4, 6, 2, 932, 2, 6, 30, 157, 4, 44, 2, 18, 4, 20, 2, 1500, 2, 6

Offset: 0

Gus Wiseman, Apr 03 2025

nonn

These are strict twice-partitions of weight n and type PRR.

			The a(1) = 1 through a(8) = 10 twice-partitions:
  (1)  (2)   (3)    (4)      (5)      (6)       (7)        (8)
       (11)  (111)  (22)     (11111)  (33)      (1111111)  (44)
                    (1111)            (222)                (2222)
                    (11)(2)           (111111)             (22)(4)
                    (2)(11)           (111)(3)             (4)(22)
                                      (3)(111)             (1111)(4)
                                                           (4)(1111)
                                                           (11111111)
                                                           (1111)(22)
                                                           (22)(1111)

Gus Wiseman, Sequences enumerating triangles of integer partitions

For distinct instead of equal block-sums we have A279786.
This is the strict case of A279789.
The orderless version is A304442, see A353833, A381995, A381871.
Multiset partitions of this type are ranked by A326534 /\ A355743 /\ A005117.
Partitions with no partition of this type are counted by A382076, strict case of A381993.
Normal multiset partitions of this type are counted by the strict case of A382204.
A006171 counts multiset partitions into constant blocks of integer partitions of n.
A050361 counts factorizations into distinct prime powers, see A381715.
A317141 counts coarsenings of prime indices, refinements A300383.
Cf. A000005, A000040, A000688, A018818, A047966, A063834, A260685, A279784, A356065, A381453, A381455.

Table[If[n==0,1,Sum[Binomial[Length[Divisors[n/d]],d]*d!,{d,Divisors[n]}]],{n,0,100}]

a(n) = Sum_{d|n} binomial(A000005(n/d),d) * d!

cp's OEIS Frontend

A050361 Number of factorizations into distinct prime powers greater than 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A023893 Number of partitions of n into prime power parts (1 included); number of nonisomorphic Abelian subgroups of symmetric group S_n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A023894 Number of partitions of n into prime power parts (1 excluded).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A355743 Numbers whose prime indices are all prime-powers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A085970 Number of integers ranging from 2 to n that are not prime-powers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

Extensions

A356068 Number of integers ranging from 1 to n that are not prime-powers (1 is not a prime-power).

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Formula

A354911 Number of factorizations of n into relatively prime prime-powers.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs