A076610 -id:A076610 - cp's OEIS frontend

A003961 Completely multiplicative with a(prime(k)) = prime(k+1).

1, 3, 5, 9, 7, 15, 11, 27, 25, 21, 13, 45, 17, 33, 35, 81, 19, 75, 23, 63, 55, 39, 29, 135, 49, 51, 125, 99, 31, 105, 37, 243, 65, 57, 77, 225, 41, 69, 85, 189, 43, 165, 47, 117, 175, 87, 53, 405, 121, 147, 95, 153, 59, 375, 91, 297, 115, 93, 61, 315, 67, 111, 275, 729, 119

Offset: 1

Marc LeBrun

Meyers (see Guy reference) conjectures that for all r >= 1, the least odd number not in the set {a(i): i < prime(r)} is prime(r+1). - N. J. A. Sloane, Jan 08 2021
Meyers' conjecture would be refuted if and only if for some r there were such a large gap between prime(r) and prime(r+1) that there existed a composite c for which prime(r) < c < a(c) < prime(r+1), in which case (by Bertrand's postulate) c would necessarily be a term of A246281. - Antti Karttunen, Mar 29 2021
a(n) is odd for all n and for each odd m there exists a k with a(k) = m (see A064216). a(n) > n for n > 1: bijection between the odd and all numbers. - Reinhard Zumkeller, Sep 26 2001
a(n) and n have the same number of distinct primes with (A001222) and without multiplicity (A001221). - Michel Marcus, Jun 13 2014
From Antti Karttunen, Nov 01 2019: (Start)
More generally, a(n) has the same prime signature as n, A046523(a(n)) = A046523(n). Also A246277(a(n)) = A246277(n) and A287170(a(n)) = A287170(n).
Many permutations and other sequences that employ prime factorization of n to encode either polynomials, partitions (via Heinz numbers) or multisets in general can be easily defined by using this sequence as one of their constituent functions. See the last line in the Crossrefs section for examples.
(End)

			a(12) = a(2^2 * 3) = a(prime(1)^2 * prime(2)) = prime(2)^2 * prime(3) = 3^2 * 5 = 45.
a(A002110(n)) = A002110(n + 1) / 2.

Richard K. Guy, editor, Problems From Western Number Theory Conferences, Labor Day, 1983, Problem 367 (Proposed by Leroy F. Meyers, The Ohio State U.).

Indranil Ghosh, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Index entries for sequences computed from indices in prime factorization.
Index entries for sequences related to Heinz numbers.

See A045965 for another version.
Row 1 of table A242378 (which gives the "k-th powers" of this sequence), row 3 of A297845 and of A306697. See also arrays A066117, A246278, A255483, A308503, A329050.
Cf. A064989 (a left inverse), A064216, A000040, A002110, A000265, A027746, A046523, A048673 (= (a(n)+1)/2), A108228 (= (a(n)-1)/2), A191002 (= a(n)*n), A252748 (= a(n)-2n), A286385 (= a(n)-sigma(n)), A283980 (= a(n)*A006519(n)), A341529 (= a(n)*sigma(n)), A326042, A049084, A001221, A001222, A122111, A225546, A260443, A245606, A244319, A246269 (= A065338(a(n))), A322361 (= gcd(n, a(n))), A305293.
Cf. A191555, A252738.
Cf. A249734, A249735 (bisections).
Cf. A246261 (a(n) is of the form 4k+1), A246263 (of the form 4k+3), A246271, A246272, A246259, A246281 (n such that a(n) < 2n), A246282 (n such that a(n) > 2n), A252742.
Cf. A275717 (a(n) > a(n-1)), A275718 (a(n) < a(n-1)).
Cf. A003972 (Möbius transform), A003973 (Inverse Möbius transform), A318321.
Cf. A300841, A305421, A322991, A250469, A269379 for analogous shift-operators in other factorization and quasi-factorization systems.
Cf. also following permutations and other sequences that can be defined with the help of this sequence: A005940, A163511, A122111, A260443, A206296, A265408, A265750, A275733, A275735, A297845, A091202 & A091203, A250245 & A250246, A302023 & A302024, A302025 & A302026.
A version for partition numbers is A003964, strict A357853.
A permutation of A005408.
Applying the same transformation again gives A357852.
Other multiplicative sequences: A064988, A357977, A357978, A357980, A357983.
A056239 adds up prime indices, row-sums of A112798.
Cf. A000720, A076610, A296150.

a003961 1 = 1
a003961 n = product $ map (a000040 . (+ 1) . a049084) $ a027746_row n
-- Reinhard Zumkeller, Apr 09 2012, Oct 09 2011
(MIT/GNU Scheme, with Aubrey Jaffer's SLIB Scheme library)
(require 'factor)
(define (A003961 n) (apply * (map A000040 (map 1+ (map A049084 (factor n))))))
;; Antti Karttunen, May 20 2014

a:= n-> mul(nextprime(i[1])^i[2], i=ifactors(n)[2]):
seq(a(n), n=1..80);  # Alois P. Heinz, Sep 13 2017

a[p_?PrimeQ] := a[p] = Prime[ PrimePi[p] + 1]; a[1] = 1; a[n_] := a[n] = Times @@ (a[#1]^#2& @@@ FactorInteger[n]); Table[a[n], {n, 1, 65}] (* Jean-François Alcover, Dec 01 2011, updated Sep 20 2019 *)
Table[Times @@ Map[#1^#2 & @@ # &, FactorInteger[n] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[n == 1], {n, 65}] (* Michael De Vlieger, Mar 24 2017 *)

a(n)=local(f); if(n<1,0,f=factor(n); prod(k=1,matsize(f)[1],nextprime(1+f[k,1])^f[k,2]))

a(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ Michel Marcus, May 17 2014

use ntheory ":all";  sub a003961 { vecprod(map { next_prime($) } factor(shift)); }  # _Dana Jacobsen, Mar 06 2016

from sympy import factorint, prime, primepi, prod
def a(n):
    f=factorint(n)
    return 1 if n==1 else prod(prime(primepi(i) + 1)**f[i] for i in f)
[a(n) for n in range(1, 11)] # Indranil Ghosh, May 13 2017

If n = Product p(k)^e(k) then a(n) = Product p(k+1)^e(k).
Multiplicative with a(p^e) = A000040(A000720(p)+1)^e. - David W. Wilson, Aug 01 2001
a(n) = Product_{k=1..A001221(n)} A000040(A049084(A027748(n,k))+1)^A124010(n,k). - Reinhard Zumkeller, Oct 09 2011 [Corrected by Peter Munn, Nov 11 2019]
A064989(a(n)) = n and a(A064989(n)) = A000265(n). - Antti Karttunen, May 20 2014 & Nov 01 2019
A001221(a(n)) = A001221(n) and A001222(a(n)) = A001222(n). - Michel Marcus, Jun 13 2014
From Peter Munn, Oct 31 2019: (Start)
a(n) = A225546((A225546(n))^2).
a(A225546(n)) = A225546(n^2).
(End)
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{p prime} ((p^2-p)/(p^2-nextprime(p))) = 2.06399637... . - Amiram Eldar, Nov 18 2022

A355739 Number of ways to choose a sequence of all different divisors, one of each prime index of n (with multiplicity).

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Jul 18 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 a(49) = 6 ways are: (1,2), (1,4), (2,1), (2,4), (4,1), (4,2).
The a(182) = 5 ways are: (1,2,3), (1,2,6), (1,4,2), (1,4,3), (1,4,6).
The a(546) = 2 ways are: (1,2,4,3), (1,2,4,6).

Wikipedia, Cartesian product.

This is the strict version of A355731, firsts A355732.
For relatively prime instead of strict we have A355737, firsts A355738.
Positions of 0's are A355740.
A000005 counts divisors.
A001221 counts distinct prime factors, with sum A001414.
A001222 counts prime factors with multiplicity.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798.
A120383 lists numbers divisible by all of their prime indices.
A289508 gives GCD of prime indices, positions of 1's A289509.
Cf. A000720, A076610, A302796, A355535, A355537, A355733, A355735, A355741, A355742, A355744, A355748.

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

A355741 Number of ways to choose a sequence of prime factors, one of each prime index of n.

Original entry on oeis.org

1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 2, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 2, 0, 1, 0, 1, 0, 1, 0, 2, 0, 2, 0, 1, 0, 2, 0, 1, 0, 2, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 2, 0, 1, 0, 2, 0, 1, 0, 1, 0, 2, 0, 2, 0, 1, 0, 1, 0, 2, 0, 1, 0, 1, 0, 1, 0, 2

Offset: 1

Gus Wiseman, Jul 18 2022

First differs from A355744 at a(169) = 4, A355744(169) = 3.

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 prime indices of 1131 are {2,6,10}, and the a(1131) = 4 choices are: {2,2,2}, {2,2,5}, {2,3,2}, {2,3,5}.

Wikipedia, Cartesian product.

Positions of 0's are A299174.
The version for all divisors is A355731, firsts A355732.
Choosing prime-power divisors gives A355742.
Positions of 1's are A355743.
Counting multisets instead of sequences gives A355744.
The weakly increasing case is A355745, all divisors A355735.
A001414 adds up distinct prime factors, counted by A001221.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A289509 lists numbers with relatively prime prime indices.
A324850 lists numbers divisible by the product of their prime indices.
Cf. A000720, A061395, A076610, A120383, A335433, A355733, A355737, A355739, A355746, A355747.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Times@@PrimeNu/@primeMS[n],{n,100}]

Totally multiplicative with a(prime(k)) = A001221(k).

A355529 Numbers of which it is not possible to choose a different prime factor of each prime index (with multiplicity).

Original entry on oeis.org

2, 4, 6, 8, 9, 10, 12, 14, 16, 18, 20, 21, 22, 24, 25, 26, 27, 28, 30, 32, 34, 36, 38, 40, 42, 44, 45, 46, 48, 49, 50, 52, 54, 56, 57, 58, 60, 62, 63, 64, 66, 68, 70, 72, 74, 75, 76, 78, 80, 81, 82, 84, 86, 88, 90, 92, 94, 96, 98, 99, 100, 102, 104, 105, 106

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.

Includes all even numbers.

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

The odd case is A355535.
The case of all divisors (not just primes) is A355740, zeros of A355739.
These choices are variously counted by A355741, A355744, A355745.
A001414 adds up distinct prime divisors, counted by A001221.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A120383 lists numbers divisible by all of their prime indices.
A324850 lists numbers divisible by the product of their prime indices.
A355731 counts choices of a divisor of each prime index, firsts A355732.
Cf. A000720, A076610, A318979, A335433, A335448, A355733.

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

A355740 Numbers of which it is not possible to choose a different divisor of each prime index.

Original entry on oeis.org

4, 8, 12, 16, 18, 20, 24, 27, 28, 32, 36, 40, 44, 48, 50, 52, 54, 56, 60, 64, 68, 72, 76, 80, 81, 84, 88, 90, 92, 96, 100, 104, 108, 112, 116, 120, 124, 125, 126, 128, 132, 135, 136, 140, 144, 148, 150, 152, 156, 160, 162, 164, 168, 172, 176, 180, 184, 188

Offset: 1

Gus Wiseman, Jul 22 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.

By Hall's marriage theorem, k is a term if and only if there is a sub-multiset S of the prime indices of k such that fewer than |S| numbers are divisors of a member of S. Equivalently, k is divisible by a member of A370348. - Robert Israel, Feb 15 2024

			The terms together with their prime indices begin:
    4: {1,1}
    8: {1,1,1}
   12: {1,1,2}
   16: {1,1,1,1}
   18: {1,2,2}
   20: {1,1,3}
   24: {1,1,1,2}
   27: {2,2,2}
   28: {1,1,4}
   32: {1,1,1,1,1}
   36: {1,1,2,2}
   40: {1,1,1,3}
   44: {1,1,5}
   48: {1,1,1,1,2}
For example, the choices of a divisor of each prime index of 90 are: (1,1,1,1), (1,1,1,3), (1,1,2,1), (1,1,2,3), (1,2,1,1), (1,2,1,3), (1,2,2,1), (1,2,2,3). But none of these has all distinct elements, so 90 is in the sequence.

Robert Israel, Table of n, a(n) for n = 1..10000
Wikipedia, Cartesian product.

Positions of 0's in A355739.
The case of just prime factors (not all divisors) is A355529, odd A355535.
The unordered case is counted by A355733, firsts A355734.
A000005 counts divisors.
A001414 adds up distinct prime divisors, counted by A001221.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A120383 lists numbers divisible by all of their prime indices.
A324850 lists numbers divisible by the product of their prime indices.
A355731 counts choices of a divisor of each prime index, firsts A355732.
A355741 chooses prime factors of prime indices, variations A355744, A355745.
Cf. A000720, A076610, A335433, A335448, A340827, A355737, A355749, A370348.

filter:= proc(n) uses numtheory, GraphTheory; local B, S, F, D, E, G, t, d;
  F:= ifactors(n)[2];
  F:= map(t -> [pi(t[1]), t[2]], F);
  D:= `union`(seq(divisors(t[1]), t = F));
  F:= map(proc(t) local i; seq([t[1], i], i=1..t[2]) end proc, F);
  if nops(D) < nops(F) then return false fi;
  E:= {seq(seq({t, d}, d=divisors(t[1])), t = F)};
  S:= map(t -> convert(t, name), [op(F), op(D)]);
  E:= map(e -> map(convert, e, name), E);
  G:= Graph(S, E);
  B:= BipartiteMatching(G);
  B[1] = nops(F);
end proc:
remove(filter, [$1..200]); # Robert Israel, Feb 15 2024

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

We have A001221(a(n)) >= A303975(a(n)).

A302569 Numbers that are either prime or whose prime indices are pairwise coprime. Heinz numbers of integer partitions with pairwise coprime parts.

Original entry on oeis.org

2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 16, 17, 19, 20, 22, 23, 24, 26, 28, 29, 30, 31, 32, 33, 34, 35, 37, 38, 40, 41, 43, 44, 46, 47, 48, 51, 52, 53, 55, 56, 58, 59, 60, 61, 62, 64, 66, 67, 68, 69, 70, 71, 73, 74, 76, 77, 79, 80, 82, 83, 85, 86, 88, 89

Offset: 1

Gus Wiseman, Apr 10 2018

nonn

A prime index of n is a number m such that prime(m) divides n.

The Heinz number of an integer partition (y_1,..,y_k) is prime(y_1)*..*prime(y_k).

			Entry A302242 describes a correspondence between positive integers and multiset multisystems. In this case it gives the following sequence of multiset systems.
02: {{}}
03: {{1}}
04: {{},{}}
05: {{2}}
06: {{},{1}}
07: {{1,1}}
08: {{},{},{}}
10: {{},{2}}
11: {{3}}
12: {{},{},{1}}
13: {{1,2}}
14: {{},{1,1}}
15: {{1},{2}}
16: {{},{},{},{}}
17: {{4}}
19: {{1,1,1}}
20: {{},{},{2}}
22: {{},{3}}
23: {{2,2}}
24: {{},{},{},{1}}
26: {{},{1,2}}
28: {{},{},{1,1}}
29: {{1,3}}
30: {{},{1},{2}}
31: {{5}}
32: {{},{},{},{},{}}

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000

Subsequence of A122132.

Cf. A000961, A001222, A005117, A007359, A007716, A051424, A056239, A076610, A101268, A275024, A302505, A302568.

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[200],Or[PrimeQ[#],CoprimeQ@@primeMS[#]]&]

is(n)=if(n<9, return(n>1)); n>>=valuation(n,2); if(n<9, return(1)); my(f=factor(n)); if(vecmax(f[,2])>1, return(0)); if(#f~==1, return(1)); my(v=apply(primepi, f[,1]),P=vecprod(v)); for(i=1,#v, if(gcd(v[i],P/v[i])>1, return(0))); 1 \\ Charles R Greathouse IV, Nov 11 2021

A355731 Number of ways to choose a sequence of divisors, one of each element of the multiset of prime indices of n (row n of A112798).

Original entry on oeis.org

1, 1, 2, 1, 2, 2, 3, 1, 4, 2, 2, 2, 4, 3, 4, 1, 2, 4, 4, 2, 6, 2, 3, 2, 4, 4, 8, 3, 4, 4, 2, 1, 4, 2, 6, 4, 6, 4, 8, 2, 2, 6, 4, 2, 8, 3, 4, 2, 9, 4, 4, 4, 5, 8, 4, 3, 8, 4, 2, 4, 6, 2, 12, 1, 8, 4, 2, 2, 6, 6, 6, 4, 4, 6, 8, 4, 6, 8, 4, 2, 16, 2, 2, 6, 4, 4

Offset: 1

Gus Wiseman, Jul 16 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 a(15) = 4 choices are: (1,1), (1,3), (2,1), (2,3).
The a(18) = 4 choices are: (1,1,1), (1,1,2), (1,2,1), (1,2,2).

Wikipedia, Cartesian product.

Positions of 1's are A000079.
Dominated by A003963 (cf. A049820), with equality at A003586.
Positions of first appearances are A355732.
Counting distinct sequences after sorting gives A355733, firsts A355734.
Requiring the result to be weakly increasing gives A355735, firsts A355736.
Requiring the result to be relatively prime gives A355737, firsts A355738.
Requiring the choices to be distinct gives A355739, zeros A355740.
For prime divisors A355741, prime-powers A355742, weakly increasing A355745.
Choosing divisors of each of 1..n and resorting gives A355747.
An ordered version (using standard order compositions) is A355748.
A000005 counts divisors.
A001414 adds up distinct prime divisors, counted by A001221.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A120383 lists numbers divisible by all of their prime indices.
A289509 lists numbers with relatively prime prime indices.
A324850 lists numbers divisible by the product of their prime indices.
A340852 lists numbers that can be factored into divisors of bigomega.
Cf. A000720, A061395, A076610, A316524, A326841, A335433, A335448, A340606, A344616, A345957, A355744, A355746.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Times@@Length/@Divisors/@primeMS[n],{n,100}]

a(n) = Product_{k=1..A001222(n)} A000005(A112798(n,k)).

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

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, 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

A320628 Products of primes of nonprime index.

Original entry on oeis.org

1, 2, 4, 7, 8, 13, 14, 16, 19, 23, 26, 28, 29, 32, 37, 38, 43, 46, 47, 49, 52, 53, 56, 58, 61, 64, 71, 73, 74, 76, 79, 86, 89, 91, 92, 94, 97, 98, 101, 103, 104, 106, 107, 112, 113, 116, 122, 128, 131, 133, 137, 139, 142, 146, 148, 149, 151, 152, 158, 161, 163

Offset: 1

Gus Wiseman, Oct 18 2018

nonn

The index of a prime number n is the number m such that n is the m-th prime.

The asymptotic density of this sequence is Product_{p in A006450} (1 - 1/p) = 1/(Sum_{n>=1} 1/A076610(n)) < 1/3. - Amiram Eldar, Feb 02 2021

			The sequence of terms begins:
   1 = 1
   2 = prime(1)
   4 = prime(1)^2
   7 = prime(4)
   8 = prime(1)^3
  13 = prime(6)
  14 = prime(1)*prime(4)
  16 = prime(1)^4
  19 = prime(8)
  23 = prime(9)
  26 = prime(1)*prime(6)
  28 = prime(1)^2*prime(4)
  29 = prime(10)
  32 = prime(1)^5
  37 = prime(12)
  38 = prime(1)*prime(8)
  43 = prime(14)
  46 = prime(1)*prime(9)
  47 = prime(15)
  49 = prime(4)^2
  52 = prime(1)^2*prime(6)
  53 = prime(16)
  56 = prime(1)^3*prime(4)
  58 = prime(1)*prime(10)
  61 = prime(18)
  64 = prime(1)^6
  71 = prime(20)
  73 = prime(21)
  74 = prime(1)*prime(12)
  76 = prime(1)^2*prime(8)
  79 = prime(22)
  86 = prime(1)*prime(14)
  89 = prime(24)
  91 = prime(4)*prime(6)
  92 = prime(1)^2*prime(9)
  94 = prime(1)*prime(15)
  97 = prime(25)
  98 = prime(1)*prime(4)^2

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A018252, A056239, A290103, A302242, A302478, A320533, A320629, A320630, A320631, A320633.
Complement of A331386.
Positions of zeros in A257994.
Primes of prime index are A006450.
Primes of nonprime index are A007821.
Products of primes of prime index are A076610.
Products of primes of nonprime index are this sequence.
The number of prime prime indices is given by A257994.
The number of nonprime prime indices is given by A330944.
Cf. A000040, A000720, A001222, A112798, A295665.

Select[Range[100],And@@Not/@PrimeQ/@PrimePi/@First/@FactorInteger[#]&]

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))))

cp's OEIS Frontend

A003961 Completely multiplicative with a(prime(k)) = prime(k+1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Perl

Python

Formula

A355739 Number of ways to choose a sequence of all different divisors, one of each prime index of n (with multiplicity).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A355741 Number of ways to choose a sequence of prime factors, one of each prime index of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A355529 Numbers of which it is not possible to choose a different prime factor of each prime index (with multiplicity).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A355740 Numbers of which it is not possible to choose a different divisor of each prime index.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A302569 Numbers that are either prime or whose prime indices are pairwise coprime. Heinz numbers of integer partitions with pairwise coprime parts.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

A355731 Number of ways to choose a sequence of divisors, one of each element of the multiset of prime indices of n (row n of A112798).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs