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 41-50 of 80 results. Next

A317624 Number of integer partitions of n where all parts are > 1 and whose LCM is n.

Original entry on oeis.org

0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 3, 1, 5, 1, 1, 1, 17, 1, 1, 1, 7, 1, 60, 1, 1, 1, 1, 1, 76, 1, 1, 1, 55, 1, 105, 1, 11, 10, 1, 1, 187, 1, 6, 1, 13, 1, 30, 1, 111, 1, 1, 1, 5043, 1, 1, 15, 1, 1, 230, 1, 17, 1, 242, 1, 4173, 1, 1, 12, 19, 1
Offset: 0

Views

Author

Gus Wiseman, Aug 01 2018

Keywords

Examples

			The a(20) = 5 partitions are (20), (10,4,4,2), (10,4,2,2,2), (5,5,4,4,2), (5,5,4,2,2,2).
The a(45) = 10 partitions:
  (45),
  (15,15,9,3,3), (15,9,9,9,3),
  (15,9,9,3,3,3,3), (15,9,5,5,5,3,3), (9,9,9,5,5,5,3),
  (15,9,3,3,3,3,3,3,3), (9,9,5,5,5,3,3,3,3), (9,5,5,5,5,5,5,3,3),
  (9,5,5,5,3,3,3,3,3,3,3).
From _David A. Corneth_, Sep 08 2018: (Start)
Let sum(t) denote the sum of elements of a tuple t. The tuples t with distinct divisors of 45 that have lcm(t) = 45 and sum(t) <= 45 are {(45) and (3, 9, 15), (3, 5, 9, 15), (3, 5, 9), (5, 9), (9, 15), (5, 9, 15)}. For each such tuple t, find the number of partitions of 45 - s(t) into distinct parts of t.
For the tuple (45), there is 1 partition of 45 - 45 = 0 into parts with 45. That is: {()}.
For the tuple (3, 9, 15), there are 4 partitions of 45 - (3 + 9 + 15) = 18 into parts with 3, 9 and 15. They are {(3, 15), (9, 9), (3, 3, 3, 9), (3, 3, 3, 3, 3, 3)}.
For the tuple (3, 5, 9), there are 4 partitions of 45 - (3 + 5 + 9) = 28 into parts with 3, 5 and 9; they are {(5, 5, 9, 9), (3, 3, 3, 5, 5, 9), (3, 5, 5, 5, 5, 5), (3, 3, 3, 3, 3, 3, 5, 5)}.
For the tuple (3, 5, 9, 15), there is 1 partition of 45 - (3 + 5 + 9 + 15) = 13 into parts with 3, 5, 9 and 15. That is (3, 5, 5).
The other tuples, (5, 9), (9, 15), and (5, 9, 15); they give no extra tuples. That's because there is no solution to the Diophantine equation for 5x + 9y = 45 - (5 + 9), corresponding to the tuple (5, 9) with nonnegative x, y.
That also excludes (9, 15); if there is a solution for that, there would also be a solution for (5, 9). This could whittle down the number of seeds even further. Similarly, (5, 9, 15) gives no solution.
Therefore a(45) = 1 + 4 + 4 + 1 = 10.
(End)
In general, there are A318670(n) (<= A069626(n)) such seed sets of divisors where to start extending the partition from. (See the second PARI program which uses subroutine toplevel_starting_sets.) - _Antti Karttunen_, Sep 08 2018
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],And[Min@@#>=2,LCM@@#==n]&]],{n,30}]
  • PARI
    strong_divisors_reversed(n) = vecsort(select(x -> (x>1), divisors(n)), , 4);
    partitions_into_lcm(orgn,n,parts,from=1,m=1) = if(!n,(m==orgn),my(k = #parts, s=0); for(i=from,k,if(parts[i]<=n, s += partitions_into_lcm(orgn,n-parts[i],parts,i,lcm(m,parts[i])))); (s));
    A317624(n) = if(n<=1,0,partitions_into_lcm(n,n,strong_divisors_reversed(n))); \\ Antti Karttunen, Sep 07 2018
    
  • PARI
    strong_divisors_reversed(n) = vecsort(select(x -> (x>1), divisors(n)), , 4);
    partitions_into(n,parts,from=1) = if(!n,1, if(#parts==from, (0==(n%parts[from])), my(s=0); for(i=from,#parts,if(parts[i]<=n, s += partitions_into(n-parts[i],parts,i))); (s)));
    toplevel_starting_sets(orgn,n,parts,from=1,ss=List([])) = { my(k = #parts, s=0, newss); if(lcm(Vec(ss))==orgn,s += partitions_into(n,ss)); for(i=from,k,if(parts[i]<=n, newss = List(ss); listput(newss,parts[i]); s += toplevel_starting_sets(orgn,n-parts[i],parts,i+1,newss))); (s) };
    A317624(n) = if(n<=1,0,toplevel_starting_sets(n,n,strong_divisors_reversed(n))); \\ Antti Karttunen, Sep 08-10 2018

A319320 Number of integer partitions of n such that every distinct submultiset has a different LCM.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 3, 2, 3, 4, 5, 4, 6, 7, 7, 9, 11, 12, 12, 15, 17, 20, 22, 24, 25, 31, 35, 39, 40, 48, 51, 55, 64, 73, 77, 85, 92, 104, 115, 126, 136, 147, 157, 176, 198, 211, 234, 246, 269, 294, 326, 350, 375, 403, 443, 475, 526, 560, 600, 650
Offset: 1

Views

Author

Gus Wiseman, Sep 17 2018

Keywords

Comments

Note that such partitions are necessarily strict.

Examples

			The a(19) = 12 partitions:
  (19),
  (10,9), (11,8), (12,7), (13,6), (14,5), (15,4), (16,3), (17,2),
  (8,6,5), (11,5,3),
  (7,5,4,3).
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&UnsameQ@@LCM@@@Union[Rest[Subsets[#]]]&]],{n,30}]

A319837 Numbers whose distinct prime indices are pairwise indivisible and whose own prime indices span an initial interval of positive integers.

Original entry on oeis.org

1, 2, 3, 4, 7, 8, 9, 13, 15, 16, 19, 27, 32, 35, 37, 45, 49, 53, 61, 64, 69, 75, 81, 89, 91, 95, 113, 128, 131, 135, 141, 143, 145, 151, 161, 165, 169, 175, 207, 223, 225, 243, 245, 247, 251, 256, 265, 281, 299, 309, 311, 329, 343, 355, 359, 361, 375, 377, 385
Offset: 1

Views

Author

Gus Wiseman, Dec 16 2018

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. The multiset multisystem with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}. This sequence lists all MM-numbers of not necessarily strict antichains of multisets spanning an initial interval of positive integers.

Examples

			The sequence of multisystems whose MM-numbers belong to the sequence begins:
   1: {}
   2: {{}}
   3: {{1}}
   4: {{},{}}
   7: {{1,1}}
   8: {{},{},{}}
   9: {{1},{1}}
  13: {{1,2}}
  15: {{1},{2}}
  16: {{},{},{},{}}
  19: {{1,1,1}}
  27: {{1},{1},{1}}
  32: {{},{},{},{},{}}
  35: {{2},{1,1}}
  37: {{1,1,2}}
  45: {{1},{1},{2}}
  49: {{1,1},{1,1}}
  53: {{1,1,1,1}}
  61: {{1,2,2}}
  64: {{},{},{},{},{},{}}
  69: {{1},{2,2}}
  75: {{1},{2},{2}}
  81: {{1},{1},{1},{1}}
  89: {{1,1,1,2}}
  91: {{1,1},{1,2}}
  95: {{2},{1,1,1}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    normQ[sys_]:=Or[Length[sys]==0,Union@@sys==Range[Max@@Max@@sys]];
    stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
    Select[Range[200],And[normQ[primeMS/@primeMS[#]],stableQ[primeMS[#],Divisible]]&]

A320275 Numbers whose distinct prime indices are pairwise indivisible and whose own prime indices are connected and span an initial interval of positive integers.

Original entry on oeis.org

2, 3, 7, 9, 13, 19, 27, 37, 49, 53, 61, 81, 89, 91, 113, 131, 151, 169, 223, 243, 247, 251, 281, 299, 311, 343, 359, 361, 377, 427, 463, 503, 593, 611, 637, 659, 689, 703, 719, 729, 791, 827, 851, 863, 923, 953, 1069, 1073, 1159, 1163, 1183, 1291, 1321, 1339
Offset: 1

Views

Author

Gus Wiseman, Dec 16 2018

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. The multiset multisystem with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}. This sequence lists all MM-numbers of not necessarily strict connected antichains of multisets spanning an initial interval of positive integers.

Examples

			The sequence of multisystems whose MM-numbers belong to the sequence begins:
    2: {{}}
    3: {{1}}
    7: {{1,1}}
    9: {{1},{1}}
   13: {{1,2}}
   19: {{1,1,1}}
   27: {{1},{1},{1}}
   37: {{1,1,2}}
   49: {{1,1},{1,1}}
   53: {{1,1,1,1}}
   61: {{1,2,2}}
   81: {{1},{1},{1},{1}}
   89: {{1,1,1,2}}
   91: {{1,1},{1,2}}
  113: {{1,2,3}}
  131: {{1,1,1,1,1}}
  151: {{1,1,2,2}}
  169: {{1,2},{1,2}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    normQ[sys_]:=Or[Length[sys]==0,Union@@sys==Range[Max@@Max@@sys]];
    zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[Less@@#,GCD@@s[[#]]]>1&]},If[c=={},s,zsm[Sort[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
    stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
    Select[Range[200],And[normQ[primeMS/@primeMS[#]],stableQ[primeMS[#],Divisible],Length[zsm[primeMS[#]]]==1]&]

A322306 Number of connected divisors of n. Number of connected submultisets of the n-th multiset multisystem (A302242).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 03 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. A positive integer is connected if its prime indices are connected (see A305078).

Examples

			The a(1365) = 12 divisors are 3, 5, 7, 13, 21, 39, 65, 91, 195, 273, 455, 1365. These correspond to the following connected submultisets of {{1},{2},{1,1},{1,2}}.
     3: {{1}}
     5: {{2}}
     7: {{1,1}}
    13: {{1,2}}
    21: {{1},{1,1}}
    39: {{1},{1,2}}
    65: {{2},{1,2}}
    91: {{1,1},{1,2}}
   195: {{1},{2},{1,2}}
   273: {{1},{1,1},{1,2}}
   455: {{2},{1,1},{1,2}}
  1365: {{1},{2},{1,1},{1,2}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[Less@@#,GCD@@s[[#]]]>1&]},If[c=={},s,zsm[Union[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
    Table[Length[Select[Union[Subsets[primeMS[n]]],Length[zsm[#]]==1&]],{n,50}]

A322307 Number of multisets in the swell of the n-th multiset multisystem.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 3, 1, 1, 2, 2, 2, 2, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 3, 1, 2, 2, 1, 2, 3, 1, 2, 2, 3, 1, 2, 1, 2, 2, 2, 2, 3, 1, 2, 1, 2, 1, 3, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Dec 03 2018

Keywords

Comments

First differs from A001221 at a(91) = 3, A001221(91) = 2.
The swell of a multiset partition is the set of possible joins of its connected submultisets, where the multiplicity of a vertex in the join of a set of multisets is the maximum multiplicity of the same vertex among the parts. For example the swell of {{1,1},{1,2},{2,2}} is:
{1,1}
{1,2}
{2,2}
{1,1,2}
{1,2,2}
{1,1,2,2}

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    zwell[y_]:=Union[y,Join@@Cases[Subsets[Union[y],{2}],{x_,z_}?(GCD@@#>1&):>zwell[Sort[Append[Fold[DeleteCases[#1,#2,{1},1]&,y,{x,z}],LCM[x,z]]]]]];
    Table[Length[zwell[primeMS[n]]],{n,100}]

A327779 Number of integer partitions of n whose LCM is greater than n.

Original entry on oeis.org

1, 0, 0, 0, 0, 1, 0, 2, 3, 7, 9, 18, 16, 31, 42, 61, 87, 133, 169, 246, 302, 411, 545, 738, 874, 1167, 1497, 1945, 2421, 3110, 3498, 4476, 5615, 7061, 8777, 10925, 12957, 16036, 19644, 24061, 28858, 35177, 41572, 50424, 60643, 72953, 87499, 104893, 123821, 147776
Offset: 0

Views

Author

Gus Wiseman, Sep 25 2019

Keywords

Examples

			The a(5) = 1 through a(12) = 16 partitions (empty columns not shown):
  (32)  (43)  (53)   (54)    (64)     (65)      (75)
        (52)  (431)  (72)    (73)     (74)      (543)
              (521)  (432)   (433)    (83)      (651)
                     (522)   (532)    (92)      (732)
                     (531)   (541)    (443)     (741)
                     (4311)  (721)    (533)     (831)
                     (5211)  (4321)   (542)     (921)
                             (5311)   (641)     (5322)
                             (43111)  (722)     (5331)
                                      (731)     (5421)
                                      (4322)    (7221)
                                      (4331)    (7311)
                                      (5321)    (53211)
                                      (5411)    (54111)
                                      (7211)    (72111)
                                      (43211)   (531111)
                                      (53111)
                                      (431111)
		

Crossrefs

The Heinz numbers of these partitions are given by A327784.
Partitions whose LCM is a multiple of their sum are A327778.
Partitions whose LCM is equal to their sum are A074761.
Partitions whose LCM is less than their sum are A327781.

Programs

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

A327781 Number of integer partitions of n whose LCM is less than n.

Original entry on oeis.org

0, 0, 1, 2, 4, 5, 9, 12, 18, 22, 30, 37, 52, 69, 89, 110, 143, 163, 204, 243, 298, 374, 451, 516, 620, 790, 932, 1064, 1243, 1454, 1699, 2365, 2733, 3071, 3524, 3945, 4526, 5600, 6361, 7111, 8057, 9405, 10621, 12836, 14395, 16066, 18047, 19860, 22143, 25748
Offset: 0

Views

Author

Gus Wiseman, Sep 25 2019

Keywords

Examples

			The a(2) = 1 through a(8) = 18 partitions:
  (11)  (21)   (22)    (41)     (33)      (61)       (44)
        (111)  (31)    (221)    (42)      (322)      (62)
               (211)   (311)    (51)      (331)      (71)
               (1111)  (2111)   (222)     (421)      (332)
                       (11111)  (411)     (511)      (422)
                                (2211)    (2221)     (611)
                                (3111)    (3211)     (2222)
                                (21111)   (4111)     (3221)
                                (111111)  (22111)    (3311)
                                          (31111)    (4211)
                                          (211111)   (5111)
                                          (1111111)  (22211)
                                                     (32111)
                                                     (41111)
                                                     (221111)
                                                     (311111)
                                                     (2111111)
                                                     (11111111)
		

Crossrefs

The Heinz numbers of these partitions are given by A327776.
Partitions whose LCM is equal to their sum are A074761.
Partitions whose LCM is greater than their sum are A327779.

Programs

  • Maple
    a:= proc(m) option remember; local b; b:=
          proc(n, i, l) option remember; `if`(n=0, 1,
           `if`(i>1, b(n, i-1, l), 0) +(h-> `if`(h0, b(m$2, 1), 0)
        end:
    seq(a(n), n=0..70);  # Alois P. Heinz, Oct 10 2019
  • Mathematica
    Table[Length[Select[IntegerPartitions[n],LCM@@#1, b[n, i - 1, l], 0] + Function[h, If[h0, b[m, m, 1], 0]];
    a /@ Range[0, 70] (* Jean-François Alcover, May 18 2021, after Alois P. Heinz *)
  • PARI
    b(m,n)={my(d=divisors(m)); polcoef(1/prod(i=1, #d, 1 - x^d[i] + O(x*x^n)), n)}
    a(n)={sum(m=1, n-1, b(m, n)*sum(i=1, (n-1)\m, moebius(i)))} \\ Andrew Howroyd, Oct 09 2019

A327783 Heinz numbers of integer partitions whose LCM is a multiple of their sum.

Original entry on oeis.org

2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 30, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 154, 157, 163, 165, 167, 173, 179, 181, 190, 191, 193, 197, 198, 199, 211, 223, 227, 229, 233, 239, 241
Offset: 1

Views

Author

Gus Wiseman, Sep 25 2019

Keywords

Comments

First differs from A319333 in having 154.
First nonsquarefree term is 198.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The sequence of terms together with their prime indices begins:
    2: {1}
    3: {2}
    5: {3}
    7: {4}
   11: {5}
   13: {6}
   17: {7}
   19: {8}
   23: {9}
   29: {10}
   30: {1,2,3}
   31: {11}
   37: {12}
   41: {13}
   43: {14}
   47: {15}
   53: {16}
   59: {17}
   61: {18}
   67: {19}
		

Crossrefs

The enumeration of these partitions by sum is A327778.
Heinz numbers of partitions whose LCM is twice their sum are A327775.
Heinz numbers of partitions whose LCM is less than their sum are A327776.
Heinz numbers of partitions whose LCM is greater than their sum are A327784.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[2,100],Divisible[LCM@@primeMS[#],Total[primeMS[#]]]&]

Formula

A056239(a(k)) | A290103(a(k)).

A060766 Least common multiple of differences between consecutive divisors of n (ordered by size).

Original entry on oeis.org

1, 2, 2, 4, 3, 6, 4, 6, 15, 10, 6, 12, 35, 10, 8, 16, 9, 18, 10, 28, 99, 22, 12, 20, 143, 18, 42, 28, 60, 30, 16, 88, 255, 28, 18, 36, 323, 130, 60, 40, 21, 42, 154, 60, 483, 46, 24, 42, 75, 238, 234, 52, 27, 132, 84, 304, 783, 58, 60, 60, 899, 84, 32, 104, 165, 66, 442
Offset: 2

Views

Author

Labos Elemer, Apr 24 2001

Keywords

Examples

			For n=98, divisors={1,2,7,14,49,98}; differences={1,5,7,35,49}; a(98) = LCM of differences = 245.
		

Crossrefs

The GCD version appears to be A258409.
The LCM of the prime indices of n is A290103(n).
The differences between consecutive divisors of n are row n of A193829.

Programs

  • Mathematica
    a[n_ ] := LCM@@(Drop[d=Divisors[n], 1]-Drop[d, -1])
    Table[LCM@@Differences[Divisors[n]],{n,2,70}] (* Harvey P. Dale, Oct 08 2012 *)

Formula

a(n) = A290103(A328023(n)). - Gus Wiseman, Oct 16 2019

Extensions

Edited by Dean Hickerson, Jan 22 2002
Previous Showing 41-50 of 80 results. Next