cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 71-80 of 114 results. Next

A372590 Numbers whose binary weight (A000120) plus bigomega (A001222) is odd.

Original entry on oeis.org

1, 3, 4, 5, 12, 14, 16, 17, 18, 20, 21, 22, 23, 25, 26, 27, 29, 30, 35, 38, 43, 45, 48, 49, 53, 55, 56, 62, 63, 64, 66, 68, 69, 71, 72, 74, 75, 78, 80, 81, 82, 83, 84, 87, 88, 89, 91, 92, 93, 94, 99, 100, 101, 102, 104, 105, 108, 113, 114, 115, 116, 118, 120
Offset: 1

Views

Author

Gus Wiseman, May 14 2024

Keywords

Comments

The even version is A372591.

Examples

			The terms (center), their binary indices (left), and their weakly decreasing prime indices (right) begin:
        {1}   1  ()
      {1,2}   3  (2)
        {3}   4  (1,1)
      {1,3}   5  (3)
      {3,4}  12  (2,1,1)
    {2,3,4}  14  (4,1)
        {5}  16  (1,1,1,1)
      {1,5}  17  (7)
      {2,5}  18  (2,2,1)
      {3,5}  20  (3,1,1)
    {1,3,5}  21  (4,2)
    {2,3,5}  22  (5,1)
  {1,2,3,5}  23  (9)
    {1,4,5}  25  (3,3)
    {2,4,5}  26  (6,1)
  {1,2,4,5}  27  (2,2,2)
  {1,3,4,5}  29  (10)
  {2,3,4,5}  30  (3,2,1)
    {1,2,6}  35  (4,3)
    {2,3,6}  38  (8,1)
  {1,2,4,6}  43  (14)
  {1,3,4,6}  45  (3,2,2)
		

Crossrefs

For sum (A372428, zeros A372427) we have A372586, complement A372587.
For minimum (A372437) we have A372439, complement A372440.
Positions of odd terms in A372441, zeros A071814.
For maximum (A372442, zeros A372436) we have A372588, complement A372589.
The complement is A372591.
For just binary indices:
- length: A000069, complement A001969
- sum: A158705, complement A158704
- minimum: A003159, complement A036554
- maximum: A053738, complement A053754
For just prime indices:
- length: A026424 (count A027193), complement A028260 (count A027187)
- sum: A300063 (count A058695), complement A300061 (count A058696)
- minimum: A340932 (count A026804), complement A340933 (count A026805)
- maximum: A244991 (count A027193), complement A244990 (count A027187)
A005408 lists odd numbers.
A019565 gives Heinz number of binary indices, adjoint A048675.
A029837 gives greatest binary index, least A001511.
A031368 lists odd-indexed primes, even A031215.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A070939 gives length of binary expansion.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.

Programs

  • Mathematica
    Select[Range[100],OddQ[DigitCount[#,2,1]+PrimeOmega[#]]&]

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

Views

Author

Gus Wiseman, Nov 08 2020

Keywords

Comments

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

Examples

			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}
		

Crossrefs

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.
Subsequence of A332820.

Programs

  • Mathematica
    Select[Range[1000],PrimeOmega[#]==3&&OddQ[Times@@(1+PrimePi/@First/@FactorInteger[#])]&]
  • PARI
    isok(m) = my(f=factor(m)); (bigomega(f)==3) && (#select(x->(x%2), apply(primepi, f[,1]~)) == 0); \\ Michel Marcus, Nov 10 2020
    
  • Python
    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

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

Views

Author

Gus Wiseman, Mar 18 2022

Keywords

Comments

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.

Examples

			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
		

Crossrefs

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.

Programs

  • Mathematica
    Select[Range[1000],#==1||And@@EvenQ/@PrimePi/@First/@FactorInteger[#]&&And@@EvenQ/@Last/@FactorInteger[#]&]
  • Python
    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

Formula

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

Views

Author

Gus Wiseman, Mar 24 2022

Keywords

Comments

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.

Examples

			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.
		

Crossrefs

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.

Programs

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

Formula

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

A366850 Number of integer partitions of n whose odd parts are relatively prime.

Original entry on oeis.org

0, 1, 1, 2, 3, 5, 7, 11, 16, 22, 32, 43, 60, 80, 110, 140, 194, 244, 327, 410, 544, 670, 883, 1081, 1401, 1708, 2195, 2651, 3382, 4069, 5129, 6157, 7708, 9194, 11438, 13599, 16788, 19911, 24432, 28858, 35229, 41507, 50359, 59201, 71489, 83776, 100731, 117784
Offset: 0

Views

Author

Gus Wiseman, Oct 28 2023

Keywords

Examples

			The a(1) = 1 through a(8) = 16 partitions:
  (1)  (11)  (21)   (31)    (41)     (51)      (61)       (53)
             (111)  (211)   (221)    (321)     (331)      (71)
                    (1111)  (311)    (411)     (421)      (431)
                            (2111)   (2211)    (511)      (521)
                            (11111)  (3111)    (2221)     (611)
                                     (21111)   (3211)     (3221)
                                     (111111)  (4111)     (3311)
                                               (22111)    (4211)
                                               (31111)    (5111)
                                               (211111)   (22211)
                                               (1111111)  (32111)
                                                          (41111)
                                                          (221111)
                                                          (311111)
                                                          (2111111)
                                                          (11111111)
		

Crossrefs

For all parts (not just odd) we have A000837, complement A018783.
The complement is counted by A366842.
These partitions have ranks A366846.
A000041 counts integer partitions, strict A000009 (also into odds).
A000740 counts relatively prime compositions.
A078374 counts relatively prime strict partitions.
A113685 counts partitions by sum of odd parts, rank statistic A366528.
A168532 counts partitions by gcd.
A239261 counts partitions with (sum of odd parts) = (sum of even parts).

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],GCD@@Select[#,OddQ]==1&]],{n,0,30}]

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

Views

Author

Gus Wiseman, Nov 08 2020

Keywords

Comments

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

Examples

			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}
		

Crossrefs

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).
Subsequence of A332820.

Programs

  • Maple
    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
  • Mathematica
    Select[Range[100],PrimeOmega[#]==3&&OddQ[Times@@PrimePi/@First/@FactorInteger[#]]&]
  • PARI
    isok(m) = my(f=factor(m)); (bigomega(f)==3) && (#select(x->!(x%2), apply(primepi, f[,1]~)) == 0); \\ Michel Marcus, Nov 10 2020
    
  • Python
    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
    
  • Python
    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

Views

Author

Gus Wiseman, Nov 08 2020

Keywords

Comments

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.

Examples

			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}
		

Crossrefs

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.

Programs

  • Mathematica
    Select[Range[1000],SquareFreeQ[#]&&PrimeOmega[#]==3&&OddQ[Times@@(1+PrimePi/@First/@FactorInteger[#])]&]
  • PARI
    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
    
  • Python
    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

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

Views

Author

Gus Wiseman, Mar 11 2022

Keywords

Comments

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

Examples

			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
		

Crossrefs

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.

Programs

  • Mathematica
    Select[Range[100],And@@EvenQ/@PrimePi/@First/@FactorInteger[#]&&And@@OddQ/@Last/@FactorInteger[#]&]
  • Python
    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

Formula

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

A356930 Numbers whose prime indices have all odd prime indices. MM-numbers of finite multisets of finite multisets of odd numbers.

Original entry on oeis.org

1, 2, 3, 4, 6, 7, 8, 9, 11, 12, 14, 16, 18, 19, 21, 22, 24, 27, 28, 29, 31, 32, 33, 36, 38, 42, 44, 48, 49, 53, 54, 56, 57, 58, 59, 62, 63, 64, 66, 71, 72, 76, 77, 79, 81, 83, 84, 87, 88, 93, 96, 97, 98, 99, 106, 108, 112, 114, 116, 118, 121, 124, 126, 127
Offset: 1

Views

Author

Gus Wiseman, Sep 11 2022

Keywords

Comments

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 multiset of multisets with MM-number n to be formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. The combined size of this multiset of multisets is A302242(n). For example, the prime indices of 78 are {1,2,6}, so the multiset of multisets with MM-number 78 is {{},{1},{1,2}}.

Examples

			The initial terms and corresponding multisets of multisets:
   1: {}
   2: {{}}
   3: {{1}}
   4: {{},{}}
   6: {{},{1}}
   7: {{1,1}}
   8: {{},{},{}}
   9: {{1},{1}}
  11: {{3}}
  12: {{},{},{1}}
  14: {{},{1,1}}
  16: {{},{},{},{}}
  18: {{},{1},{1}}
  19: {{1,1,1}}
  21: {{1},{1,1}}
  22: {{},{3}}
  24: {{},{},{},{1}}
  27: {{1},{1},{1}}
  28: {{},{},{1,1}}
  29: {{1,3}}
  31: {{5}}
  32: {{},{},{},{},{}}
		

Crossrefs

Multisets of odd numbers are counted by A000009, ranked by A066208.
Factorizations of this type are counted by A356931.
The version for odd lengths instead of parts is A356935, ranked by A089259.
Other conditions: A302478, A302492, A356939, A356940, A356944, A356955.
A000041 counts integer partitions, strict A000009.
A000688 counts factorizations into prime powers.
A001055 counts factorizations.
A001221 counts prime divisors, sum A001414.
A001222 counts prime factors with multiplicity.
A056239 adds up prime indices, row sums of A112798.

Programs

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

A366846 Numbers whose odd prime indices are relatively prime.

Original entry on oeis.org

2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42, 44, 46, 48, 50, 52, 54, 55, 56, 58, 60, 62, 64, 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 85, 86, 88, 90, 92, 94, 96, 98, 100, 102, 104, 106, 108, 110, 112, 114, 116, 118, 120, 122
Offset: 1

Views

Author

Gus Wiseman, Oct 29 2023

Keywords

Comments

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.

Examples

			The odd prime indices of 115 are {3,9}, and these are not relatively prime, so 115 is not in the sequence.
The odd prime indices of 825 are {3,3,5}, and these are relatively prime, so 825 is in the sequence.
		

Crossrefs

Including even indices gives A289509, ones of A289508, counted by A000837.
The complement when including even indices is A318978, counted by A018783.
The nonzero complement ranks the partitions counted by A366842.
The version for halved even indices is A366847.
The odd case is A366848.
The partitions with these Heinz numbers are counted by A366850.
A000041 counts integer partitions, strict A000009 (also into odds).
A112798 lists prime indices, length A001222, sum A056239.
A257992 counts even prime indices, odd A257991.
A366528 adds up odd prime indices, partition triangle A113685.
A366531 = 2*A366533 adds up even prime indices, triangle A113686/A174713.

Programs

  • Mathematica
    Select[Range[100], GCD@@Select[PrimePi/@First/@FactorInteger[#], OddQ]==1&]
Previous Showing 71-80 of 114 results. Next