A325131 -id:A325131 - cp's OEIS frontend

A109297 Primal codes of finite permutations on positive integers.

1, 2, 9, 12, 18, 40, 112, 125, 250, 352, 360, 540, 600, 675, 832, 1008, 1125, 1350, 1500, 2176, 2250, 2268, 2352, 2401, 3168, 3969, 4802, 4864, 7488, 7938, 10692, 11616, 11776, 14000, 19584, 21609, 27440, 28812, 29403, 29696, 32448, 35000, 37908, 43218, 43776

Offset: 1

Jon Awbrey, Jul 08 2005

nonn

A finite permutation is a bijective mapping from a finite set to itself, counting the empty mapping as a permutation of the empty set.

Also Heinz numbers of integer partitions where the set of distinct parts is equal to the set of distinct multiplicities. These partitions are counted by A114640. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Apr 02 2019

			Writing (prime(i))^j as i:j, we have the following table:
Primal Codes of Finite Permutations on Positive Integers
` ` ` 1 = { }
` ` ` 2 = 1:1
` ` ` 9 = 2:2
` ` `12 = 1:2 2:1
` ` `18 = 1:1 2:2
` ` `40 = 1:3 3:1
` ` 112 = 1:4 4:1
` ` 125 = 3:3
` ` 250 = 1:1 3:3
` ` 352 = 1:5 5:1
` ` 360 = 1:3 2:2 3:1
` ` 540 = 1:2 2:3 3:1
` ` 600 = 1:3 2:1 3:2
` ` 675 = 2:3 3:2
` ` 832 = 1:6 6:1
` `1008 = 1:4 2:2 4:1
` `1125 = 2:2 3:3
` `1350 = 1:1 2:3 3:2
` `1500 = 1:2 2:1 3:3
` `2176 = 1:7 7:1
` `2250 = 1:1 2:2 3:3

Amiram Eldar, Table of n, a(n) for n = 1..300 (terms 1..100 from Alois P. Heinz)
Jon Awbrey, Riffs and Rotes.

Subsequence of A130091.
Cf. A076954, A106177, A108352, A108371, A109298, A109301, A056239, A112798, A114640, A118914.
Cf. A324524, A324525, A324571, A325127, A325128, A325130, A325131.

a:= proc(n) option remember; local k; for k from 1+`if`(n=1, 0,
      a(n-1)) while (l-> sort(map(i-> i[2], l)) <> sort(map(
      i-> numtheory[pi](i[1]), l)))(ifactors(k)[2]) do od; k
    end:
seq(a(n), n=1..45);  # Alois P. Heinz, Mar 08 2019

Select[Range[1000],#==1||Union[PrimePi/@First/@FactorInteger[#]]==Union[Last/@FactorInteger[#]]&] (* Gus Wiseman, Apr 02 2019 *)

is(n) = {my(f = factor(n), p = f[,1], e = vecsort(f[,2])); for(i=1, #p, if(primepi(p[i]) != e[i], return(0))); 1}; \\ Amiram Eldar, Jul 30 2022

More terms from Franklin T. Adams-Watters, Dec 19 2005

Offset set to 1 by Alois P. Heinz, Mar 08 2019

A114639 Number of partitions of n such that the set of parts and the set of multiplicities of parts are disjoint.

Original entry on oeis.org

1, 0, 2, 2, 2, 3, 5, 4, 7, 7, 13, 16, 19, 23, 33, 34, 44, 58, 63, 80, 101, 112, 139, 171, 196, 234, 288, 328, 394, 478, 545, 658, 777, 881, 1050, 1236, 1414, 1666, 1936, 2216, 2592, 3018, 3428, 3992, 4604, 5243, 6069, 6986, 7951, 9139, 10447, 11892, 13625

Offset: 0

Vladeta Jovovic, Feb 18 2006

nonn

The Heinz numbers of these partitions are given by A325131. - Gus Wiseman, Apr 02 2019

			From _Gus Wiseman_, Apr 02 2019: (Start)
The a(2) = 2 through a(9) = 7 partitions:
  (2)   (3)    (4)     (5)      (6)       (7)        (8)         (9)
  (11)  (111)  (1111)  (32)     (33)      (43)       (44)        (54)
                       (11111)  (42)      (52)       (53)        (63)
                                (222)     (1111111)  (62)        (72)
                                (111111)             (2222)      (432)
                                                     (3311)      (222111)
                                                     (11111111)  (111111111)
(End)

Alois P. Heinz, Table of n, a(n) for n = 0..100

Cf. A052335, A087153, A114640, A115584, A117144, A276429, A324572, A325130, A325131, A336032.

b:= proc(n, i, p, m) option remember; `if`(n=0, 1,
      `if`(i<1, 0, b(n, i-1, p, select(x-> x x<=n-i*j, p union {i}),
         select(x-> x b(n$2, {}$2):
seq(a(n), n=0..40);  # Alois P. Heinz, Aug 09 2016

b[n_, i_, p_, m_] := b[n, i, p, m] = If[n == 0, 1, If[i<1, 0, b[n, i-1, p, Select[m, #Jean-François Alcover, Feb 05 2017, after Alois P. Heinz *)
Table[Length[Select[IntegerPartitions[n],Intersection[#,Length/@Split[#]]=={}&]],{n,0,30}] (* Gus Wiseman, Apr 02 2019 *)

a(0)=1 prepended and more terms from Alois P. Heinz, Aug 09 2016

A276429 Number of partitions of n containing no part i of multiplicity i.

Original entry on oeis.org

1, 0, 2, 2, 3, 5, 8, 9, 16, 19, 29, 36, 53, 65, 92, 115, 154, 195, 257, 318, 419, 516, 663, 821, 1039, 1277, 1606, 1963, 2441, 2978, 3675, 4454, 5469, 6603, 8043, 9688, 11732, 14066, 16963, 20260, 24310, 28953, 34586, 41047, 48857, 57802, 68528, 80862, 95534, 112388, 132391

Offset: 0

Emeric Deutsch, Sep 19 2016

nonn

The Heinz numbers of these partitions are given by A325130. - Gus Wiseman, Apr 02 2019

			a(4) = 3 because we have [1,1,1,1], [1,1,2], and [4]; the partitions [1,3], [2,2] do not qualify.
From _Gus Wiseman_, Apr 02 2019: (Start)
The a(2) = 2 through a(7) = 9 partitions:
  (2)   (3)    (4)     (5)      (6)       (7)
  (11)  (111)  (211)   (32)     (33)      (43)
               (1111)  (311)    (42)      (52)
                       (2111)   (222)     (511)
                       (11111)  (411)     (3211)
                                (3111)    (4111)
                                (21111)   (31111)
                                (111111)  (211111)
                                          (1111111)
(End)

Vaclav Kotesovec, Table of n, a(n) for n = 0..12782 (terms 0..5000 from Alois P. Heinz)

Cf. A052335, A087153, A114639, A115584, A117144, A276427, A324572, A325130, A325131, A336269.

g := product(1/(1-x^i)-x^(i^2), i = 1 .. 100): gser := series(g, x = 0, 53): seq(coeff(gser, x, n), n = 0 .. 50);
# second Maple program:
b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
      add(`if`(i=j, 0, b(n-i*j, i-1)), j=0..n/i)))
    end:
a:= n-> b(n$2):
seq(a(n), n=0..60);  # Alois P. Heinz, Sep 19 2016

b[n_, i_] := b[n, i] = Expand[If[n == 0, 1, If[i < 1, 0, Sum[If[i == j, x, 1]*b[n - i*j, i - 1], {j, 0, n/i}]]]]; T[n_] := Function[p, Table[Coefficient[p, x, i], {i, 0, Exponent[p, x]}]][b[n, n]]; Table[T[n][[1]], {n, 0, 60}] (* Jean-François Alcover, Nov 28 2016 after Alois P. Heinz's Maple code for A276427 *)
Table[Length[Select[IntegerPartitions[n],And@@Table[Count[#,i]!=i,{i,Union[#]}]&]],{n,0,30}] (* Gus Wiseman, Apr 02 2019 *)

a(n) = A276427(n,0).

G.f.: g(x) = Product_{i>=1} (1/(1-x^i) - x^{i^2}).

A353393 Positive integers m > 1 that are prime or whose prime shadow A181819(m) is a divisor of m that is already in the sequence.

Original entry on oeis.org

2, 3, 5, 7, 9, 11, 13, 17, 19, 23, 29, 31, 36, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 125, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 225, 227, 229, 233, 239, 241, 251

Offset: 1

Gus Wiseman, May 15 2022

nonn

We define the prime shadow A181819(n) to be the product of primes indexed by the exponents in the prime factorization of n. For example, 90 = prime(1)*prime(2)^2*prime(3) has prime shadow prime(1)*prime(2)*prime(1) = 12.

			The terms together with their prime indices begin:
    2: {1}
    3: {2}
    5: {3}
    7: {4}
    9: {2,2}
   11: {5}
   13: {6}
   17: {7}
   19: {8}
   23: {9}
   29: {10}
   31: {11}
   36: {1,1,2,2}

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

The first term that is not a prime power A000961 is 36.
The first term that is not a prime or a perfect power A001597 is 1260. - Corrected by Robert Israel, Mar 10 2025
The non-recursive version is A325755, counted by A325702.
Removing all primes gives A353389.
These partitions are counted by A353426.
The version for compositions is A353431.
A001222 counts prime factors with multiplicity, distinct A001221.
A003963 gives product of prime indices.
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime signature, sorted A118914.
A130091 lists numbers with all distinct prime exponents, counted by A098859.
A181819 gives prime shadow, with an inverse A181821.
A325131 lists numbers relatively prime to their prime shadow.
Cf. A000005, A000040, A047966, A182850, A316413, A316428, A325756, A353394, A353395, A353399.

pshadow:= proc(n) local F,i;
  F:= ifactors(n)[2];
  mul(ithprime(i),i=F[..,2])
end proc:
filter:= proc(n) local s;
  if isprime(n) then return true fi;
  s:= pshadow(n);
  n mod s = 0 and member(s,R)
end proc:
R:= {}:
for i from 2 to 2000 do if filter(i) then R:= R union {i} fi od:
sort(convert(R,list)); # Robert Israel, Mar 10 2025

red[n_]:=If[n==1,1,Times@@Prime/@Last/@FactorInteger[n]];
suQ[n_]:=PrimeQ[n]||Divisible[n,red[n]]&&suQ[red[n]];
Select[Range[2,200],suQ[#]&]

Equals A353389 U A000040.

A353394 Product of prime shadows of prime indices of n (with multiplicity).

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, May 17 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.

We define the prime shadow A181819(n) to be the product of primes indexed by the exponents in the prime factorization of n. For example, 90 = prime(1)*prime(2)^2*prime(3) has prime shadow prime(1)*prime(2)*prime(1) = 12.

			We have 42 = prime(1)*prime(2)*prime(4), so a(42) = 1*2*3 = 6.

Positions of first appearances are A353397.
A001222 counts prime factors with multiplicity, distinct A001221.
A003963 gives product of prime indices.
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime signature, sorted A118914, product A005361.
A181819 gives prime shadow, with an inverse A181821.
A324850 lists numbers divisible by the product of their prime indices.
A325131 lists numbers relatively prime to their prime shadow.
A325755 lists numbers divisible by their prime shadow, quotient also A325756, with recursion A353393.
Cf. A003586, A005117, A008480, A182850, A266477, A289509, A316428, A320325, A325702, A353389, A353395, A353399.

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

a(n) = Product_i A181819(A112798(n,i)).
Positions where a(n) = A003963(n) are A003586.
Positions where a(n) = A005361(n) are A353399, counted by A353398.
Positions where a(n) = A181819(n) are A353395, counted by A353396.

A353399 Numbers whose product of prime exponents equals the product of prime shadows of its prime indices.

Original entry on oeis.org

1, 2, 12, 20, 36, 44, 56, 68, 100, 124, 164, 184, 208, 236, 240, 268, 332, 436, 464, 484, 508, 528, 608, 628, 688, 716, 720, 752, 764, 776, 816, 844, 880, 964, 1108, 1132, 1156, 1168, 1200, 1264, 1296, 1324, 1344, 1360, 1412, 1468, 1488, 1584, 1604, 1616, 1724

Offset: 1

Gus Wiseman, May 17 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.

We define the prime shadow A181819(n) to be the product of primes indexed by the exponents in the prime factorization of n. For example, 90 = prime(1)*prime(2)^2*prime(3) has prime shadow prime(1)*prime(2)*prime(1) = 12.

			The terms together with their prime indices begin:
     1: {}
     2: {1}
    12: {1,1,2}
    20: {1,1,3}
    36: {1,1,2,2}
    44: {1,1,5}
    56: {1,1,1,4}
    68: {1,1,7}
   100: {1,1,3,3}
   124: {1,1,11}
   164: {1,1,13}
   184: {1,1,1,9}
   208: {1,1,1,1,6}
   236: {1,1,17}
   240: {1,1,1,1,2,3}

Product of prime indices is A003963, counted by A339095.
The LHS (product of exponents) is A005361, counted by A266477.
The RHS (product of shadows) is A353394, first appearances A353397.
A related comparison is A353395, counted by A353396.
The partitions are counted by A353398.
Taking indices instead of exponents on the LHS gives A353503.
A001222 counts prime factors with multiplicity, distinct A001221.
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime signature, sorted A118914.
A130091 lists numbers with distinct prime exponents, counted by A098859.
A181819 gives prime shadow, with an inverse A181821.
A325131 lists numbers relatively prime to their prime shadow.
Numbers divisible by their prime shadow:
- counted by A325702
- listed by A325755
- co-recursive version A325756
- nonprime recursive version A353389
- recursive version A353393
- recursive version counted by A353426
Cf. A000720, A003586, A005117, A143773, A182850, A316428, A320325, A324850.

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

A005361(a(n)) = A353394(a(n)).

A353503 Numbers whose product of prime indices equals their product of prime exponents (prime signature).

Original entry on oeis.org

1, 2, 12, 36, 40, 112, 352, 832, 960, 1296, 2176, 2880, 4864, 5376, 11776, 12544, 16128, 29696, 33792, 34560, 38400, 63488, 64000, 101376, 115200, 143360, 151552, 159744, 335872, 479232, 704512, 835584, 1540096, 1658880, 1802240

Offset: 1

Gus Wiseman, May 17 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. A number's prime signature (row n A124010) is the sequence of positive exponents in its prime factorization.

			The terms together with their prime indices begin:
     1: {}
     2: {1}
    12: {1,1,2}
    36: {1,1,2,2}
    40: {1,1,1,3}
   112: {1,1,1,1,4}
   352: {1,1,1,1,1,5}
   832: {1,1,1,1,1,1,6}
   960: {1,1,1,1,1,1,2,3}
  1296: {1,1,1,1,2,2,2,2}
  2176: {1,1,1,1,1,1,1,7}
  2880: {1,1,1,1,1,1,2,2,3}
  4864: {1,1,1,1,1,1,1,1,8}
  5376: {1,1,1,1,1,1,1,1,2,4}

For shadows instead of exponents we get A003586, counted by A008619.
The LHS (product of prime indices) is A003963, counted by A339095.
The RHS (product of prime exponents) is A005361, counted by A266477.
The version for shadows instead of indices is A353399, counted by A353398.
These partitions are counted by A353506.
A001222 counts prime factors with multiplicity, distinct A001221.
A056239 adds up prime indices, row sums of A112798 and A296150.
A130091 lists numbers with distinct prime exponents, counted by A098859.
A124010 gives prime signature, sorted A118914.
A181819 gives prime shadow, with an inverse A181821.
A353394 gives product of shadows of prime indices, firsts A353397.
Cf. A000720, A008480, A085629, A097318, A109297, A304678, A318871, A320325, A325131, A325755, A353500, A353507.

Select[Range[1000],Times@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>PrimePi[p]^k]==Times@@Last/@FactorInteger[#]&]

from itertools import count, islice
from math import prod
from sympy import primepi, factorint
def A353503_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda n: n == 1 or prod((f:=factorint(n)).values()) == prod(primepi(p)**e for p,e in f.items()), count(max(startvalue,1)))
A353503_list = list(islice(A353503_gen(),20)) # Chai Wah Wu, May 20 2022

A003963(a(n)) = A005361(a(n)).

A054744 p-full numbers: numbers such that if any prime p divides it, then so does p^p.

Original entry on oeis.org

1, 4, 8, 16, 27, 32, 64, 81, 108, 128, 216, 243, 256, 324, 432, 512, 648, 729, 864, 972, 1024, 1296, 1728, 1944, 2048, 2187, 2592, 2916, 3125, 3456, 3888, 4096, 5184, 5832, 6561, 6912, 7776, 8192, 8748, 10368, 11664, 12500, 13824, 15552, 15625, 16384

Offset: 1

James Sellers, Apr 22 2000

A027748(a(n),k) <= A124010(a(n),k), 1<=k<=A001221(a(n)). [Reinhard Zumkeller, Apr 28 2012]

Heinz numbers of integer partitions where the multiplicity of each part k is at least prime(k). These partitions are counted by A325132. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Apr 02 2019

			8 is an element because 8 = 2^3 and 2<=3, while 25 is not an element because 25 = 5^2 and 5>2.
From _Gus Wiseman_, Apr 02 2019: (Start)
The sequence of terms together with their prime indices begins:
    1: {}
    4: {1,1}
    8: {1,1,1}
   16: {1,1,1,1}
   27: {2,2,2}
   32: {1,1,1,1,1}
   64: {1,1,1,1,1,1}
   81: {2,2,2,2}
  108: {1,1,2,2,2}
  128: {1,1,1,1,1,1,1}
  216: {1,1,1,2,2,2}
  243: {2,2,2,2,2}
  256: {1,1,1,1,1,1,1,1}
  324: {1,1,2,2,2,2}
  432: {1,1,1,1,2,2,2}
  512: {1,1,1,1,1,1,1,1,1}
  648: {1,1,1,2,2,2,2}
  729: {2,2,2,2,2,2}
  864: {1,1,1,1,1,2,2,2}
  972: {1,1,2,2,2,2,2}
(End)

Donovan Johnson, Table of n, a(n) for n = 1..10000

Cf. A048103, A001694, A036966.
Cf. A100716.
Cf. A056239, A109298, A112798, A115584, A117144, A124010, A276078, A324525, A324571, A325127, A325128, A325130, A325131, A325132.

a054744 n = a054744_list !! (n-1)
a054744_list = filter (\x -> and $
   zipWith (<=) (a027748_row x) (map toInteger $ a124010_row x)) [1..]
-- Reinhard Zumkeller, Apr 28 2012

Select[Range[1000],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>k>=p]&] (* Gus Wiseman, Apr 02 2019 *)

If n = Product p_i^e_i then p_i<=e_i for all i.

Sum_{n>=1} 1/a(n) = Product_{p prime} 1 + 1/(p^(p-1)*(p-1)) = 1.58396891058853238595.... - Amiram Eldar, Oct 24 2020

A353397 Replace prime(k) with prime(2^k) in the prime factorization of n.

Original entry on oeis.org

1, 3, 7, 9, 19, 21, 53, 27, 49, 57, 131, 63, 311, 159, 133, 81, 719, 147, 1619, 171, 371, 393, 3671, 189, 361, 933, 343, 477, 8161, 399, 17863, 243, 917, 2157, 1007, 441, 38873, 4857, 2177, 513, 84017, 1113, 180503, 1179, 931, 11013, 386093, 567, 2809, 1083

Offset: 1

Gus Wiseman, May 17 2022

			The terms together with their prime indices begin:
      1: {}
      3: {2}
      7: {4}
      9: {2,2}
     19: {8}
     21: {2,4}
     53: {16}
     27: {2,2,2}
     49: {4,4}
     57: {2,8}
    131: {32}
     63: {2,2,4}

Amiram Eldar, Table of n, a(n) for n = 1..397 (calculated using the b-file at A033844)

These are the positions of first appearances in A353394.
A001222 counts prime factors with multiplicity, distinct A001221.
A003963 gives product of prime indices, counted by A339095.
A033844 lists primes indexed by powers of 2.
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime signature, sorted A118914, product A005361.
A181819 gives prime shadow, firsts A181821, relatively prime A325131.
Equivalent sequence with prime(2*k) instead of prime(2^k): A297002.
Cf. A003586, A005117, A130091, A182850, A289509, A324850, A325755, A353393, A353395, A353399.

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

a(n) = my(f=factor(n)); for(k=1, #f~, f[k,1] = prime(2^primepi(f[k,1]))); factorback(f); \\ Michel Marcus, May 20 2022

from math import prod
from sympy import prime, primepi, factorint
def A353397(n): return prod(prime(2**primepi(p))**e for p, e in factorint(n).items()) # Chai Wah Wu, May 20 2022

If n = prime(e_1)...prime(e_k), then a(n) = prime(2^(e_1))...prime(2^(e_k)).

Sum_{n>=1} 1/a(n) = 1/Product_{k>=1} (1 - 1/prime(2^k)) = 1.90812936178871496289... . - Amiram Eldar, Dec 09 2022

A325130 Numbers in whose prime factorization the exponent of prime(k) is not equal to k for any prime index k.

Original entry on oeis.org

1, 3, 4, 5, 7, 8, 11, 12, 13, 15, 16, 17, 19, 20, 21, 23, 24, 25, 27, 28, 29, 31, 32, 33, 35, 37, 39, 40, 41, 43, 44, 47, 48, 49, 51, 52, 53, 55, 56, 57, 59, 60, 61, 64, 65, 67, 68, 69, 71, 73, 75, 76, 77, 79, 80, 81, 83, 84, 85, 87, 88, 89, 91, 92, 93, 95, 96

Offset: 1

Gus Wiseman, Apr 01 2019

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 Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are Heinz numbers of the integer partitions counted by A276429.
The asymptotic density of this sequence is Product_{k>=1} (1 - 1/prime(k)^k + 1/prime(k)^(k+1)) = 0.68974964705635552968... - Amiram Eldar, Jan 09 2021

			The sequence of terms together with their prime indices begins:
   1: {}
   3: {2}
   4: {1,1}
   5: {3}
   7: {4}
   8: {1,1,1}
  11: {5}
  12: {1,1,2}
  13: {6}
  15: {2,3}
  16: {1,1,1,1}
  17: {7}
  19: {8}
  20: {1,1,3}
  21: {2,4}
  23: {9}
  24: {1,1,1,2}
  25: {3,3}
  27: {2,2,2}
  28: {1,1,4}

Alois P. Heinz, Table of n, a(n) for n = 1..10000

Cf. A056239, A087153, A112798, A124010, A276078, A276429.
Cf. A324525, A324571, A325127, A325128, A325130, A325131.
Complement of A276936.

q:= n-> andmap(i-> numtheory[pi](i[1])<>i[2], ifactors(n)[2]):
a:= proc(n) option remember; local k; for k from 1+
     `if`(n=1, 0, a(n-1)) while not q(k) do od; k
    end:
seq(a(n), n=1..80);  # Alois P. Heinz, Oct 28 2019

Select[Range[100],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>k!=PrimePi[p]]&]

cp's OEIS Frontend

A109297 Primal codes of finite permutations on positive integers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Extensions

A114639 Number of partitions of n such that the set of parts and the set of multiplicities of parts are disjoint.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Extensions

A276429 Number of partitions of n containing no part i of multiplicity i.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A353393 Positive integers m > 1 that are prime or whose prime shadow A181819(m) is a divisor of m that is already in the sequence.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A353394 Product of prime shadows of prime indices of n (with multiplicity).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A353399 Numbers whose product of prime exponents equals the product of prime shadows of its prime indices.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A353503 Numbers whose product of prime indices equals their product of prime exponents (prime signature).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Python