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 33 results. Next

A323093 Number of integer partitions of n where no part is 2^k times any other part, for any k > 0.

Original entry on oeis.org

1, 1, 2, 2, 4, 4, 6, 9, 12, 13, 18, 23, 29, 37, 49, 55, 71, 84, 104, 126, 153, 185, 221, 261, 317, 375, 446, 523, 623, 721, 854, 994, 1168, 1357, 1579, 1833, 2126, 2455, 2843, 3270, 3766, 4320, 4980, 5687, 6521, 7444, 8498, 9684, 11039, 12540, 14262
Offset: 0

Views

Author

Gus Wiseman, Jan 04 2019

Keywords

Examples

			The a(1) = 1 through a(8) = 12 integer partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (111)  (22)    (32)     (33)      (43)       (44)
                    (31)    (311)    (51)      (52)       (53)
                    (1111)  (11111)  (222)     (61)       (62)
                                     (3111)    (322)      (71)
                                     (111111)  (331)      (332)
                                               (511)      (611)
                                               (31111)    (2222)
                                               (1111111)  (3311)
                                                          (5111)
                                                          (311111)
                                                          (11111111)
		

Crossrefs

Programs

  • Mathematica
    stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
    Table[Length[Select[IntegerPartitions[n],stableQ[#,IntegerQ[Log[2,#1/#2]]&]&]],{n,30}]

A350841 Heinz numbers of integer partitions with a difference < -1 and a conjugate difference < -1.

Original entry on oeis.org

20, 28, 40, 44, 52, 56, 63, 68, 76, 80, 84, 88, 92, 99, 100, 104, 112, 116, 117, 124, 126, 132, 136, 140, 148, 152, 153, 156, 160, 164, 168, 171, 172, 176, 184, 188, 189, 196, 198, 200, 204, 207, 208, 212, 220, 224, 228, 232, 234, 236, 244, 248, 252, 260, 261
Offset: 1

Views

Author

Gus Wiseman, Jan 26 2022

Keywords

Comments

We define a difference of a partition to be a difference of two adjacent parts.

Examples

			The terms together with their prime indices begin:
   20: (3,1,1)
   28: (4,1,1)
   40: (3,1,1,1)
   44: (5,1,1)
   52: (6,1,1)
   56: (4,1,1,1)
   63: (4,2,2)
   68: (7,1,1)
   76: (8,1,1)
   80: (3,1,1,1,1)
   84: (4,2,1,1)
   88: (5,1,1,1)
   92: (9,1,1)
   99: (5,2,2)
		

Crossrefs

Heinz number rankings are in parentheses below.
Taking just one condition gives (A073492) and (A065201), counted by A239955.
These partitions are counted by A350839.
A000041 = integer partitions, strict A000009.
A034296 = partitions with no gaps (A073491), strict A001227 (A073485).
A090858 = partitions with a single gap of size 1 (A325284).
A116931 = partitions with no successions (A319630), strict A003114.
A116932 = partitions with no successions or gaps of size 1, strict A025157.
A350842 = partitions with no gaps of size 1, strict A350844, sets A005314.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    Select[Range[100],(Min@@Differences[Reverse[primeMS[#]]]<-1)&&(Min@@Differences[conj[primeMS[#]]]<-1)&]

A305150 Number of factorizations of n into distinct, pairwise indivisible factors greater than 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 3, 1, 2, 1, 2, 1, 5, 1, 1, 2, 2, 2, 2, 1, 2, 2, 3, 1, 5, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 1, 3, 2, 3, 2, 2, 1, 6, 1, 2, 2, 1, 2, 5, 1, 2, 2, 5, 1, 3, 1, 2, 2, 2, 2, 5, 1, 3, 1, 2, 1, 6, 2, 2, 2, 3, 1, 6, 2, 2, 2, 2, 2, 4, 1, 2, 2, 2, 1, 5, 1, 3, 5
Offset: 1

Views

Author

Gus Wiseman, May 26 2018

Keywords

Examples

			The a(60) = 6 factorizations are (3 * 4 * 5), (3 * 20), (4 * 15), (5 * 12), (6 * 10), (60). Missing from this list are (2 * 3 * 10), (2 * 5 * 6), (2 * 30).
		

Crossrefs

Programs

  • Mathematica
    facs[n_] := If[n <= 1, {{}}, Join@@Table[Map[Prepend[#, d] &, Select[facs[n/d], Min@@ # >= d &]], {d, Rest[Divisors[n]]}]]; Table[Length[Select[facs[n], UnsameQ@@ # && Select[Tuples[Union[#], 2], UnsameQ@@ # && Divisible@@ # &] == {} &]], {n, 100}]
  • PARI
    A305150(n, m=n, facs=List([])) = if(1==n, 1, my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m)&&factorback(apply(x -> (x%d),Vec(facs))), newfacs = List(facs); listput(newfacs,d); s += A305150(n/d, d-1, newfacs))); (s)); \\ Antti Karttunen, Dec 06 2018

Formula

a(n) <= A045778(n) <= A001055(n). - Antti Karttunen, Dec 06 2018

Extensions

More terms from Antti Karttunen, Dec 06 2018

A326077 Number of maximal primitive subsets of {1..n}.

Original entry on oeis.org

1, 1, 2, 2, 3, 3, 4, 4, 6, 7, 11, 11, 13, 13, 23, 24, 36, 36, 48, 48, 64, 66, 126, 126, 150, 151, 295, 363, 507, 507, 595, 595, 895, 903, 1787, 1788, 2076, 2076, 4132, 4148, 5396, 5396, 6644, 6644, 9740, 11172, 22300, 22300, 26140, 26141, 40733, 40773, 60333, 60333, 80781, 80783
Offset: 0

Views

Author

Gus Wiseman, Jun 05 2019

Keywords

Comments

a(n) is the number of maximal primitive subsets of {1, ..., n}. Here primitive means that no element of the subset divides any other and maximal means that no element can be added to the subset while maintaining the property of being pairwise indivisible. - Nathan McNew, Aug 10 2020

Examples

			The a(0) = 1 through a(9) = 7 sets:
  {}  {1}  {1}  {1}   {1}   {1}    {1}    {1}     {1}     {1}
           {2}  {23}  {23}  {235}  {235}  {2357}  {2357}  {2357}
                      {34}  {345}  {345}  {3457}  {3457}  {2579}
                                   {456}  {4567}  {3578}  {3457}
                                                  {4567}  {3578}
                                                  {5678}  {45679}
                                                          {56789}
		

Crossrefs

Programs

  • Mathematica
    stableQ[u_, Q_]:=!Apply[Or, Outer[#1=!=#2&&Q[#1, #2]&, u, u, 1], {0, 1}];
    fasmax[y_]:=Complement[y, Union@@(Most[Subsets[#]]&/@y)];
    Table[Length[fasmax[Select[Subsets[Range[n]],stableQ[#,Divisible]&]]],{n,0,10}]
  • PARI
    divset(n)={sumdiv(n, d, if(dif(k>#p, ismax(b), my(f=!bitand(p[k], b)); if(!f || bittest(d, k), self()(k+1, b)) + if(f, self()(k+1, b+(1<Andrew Howroyd, Aug 30 2019

Extensions

Terms a(19) to a(55) from Andrew Howroyd, Aug 30 2019
Name edited by Nathan McNew, Aug 10 2020

A350837 Number of integer partitions of n with no adjacent parts of quotient 2.

Original entry on oeis.org

1, 1, 2, 2, 4, 5, 7, 10, 14, 18, 24, 31, 41, 53, 70, 87, 112, 140, 178, 221, 277, 344, 428, 526, 648, 792, 971, 1180, 1436, 1738, 2103, 2533, 3049, 3660, 4387, 5242, 6259, 7450, 8860, 10511, 12453, 14723, 17387, 20489, 24121, 28343, 33269, 38982, 45632, 53327
Offset: 0

Views

Author

Gus Wiseman, Jan 18 2022

Keywords

Comments

The first of these partitions that is not double-free (see A323092 for definition) is (4,3,2).

Examples

			The a(1) = 1 through a(7) = 10 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)
       (11)  (111)  (22)    (32)     (33)      (43)
                    (31)    (41)     (51)      (52)
                    (1111)  (311)    (222)     (61)
                            (11111)  (411)     (322)
                                     (3111)    (331)
                                     (111111)  (511)
                                               (4111)
                                               (31111)
                                               (1111111)
		

Crossrefs

The version with quotients >= 2 is A000929, sets A018819.
<= 2 is A342094, ranked by A342191.
< 2 is A342096, sets A045690, strict A342097.
> 2 is A342098, sets A040039.
The sets version (subsets of prescribed maximum) is A045691.
These partitions are ranked by A350838.
The strict case is A350840.
A version for differences is A350842, strict A350844.
The complement is counted by A350846, ranked by A350845.
A000041 = integer partitions.
A116931 = partitions with no successions, ranked by A319630.
A116932 = partitions with differences != 1 or 2, strict A025157.
A323092 = double-free partitions, ranked by A320340.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n], FreeQ[Divide@@@Partition[#,2,1],2]&]],{n,0,15}]

A371178 Number of integer partitions of n containing all divisors of all parts.

Original entry on oeis.org

1, 1, 1, 2, 3, 4, 6, 9, 12, 16, 21, 28, 37, 48, 62, 80, 101, 127, 162, 202, 252, 312, 386, 475, 585, 713, 869, 1056, 1278, 1541, 1859, 2232, 2675, 3196, 3811, 4534, 5386, 6379, 7547, 8908, 10497, 12345, 14501, 16999, 19897, 23253, 27135, 31618, 36796, 42756
Offset: 0

Views

Author

Gus Wiseman, Mar 17 2024

Keywords

Comments

The Heinz numbers of these partitions are given by A371177.
Also partitions such that the number of distinct parts is equal to the number of distinct divisors of parts.

Examples

			The partition (4,2,1,1) contains all distinct divisors {1,2,4}, so is counted under a(8).
The partition (4,4,3,2,2,2,1) contains all distinct divisors {1,2,3,4} so is counted under 4 + 4 + 3 + 2 + 2 + 2 + 1 = 18. - _David A. Corneth_, Mar 18 2024
The a(0) = 1 through a(8) = 12 partitions:
  ()  (1)  (11)  (21)   (31)    (221)    (51)      (331)      (71)
                 (111)  (211)   (311)    (321)     (421)      (521)
                        (1111)  (2111)   (2211)    (511)      (3221)
                                (11111)  (3111)    (2221)     (3311)
                                         (21111)   (3211)     (4211)
                                         (111111)  (22111)    (5111)
                                                   (31111)    (22211)
                                                   (211111)   (32111)
                                                   (1111111)  (221111)
                                                              (311111)
                                                              (2111111)
                                                              (11111111)
		

Crossrefs

The LHS is represented by A001221, distinct case of A001222.
For partitions with no divisors of parts we have A305148, ranks A316476.
The RHS is represented by A370820, for prime factors A303975.
The strict case is A371128.
Counting all parts on the LHS gives A371130, ranks A370802.
The complement is counted by A371132.
For submultisets instead of distinct parts we have A371172, ranks A371165.
These partitions have ranks A371177.
A000005 counts divisors.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],SubsetQ[#,Union@@Divisors/@#]&]],{n,0,30}]

A328677 Numbers whose distinct prime indices are relatively prime and pairwise indivisible.

Original entry on oeis.org

2, 4, 8, 15, 16, 32, 33, 35, 45, 51, 55, 64, 69, 75, 77, 85, 93, 95, 99, 119, 123, 128, 135, 141, 143, 145, 153, 155, 161, 165, 175, 177, 187, 201, 205, 207, 209, 215, 217, 219, 221, 225, 245, 249, 253, 255, 256, 265, 275, 279, 287, 291, 295, 297, 309, 323
Offset: 1

Views

Author

Gus Wiseman, Oct 30 2019

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. Stable numbers are listed in A316476.

Examples

			The sequence of terms together with their prime indices begins:
    2: {1}
    4: {1,1}
    8: {1,1,1}
   15: {2,3}
   16: {1,1,1,1}
   32: {1,1,1,1,1}
   33: {2,5}
   35: {3,4}
   45: {2,2,3}
   51: {2,7}
   55: {3,5}
   64: {1,1,1,1,1,1}
   69: {2,9}
   75: {2,3,3}
   77: {4,5}
   85: {3,7}
   93: {2,11}
   95: {3,8}
   99: {2,2,5}
  119: {4,7}
		

Crossrefs

These are the Heinz numbers of the partitions counted by A328676.
Numbers whose prime indices are relatively prime are A289509.
Partitions whose distinct parts are pairwise indivisible are A305148.
The version for binary indices (instead of prime indices) is A328671.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
    Select[Range[100],GCD@@primeMS[#]==1&&stableQ[primeMS[#],Divisible]&]

Formula

Intersection of A316476 and A289509.

A202425 Number of partitions of n into parts having pairwise common factors but no overall common factor.

Original entry on oeis.org

1, 0, 0, 0, 0, 0, 2, 0, 0, 0, 3, 0, 3, 0, 0, 1, 6, 0, 5, 0, 2, 2, 9, 0, 8, 2, 4, 3, 16, 0, 22, 5, 6, 5, 19, 2, 35, 8, 14, 6, 44, 4, 55, 13, 16, 19, 64, 6, 82, 17, 39, 31, 108, 10, 105, 40, 66, 46, 161, 14, 182, 61, 97, 72, 207, 37, 287, 85, 144, 93, 357, 59
Offset: 31

Views

Author

Alois P. Heinz, Dec 19 2011

Keywords

Examples

			a(31) = 1: [6,10,15] = [2*3,2*5,3*5].
a(37) = 2: [6,6,10,15], [10,12,15].
a(41) = 3: [6,10,10,15], [6,15,20], [6,14,21].
a(47) = 6: [6,6,10,10,15], [10,10,12,15], [6,6,15,20], [12,15,20], [6,6,14,21], [12,14,21].
a(49) = 5: [6,6,6,6,10,15], [6,6,10,12,15], [10,12,12,15], [6,10,15,18], [10,15,24].
		

Crossrefs

The version with only distinct parts compared is A328672.
The Heinz numbers of these partitions are A328868.
The strict case is A202385, which is essentially the same as A318715.
The version for non-isomorphic multiset partitions is A319759.
The version for set-systems is A326364.
Intersecting partitions are A200976.

Programs

  • Maple
    with(numtheory):
    w:= (m, h)-> mul(`if`(j>=h, 1, j), j=factorset(m)):
    b:= proc(n, i, g, s) option remember; local j, ok, si;
          if n<0 then 0
        elif n=0 then `if`(g>1, 0, 1)
        elif i<2 or member(1, s) then 0
        else ok:= evalb(i<=n);
             si:= map(x->w(x, i), s);
             for j in s while ok do ok:= igcd(i, j)>1 od;
             b(n, i-1, g, si) +`if`(ok, add(b(n-t*i, i-1, igcd(i, g),
                          si union {w(i,i)} ), t=1..iquo(n, i)), 0)
          fi
        end:
    a:= n-> b(n, n, 0, {}):
    seq(a(n), n=31..100);
  • Mathematica
    w[m_, h_] := Product[If[j >= h, 1, j], {j, FactorInteger[m][[All, 1]]}]; b[n_, i_, g_, s_] := b[n, i, g, s] = Module[{j, ok, si}, Which[n<0, 0, n == 0, If[g>1, 0, 1], i<2 || MemberQ[s, 1], 0, True, ok = (i <= n); si = w[#, i]& /@ s; Do[If[ok, ok = (GCD[i, j]>1)], {j, s}]; b[n, i-1, g, si] + If[ok, Sum[b[n-t*i, i-1, GCD[i, g], si ~Union~ {w[i, i]}], {t, 1, Quotient[n, i]}], 0]]]; a[n_] := b[n, n, 0, {}]; Table[a[n], {n, 31, 100}] (* Jean-François Alcover, Feb 16 2017, translated from Maple *)
    Table[Length[Select[IntegerPartitions[n],GCD@@#==1&&And@@(GCD[##]>1&)@@@Tuples[#,2]&]],{n,0,40}] (* Gus Wiseman, Nov 04 2019 *)

Formula

a(n > 0) = A328672(n) - 1. - Gus Wiseman, Nov 04 2019

A328671 Numbers whose binary indices are relatively prime and pairwise indivisible.

Original entry on oeis.org

1, 6, 12, 18, 20, 22, 24, 28, 48, 56, 66, 68, 70, 72, 76, 80, 82, 84, 86, 88, 92, 96, 104, 112, 120, 132, 144, 148, 176, 192, 196, 208, 212, 224, 240, 258, 264, 272, 274, 280, 296, 304, 312, 320, 322, 328, 336, 338, 344, 352, 360, 368, 376, 384, 400, 416, 432
Offset: 1

Views

Author

Gus Wiseman, Oct 29 2019

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.

Examples

			The sequence of terms together with their binary expansions and binary indices begins:
    1:         1 ~ {1}
    6:       110 ~ {2,3}
   12:      1100 ~ {3,4}
   18:     10010 ~ {2,5}
   20:     10100 ~ {3,5}
   22:     10110 ~ {2,3,5}
   24:     11000 ~ {4,5}
   28:     11100 ~ {3,4,5}
   48:    110000 ~ {5,6}
   56:    111000 ~ {4,5,6}
   66:   1000010 ~ {2,7}
   68:   1000100 ~ {3,7}
   70:   1000110 ~ {2,3,7}
   72:   1001000 ~ {4,7}
   76:   1001100 ~ {3,4,7}
   80:   1010000 ~ {5,7}
   82:   1010010 ~ {2,5,7}
   84:   1010100 ~ {3,5,7}
   86:   1010110 ~ {2,3,5,7}
   88:   1011000 ~ {4,5,7}
		

Crossrefs

The version for prime indices (instead of binary indices) is A328677.
Numbers whose binary indices are relatively prime are A291166.
Numbers whose distinct prime indices are pairwise indivisible are A316476.
BII-numbers of antichains are A326704.
Relatively prime partitions whose distinct parts are pairwise indivisible are A328676, with strict case A328678.

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
    Select[Range[100],GCD@@bpe[#]==1&&stableQ[bpe[#],Divisible]&]

Formula

Intersection of A291166 with A326704.

A328672 Number of integer partitions of n with relatively prime parts in which no two distinct parts are relatively prime.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 4, 1, 4, 1, 1, 2, 7, 1, 6, 1, 3, 3, 10, 1, 9, 3, 5, 4, 17, 1, 23, 6, 7, 6, 20, 3, 36, 9, 15, 7, 45, 5, 56, 14, 17, 20, 65, 7, 83, 18, 40
Offset: 0

Views

Author

Gus Wiseman, Oct 29 2019

Keywords

Comments

Positions of terms greater than 1 are {31, 37, 41, 43, 46, 47, 49, ...}.
A partition with no two distinct parts relatively prime is said to be intersecting.

Examples

			Examples:
  a(31) = 2:         a(46) = 2:
    (15,10,6)          (15,15,10,6)
    (1^31)             (1^46)
  a(37) = 3:         a(47) = 7:
    (15,12,10)         (20,15,12)
    (15,10,6,6)        (21,14,12)
    (1^37)             (20,15,6,6)
  a(41) = 4:           (21,14,6,6)
    (20,15,6)          (15,12,10,10)
    (21,14,6)          (15,10,10,6,6)
    (15,10,10,6)       (1^47)
    (1^41)           a(49) = 6:
  a(43) = 4:           (24,15,10)
    (18,15,10)         (18,15,10,6)
    (15,12,10,6)       (15,12,12,10)
    (15,10,6,6,6)      (15,12,10,6,6)
    (1^43)             (15,10,6,6,6,6)
                       (1^39)
		

Crossrefs

The Heinz numbers of these partitions are A328679.
The strict case is A318715.
The version for non-isomorphic multiset partitions is A319759.
Relatively prime partitions are A000837.
Intersecting partitions are A328673.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],GCD@@#==1&&And[And@@(GCD[##]>1&)@@@Subsets[Union[#],{2}]]&]],{n,0,32}]

Formula

a(n > 0) = A202425(n) + 1.
Previous Showing 11-20 of 33 results. Next