A258116 -id:A258116 - cp's OEIS frontend

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

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.

A338556 Products of three prime numbers of even index.

Original entry on oeis.org

27, 63, 117, 147, 171, 261, 273, 333, 343, 387, 399, 477, 507, 549, 609, 637, 639, 711, 741, 777, 801, 903, 909, 931, 963, 1017, 1083, 1113, 1131, 1179, 1183, 1251, 1281, 1359, 1421, 1443, 1467, 1491, 1557, 1629, 1653, 1659, 1677, 1729, 1737, 1791, 1813, 1869

Offset: 1

Gus Wiseman, Nov 08 2020

nonn

All terms are odd.

Also Heinz numbers of integer partitions with 3 parts, all of which are even. These partitions are counted by A001399.

			The sequence of terms together with their prime indices begins:
      27: {2,2,2}      637: {4,4,6}     1183: {4,6,6}
      63: {2,2,4}      639: {2,2,20}    1251: {2,2,34}
     117: {2,2,6}      711: {2,2,22}    1281: {2,4,18}
     147: {2,4,4}      741: {2,6,8}     1359: {2,2,36}
     171: {2,2,8}      777: {2,4,12}    1421: {4,4,10}
     261: {2,2,10}     801: {2,2,24}    1443: {2,6,12}
     273: {2,4,6}      903: {2,4,14}    1467: {2,2,38}
     333: {2,2,12}     909: {2,2,26}    1491: {2,4,20}
     343: {4,4,4}      931: {4,4,8}     1557: {2,2,40}
     387: {2,2,14}     963: {2,2,28}    1629: {2,2,42}
     399: {2,4,8}     1017: {2,2,30}    1653: {2,8,10}
     477: {2,2,16}    1083: {2,8,8}     1659: {2,4,22}
     507: {2,6,6}     1113: {2,4,16}    1677: {2,6,14}
     549: {2,2,18}    1131: {2,6,10}    1729: {4,6,8}
     609: {2,4,10}    1179: {2,2,32}    1737: {2,2,44}

A014612 allows all prime indices (not just even) (strict: A007304).
A066207 allows products of any length (strict: A258117).
A338471 is the version for odds instead of evens (strict: A307534).
A338557 is the strict case.
A014311 is a ranking of ordered triples (strict: A337453).
A001399(n-3) counts 3-part partitions (strict: A001399(n-6)).
A005117 lists squarefree numbers, with even case A039956.
A008284 counts partitions by sum and length (strict: A008289).
A023023 counts 3-part relatively prime partitions (strict: A101271).
A046316 lists products of exactly three odd primes (strict: A046389).
A066208 lists numbers with all odd prime indices (strict: A258116).
A075818 lists even Heinz numbers of 3-part partitions (strict: A075819).
A307719 counts 3-part pairwise coprime partitions (strict: A220377).
A285508 lists Heinz numbers of non-strict triples.
Cf. A000217, A001221, A001222, A037144, A056239, A112798, A337599, A337600.
Subsequence of A332820.

Select[Range[1000],PrimeOmega[#]==3&&OddQ[Times@@(1+PrimePi/@First/@FactorInteger[#])]&]

isok(m) = my(f=factor(m)); (bigomega(f)==3) && (#select(x->(x%2), apply(primepi, f[,1]~)) == 0); \\ Michel Marcus, Nov 10 2020

from itertools import filterfalse
from math import isqrt
from sympy import primepi, primerange, integer_nthroot
def A338556(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return int(n+x-sum((primepi(x//(k*m))>>1)-(b>>1)+1 for a,k in filterfalse(lambda x:x[0]&1,enumerate(primerange(3,integer_nthroot(x,3)[0]+1),2)) for b,m in filterfalse(lambda x:x[0]&1,enumerate(primerange(k,isqrt(x//k)+1),a))))
    return bisection(f,n,n) # Chai Wah Wu, Oct 18 2024

A238744 Irregular table read by rows: T (n, k) gives the number of primes p such that p^k divides n; table omits all zero values.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 2, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 2, 2, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 1, 1, 3, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 1, 2

Offset: 2

Matthew Vandermast, Apr 28 2014

If the prime signature of n (nonincreasing version) is viewed as a partition, row n gives the conjugate partition.

			24 = 2^3*3 is divisible by two prime numbers (2 and 3), one square of a prime (4 = 2^2), and one cube of a prime (8 = 2^3); therefore, row 24 of the table is {2,1,1}.
From _Gus Wiseman_, Mar 31 2022: (Start)
Rows begin:
     1: ()        16: (1,1,1,1)    31: (1)
     2: (1)       17: (1)          32: (1,1,1,1,1)
     3: (1)       18: (2,1)        33: (2)
     4: (1,1)     19: (1)          34: (2)
     5: (1)       20: (2,1)        35: (2)
     6: (2)       21: (2)          36: (2,2)
     7: (1)       22: (2)          37: (1)
     8: (1,1,1)   23: (1)          38: (2)
     9: (1,1)     24: (2,1,1)      39: (2)
    10: (2)       25: (1,1)        40: (2,1,1)
    11: (1)       26: (2)          41: (1)
    12: (2,1)     27: (1,1,1)      42: (3)
    13: (1)       28: (2,1)        43: (1)
    14: (2)       29: (1)          44: (2,1)
    15: (2)       30: (3)          45: (2,1)
(End)

Row lengths are A051903(n); row sums are A001222(n).
Cf. A217171.
These partitions are ranked by A238745.
For prime indices A296150 instead of exponents we get A321649, rev A321650.
A000700 counts self-conjugate partitions, ranked by A088902.
A003963 gives product of prime indices, conjugate A329382.
A008480 gives number of permutations of prime indices, conjugate A321648.
A056239 adds up prime indices, row sums of A112798.
A124010 gives prime signature, sorted A118914, length A001221.
A352486-A352490 are sets related to the fixed points of A122111.
Cf. A000701, A000720, A046682, A238747, A258116, A319005, A330644, A352491.

Table[Length/@Table[Select[Last/@FactorInteger[n],#>=k&],{k,Max@@Last/@FactorInteger[n]}],{n,2,100}] (* Gus Wiseman, Mar 31 2022 *)

Row n is identical to row A124859(n) of table A212171.

A338471 Products of three prime numbers of odd index.

Original entry on oeis.org

8, 20, 44, 50, 68, 92, 110, 124, 125, 164, 170, 188, 230, 236, 242, 268, 275, 292, 310, 332, 374, 388, 410, 412, 425, 436, 470, 506, 508, 548, 575, 578, 590, 596, 605, 628, 668, 670, 682, 716, 730, 764, 775, 782, 788, 830, 844, 902, 908, 932, 935, 964, 970

Offset: 1

Gus Wiseman, Nov 08 2020

nonn

Also Heinz numbers of integer partitions with 3 parts, all of which are odd. These partitions are counted by A001399.

			The sequence of terms together with their prime indices begins:
       8: {1,1,1}      268: {1,1,19}     575: {3,3,9}
      20: {1,1,3}      275: {3,3,5}      578: {1,7,7}
      44: {1,1,5}      292: {1,1,21}     590: {1,3,17}
      50: {1,3,3}      310: {1,3,11}     596: {1,1,35}
      68: {1,1,7}      332: {1,1,23}     605: {3,5,5}
      92: {1,1,9}      374: {1,5,7}      628: {1,1,37}
     110: {1,3,5}      388: {1,1,25}     668: {1,1,39}
     124: {1,1,11}     410: {1,3,13}     670: {1,3,19}
     125: {3,3,3}      412: {1,1,27}     682: {1,5,11}
     164: {1,1,13}     425: {3,3,7}      716: {1,1,41}
     170: {1,3,7}      436: {1,1,29}     730: {1,3,21}
     188: {1,1,15}     470: {1,3,15}     764: {1,1,43}
     230: {1,3,9}      506: {1,5,9}      775: {3,3,11}
     236: {1,1,17}     508: {1,1,31}     782: {1,7,9}
     242: {1,5,5}      548: {1,1,33}     788: {1,1,45}

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

A066208 allows products of any length (strict: A258116).
A307534 is the squarefree case.
A338469 is the restriction to odds.
A338556 is the version for evens (strict: A338557).
A000009 counts partitions into odd parts (strict: A000700).
A001399(n-3) counts 3-part partitions (strict: A001399(n-6)).
A008284 counts partitions by sum and length.
A014311 is a ranking of ordered triples (strict: A337453).
A014612 lists Heinz numbers of all triples (strict: A007304).
A023023 counts 3-part relatively prime partitions (strict: A101271).
A023023 counts 3-part relatively prime partitions (strict: A078374).
A046316 lists products of exactly three odd primes (strict: A046389).
A066207 lists numbers with all even prime indices (strict: A258117).
A075818 lists even Heinz numbers of 3-part partitions (strict: A075819).
A285508 lists Heinz numbers of non-strict triples.
A307719 counts 3-part pairwise coprime partitions (strict: A220377).
Cf. A000217, A001221, A001222, A005117, A037144, A056239, A112798.
Subsequence of A332820.

N:= 1000: # for terms <= N
R:= NULL:
for i from 1 by 2 do
  p:= ithprime(i);
  if p^3 >= N then break fi;
  for j from i by 2 do
    q:= ithprime(j);
    if p*q^2 >= N then break fi;
    for k from j by 2 do
      x:= p*q*ithprime(k);
      if x > N then break fi;
      R:= R,x;
od od od:
sort([R]); # Robert Israel, Jun 11 2025

Select[Range[100],PrimeOmega[#]==3&&OddQ[Times@@PrimePi/@First/@FactorInteger[#]]&]

isok(m) = my(f=factor(m)); (bigomega(f)==3) && (#select(x->!(x%2), apply(primepi, f[,1]~)) == 0); \\ Michel Marcus, Nov 10 2020

from sympy import primerange
from itertools import combinations_with_replacement as mc
def aupto(limit):
    pois = [p for i, p in enumerate(primerange(2, limit//4+1)) if i%2 == 0]
    return sorted(set(a*b*c for a, b, c in mc(pois, 3) if a*b*c <= limit))
print(aupto(971)) # Michael S. Branicky, Aug 20 2021

from math import isqrt
from sympy import primepi, primerange, integer_nthroot
def A338471(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return int(n+x-sum((primepi(x//(k*m))+1>>1)-(b+1>>1)+1 for a,k in filter(lambda x:x[0]&1,enumerate(primerange(integer_nthroot(x,3)[0]+1),1)) for b,m in filter(lambda x:x[0]&1,enumerate(primerange(k,isqrt(x//k)+1),a))))
    return bisection(f,n,n) # Chai Wah Wu, Oct 18 2024

A338557 Products of three distinct prime numbers of even index.

Original entry on oeis.org

273, 399, 609, 741, 777, 903, 1113, 1131, 1281, 1443, 1491, 1653, 1659, 1677, 1729, 1869, 2067, 2109, 2121, 2247, 2373, 2379, 2451, 2639, 2751, 2769, 2919, 3021, 3081, 3171, 3219, 3367, 3423, 3471, 3477, 3633, 3741, 3801, 3857, 3913, 3939, 4047, 4053, 4173

Offset: 1

Gus Wiseman, Nov 08 2020

nonn

All terms are odd.
Also sphenic numbers (A007304) with all even prime indices (A031215).
Also Heinz numbers of strict integer partitions with 3 parts, all of which are even. These partitions are counted by A001399.

			The sequence of terms together with their prime indices begins:
     273: {2,4,6}     1869: {2,4,24}    3219: {2,10,12}
     399: {2,4,8}     2067: {2,6,16}    3367: {4,6,12}
     609: {2,4,10}    2109: {2,8,12}    3423: {2,4,38}
     741: {2,6,8}     2121: {2,4,26}    3471: {2,6,24}
     777: {2,4,12}    2247: {2,4,28}    3477: {2,8,18}
     903: {2,4,14}    2373: {2,4,30}    3633: {2,4,40}
    1113: {2,4,16}    2379: {2,6,18}    3741: {2,10,14}
    1131: {2,6,10}    2451: {2,8,14}    3801: {2,4,42}
    1281: {2,4,18}    2639: {4,6,10}    3857: {4,8,10}
    1443: {2,6,12}    2751: {2,4,32}    3913: {4,6,14}
    1491: {2,4,20}    2769: {2,6,20}    3939: {2,6,26}
    1653: {2,8,10}    2919: {2,4,34}    4047: {2,8,20}
    1659: {2,4,22}    3021: {2,8,16}    4053: {2,4,44}
    1677: {2,6,14}    3081: {2,6,22}    4173: {2,6,28}
    1729: {4,6,8}     3171: {2,4,36}    4179: {2,4,46}

For the following, NNS means "not necessarily strict".
A007304 allows all prime indices (not just even) (NNS: A014612).
A046389 allows all odd primes (NNS: A046316).
A258117 allows products of any length (NNS: A066207).
A307534 is the version for odds instead of evens (NNS: A338471).
A337453 is a different ranking of ordered triples (NNS: A014311).
A338556 is the NNS version.
A001399(n-6) counts strict 3-part partitions (NNS: A001399(n-3)).
A005117 lists squarefree numbers, with even case A039956.
A078374 counts 3-part relatively prime strict partitions (NNS: A023023).
A075819 lists even Heinz numbers of strict triples (NNS: A075818).
A220377 counts 3-part pairwise coprime strict partitions (NNS: A307719).
A258116 lists squarefree numbers with all odd prime indices (NNS: A066208).
A285508 lists Heinz numbers of non-strict triples.
Cf. A000217, A001221, A001222, A037144, A056239, A112798, A337605.

Select[Range[1000],SquareFreeQ[#]&&PrimeOmega[#]==3&&OddQ[Times@@(1+PrimePi/@First/@FactorInteger[#])]&]

isok(m) = my(f=factor(m)); (bigomega(f)==3) && (omega(f)==3) && (#select(x->(x%2), apply(primepi, f[,1]~)) == 0); \\ Michel Marcus, Nov 10 2020

from itertools import filterfalse
from math import isqrt
from sympy import primepi, primerange, nextprime, integer_nthroot
def A338557(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return int(n+x-sum((primepi(x//(k*m))>>1)-(b>>1) for a,k in filterfalse(lambda x:x[0]&1,enumerate(primerange(3,integer_nthroot(x,3)[0]+1),2)) for b,m in filterfalse(lambda x:x[0]&1,enumerate(primerange(nextprime(k)+1,isqrt(x//k)+1),a+2))))
    return bisection(f,n,n) # Chai Wah Wu, Oct 18 2024

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

A352143 Numbers whose prime indices and conjugate prime indices are all odd.

Original entry on oeis.org

1, 2, 5, 8, 11, 17, 20, 23, 31, 32, 41, 44, 47, 59, 67, 68, 73, 80, 83, 92, 97, 103, 109, 124, 125, 127, 128, 137, 149, 157, 164, 167, 176, 179, 188, 191, 197, 211, 227, 233, 236, 241, 257, 268, 269, 272, 275, 277, 283, 292, 307, 313, 320, 331, 332, 347, 353

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 whose parts and conjugate parts are all odd. They are counted by A053253.

			The terms together with their prime indices begin:
   1: {}
   2: {1}
   5: {3}
   8: {1,1,1}
  11: {5}
  17: {7}
  20: {1,1,3}
  23: {9}
  31: {11}
  32: {1,1,1,1,1}
  41: {13}
  44: {1,1,5}
  47: {15}
  59: {17}
  67: {19}
  68: {1,1,7}
  73: {21}
  80: {1,1,1,1,3}

The restriction to primes is A031368.
These partitions appear to be counted by A053253.
The even version is A066207^2.
For even instead of odd conjugate parts we get A066208^2.
The first condition alone (all odd indices) is A066208, counted by A000009.
The second condition alone is A346635, counted by A000009.
A055922 counts partitions with odd multiplicities, ranked by A268335.
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.
A162642 counts odd prime exponents, even A162641.
A238745 gives the Heinz number of the conjugate prime signature.
A257991 counts odd indices, even A257992.
A258116 ranks strict partitions with all odd parts, even A258117.
A351979 = odd indices and even multiplicities, counted by A035457.
A352140 = even indices and odd multiplicities, counted by A055922 aerated.
A352141 = even indices and even multiplicities, counted by A035444.
A352142 = odd indices and odd multiplicities, counted by A117958.
Cf. A000290, A000701, A000720, A028260, A045931, A046682, A055396, A061395, A195017, A241638, A325698, A325700.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
Select[Range[100],And@@OddQ/@primeMS[#]&&And@@OddQ/@conj[primeMS[#]]&]

Intersection of A066208 and A346635.

A338469 Products of three odd prime numbers of odd index.

Original entry on oeis.org

125, 275, 425, 575, 605, 775, 935, 1025, 1175, 1265, 1331, 1445, 1475, 1675, 1705, 1825, 1955, 2057, 2075, 2255, 2425, 2575, 2585, 2635, 2645, 2725, 2783, 3175, 3179, 3245, 3425, 3485, 3565, 3685, 3725, 3751, 3925, 3995, 4015, 4175, 4301, 4475, 4565, 4715

Offset: 1

Gus Wiseman, Nov 08 2020

nonn

Also Heinz numbers of integer partitions with 3 parts, all of which are odd and > 1. These partitions are counted by A001399.

			The sequence of terms together with their prime indices begins:
     125: {3,3,3}     1825: {3,3,21}    3425: {3,3,33}
     275: {3,3,5}     1955: {3,7,9}     3485: {3,7,13}
     425: {3,3,7}     2057: {5,5,7}     3565: {3,9,11}
     575: {3,3,9}     2075: {3,3,23}    3685: {3,5,19}
     605: {3,5,5}     2255: {3,5,13}    3725: {3,3,35}
     775: {3,3,11}    2425: {3,3,25}    3751: {5,5,11}
     935: {3,5,7}     2575: {3,3,27}    3925: {3,3,37}
    1025: {3,3,13}    2585: {3,5,15}    3995: {3,7,15}
    1175: {3,3,15}    2635: {3,7,11}    4015: {3,5,21}
    1265: {3,5,9}     2645: {3,9,9}     4175: {3,3,39}
    1331: {5,5,5}     2725: {3,3,29}    4301: {5,7,9}
    1445: {3,7,7}     2783: {5,5,9}     4475: {3,3,41}
    1475: {3,3,17}    3175: {3,3,31}    4565: {3,5,23}
    1675: {3,3,19}    3179: {5,7,7}     4715: {3,9,13}
    1705: {3,5,11}    3245: {3,5,17}    4775: {3,3,43}

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

A046316 allows all primes (strict: A046389).
A338471 allows all odd primes (strict: A307534).
A338556 is the version for evens (strict: A338557).
A000009 counts partitions into odd parts (strict: A000700).
A001399(n-3) counts 3-part partitions (strict: A001399(n-6)).
A005408 lists odds (strict: A056911).
A008284 counts partitions by sum and length.
A014311 is a ranking of 3-part compositions (strict: A337453).
A014612 lists Heinz numbers of 3-part partitions (strict: A007304).
A023023 counts 3-part relatively prime partitions (strict: A101271).
A066207 lists numbers with all even prime indices (strict: A258117).
A066208 lists numbers with all odd prime indices (strict: A258116).
A075818 lists even Heinz numbers of 3-part partitions (strict: A075819).
A285508 lists Heinz numbers of non-strict 3-part partitions.
Cf. A001221, A001222, A002620, A005117, A037144, A056239, A112798.

N:= 10000: # for terms <= N
P0:= [seq(ithprime(i),i=3..numtheory:-pi(floor(N/25)),2)]:
sort(select(`<=`,[seq(seq(seq(P0[i]*P0[j]*P0[k],k=1..j),j=1..i),i=1..nops(P0))], N)); # Robert Israel, Nov 12 2020

Select[Range[1,1000,2],PrimeOmega[#]==3&&OddQ[Times@@PrimePi/@First/@FactorInteger[#]]&]

isok(m) = my(f=factor(m)); (m%2) && (bigomega(f)==3) && (#select(x->!(x%2), apply(primepi, f[,1]~)) == 0); \\ Michel Marcus, Nov 10 2020

from math import isqrt
from sympy import primepi, primerange, integer_nthroot
def A338469(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return int(n+x-sum((primepi(x//(k*m))+1>>1)-(b+1>>1)+1 for a,k in filter(lambda x:x[0]&1,enumerate(primerange(5,integer_nthroot(x,3)[0]+1),3)) for b,m in filter(lambda x:x[0]&1,enumerate(primerange(k,isqrt(x//k)+1),a))))
    return bisection(f,n,n) # Chai Wah Wu, Oct 18 2024

A316094 FDH numbers of strict integer partitions with odd parts.

Original entry on oeis.org

1, 2, 4, 7, 8, 11, 14, 16, 19, 22, 25, 28, 31, 32, 38, 41, 44, 47, 50, 53, 56, 61, 62, 64, 71, 76, 77, 79, 82, 83, 88, 94, 97, 100, 101, 103, 106, 107, 109, 112, 113, 121, 122, 124, 127, 128, 131, 133, 137, 139, 142, 149, 151, 152, 154, 157, 158, 163, 164, 166

Offset: 1

Gus Wiseman, Jun 24 2018

nonn

Also numbers n such that A305829(n) is odd.

Let f(n) = A050376(n) be the n-th Fermi-Dirac prime. The FDH number of a strict integer partition (y_1,...,y_k) is f(y_1)*...*f(y_k).

			Sequence of all integer partitions with distinct odd parts begins (), (1), (3), (5), (3,1), (7), (5,1), (9), (11), (7,1), (13), (5,3), (15), (9,1), (11,1), (17), (7,3), (19), (13,1), (21), (5,3,1), (23), (15,1), (9,3), (25), (11,3), (7,5), (27), (17,1), (29), (7,3,1), (19,1), (31).

Cf. A000700, A050376, A064547, A213925, A258116, A279065, A299755, A299757, A305829.

nn=100;
FDfactor[n_]:=If[n==1,{},Sort[Join@@Cases[FactorInteger[n],{p_,k_}:>Power[p,Cases[Position[IntegerDigits[k,2]//Reverse,1],{m_}->2^(m-1)]]]]];
FDprimeList=Array[FDfactor,nn,1,Union];FDrules=MapIndexed[(#1->#2[[1]])&,FDprimeList];
Select[Range[nn],OddQ[Times@@(FDfactor[#]/.FDrules)]&]

A363219 Twice the median of the conjugate of the integer partition with Heinz number n.

Original entry on oeis.org

0, 2, 2, 4, 2, 3, 2, 6, 4, 2, 2, 4, 2, 2, 4, 8, 2, 5, 2, 2, 3, 2, 2, 5, 4, 2, 6, 2, 2, 4, 2, 10, 2, 2, 4, 6, 2, 2, 2, 2, 2, 3, 2, 2, 6, 2, 2, 6, 4, 4, 2, 2, 2, 7, 4, 2, 2, 2, 2, 4, 2, 2, 4, 12, 3, 2, 2, 2, 2, 4, 2, 7, 2, 2, 6, 2, 4, 2, 2, 2, 8, 2, 2, 3, 2, 2

Offset: 1

Gus Wiseman, May 25 2023

nonn

The median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length). Since the denominator is always 1 or 2, the median can be represented as an integer by multiplying by 2.

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.

			The partition (4,2,1) has Heinz number 42 and conjugate (3,2,1,1) with median 3/2, so a(42) = 3.

Twice the row media of A321649 or A321650.
For mean instead of twice median we have A326839/A326840.
This is the conjugate version of A360005.
A000700 counts self-conjugate partitions, ranked by A088902 (cf. A258116).
A056239 adds up prime indices, row sums of A112798 and A296150.
A122111 is partition conjugation in terms of Heinz numbers.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.
A352491 gives n minus Heinz number of conjugate.
A363220 counts partitions with same median as conjugate.
Cf. A046682, A321648, A325040, A326567/A326568, A326848, A362617, A362618, A363223, A363224.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
Table[If[n==1,0,2*Median[conj[prix[n]]]],{n,100}]

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Formula

A338556 Products of three prime numbers of even index.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Python

A238744 Irregular table read by rows: T (n, k) gives the number of primes p such that p^k divides n; table omits all zero values.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A338471 Products of three prime numbers of odd index.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Python

A338557 Products of three distinct prime numbers of even index.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Python

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Formula

A352143 Numbers whose prime indices and conjugate prime indices are all odd.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica