A325131 -id:A325131 - cp's OEIS frontend

A353500 Numbers that are the smallest number with product of prime exponents k for some k. Sorted positions of first appearances in A005361, unsorted version A085629.

1, 4, 8, 16, 32, 64, 128, 144, 216, 288, 432, 864, 1152, 1296, 1728, 2048, 2592, 3456, 5184, 7776, 8192, 10368, 13824, 15552, 18432, 20736, 31104, 41472, 55296, 62208, 73728, 86400, 108000, 129600, 131072, 165888, 194400, 216000, 221184, 259200, 279936, 324000

Offset: 1

Gus Wiseman, May 17 2022

nonn

All terms are highly powerful (A005934), but that sequence looks only at first appearances that reach a record, and is missing 1152, 2048, 8192, etc.

			The prime exponents of 86400 are (7,3,2), and this is the first case of product 42, so 86400 is in the sequence.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Wikipedia, Highly powerful number.

These are the positions of first appearances in A005361, counted by A266477.
This is the sorted version of A085629.
The version for shadows instead of exponents is A353397, firsts in A353394.
A001222 counts prime factors with multiplicity, distinct A001221.
A003963 gives product of prime indices, counted by A339095.
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime exponents, sorted A118914.
A130091 lists numbers with distinct prime exponents, counted by A098859.
A181819 gives prime shadow, with an inverse A181821.
Cf. A070175, A097318, A116608, A162247, A182850, A304678, A325131, A325238, A353399, A353503, A353506, A353507.
Subsequence of A181800.

nn=1000;
d=Table[Times@@Last/@FactorInteger[n],{n,nn}];
Select[Range[nn],!MemberQ[Take[d,#-1],d[[#]]]&]
lps[fct_] := Module[{nf = Length[fct]}, Times @@ (Prime[Range[nf]]^Reverse[fct])]; lps[{1}] = 1; q[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, (n == 1 || AllTrue[e, # > 1 &]) && n == Min[lps /@ f[Times @@ e]]]; Select[Cases[Import["https://oeis.org/A025487/b025487.txt", "Table"], {, }][[;; , 2]], q] (* Amiram Eldar, Sep 29 2024, using the function f by T. D. Noe at A162247 *)

A353389 Create the sequence of all positive integers > 1 that are prime or whose prime shadow (A181819) is a divisor that is already in the sequence. Then remove all the primes.

Original entry on oeis.org

9, 36, 125, 225, 441, 1089, 1260, 1521, 1980, 2340, 2401, 2601, 2772, 3060, 3249, 3276, 3420, 4140, 4284, 4761, 4788, 5148, 5220, 5580, 5796, 6660, 6732, 7308, 7380, 7524, 7569, 7740, 7812, 7956, 8460, 8649, 8892, 9108, 9324, 9540, 10332, 10620, 10764, 10836

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.

Said differently, these are nonprime numbers > 1 whose prime shadow is a divisor that is either a prime number or a number already in the sequence.

			The initial terms and their prime indices:
     9: {2,2}
    36: {1,1,2,2}
   125: {3,3,3}
   225: {2,2,3,3}
   441: {2,2,4,4}
  1089: {2,2,5,5}
  1260: {1,1,2,2,3,4}
  1521: {2,2,6,6}
  1980: {1,1,2,2,3,5}

The first term that is not a perfect power A001597 is 1260.
Without the recursion we have A325755 (a superset), counted by A325702.
Before removing the primes we had A353393.
These partitions are counted by A353426 minus one.
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.
A181819 gives prime shadow, with an inverse A181821.
A182850 and A323014 give frequency depth, counted by A225485 and A325280.
A325131 lists numbers relatively prime to their prime shadow.
Cf. A000005, A000040, A005117, A316413, A316428, A324850, A353394, A353395, A353397, A353399.

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,2000],suQ[#]&&!PrimeQ[#]&]

A353398 Number of integer partitions of n where the product of multiplicities equals the product of prime shadows of the parts.

Original entry on oeis.org

1, 1, 0, 0, 1, 1, 1, 2, 1, 2, 1, 2, 6, 5, 4, 4, 6, 6, 8, 8, 13, 16, 13, 16, 18, 16, 20, 21, 27, 30, 27, 33, 41, 44, 51, 48, 58, 61, 66, 66, 74, 83, 86, 99, 102, 111, 115, 126, 137, 147, 156

Offset: 0

Gus Wiseman, May 17 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 a(8) = 1 through a(14) = 4 partitions (A = 10, B = 11):
  3311  711     61111  521111   5511      B11       A1111
        321111         3221111  9111      721111    731111
                                531111    811111    33221111
                                3321111   5221111   422111111
                                22221111  43111111
                                42111111

The LHS (product of multiplicities) is A005361, counted by A266477.
The RHS (product of prime shadows) is A353394, first appearances A353397.
A related comparison is A353396, ranked by A353395.
These partitions are ranked by A353399.
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.
A181819 gives prime shadow, with an inverse A181821.
A325131 lists numbers relatively prime to their prime shadow.
A325755 lists numbers divisible by their prime shadow, counted by A325702.
A339095 counts partitions by product (or factorizations by sum).
Cf. A002033, A003963, A143773, A182850, A239455, A324850, A325756, A353393, A353426.

red[n_]:=If[n==1,1,Times@@Prime/@Last/@FactorInteger[n]];
Table[Length[Select[IntegerPartitions[n],Times@@red/@#==Times@@Length/@Split[#]&]],{n,0,30}]

A352142 Numbers whose prime factorization has all odd indices and all odd exponents.

Original entry on oeis.org

1, 2, 5, 8, 10, 11, 17, 22, 23, 31, 32, 34, 40, 41, 46, 47, 55, 59, 62, 67, 73, 82, 83, 85, 88, 94, 97, 103, 109, 110, 115, 118, 125, 127, 128, 134, 136, 137, 146, 149, 155, 157, 160, 166, 167, 170, 179, 184, 187, 191, 194, 197, 205, 206, 211, 218, 227, 230

Offset: 1

Gus Wiseman, Mar 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, sum A056239, length A001222.
A number's prime signature is the sequence of positive exponents in its prime factorization, which is row n of A124010, length A001221, sum A001222.
These are the Heinz numbers of integer partitions with all odd parts and all odd multiplicities, counted by A117958.

			The terms together with their prime indices begin:
   1 = 1
   2 = prime(1)
   5 = prime(3)
   8 = prime(1)^3
  10 = prime(1) prime(3)
  11 = prime(5)
  17 = prime(7)
  22 = prime(1) prime(5)
  23 = prime(9)
  31 = prime(11)
  32 = prime(1)^5
  34 = prime(1) prime(7)
  40 = prime(1)^3 prime(3)

David A. Corneth, Table of n, a(n) for n = 1..10000

The restriction to primes is A031368.
The first condition alone is A066208, counted by A000009.
These partitions are counted by A117958.
The squarefree case is A258116, even A258117.
The second condition alone is A268335, counted by A055922.
The even-even version is A352141 counted by A035444.
A000290 = exponents all even, counted by A035363.
A056166 = exponents all prime, counted by A055923.
A066207 = indices all even, counted by A035363 (complement A086543).
A109297 = same indices as exponents, counted by A114640.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.
A162641 counts even prime exponents, odd A162642.
A257991 counts odd prime indices, even A257992.
A325131 = disjoint indices from exponents, counted by A114639.
A346068 = indices and exponents all prime, counted by A351982.
A351979 = odd indices with even exponents, counted by A035457.
A352140 = even indices with odd exponents, counted by A055922 aerated.
A352143 = odd indices with odd conjugate indices, counted by A053253 aerated.
Cf. A000720, A028260, A055396, A061395, A106529, A181819, A195017, A241638, A276078, A324517, A324524, A324525, A325698, A325700.

Select[Range[100],#==1||And@@OddQ/@PrimePi/@First/@FactorInteger[#]&&And@@OddQ/@Last/@FactorInteger[#]&]

from itertools import count, islice
from sympy import primepi, factorint
def A352142_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda k:all(map(lambda x:x[1]%2 and primepi(x[0])%2, factorint(k).items())),count(max(startvalue,1)))
A352142_list = list(islice(A352142_gen(),30)) # Chai Wah Wu, Mar 18 2022

Intersection of A066208 and A268335.
A257991(a(n)) = A001222(a(n)).
A162642(a(n)) = A001221(a(n)).
A257992(a(n)) = A162641(a(n)) = 0.

A352141 Numbers whose prime factorization has all even indices and all even exponents.

Original entry on oeis.org

1, 9, 49, 81, 169, 361, 441, 729, 841, 1369, 1521, 1849, 2401, 2809, 3249, 3721, 3969, 5041, 6241, 6561, 7569, 7921, 8281, 10201, 11449, 12321, 12769, 13689, 16641, 17161, 17689, 19321, 21609, 22801, 25281, 26569, 28561, 29241, 29929, 32761, 33489, 35721

Offset: 1

Gus Wiseman, Mar 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, sum A056239, length A001222.
A number's prime signature is the sequence of positive exponents in its prime factorization, which is row n of A124010, length A001221, sum A001222.
These are the Heinz numbers of partitions with all even parts and all even multiplicities, counted by A035444.

			The terms together with their prime indices begin:
     1 = 1
     9 = prime(2)^2
    49 = prime(4)^2
    81 = prime(2)^4
   169 = prime(6)^2
   361 = prime(8)^2
   441 = prime(2)^2 prime(4)^2
   729 = prime(2)^6
   841 = prime(10)^2
  1369 = prime(12)^2
  1521 = prime(2)^2 prime(6)^2
  1849 = prime(14)^2
  2401 = prime(4)^4
  2809 = prime(16)^2
  3249 = prime(2)^2 prime(8)^2
  3721 = prime(18)^2
  3969 = prime(2)^4 prime(4)^2

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

The second condition alone (all even exponents) is A000290, counted by A035363.
The restriction to primes is A031215.
These partitions are counted by A035444.
The first condition alone is A066207, counted by A035363, squarefree A258117.
A056166 = exponents all prime, counted by A055923.
A066208 = prime indices all odd, counted by A000009.
A109297 = same indices as exponents, counted by A114640.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.
A162641 counts even exponents, odd A162642.
A257991 counts odd indices, even A257992.
A325131 = disjoint indices from exponents, counted by A114639.
A346068 = indices and exponents all prime, counted by A351982.
A351979 = odd indices with even exponents, counted by A035457.
A352140 = even indices with odd exponents, counted by A055922 aerated.
A352142 = odd indices with odd exponents, counted by A117958.
Cf. A000720, A028260, A055396, A061395, A181819, A195017, A241638, A268335, A276078, A324524, A324525, A324588, A325698, A325700, A352143.

Select[Range[1000],#==1||And@@EvenQ/@PrimePi/@First/@FactorInteger[#]&&And@@EvenQ/@Last/@FactorInteger[#]&]

from itertools import count, islice
from sympy import factorint, primepi
def A352141_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda k:all(map(lambda x: not (x[1]%2 or primepi(x[0])%2), factorint(k).items())),count(max(startvalue,1)))
A352141_list = list(islice(A352141_gen(),30)) # Chai Wah Wu, Mar 18 2022

Intersection of A000290 and A066207.
A257991(a(n)) = A162642(a(n)) = 0.
A257992(a(n)) = A001222(a(n)).
A162641(a(n)) = A001221(a(n)).
Sum_{n>=1} 1/a(n) = 1/Product_{k>=1} (1 - 1/prime(2*k)^2) = 1.163719... . - Amiram Eldar, Sep 19 2022

A352492 Powerful numbers whose prime indices are all prime numbers.

Original entry on oeis.org

1, 9, 25, 27, 81, 121, 125, 225, 243, 289, 625, 675, 729, 961, 1089, 1125, 1331, 1681, 2025, 2187, 2601, 3025, 3125, 3267, 3375, 3481, 4489, 4913, 5625, 6075, 6561, 6889, 7225, 7803, 8649, 9801, 10125, 11881, 11979, 14641, 15125, 15129, 15625, 16129, 16875

Offset: 1

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

			The terms together with their prime indices (not prime factors) begin:
    1: {}
    9: {2,2}
   25: {3,3}
   27: {2,2,2}
   81: {2,2,2,2}
  121: {5,5}
  125: {3,3,3}
  225: {2,2,3,3}
  243: {2,2,2,2,2}
  289: {7,7}
  625: {3,3,3,3}
  675: {2,2,2,3,3}
  729: {2,2,2,2,2,2}
  961: {11,11}
For example, 675 = prime(2)^3 prime(3)^2 = 3^3 * 5^2.

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

Powerful numbers are A001694, counted by A007690.
The version for prime exponents instead of indices is A056166, counted by A055923.
This is the powerful case of A076610 (products of A006450), counted by A000607.
The partitions with these Heinz numbers are counted by A339218.
A000040 lists primes.
A031368 lists primes of odd index, products A066208.
A101436 counts exponents in prime factorization that are themselves prime.
A112798 lists prime indices, reverse A296150, sum A056239.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.
A053810 lists all numbers p^q with p and q prime, counted by A230595.
A257994 counts prime indices that are themselves prime, complement A330944.
Cf. A000720, A000961, A001597, A007821, A007916, A052485, A109297, A140319, A164336, A181819, A325131, A330945, A346068.

Select[Range[1000],#==1||And@@PrimeQ/@PrimePi/@First/@FactorInteger[#]&&Min@@Last/@FactorInteger[#]>1&]

Intersection of A001694 and A076610.

Sum_{n>=1} 1/a(n) = Product_{p in A006450} (1 + 1/(p*(p-1))) = 1.24410463... - Amiram Eldar, May 04 2022

A353395 Numbers k such that the prime shadow of k equals the product of prime shadows of the prime indices of k.

Original entry on oeis.org

1, 3, 5, 11, 15, 17, 26, 31, 33, 41, 51, 55, 58, 59, 67, 78, 83, 85, 86, 93, 94, 109, 123, 126, 127, 130, 146, 148, 155, 157, 158, 165, 174, 177, 179, 187, 191, 196, 201, 202, 205, 211, 241, 244, 249, 255, 258, 274, 277, 278, 282, 283, 284, 286, 290, 295, 298

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: {}         78: {1,2,6}      158: {1,22}
      3: {2}        83: {23}         165: {2,3,5}
      5: {3}        85: {3,7}        174: {1,2,10}
     11: {5}        86: {1,14}       177: {2,17}
     15: {2,3}      93: {2,11}       179: {41}
     17: {7}        94: {1,15}       187: {5,7}
     26: {1,6}     109: {29}         191: {43}
     31: {11}      123: {2,13}       196: {1,1,4,4}
     33: {2,5}     126: {1,2,2,4}    201: {2,19}
     41: {13}      127: {31}         202: {1,26}
     51: {2,7}     130: {1,3,6}      205: {3,13}
     55: {3,5}     146: {1,21}       211: {47}
     58: {1,10}    148: {1,1,12}     241: {53}
     59: {17}      155: {3,11}       244: {1,1,18}
     67: {19}      157: {37}         249: {2,23}
For example, 126 is in the sequence because its prime indices {1,2,2,4} have shadows {1,2,2,3}, with product 12, which is also the prime shadow of 126.

The prime terms are A006450.
The LHS (prime shadow) is A181819, with an inverse A181821.
The RHS (product of shadows) is A353394, first appearances A353397.
This is a ranking of the partitions counted by A353396.
Another related comparison is A353399, counted by A353398.
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.
A130091 lists numbers with distinct prime exponents, counted by A098859.
A324850 lists numbers divisible by the product of their prime indices.
Numbers divisible by their prime shadow:
- counted by A325702
- listed by A325755
- co-recursive version A325756
- nonprime recursive version A353389
- recursive version A353393, counted by A353426
Cf. A000005, A000961, A003586, A005117, A143773, A182850, A316428, A316438, A320325, A325131, A339095.

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[#]==red[#]&]

A181819(a(n)) = A353394(a(n)) = Product_i A181819(A112798(a(n),i)).

A352140 Numbers whose prime factorization has all even prime indices and all odd exponents.

Original entry on oeis.org

1, 3, 7, 13, 19, 21, 27, 29, 37, 39, 43, 53, 57, 61, 71, 79, 87, 89, 91, 101, 107, 111, 113, 129, 131, 133, 139, 151, 159, 163, 173, 181, 183, 189, 193, 199, 203, 213, 223, 229, 237, 239, 243, 247, 251, 259, 263, 267, 271, 273, 281, 293, 301, 303, 311, 317

Offset: 1

Gus Wiseman, Mar 11 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, sum A056239, length A001222.
A number's prime signature is the sequence of positive exponents in its prime factorization, which is row n of A124010, length A001221, sum A001222.
Also Heinz numbers of integer partitions with all even parts and all odd multiplicities, counted by A055922 aerated.
All terms are odd. - Michael S. Branicky, Mar 12 2022

			The terms together with their prime indices begin:
      1 = 1
      3 = prime(2)^1
      7 = prime(4)^1
     13 = prime(6)^1
     19 = prime(8)^1
     21 = prime(4)^1 prime(2)^1
     27 = prime(2)^3
     29 = prime(10)^1
     37 = prime(12)^1
     39 = prime(6)^1 prime(2)^1
     43 = prime(14)^1
     53 = prime(16)^1
     57 = prime(8)^1 prime(2)^1
     61 = prime(18)^1
     71 = prime(20)^1

The restriction to primes is A031215.
These partitions are counted by A055922 (aerated).
The first condition alone is A066207, counted by A035363.
The squarefree case is A258117.
The second condition alone is A268335, counted by A055922.
A056166 = exponents all prime, counted by A055923.
A066208 = prime indices all odd, counted by A000009.
A109297 = same indices as exponents, counted by A114640.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.
A162641 counts even prime exponents, odd A162642.
A257991 counts odd prime indices, even A257992.
A325131 = disjoint indices from exponents, counted by A114639.
A346068 = indices and exponents all prime, counted by A351982.
A351979 = odd indices with even exponents, counted by A035457.
A352141 = even indices with even exponents, counted by A035444.
A352142 = odd indices with odd exponents, counted by A117958.
Cf. A000720, A028260, A055396, A061395, A181819, A195017, A241638, A276078, A324517, A324524, A324525, A325698.

Select[Range[100],And@@EvenQ/@PrimePi/@First/@FactorInteger[#]&&And@@OddQ/@Last/@FactorInteger[#]&]

from sympy import factorint, primepi
def ok(n):
    if n%2 == 0: return False
    return all(primepi(p)%2==0 and e%2==1 for p, e in factorint(n).items())
print([k for k in range(318) if ok(k)]) # Michael S. Branicky, Mar 12 2022

Intersection of A066207 and A268335.
A257991(a(n)) = A162641(a(n)) = 0.
A162642(a(n)) = A001221(a(n)).
A257992(a(n)) = A001222(a(n)).

A351979 Numbers whose prime factorization has all odd prime indices and all even prime exponents.

Original entry on oeis.org

1, 4, 16, 25, 64, 100, 121, 256, 289, 400, 484, 529, 625, 961, 1024, 1156, 1600, 1681, 1936, 2116, 2209, 2500, 3025, 3481, 3844, 4096, 4489, 4624, 5329, 6400, 6724, 6889, 7225, 7744, 8464, 8836, 9409, 10000, 10609, 11881, 12100, 13225, 13924, 14641, 15376

Offset: 1

Gus Wiseman, Mar 11 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, sum A056239, length A001222.
A number's prime signature is the sequence of positive exponents in its prime factorization, which is row n of A124010, length A001221, sum A001222.
Also Heinz numbers of integer partitions with all odd parts and all even multiplicities, counted by A035457 (see Emeric Deutsch's comment there).

			The terms together with their prime indices begin:
     1: 1
     4: prime(1)^2
    16: prime(1)^4
    25: prime(3)^2
    64: prime(1)^6
   100: prime(1)^2 prime(3)^2
   121: prime(5)^2
   256: prime(1)^8
   289: prime(7)^2
   400: prime(1)^4 prime(3)^2
   484: prime(1)^2 prime(5)^2
   529: prime(9)^2
   625: prime(3)^4
   961: prime(11)^2
  1024: prime(1)^10
  1156: prime(1)^2 prime(7)^2
  1600: prime(1)^6 prime(3)^2
  1681: prime(13)^2
  1936: prime(1)^4 prime(5)^2

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

The second condition alone (exponents all even) is A000290, counted by A035363.
The distinct prime factors of terms all come from A031368.
These partitions are counted by A035457 or A000009 aerated.
The first condition alone (indices all odd) is A066208, counted by A000009.
The squarefree square roots are A258116, even A258117.
A056166 = exponents all prime, counted by A055923.
A066207 = indices all even, counted by complement of A086543.
A076610 = indices all prime, counted by A000607.
A109297 = same indices as exponents, counted by A114640.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.
A162641 counts even exponents, odd A162642.
A257991 counts odd indices, even A257992.
A268335 = exponents all odd, counted by A055922.
A325131 = disjoint indices from exponents, counted by A114639.
A346068 = indices and exponents all prime, counted by A351982.
A352140 = even indices with odd exponents, counted by A055922 (aerated).
A352141 = even indices with even exponents, counted by A035444.
A352142 = odd indices and odd multiplicities, counted by A117958.
Cf. A000720, A028260, A045931, A055396, A061395, A106529, A181819, A195017, A276078, A324588, A325698, A325700.

Select[Range[1000],#==1||And@@OddQ/@PrimePi/@First/@FactorInteger[#]&&And@@EvenQ/@Last/@FactorInteger[#]&]

from sympy import factorint, primepi
def ok(n):
    return all(primepi(p)%2==1 and e%2==0 for p, e in factorint(n).items())
print([k for k in range(15500) if ok(k)]) # Michael S. Branicky, Mar 12 2022

Squares of elements of A066208.
Intersection of A066208 and A000290.
A257991(a(n)) = A001222(a(n)).
A162641(a(n)) = A001221(a(n)).
A162642(a(n)) = A257992(a(n)) = 0.
Sum_{n>=1} 1/a(n) = 1/Product_{k>=1} (1 - 1/prime(2*k-1)^2) = 1.4135142... . - Amiram Eldar, Sep 19 2022

A328830 The second prime shadow of n: a(1) = 1; for n > 1, a(n) = a(A003557(n)) * prime(A056169(n)) when A056169(n) > 0, otherwise a(n) = a(A003557(n)).

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Oct 29 2019

nonn

a(n) depends only on prime signature of n (cf. A025487).

			For n = 30 = 2 * 3 * 5, there are three unitary prime factors, while A003557(30) = 1, which terminates the recursion, thus a(30) = prime(3) = 5.
For n = 60060 = 2^2 * 3 * 5 * 7 * 11 * 13, there are 5 unitary prime factors, while in A003557(60060) = 2 there is only one, thus a(60060) = prime(5) * prime(1) = 11 * 2 = 22.
The number 1260 = 2^2*3^2*5*7 has prime exponents (2,2,1,1) so its prime shadow is prime(2)*prime(2)*prime(1)*prime(1) = 36.  Next, 36 = 2^2*3^2 has prime exponents (2,2) so its prime shadow is prime(2)*prime(2) = 9. In fact, the term a(1260) = 9 is the first appearance of 9 in the sequence. - _Gus Wiseman_, Apr 28 2022

Antti Karttunen, Table of n, a(n) for n = 1..16384
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000
Index entries for sequences computed from exponents in factorization of n

Column 2 of A353510.
Cf. A001221, A001222, A003557, A008578, A056169, A071625, A323022.
Differs from A182860 for the first time at a(30) = 5, while A182860(30) = 4.
Cf. A182863 for the first appearances.
A005361 gives product of prime exponents.
A112798 gives prime indices, sum A056239.
A124010 gives prime signature, sorted A118914.
A181819 gives prime shadow, with an inverse A181821.
A325131 lists numbers relatively prime to their prime shadow.
A325755 lists numbers divisible by their prime shadow.
Cf. A109297, A182850, A353393, A353507.

A003557(n) = { my(f=factor(n)); for (i=1, #f~, f[i, 2] = f[i, 2]-1); factorback(f); }; \\ From A003557
A056169(n) = { my(f=factor(n)[, 2]); sum(i=1, #f, f[i]==1); }; \\ From A056169
A328830(n) = if(1==n,n, my(u=A056169(n)); if(0==u,1,prime(u)) * A328830(A003557(n)));

a(1) = 1; for n > 1, a(n) = A008578(1+A056169(n)) * a(A003557(n)).
A001221(a(n)) = A323022(n).
A001222(a(n)) = A071625(n).
a(n) = A181819(A181819(n)). - Gus Wiseman, Apr 27 2022

Added Gus Wiseman's new name to the front of the definition. - Antti Karttunen, Apr 27 2022

cp's OEIS Frontend

A353500 Numbers that are the smallest number with product of prime exponents k for some k. Sorted positions of first appearances in A005361, unsorted version A085629.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A353389 Create the sequence of all positive integers > 1 that are prime or whose prime shadow (A181819) is a divisor that is already in the sequence. Then remove all the primes.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A353398 Number of integer partitions of n where the product of multiplicities equals the product of prime shadows of the parts.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A352142 Numbers whose prime factorization has all odd indices and all odd exponents.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Formula

A352141 Numbers whose prime factorization has all even indices and all even exponents.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Formula

A352492 Powerful numbers whose prime indices are all prime numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A353395 Numbers k such that the prime shadow of k equals the product of prime shadows of the prime indices of k.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A352140 Numbers whose prime factorization has all even prime indices and all odd exponents.

Views

Author

Keywords

Comments

Examples