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 11-20 of 35 results. Next

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

A303140 Number of strict integer partitions of n with at least two but not all parts having a common divisor greater than 1.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 0, 2, 2, 4, 2, 8, 7, 14, 14, 21, 18, 33, 32, 50, 54, 72, 67, 103, 110, 145, 155, 201, 196, 271, 293, 372, 400, 493, 512, 647, 704, 858, 924, 1115, 1167, 1436, 1560, 1854, 2022, 2368, 2510, 3005, 3255, 3804, 4144, 4792, 5116, 5989, 6514, 7486
Offset: 1

Views

Author

Gus Wiseman, Apr 19 2018

Keywords

Examples

			The a(14) = 7 partitions are (932), (8321), (7421), (653), (6521), (6431), (5432).
		

Crossrefs

Programs

  • Mathematica
    Table[Select[IntegerPartitions[n],UnsameQ@@#&&!CoprimeQ@@#&&GCD@@#===1&]//Length,{n,20}]

A355535 Odd numbers of which it is not possible to choose a different prime factor of each prime index.

Original entry on oeis.org

9, 21, 25, 27, 45, 49, 57, 63, 75, 81, 99, 105, 115, 117, 121, 125, 133, 135, 147, 153, 159, 171, 175, 189, 195, 207, 225, 231, 243, 245, 261, 273, 275, 279, 285, 289, 297, 315, 325, 333, 343, 345, 351, 357, 361, 363, 369, 371, 375, 387, 393, 399, 405, 423
Offset: 1

Views

Author

Gus Wiseman, Jul 22 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 begin:
    9: {2,2}
   21: {2,4}
   25: {3,3}
   27: {2,2,2}
   45: {2,2,3}
   49: {4,4}
   57: {2,8}
   63: {2,2,4}
   75: {2,3,3}
   81: {2,2,2,2}
   99: {2,2,5}
  105: {2,3,4}
For example, the prime indices of 897 are {2,6,9}, of which we can choose prime factors in two ways: (2,2,3) or (2,3,3); but neither of these has all distinct elements, so 897 is in the sequence.
		

Crossrefs

Including evens gives A355529.
The version for all divisors including evens is A355740, zeros of A355739.
Choices of a prime factor of each prime index: A355741, unordered A355744.
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.

Programs

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

A303282 Numbers whose prime indices have no common divisor other than 1 but are not pairwise coprime.

Original entry on oeis.org

18, 36, 42, 45, 50, 54, 72, 75, 78, 84, 90, 98, 99, 100, 105, 108, 114, 126, 130, 135, 144, 150, 153, 156, 162, 168, 174, 175, 180, 182, 195, 196, 198, 200, 207, 210, 216, 222, 225, 228, 230, 231, 234, 242, 245, 250, 252, 258, 260, 266, 270, 275, 279, 285, 288
Offset: 1

Views

Author

Gus Wiseman, Apr 20 2018

Keywords

Comments

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.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The sequence of integer partitions whose Heinz numbers belong to this sequence begins (221), (2211), (421), (322), (331), (2221), (22111), (332), (621), (4211), (3221), (441), (522), (3311), (432), (22211).
		

Crossrefs

Programs

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

A304712 Number of integer partitions of n whose parts are all equal or whose distinct parts are pairwise coprime.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 10, 14, 19, 25, 32, 43, 54, 70, 86, 105, 130, 162, 196, 240, 286, 339, 405, 485, 573, 674, 790, 922, 1072, 1252, 1456, 1685, 1939, 2226, 2557, 2923, 3349, 3822, 4347, 4931, 5593, 6335, 7170, 8092, 9105, 10233, 11495, 12903, 14458, 16169, 18063
Offset: 0

Views

Author

Gus Wiseman, May 17 2018

Keywords

Comments

Two parts are coprime if they have no common divisor greater than 1.

Examples

			The a(6) = 10 partitions whose parts are all equal or whose distinct parts are pairwise coprime are (6), (51), (411), (33), (321), (3111), (222), (2211), (21111), (111111).
		

Crossrefs

Programs

  • Maple
    g:= proc(n, i, s) `if`(n=0, 1, `if`(i<1, 0,
          b(n, i, select(x-> x<=i, s))))
        end:
    b:= proc(n, i, s) option remember; g(n, i-1, s)+(f->
         `if`(f intersect s={}, add(g(n-i*j, i-1, s union f)
            , j=1..n/i), 0))(numtheory[factorset](i))
        end:
    a:= n-> g(n$2, {}):
    seq(a(n), n=0..60);  # Alois P. Heinz, May 17 2018
  • Mathematica
    Table[Select[IntegerPartitions[n],Or[SameQ@@#,CoprimeQ@@Union[#]]&]//Length,{n,20}]
    (* Second program: *)
    g[n_, i_, s_] := If[n == 0, 1, If[i < 1, 0, b[n, i, Select[s, # <= i &]]]];
    b[n_, i_, s_] := b[n, i, s] = g[n, i - 1, s] + Function[f,
         If[f ~Intersection~ s == {}, Sum[g[n - i*j, i - 1, s ~Union~ f],
         {j, 1, n/i}], 0]][FactorInteger[i][[All, 1]]];
    a[n_] := g[n, n, {}];
    a /@ Range[0, 60] (* Jean-François Alcover, May 10 2021, after Alois P. Heinz *)

A303139 Number of integer partitions of n with at least two but not all parts having a common divisor greater than 1.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 5, 6, 13, 17, 33, 37, 68, 82, 125, 159, 237, 278, 409, 491, 674, 830, 1121, 1329, 1781, 2144, 2770, 3345, 4299, 5086, 6507, 7752, 9687, 11571, 14378, 16985, 21039, 24876, 30379, 35924, 43734, 51320, 62238, 73068, 87747, 103021, 123347, 143955
Offset: 1

Views

Author

Gus Wiseman, Apr 19 2018

Keywords

Examples

			The a(7) = 5 partitions are (421), (331), (322), (2221), (22111).
		

Crossrefs

Programs

  • Mathematica
    Table[Select[IntegerPartitions[n],!CoprimeQ@@#&&GCD@@#===1&]//Length,{n,30}]

A303280 Number of strict integer partitions of n whose parts have a common divisor other than 1.

Original entry on oeis.org

0, 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, 1, 32, 13, 38, 7, 57, 1, 54, 19, 68, 1, 95, 1, 90, 33, 104, 1, 148, 5, 149, 39, 166, 1, 230, 14, 226, 55, 256, 1, 360, 1, 340, 82, 390, 20, 527, 1, 513, 105, 609, 1
Offset: 1

Views

Author

Gus Wiseman, Apr 20 2018

Keywords

Examples

			The a(18) = 10 strict partitions are (18), (10,8), (12,6), (14,4), (15,3), (16,2), (8,6,4), (9,6,3), (10,6,2), (12,4,2).
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    b:= proc(n) option remember; `if`(n=0, 1, add(add(
          `if`(d::odd, d, 0), d=divisors(j))*b(n-j), j=1..n)/n)
        end:
    a:= n-> -add(mobius(d)*b(n/d), d=divisors(n) minus {1}):
    seq(a(n), n=1..100);  # Alois P. Heinz, Apr 23 2018
  • Mathematica
    Table[-Sum[MoebiusMu[d]*PartitionsQ[n/d],{d,Rest[Divisors[n]]}],{n,100}]

Formula

a(n) = -Sum_{d|n, d > 1} mu(d) * A000009(n/d).

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

Original entry on oeis.org

2, 17719, 40807, 43381, 50431, 74269, 83143, 101543, 105703, 116143, 121307, 123469, 139919, 140699, 142883, 171613, 181831, 185803, 191479, 203557, 205813, 211381, 213239, 215267, 219271, 246703, 249587, 249899, 279371, 286897, 289007, 296993, 300847, 303949
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).

Examples

			The sequence of strict integer partitions with Heinz numbers in the sequence begins: (1), (15,10,6), (21,14,6), (20,15,6), (15,12,10), (45,10,6), (18,15,10).
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100000],With[{m=PrimePi/@FactorInteger[#][[All,1]]},And[SquareFreeQ[#],GCD@@m==1,And@@(GCD[##]>1&)@@@Select[Tuples[m,2],Less@@#&]]]&]

A303283 Squarefree numbers whose prime indices have no common divisor other than 1 but are not pairwise coprime.

Original entry on oeis.org

42, 78, 105, 114, 130, 174, 182, 195, 210, 222, 230, 231, 258, 266, 285, 318, 345, 357, 366, 370, 390, 406, 426, 429, 435, 455, 462, 470, 474, 483, 494, 518, 534, 546, 555, 570, 598, 602, 606, 610, 627, 638, 642, 645, 651, 663, 665, 678, 690, 705, 714, 715
Offset: 1

Views

Author

Gus Wiseman, Apr 20 2018

Keywords

Comments

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.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The sequence of strict integer partitions whose Heinz numbers belong to this sequence begins (4,2,1), (6,2,1), (4,3,2), (8,2,1), (6,3,1), (10,2,1), (6,4,1), (6,3,2), (4,3,2,1), (12,2,1), (9,3,1), (5,4,2), (14,2,1), (8,4,1), (8,3,2), (16,2,1), (9,3,2), (7,4,2), (18,2,1), (12,3,1), (6,3,2,1).
		

Crossrefs

Programs

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

A316467 Matula-Goebel numbers of locally stable rooted identity trees, meaning no branch is a subset of any other branch of the same root.

Original entry on oeis.org

1, 2, 3, 5, 11, 15, 31, 33, 47, 55, 93, 127, 137, 141, 155, 165, 211, 257, 341, 381, 411, 465, 487, 633, 635, 709, 771, 773, 811, 907, 977, 1023, 1055, 1285, 1297, 1397, 1457, 1461, 1507, 1621, 1705, 1905, 2127, 2293, 2319, 2321, 2433, 2621, 2721, 2833, 2931
Offset: 1

Views

Author

Gus Wiseman, Jul 04 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. A number belongs to this sequence iff it is squarefree, its distinct prime indices are pairwise indivisible, and its prime indices also belong to this sequence.

Examples

			165 = prime(2)*prime(3)*prime(5) belongs to the sequence because it is squarefree, the indices {2,3,5} are pairwise indivisible, and each of them already belongs to the sequence.
Sequence of locally stable rooted identity trees preceded by their Matula-Goebel numbers begins:
    1: o
    2: (o)
    3: ((o))
    5: (((o)))
   11: ((((o))))
   15: ((o)((o)))
   31: (((((o)))))
   33: ((o)(((o))))
   47: (((o)((o))))
   55: (((o))(((o))))
   93: ((o)((((o)))))
  127: ((((((o))))))
  137: (((o)(((o)))))
  141: ((o)((o)((o))))
  155: (((o))((((o)))))
  165: ((o)((o))(((o))))
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    ain[n_]:=And[Select[Tuples[primeMS[n],2],UnsameQ@@#&&Divisible@@#&]=={},SquareFreeQ[n],And@@ain/@primeMS[n]];
    Select[Range[100],ain]
Previous Showing 11-20 of 35 results. Next