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.

Showing 1-10 of 23 results. Next

A051424 Number of partitions of n into pairwise relatively prime parts.

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 7, 10, 12, 15, 18, 23, 27, 33, 38, 43, 51, 60, 70, 81, 92, 102, 116, 134, 153, 171, 191, 211, 236, 266, 301, 335, 367, 399, 442, 485, 542, 598, 649, 704, 771, 849, 936, 1023, 1103, 1185, 1282, 1407, 1535, 1662, 1790, 1917, 2063, 2245, 2436
Offset: 0

Views

Author

Keywords

Examples

			a(4) = 4 since all partitions of 4 consist of relatively prime numbers except 2+2.
The a(6) = 7 partitions with pairwise coprime parts: (111111), (21111), (3111), (321), (411), (51), (6). - _Gus Wiseman_, Apr 14 2018
		

Crossrefs

Number of partitions of n into relatively prime parts = A000837.
Row sums of A282749.

Programs

  • Haskell
    a051424 = length . filter f . partitions where
       f [] = True
       f (p:ps) = (all (== 1) $ map (gcd p) ps) && f ps
       partitions n = ps 1 n where
         ps x 0 = [[]]
         ps x y = [t:ts | t <- [x..y], ts <- ps t (y - t)]
    -- Reinhard Zumkeller, Dec 16 2013
  • Maple
    with(numtheory):
    b:= proc(n, i, s) option remember; local f;
          if n=0 or i=1 then 1
        elif i<2 then 0
        else f:= factorset(i);
             b(n, i-1, select(x->is(xis(x b(n, n, {}):
    seq(a(n), n=0..80);  # Alois P. Heinz, Mar 14 2012
  • Mathematica
    b[n_, i_, s_] := b[n, i, s] = Module[{f}, If[n == 0 || i == 1, 1, If[i < 2, 0, f = FactorInteger[i][[All, 1]]; b[n, i-1, Select[s, # < i &]] + If[i <= n && f ~Intersection~ s == {}, b[n-i, i-1, Select[s ~Union~ f, # < i &]], 0]]]]; a[n_] := b[n, n, {}]; Table[a[n], {n, 0, 54}] (* Jean-François Alcover, Oct 03 2013, translated from Maple, after Alois P. Heinz *)

Formula

log a(n) ~ (2*Pi/sqrt(6)) sqrt(n/log n). - Eric M. Schmidt, Jul 04 2013
Apparently no formula or recurrence is known. - N. J. A. Sloane, Mar 05 2017

Extensions

More precise definition from Vladeta Jovovic, Dec 11 2004

A007359 Number of partitions of n into pairwise coprime parts that are >= 2.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 1, 3, 2, 3, 3, 5, 4, 6, 5, 5, 8, 9, 10, 11, 11, 10, 14, 18, 19, 18, 20, 20, 25, 30, 35, 34, 32, 32, 43, 43, 57, 56, 51, 55, 67, 78, 87, 87, 80, 82, 97, 125, 128, 127, 128, 127, 146, 182, 191, 185, 184, 193, 213, 263, 290, 279, 258, 271, 312, 354, 404, 402
Offset: 0

Views

Author

N. J. A. Sloane and Mira Bernstein, following a suggestion from Marc LeBrun, Apr 28 1994

Keywords

Comments

This sequence is of interest for group theory. The partitions counted by a(n) correspond to conjugacy classes of optimal order of the symmetric group of n elements: they have no fixed point, their order is the direct product of their cycle lengths and they are not contained in a subgroup of Sym_p for p < n. A123131 gives the maximum order (LCM) reachable by these partitions.

Examples

			The a(17) = 9 strict partitions into pairwise coprime parts that are greater than 1 are (17), (15,2), (14,3), (13,4), (12,5), (11,6), (10,7), (9,8), (7,5,3,2). - _Gus Wiseman_, Apr 14 2018
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Maple
    with(numtheory):
    b:= proc(n, i, s) option remember; local f;
          if n=0 then 1
        elif i<2 then 0
        else f:= factorset(i);
             b(n, i-1, select(x-> is(x is(x b(n, n, {}):
    seq(a(n), n=0..80);  # Alois P. Heinz, Mar 14 2012
  • Mathematica
    b[n_, i_, s_] := b[n, i, s] = Module[{f}, If[n == 0 || i == 1, 1, If[i<2, 0, f = FactorInteger[i][[All, 1]]; b[n, i-1, Select[s, #Jean-François Alcover, Feb 17 2014, after Alois P. Heinz *)
    Table[Length[Select[IntegerPartitions[n],FreeQ[#,1]&&(Length[#]===1||CoprimeQ@@#)&]],{n,20}] (* Gus Wiseman, Apr 14 2018 *)

Formula

a(n) = A051424(n) - A051424(n-1). - Vladeta Jovovic, Dec 11 2004

Extensions

More precise definition from Vladeta Jovovic, Dec 11 2004
More terms from Pab Ter (pabrlos2(AT)yahoo.com), Nov 13 2005

A305713 Number of strict integer partitions of n into pairwise coprime parts.

Original entry on oeis.org

1, 0, 1, 1, 2, 2, 3, 4, 4, 5, 7, 8, 9, 10, 9, 12, 16, 18, 20, 21, 20, 23, 31, 36, 36, 37, 39, 44, 54, 64, 68, 65, 63, 74, 85, 99, 112, 106, 105, 121, 144, 164, 173, 166, 161, 178, 221, 252, 254, 254, 254, 272, 327, 372, 375, 368, 376, 405, 475, 552, 568, 536
Offset: 1

Views

Author

Gus Wiseman, Sep 01 2018

Keywords

Examples

			The a(13) = 9 strict partitions are (7,6), (8,5), (9,4), (10,3), (11,2), (12,1), (7,5,1), (5,4,3,1), (7,3,2,1).
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&CoprimeQ@@#&]],{n,30}]

A220377 Number of partitions of n into three distinct and mutually relatively prime parts.

Original entry on oeis.org

1, 0, 2, 1, 3, 1, 6, 1, 7, 3, 7, 3, 14, 3, 15, 6, 14, 6, 25, 6, 22, 10, 25, 9, 42, 8, 34, 15, 37, 15, 53, 13, 48, 22, 53, 17, 78, 17, 65, 30, 63, 24, 99, 24, 88, 35, 84, 30, 126, 34, 103, 45, 103, 38, 166, 35, 124, 57, 128, 51, 184, 44, 150, 67, 172, 52, 218
Offset: 6

Views

Author

Carl Najafi, Dec 13 2012

Keywords

Comments

The Heinz numbers of these partitions are the intersection of A005117 (strict), A014612 (triples), and A302696 (coprime). - Gus Wiseman, Oct 14 2020

Examples

			For n=10 we have three such partitions: 1+2+7, 1+4+5 and 2+3+5.
From _Gus Wiseman_, Oct 14 2020: (Start)
The a(6) = 1 through a(20) = 15 triples (empty column indicated by dot, A..H = 10..17):
321  .  431  531  532  731  543  751  743  753  754  971  765  B53  875
        521       541       651       752  951  853  B51  873  B71  974
                  721       732       761  B31  871  D31  954  D51  A73
                            741       851       952       972       A91
                            831       941       B32       981       B54
                            921       A31       B41       A71       B72
                                      B21       D21       B43       B81
                                                          B52       C71
                                                          B61       D43
                                                          C51       D52
                                                          D32       D61
                                                          D41       E51
                                                          E31       F41
                                                          F21       G31
                                                                    H21
(End)
		

Crossrefs

A023022 is the 2-part version.
A101271 is the relative prime instead of pairwise coprime version.
A220377*6 is the ordered version.
A305713 counts these partitions of any length, with Heinz numbers A302797.
A307719 is the non-strict version.
A337461 is the non-strict ordered version.
A337563 is the case with no 1's.
A337605 is the pairwise non-coprime instead of pairwise coprime version.
A001399(n-6) counts strict 3-part partitions, with Heinz numbers A007304.
A008284 counts partitions by sum and length, with strict case A008289.
A318717 counts pairwise non-coprime strict partitions.
A326675 ranks pairwise coprime sets.
A327516 counts pairwise coprime partitions.
A337601 counts 3-part partitions whose distinct parts are pairwise coprime.

Programs

  • Mathematica
    Table[Length@Select[ IntegerPartitions[ n, {3}], #[[1]] != #[[2]] != #[[3]] && GCD[#[[1]], #[[2]]] == 1 && GCD[#[[1]], #[[3]]] == 1 && GCD[#[[2]], #[[3]]] == 1 &], {n, 6, 100}]
    Table[Count[IntegerPartitions[n,{3}],?(CoprimeQ@@#&&Length[ Union[#]] == 3&)],{n,6,100}] (* _Harvey P. Dale, May 22 2020 *)
  • PARI
    a(n)=my(P=partitions(n));sum(i=1,#P,#P[i]==3&&P[i][1]Charles R Greathouse IV, Dec 14 2012

Formula

a(n > 2) = A307719(n) - 1. - Gus Wiseman, Oct 15 2020

A302796 Squarefree numbers whose prime indices are relatively prime. Nonprime Heinz numbers of strict integer partitions with relatively prime parts.

Original entry on oeis.org

1, 2, 6, 10, 14, 15, 22, 26, 30, 33, 34, 35, 38, 42, 46, 51, 55, 58, 62, 66, 69, 70, 74, 77, 78, 82, 85, 86, 93, 94, 95, 102, 105, 106, 110, 114, 118, 119, 122, 123, 130, 134, 138, 141, 142, 143, 145, 146, 154, 155, 158, 161, 165, 166, 170, 174, 177, 178, 182
Offset: 1

Views

Author

Gus Wiseman, Apr 13 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. Two or more numbers are relatively prime if they have no common divisor other than 1. A single number is not considered relatively prime unless it is equal to 1.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			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}
42 : {1,2,4}
46 : {1,9}
51 : {2,7}
55 : {3,5}
58 : {1,10}
62 : {1,11}
66 : {1,2,5}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],Or[#===1,SquareFreeQ[#]&&GCD@@PrimePi/@FactorInteger[#][[All,1]]===1]&]
  • PARI
    isok(n) = {if (n == 1, return (1)); if (issquarefree(n), my(f = factor(n)); return (gcd(vector(#f~, k, primepi(f[k,1]))) == 1););} \\ Michel Marcus, Apr 13 2018

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

Views

Author

Gus Wiseman, May 17 2018

Keywords

Comments

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.

Examples

			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).
		

Crossrefs

Programs

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

A318717 Number of strict integer partitions of n in which no two parts are relatively prime.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 2, 2, 3, 1, 5, 1, 5, 4, 6, 1, 10, 1, 11, 6, 12, 1, 19, 3, 18, 8, 23, 1, 36, 2, 32, 13, 38, 7, 57, 2, 54, 19, 68, 3, 95, 3, 90, 33, 104, 3, 148, 7, 149, 40, 166, 5, 230, 17, 226, 56, 256, 6, 360, 9, 340, 84, 390, 25, 527, 11, 513, 109
Offset: 0

Views

Author

Gus Wiseman, Sep 02 2018

Keywords

Examples

			The a(20) = 11 partitions:
  (20),
  (12,8), (14,6), (15,5), (16,4), (18,2),
  (10,6,4), (10,8,2), (12,6,2), (14,4,2),
  (8,6,4,2).
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],And[UnsameQ@@#,And@@(GCD[##]>1&)@@@Select[Tuples[#,2],Less@@#&]]&]],{n,30}]

Extensions

a(51)-a(69) from Alois P. Heinz, Sep 02 2018

A304709 Number of integer partitions of n whose distinct parts are pairwise coprime.

Original entry on oeis.org

1, 1, 2, 3, 6, 7, 13, 16, 23, 29, 42, 49, 69, 83, 102, 126, 161, 191, 239, 281, 336, 402, 484, 566, 672, 787, 919, 1067, 1251, 1449, 1684, 1934, 2223, 2554, 2920, 3341, 3821, 4344, 4928, 5586, 6334, 7163, 8091, 9100, 10228, 11492, 12902, 14449, 16167, 18058
Offset: 1

Views

Author

Gus Wiseman, May 17 2018

Keywords

Comments

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.

Examples

			The a(6) = 7 integer partitions of 6 whose distinct parts are pairwise coprime are (51), (411), (321), (3111), (2211), (21111), (111111).
		

Crossrefs

Programs

  • Mathematica
    Table[Select[IntegerPartitions[n],CoprimeQ@@Union[#]&]//Length,{n,20}]
  • PARI
    lista(nn)={local(Cache=Map());
      my(excl=vector(nn, n, sum(i=1, n-1, if(gcd(i,n)>1, 2^(n-i)))));
      my(c(n, m, b)=
         if(n==0, 1,
            while(m>n || bittest(b,0), m--; b>>=1);
            my(hk=[n, m, b], z);
            if(!mapisdefined(Cache, hk, &z),
              z = if(m, self()(n, m-1, b>>1) + self()(n-m, m, bitor(b, excl[m])), 0);
              mapput(Cache, hk, z)); z));
      my(a(n)=c(n, n, 0) + 1 - numdiv(n));
      for(n=1, nn, print1(a(n), ", "))
    } \\ Andrew Howroyd, Nov 02 2019

Formula

a(n) = A304712(n) + 1 - A000005(n). - Andrew Howroyd, Nov 02 2019

A318719 Heinz numbers of strict integer partitions in which no two parts are relatively prime.

Original entry on oeis.org

1, 2, 3, 5, 7, 11, 13, 17, 19, 21, 23, 29, 31, 37, 39, 41, 43, 47, 53, 57, 59, 61, 65, 67, 71, 73, 79, 83, 87, 89, 91, 97, 101, 103, 107, 109, 111, 113, 115, 127, 129, 131, 133, 137, 139, 149, 151, 157, 159, 163, 167, 173, 179, 181, 183, 185, 191, 193, 197
Offset: 1

Views

Author

Gus Wiseman, Sep 02 2018

Keywords

Comments

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

Crossrefs

Programs

  • Mathematica
    Select[Range[200],And[SquareFreeQ[#],And@@(GCD[##]>1&)@@@Select[Tuples[PrimePi/@FactorInteger[#][[All,1]],2],Less@@#&]]&]

A302568 Odd numbers that are either prime or whose prime indices are pairwise coprime.

Original entry on oeis.org

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, 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, 177, 179
Offset: 1

Views

Author

Gus Wiseman, Apr 10 2018

Keywords

Comments

Also Heinz numbers of partitions with pairwise coprime parts all greater than 1 (A007359), where singletons are considered coprime. The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

Examples

			The sequence of terms together with their prime indices begins:
      3: {2}       43: {14}      89: {24}      141: {2,15}
      5: {3}       47: {15}      93: {2,11}    143: {5,6}
      7: {4}       51: {2,7}     95: {3,8}     145: {3,10}
     11: {5}       53: {16}      97: {25}      149: {35}
     13: {6}       55: {3,5}    101: {26}      151: {36}
     15: {2,3}     59: {17}     103: {27}      155: {3,11}
     17: {7}       61: {18}     107: {28}      157: {37}
     19: {8}       67: {19}     109: {29}      161: {4,9}
     23: {9}       69: {2,9}    113: {30}      163: {38}
     29: {10}      71: {20}     119: {4,7}     165: {2,3,5}
     31: {11}      73: {21}     123: {2,13}    167: {39}
     33: {2,5}     77: {4,5}    127: {31}      173: {40}
     35: {3,4}     79: {22}     131: {32}      177: {2,17}
     37: {12}      83: {23}     137: {33}      179: {41}
     41: {13}      85: {3,7}    139: {34}      181: {42}
Entry A302242 describes a correspondence between positive integers and multiset multisystems. In this case it gives the following sequence of multiset systems.
03: {{1}}
05: {{2}}
07: {{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}}
37: {{1,1,2}}
41: {{6}}
43: {{1,4}}
47: {{2,3}}
51: {{1},{4}}
53: {{1,1,1,1}}
		

Crossrefs

A005117 is a superset.
A007359 counts partitions with these Heinz numbers.
A302569 allows evens, with squarefree version A302798.
A337694 is the pairwise non-coprime instead of pairwise coprime version.
A337984 does not include the primes.
A305713 counts pairwise coprime strict partitions.
A327516 counts pairwise coprime partitions, ranked by A302696.
A337462 counts pairwise coprime compositions, ranked by A333227.
A337561 counts pairwise coprime strict compositions.
A337667 counts pairwise non-coprime compositions, ranked by A337666.
A337697 counts pairwise coprime compositions with no 1's.

Programs

  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[1,400,2],Or[PrimeQ[#],CoprimeQ@@primeMS[#]]&]

Formula

Equals A065091 \/ A337984.
Equals A302569 /\ A005408.

Extensions

Extended by Gus Wiseman, Oct 29 2020
Showing 1-10 of 23 results. Next