A302796 -id:A302796 - cp's OEIS frontend

A289509 Numbers k such that the gcd of the indices j for which the j-th prime prime(j) divides k is 1.

2, 4, 6, 8, 10, 12, 14, 15, 16, 18, 20, 22, 24, 26, 28, 30, 32, 33, 34, 35, 36, 38, 40, 42, 44, 45, 46, 48, 50, 51, 52, 54, 55, 56, 58, 60, 62, 64, 66, 68, 69, 70, 72, 74, 75, 76, 77, 78, 80, 82, 84, 85, 86, 88, 90, 92, 93, 94, 95, 96, 98, 99, 100, 102, 104

Offset: 1

Christopher J. Smyth, Jul 11 2017

nonn

Any integer k in the sequence encodes (by 'Heinz encoding' cf. A056239) a multiset of integers whose gcd is 1, namely the multiset containing r_j copies of j if k factors as Product_j prime(j)^{r_j} with gcd_j j = 1.
Clearly the sequence contains all even numbers and no odd primes or odd prime powers. It also clearly contains all numbers that are divisible by consecutive primes.
The sequence is the list of those k such that A289508(k) = 1.
It is also the list of those k such that A289506(k) = A289507(k).
Heinz numbers of integer partitions with relatively prime parts, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Apr 13 2018

			6 is a term because 6 = p_1*p_2 and gcd(1,2) = 1.
From _Gus Wiseman_, Apr 13 2018: (Start)
Sequence of integer partitions with relatively prime parts begins:
02 : (1)
04 : (11)
06 : (21)
08 : (111)
10 : (31)
12 : (211)
14 : (41)
15 : (32)
16 : (1111)
18 : (221)
20 : (311)
22 : (51)
24 : (2111)
26 : (61)
28 : (411)
30 : (321)
32 : (11111)
33 : (52)
34 : (71)
35 : (43)
36 : (2211)
38 : (81)
40 : (3111)
(End)

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

Cf. A001222, A007359, A051424, A056239, A289506, A289507, A289508, A296150, A302696, A302697, A302698, A302796.

p:=1:for ind to 10000 do p:=nextprime(p);primeindex[p]:=ind;od:
out:=[]:for n from 2 to 100 do m:=[];f:=ifactors(n)[2];g:=0;
for k to nops(f) do mk:=primeindex[f[k][1]];m:=[op(m),mk];
g:=gcd(g,mk);od; if g=1 then out:=[op(out),n];fi;od:out;

Select[Range[200],GCD@@PrimePi/@FactorInteger[#][[All,1]]===1&] (* Gus Wiseman, Apr 13 2018 *)

isok(n) = my(f=factor(n)); gcd(apply(x->primepi(x), f[,1])) == 1; \\ Michel Marcus, Jul 19 2017

from sympy import gcd, primepi, primefactors
def ok(n): return gcd([primepi(p) for p in primefactors(n)]) == 1
print([n for n in range(1, 151) if ok(n)]) # Indranil Ghosh, Aug 06 2017

A072774 Powers of squarefree numbers.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 46, 47, 49, 51, 53, 55, 57, 58, 59, 61, 62, 64, 65, 66, 67, 69, 70, 71, 73, 74, 77, 78, 79, 81, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97

Offset: 1

Reinhard Zumkeller, Jul 10 2002

nonn

Essentially the same as A062770. - R. J. Mathar, Sep 25 2008
Numbers m such that in canonical prime factorization all prime exponents are identical: A124010(m,k) = A124010(m,1) for k = 2..A000005(m). - Reinhard Zumkeller, Apr 06 2014
Heinz numbers of uniform partitions. An integer partition is uniform if all parts appear with the same multiplicity. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Apr 16 2018

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Complement of A059404.
Cf. A072775, A072776, A072777 (subsequence), A005117, A072778, A124010, A329332 (tabular arrangement), A384667 (characteristic function).
A subsequence of A242414.
Cf. A000005, A000009, A000837, A007916, A013661, A047966, A052409, A052410, A072774, A078374, A289023, A289509, A300486, A302491, A302796, A302979.

import Data.Map (empty, findMin, deleteMin, insert)
import qualified Data.Map.Lazy as Map (null)
a072774 n = a072774_list !! (n-1)
(a072774_list, a072775_list, a072776_list) = unzip3 $
   (1, 1, 1) : f (tail a005117_list) empty where
   f vs'@(v:vs) m
    | Map.null m || xx > v = (v, v, 1) :
                             f vs (insert (v^2) (v, 2) m)
    | otherwise = (xx, bx, ex) :
                  f vs' (insert (bx*xx) (bx, ex+1) $ deleteMin m)
    where (xx, (bx, ex)) = findMin m
-- Reinhard Zumkeller, Apr 06 2014

isA := n -> n=1 or is(1 = nops({seq(p[2], p in ifactors(n)[2])})):
select(isA, [seq(1..97)]);  # Peter Luschny, Jun 10 2025

Select[Range[100], Length[Union[FactorInteger[#][[All, 2]]]] == 1 &] (* Geoffrey Critzer, Mar 30 2015 *)

is(n)=ispower(n,,&n); issquarefree(n) \\ Charles R Greathouse IV, Oct 16 2015

from math import isqrt
from sympy import mobius, integer_nthroot
def A072774(n):
    def g(x): return int(sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1)))-1
    def f(x): return n-2+x-sum(g(integer_nthroot(x,k)[0]) for k in range(1,x.bit_length()))
    kmin, kmax = 1,2
    while f(kmax) >= kmax:
        kmax <<= 1
    while True:
        kmid = kmax+kmin>>1
        if f(kmid) < kmid:
            kmax = kmid
        else:
            kmin = kmid
        if kmax-kmin <= 1:
            break
    return kmax # Chai Wah Wu, Aug 19 2024

a(n) = A072775(n)^A072776(n).
Sum_{n>=1} 1/a(n)^s = 1 + Sum_{k>=1} (zeta(k*s)/zeta(2*k*s)-1) for s > 1. - Amiram Eldar, Mar 20 2025
a(n)/n ~ Pi^2/6 (A013661). - Friedjof Tellkamp, Jun 09 2025

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

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Jul 18 2022

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

			The a(49) = 6 ways are: (1,2), (1,4), (2,1), (2,4), (4,1), (4,2).
The a(182) = 5 ways are: (1,2,3), (1,2,6), (1,4,2), (1,4,3), (1,4,6).
The a(546) = 2 ways are: (1,2,4,3), (1,2,4,6).

Wikipedia, Cartesian product.

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

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

A078374 Number of partitions of n into distinct and relatively prime parts.

Original entry on oeis.org

1, 0, 1, 1, 2, 2, 4, 4, 6, 7, 11, 10, 17, 17, 23, 26, 37, 36, 53, 53, 70, 77, 103, 103, 139, 147, 184, 199, 255, 260, 339, 358, 435, 474, 578, 611, 759, 810, 963, 1045, 1259, 1331, 1609, 1726, 2015, 2200, 2589, 2762, 3259, 3509, 4058, 4416, 5119, 5488, 6364, 6882

Offset: 1

Vladeta Jovovic, Dec 24 2002

nonn

The Heinz numbers of these partitions are given by A302796, which is the intersection of A005117 (strict) and A289509 (relatively prime). - Gus Wiseman, Oct 18 2020

			From _Gus Wiseman_, Oct 18 2020: (Start)
The a(1) = 1 through a(13) = 17 partitions (empty column indicated by dot, A = 10, B = 11, C = 12):
  1   .  21   31   32   51    43    53    54    73     65     75     76
                   41   321   52    71    72    91     74     B1     85
                              61    431   81    532    83     543    94
                              421   521   432   541    92     651    A3
                                          531   631    A1     732    B2
                                          621   721    542    741    C1
                                                4321   632    831    643
                                                       641    921    652
                                                       731    5421   742
                                                       821    6321   751
                                                       5321          832
                                                                     841
                                                                     931
                                                                     A21
                                                                     5431
                                                                     6421
                                                                     7321
(End)

Seiichi Manyama, Table of n, a(n) for n = 1..10000
N. J. A. Sloane, Transforms

Cf. A047966.
A000837 is the not necessarily strict version.
A302796 gives the Heinz numbers of these partitions.
A305713 is the pairwise coprime instead of relatively prime version.
A332004 is the ordered version.
A337452 is the case without 1's.
A000009 counts strict partitions.
A000740 counts relatively prime compositions.
Cf. A007359, A101268, A289508, A289509, A291166, A298748, A337451, A337485, A337451, A337561, A337563.

Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&GCD@@#==1&]],{n,15}] (* Gus Wiseman, Oct 18 2020 *)

Moebius transform of A000009.

G.f.: 1 + Sum_{n>=1} a(n)*x^n/(1 - x^n) = Product_{n>=1} (1 + x^n). - Ilya Gutkovskiy, Apr 26 2017

A304713 Squarefree numbers whose prime indices are pairwise indivisible. Heinz numbers of strict integer partitions with pairwise indivisible parts.

Original entry on oeis.org

1, 2, 3, 5, 7, 11, 13, 15, 17, 19, 23, 29, 31, 33, 35, 37, 41, 43, 47, 51, 53, 55, 59, 61, 67, 69, 71, 73, 77, 79, 83, 85, 89, 91, 93, 95, 97, 101, 103, 107, 109, 113, 119, 123, 127, 131, 137, 139, 141, 143, 145, 149, 151, 155, 157, 161, 163, 165, 167, 173

Offset: 1

Gus Wiseman, May 17 2018

nonn

A prime index of n is a number m such that prime(m) divides n. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

			Sequence of entries together with their corresponding multiset multisystems (see A302242) begins:
1:  {}
2:  {{}}
3:  {{1}}
5:  {{2}}
7:  {{1,1}}
11: {{3}}
13: {{1,2}}
15: {{1},{2}}
17: {{4}}
19: {{1,1,1}}
23: {{2,2}}
29: {{1,3}}
31: {{5}}
33: {{1},{3}}
35: {{2},{1,1}}

Cf. A000009, A005117, A006126, A056239, A073576, A285572, A285573, A293606, A293993, A302696, A302796, A303362, A303365, A304711.

Select[Range[300],SquareFreeQ[#]&&Select[Tuples[PrimePi/@First/@FactorInteger[#],2],UnsameQ@@#&&Divisible@@#&]==={}&]

A304711 Heinz numbers of integer partitions whose distinct parts are pairwise coprime.

Original entry on oeis.org

2, 4, 6, 8, 10, 12, 14, 15, 16, 18, 20, 22, 24, 26, 28, 30, 32, 33, 34, 35, 36, 38, 40, 44, 45, 46, 48, 50, 51, 52, 54, 55, 56, 58, 60, 62, 64, 66, 68, 69, 70, 72, 74, 75, 76, 77, 80, 82, 85, 86, 88, 90, 92, 93, 94, 95, 96, 98, 99, 100, 102, 104, 106, 108, 110

Offset: 1

Gus Wiseman, May 17 2018

nonn

Two parts are coprime if they have no common divisor greater than 1. For partitions of length 1 note that (1) is coprime but (x) is not coprime for x > 1.

First differs from A289509 at a(24) = 44, A289509(24) = 42.

			Sequence of all partitions whose distinct parts are pairwise coprime begins (1), (11), (21), (111), (31), (211), (41), (32), (1111), (221), (311), (51), (2111), (61), (411), (321), (11111), (52), (71), (43), (2211), (81), (3111).

Cf. A000837, A007359, A018783, A051424, A056239, A078374, A101268, A289508, A289509, A298748, A300486, A302569, A302696, A302698, A302796, A302797, A304709.

Select[Range[200],CoprimeQ@@PrimePi/@FactorInteger[#][[All,1]]&]

A355737 Number of ways to choose a sequence of divisors, one of each prime index of n (with multiplicity), such that the result has no common divisor > 1.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 2, 1, 3, 4, 1, 1, 4, 1, 2, 4, 2, 1, 2, 3, 4, 7, 3, 1, 4, 1, 1, 4, 2, 6, 4, 1, 4, 6, 2, 1, 6, 1, 2, 8, 3, 1, 2, 5, 4, 4, 4, 1, 8, 4, 3, 5, 4, 1, 4, 1, 2, 10, 1, 6, 4, 1, 2, 6, 6, 1, 4, 1, 6, 8, 4, 6, 8, 1, 2, 15, 2, 1, 6, 4, 4

Offset: 1

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

			The a(2) = 1 through a(18) = 4 choices:
  1  1  11  1  11  1  111  11  11  1  111  1  11  11  1111  1  111
               12          12  13     112     12  13           112
                           21                 14  21           121
                                                  23           122

Wikipedia, Coprime integers.

Dominated by A355731, firsts A355732, primes A355741, prime-powers A355742.
For weakly increasing instead of coprime we have A355735, primes A355745.
Positions of first appearances are A355738.
For strict instead of coprime we have A355739, zeros A355740.
A000005 counts divisors.
A001221 counts distinct prime factors, with sum A001414.
A001222 counts prime factors with multiplicity.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798.
A120383 lists numbers divisible by all of their prime indices.
A289508 gives GCD of prime indices.
A289509 ranks relatively prime partitions, odd A302697, squarefree A302796.
A324850 lists numbers divisible by the product of their prime indices.
Cf. A000720, A007359, A051424, A076610, A289507, A296150, A302696, A302698, A355733, A355744, A355748.

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

A302797 Squarefree numbers whose prime indices are pairwise coprime. Heinz numbers of strict integer partitions with pairwise coprime parts.

Original entry on oeis.org

1, 2, 6, 10, 14, 15, 22, 26, 30, 33, 34, 35, 38, 46, 51, 55, 58, 62, 66, 69, 70, 74, 77, 82, 85, 86, 93, 94, 95, 102, 106, 110, 118, 119, 122, 123, 134, 138, 141, 142, 143, 145, 146, 154, 155, 158, 161, 165, 166, 170, 177, 178, 186, 187, 190, 194, 201, 202

Offset: 1

Gus Wiseman, Apr 13 2018

nonn

A prime index of n is a number m such that prime(m) divides n. Two or more numbers are coprime if no pair of them has a common divisor other than 1. A single number is not considered coprime unless it is equal to 1.

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

			Sequence of terms together with their sets of prime indices begins:
01 : {}
02 : {1}
06 : {1,2}
10 : {1,3}
14 : {1,4}
15 : {2,3}
22 : {1,5}
26 : {1,6}
30 : {1,2,3}
33 : {2,5}
34 : {1,7}
35 : {3,4}
38 : {1,8}
46 : {1,9}
51 : {2,7}
55 : {3,5}
58 : {1,10}
62 : {1,11}
66 : {1,2,5}
69 : {2,9}
70 : {1,3,4}

Cf. A001222, A003963, A005117, A007359, A051424, A056239, A275024, A289509, A302242, A302505, A302696, A302697, A302698, A302796, A302798.

Select[Range[100],Or[#===1,SquareFreeQ[#]&&CoprimeQ@@PrimePi/@FactorInteger[#][[All,1]]]&]

A101271 Number of partitions of n into 3 distinct and relatively prime parts.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 6, 8, 9, 12, 12, 16, 15, 21, 20, 26, 25, 33, 28, 40, 36, 45, 42, 56, 44, 65, 56, 70, 64, 84, 66, 96, 81, 100, 88, 120, 90, 133, 110, 132, 121, 161, 120, 175, 140, 176, 156, 208, 153, 220, 180, 222, 196, 261, 184, 280, 225, 270, 240, 312, 230, 341, 272

Offset: 6

Vladeta Jovovic, Dec 19 2004

The Heinz numbers of these partitions are the intersection of A289509 (relatively prime), A005117 (strict), and A014612 (triple). - Gus Wiseman, Oct 15 2020

			For n=10 we have 4 such partitions: 1+2+7, 1+3+6, 1+4+5 and 2+3+5.
From _Gus Wiseman_, Oct 13 2020: (Start)
The a(6) = 1 through a(18) = 15 triples (A..F = 10..15):
  321  421  431  432  532  542  543  643  653  654  754  764  765
            521  531  541  632  651  652  743  753  763  854  873
                 621  631  641  732  742  752  762  853  863  954
                      721  731  741  751  761  843  871  872  972
                           821  831  832  851  852  943  953  981
                                921  841  932  861  952  962  A53
                                     931  941  942  961  971  A71
                                     A21  A31  951  A51  A43  B43
                                          B21  A32  B32  A52  B52
                                               A41  B41  A61  B61
                                               B31  C31  B42  C51
                                               C21  D21  B51  D32
                                                         C32  D41
                                                         C41  E31
                                                         D31  F21
                                                         E21
(End)

Fausto A. C. Cariboni, Table of n, a(n) for n = 6..10000

Cf. A023024-A023030, A000742-A000743, A023031-A023035.
A000741 is the ordered non-strict version.
A001399(n-6) does not require relative primality.
A023022 counts pairs instead of triples.
A023023 is the not necessarily strict version.
A078374 counts these partitions of any length, with Heinz numbers A302796.
A101271*6 is the ordered version.
A220377 is the pairwise coprime instead of relatively prime version.
A284825 counts the case that is pairwise non-coprime also.
A337605 is the pairwise non-coprime instead of relatively prime version.
A008289 counts strict partitions by sum and length.
A007304 gives the Heinz numbers of 3-part strict partitions.
A307719 counts 3-part pairwise coprime partitions.
A337601 counts 3-part partitions whose distinct parts are pairwise coprime.
Cf. A000010, A000217, A000837, A007360, A014612, A055684, A289509, A332004, A337452, A337563.

m:=3: with(numtheory): g:=sum(mobius(k)*x^(m*(m+1)/2*k)/Product(1-x^(i*k),i=1..m),k=1..20): gser:=series(g,x=0,80): seq(coeff(gser,x^n),n=6..77); # Emeric Deutsch, May 31 2005

Table[Length[Select[IntegerPartitions[n,{3}],UnsameQ@@#&&GCD@@#==1&]],{n,6,50}] (* Gus Wiseman, Oct 13 2020 *)

G.f. for the number of partitions of n into m distinct and relatively prime parts is Sum(moebius(k)*x^(m*(m+1)/2*k)/Product(1-x^(i*k), i=1..m), k=1..infinity).

More terms from Emeric Deutsch, May 31 2005

A316470 Matula-Goebel numbers of unlabeled rooted RPMG-trees, meaning the Matula-Goebel numbers of the branches of any non-leaf node are relatively prime.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 14, 16, 18, 24, 26, 28, 32, 36, 38, 42, 48, 52, 54, 56, 64, 72, 74, 76, 78, 84, 86, 96, 98, 104, 106, 108, 112, 114, 122, 126, 128, 144, 148, 152, 156, 162, 168, 172, 178, 182, 192, 196, 202, 208, 212, 214, 216, 222, 224, 228, 234, 244, 252

Offset: 1

Gus Wiseman, Jul 04 2018

nonn

A prime index of n is a number m such that prime(m) divides n. A number is in the sequence iff it is 1 or its prime indices are relatively prime and already belong to the sequence.

			The sequence of all RPMG-trees preceded by their Matula-Goebel numbers begins:
   1: o
   2: (o)
   4: (oo)
   6: (o(o))
   8: (ooo)
  12: (oo(o))
  14: (o(oo))
  16: (oooo)
  18: (o(o)(o))
  24: (ooo(o))
  26: (o(o(o)))
  28: (oo(oo))
  32: (ooooo)
  36: (oo(o)(o))
  38: (o(ooo))
  42: (o(o)(oo))

Cf. A000081, A000837, A007097, A289509, A302796, A316468, A316469, A316473, A316475, A316495.

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[1000],Or[#==1,And[GCD@@primeMS[#]==1,And@@#0/@primeMS[#]]]&]

cp's OEIS Frontend

A289509 Numbers k such that the gcd of the indices j for which the j-th prime prime(j) divides k is 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

A072774 Powers of squarefree numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A078374 Number of partitions of n into distinct and relatively prime parts.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A304713 Squarefree numbers whose prime indices are pairwise indivisible. Heinz numbers of strict integer partitions with pairwise indivisible parts.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A304711 Heinz numbers of integer partitions whose distinct parts are pairwise coprime.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A355737 Number of ways to choose a sequence of divisors, one of each prime index of n (with multiplicity), such that the result has no common divisor > 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A302797 Squarefree numbers whose prime indices are pairwise coprime. Heinz numbers of strict integer partitions with pairwise coprime parts.

Views

Author