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-7 of 7 results.

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}]

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&]

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

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&]

A319301 Sum of GCDs of strict integer partitions of n.

Original entry on oeis.org

1, 2, 4, 5, 7, 10, 11, 14, 18, 21, 22, 33, 30, 39, 49, 54, 54, 78, 72, 100, 110, 121, 126, 181, 174, 207, 238, 284, 284, 389, 370, 466, 512, 582, 647, 806, 796, 954, 1066, 1265, 1300, 1616, 1652, 1979, 2192, 2452, 2636, 3202, 3336, 3892, 4237, 4843, 5172, 6090
Offset: 1

Views

Author

Gus Wiseman, Sep 16 2018

Keywords

Crossrefs

Programs

  • Maple
    b:= proc(n, i, r) option remember; `if`(i*(i+1)/2 `if`(i b(n$2, 0):
    seq(a(n), n=1..61);  # Alois P. Heinz, Mar 17 2019
  • Mathematica
    Table[Sum[GCD@@ptn,{ptn,Select[IntegerPartitions[n],UnsameQ@@#&]}],{n,30}]
    (* Second program: *)
    b[n_, i_, r_] := b[n, i, r] = If[i(i+1)/2 < n, 0,
         With[{t = GCD[i, r]}, If[i < n, b[n - i, Min[i - 1, n - i], t], 0] +
         If[i == n, t, 0] + b[n, i - 1, r]]];
    a[n_] := b[n, n, 0];
    Array[a, 61] (* Jean-François Alcover, May 20 2021, after Alois P. Heinz *)

Formula

From Richard L. Ollerton, May 06 2021: (Start)
a(n) = Sum_{d|n} A000010(n/d)*A000009(d).
a(n) = Sum_{k=1..n} A000009(gcd(n,k)).
a(n) = Sum_{k=1..n} A000009(n/gcd(n,k))*A000010(gcd(n,k))/A000010(n/gcd(n,k)). (End)

A319300 Irregular triangle where T(n,k) is the number of strict integer partitions of n with GCD equal to the k-th divisor of n.

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 0, 1, 2, 1, 2, 1, 0, 1, 4, 1, 4, 1, 0, 1, 6, 1, 1, 7, 2, 0, 1, 11, 1, 10, 2, 1, 1, 0, 1, 17, 1, 17, 4, 0, 1, 23, 2, 1, 1, 26, 4, 1, 0, 1, 37, 1, 36, 6, 2, 1, 0, 1, 53, 1, 53, 7, 2, 1, 0, 1, 70, 4, 1, 1, 77, 11, 0, 1, 103, 1, 103, 10, 4, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Sep 16 2018

Keywords

Examples

			Triangle begins:
   1
   0  1
   1  1
   1  0  1
   2  1
   2  1  0  1
   4  1
   4  1  0  1
   6  1  1
   7  2  0  1
  11  1
  10  2  1  1  0  1
  17  1
  17  4  0  1
  23  2  1  1
  26  4  1  0  1
  37  1
  36  6  2  1  0  1
  53  1
  53  7  2  1  0  1
		

Crossrefs

A regular version is A303138. Row lengths are A000005. Row sums are A000009. First column is A078374.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],And[UnsameQ@@#,GCD@@#==k]&]],{n,20},{k,Divisors[n]}]

Formula

T(n,k) = A078374(n/A027750(n,k)).
Showing 1-7 of 7 results.