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

A305657 Record values of A305566.

Original entry on oeis.org

0, 2, 10, 44, 84, 122, 184, 1590, 26508, 119304, 428568, 7911144, 8275066, 510582192, 2122131332, 34308911168, 543326090824, 140627994008752, 140730946295934, 576036971191781024, 9222527848311467840, 9222949803150423996, 2359453722081533041216, 2361155375874259181576, 2361182396716073890056, 604433511573276736205056, 604435239966141944446584, 9664324052998090353961088, 618962932654137056803769392, 39614074145912329143754325568, 633821673403520815917635373056, 162257419909717745759000144646368, 664613879048405558070320228193928512, 170141183222784760671064904320665465984, 11150370605423278283522757581632660569339264
Offset: 1

Views

Author

Robert Israel, Jun 07 2018

Keywords

Comments

Numbers k such that for some j, k = A305566(j) and k > A305566(i) for i < j.
The corresponding record locations appear to be the members of A097212 except 2 and 4.

Crossrefs

Programs

  • Maple
    # with f as in A305566
    recs:= 0: count:= 0: m:= 0:
    Agenda:= heap[new]((s,t) -> s[2]>t[2], [[1],2]);
    while count < 34 do
      T:= heap[extract](Agenda);
      v:= f(T[2]);
      if v > m then
        count:= count+1;
        recs:= recs, v;
        m:= v;
      fi;
      L:= T[1];
      heap[insert]([[op(L),1],T[2]*ithprime(nops(L)+1)],Agenda);
      heap[insert]([L+[1,0$(nops(L)-1)], 2*T[2]],Agenda);
      for j in select(i -> L[i]
    				

A074761 Number of partitions of n of order n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 9, 1, 4, 5, 1, 1, 12, 1, 27, 7, 6, 1, 81, 1, 7, 1, 54, 1, 407, 1, 1, 11, 9, 13, 494, 1, 10, 13, 423, 1, 981, 1, 137, 115, 12, 1, 1309, 1, 59, 17, 193, 1, 240, 21, 1207, 19, 15, 1, 47274, 1, 16, 239, 1, 25, 3284, 1, 333, 23, 3731, 1, 42109, 1, 19
Offset: 1

Views

Author

Vladeta Jovovic, Sep 28 2002

Keywords

Comments

Order of partition is lcm of its parts.
a(n) is the number of conjugacy classes of the symmetric group S_n such that a representative of the class has order n. Here order means the order of an element of a group. Note that a(n) = 1 if and only if n is a prime power. - W. Edwin Clark, Aug 05 2014

Examples

			The a(15) = 5 partitions are (15), (5,3,3,3,1), (5,5,3,1,1), (5,3,3,1,1,1,1), (5,3,1,1,1,1,1,1,1). - _Gus Wiseman_, Aug 01 2018
		

Crossrefs

Programs

  • Maple
    A:= proc(n)
          uses numtheory;
          local S;
        S:= add(mobius(n/i)*1/mul(1-x^j,j=divisors(i)),i=divisors(n));
        coeff(series(S,x,n+1),x,n);
    end proc:
    seq(A(n),n=1..100); # Robert Israel, Aug 06 2014
  • Mathematica
    a[n_] := With[{s = Sum[MoebiusMu[n/i]*1/Product[1-x^j, {j, Divisors[i]}], {i, Divisors[n]}]}, SeriesCoefficient[s, {x, 0, n}]]; Array[a, 80] (* Jean-François Alcover, Feb 29 2016 *)
    Table[Length[Select[IntegerPartitions[n],LCM@@#==n&]],{n,50}] (* Gus Wiseman, Aug 01 2018 *)
  • PARI
    pr(k, x)={my(t=1); fordiv(k, d, t *= (1-x^d) ); return(t); }
    a(n) =
    {
        my( x = 'x+O('x^(n+1)) );
        polcoeff( Pol( sumdiv(n, i, moebius(n/i) / pr(i, x) ) ), n );
    }
    vector(66, n, a(n) )
    \\ Joerg Arndt, Aug 06 2014

Formula

Coefficient of x^n in expansion of Sum_{i divides n} A008683(n/i)*1/Product_{j divides i} (1-x^j).

A305563 Number of reducible integer partitions of n.

Original entry on oeis.org

1, 2, 3, 4, 7, 7, 15, 16, 27, 30, 56, 56, 100, 105, 157, 188, 287, 303, 470, 524, 724, 850, 1197, 1339, 1856, 2135, 2814, 3305, 4360, 4951, 6532, 7561, 9563, 11195, 14165, 16328, 20631, 23866, 29471, 34320, 42336, 48672, 59872, 69139, 83625, 96911, 117153
Offset: 1

Views

Author

Gus Wiseman, Jun 05 2018

Keywords

Comments

A multiset m whose distinct elements are m_1, m_2, ..., m_k with multiplicities y_1, y_2, ..., y_k is reducible if either m is of size 1 or gcd(m_1, ..., m_k) = 1 and the multiset {y_1, ..., y_k} is also reducible.

Examples

			The a(6) = 7 reducible integer partitions are (6), (51), (411), (321), (3111), (21111), (111111). Missing from this list are (42), (33), (222), (2211).
		

Crossrefs

Programs

  • Mathematica
    ptnredQ[y_]:=Or[Length[y]==1,And[GCD@@y==1,ptnredQ[Sort[Length/@Split[y],Greater]]]];
    Table[Length[Select[IntegerPartitions[n],ptnredQ]],{n,20}]

A074971 Number of partitions of n into distinct parts of order n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 1, 4, 1, 1, 1, 2, 1, 4, 1, 1, 1, 1, 1, 6, 1, 1, 1, 3, 1, 3, 1, 1, 1, 1, 1, 6, 1, 1, 1, 1, 1, 4, 1, 2, 1, 1, 1, 32, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 25, 1, 1, 1, 1, 1, 2, 1, 4, 1, 1, 1, 24, 1, 1, 1, 2, 1, 24, 1, 1, 1, 1, 1, 12, 1, 1, 1, 3, 1, 2
Offset: 1

Views

Author

Vladeta Jovovic, Oct 05 2002

Keywords

Comments

Order of partition is lcm of its parts.

Examples

			The a(36) = 6 partitions are (36), (18,12,6), (18,12,4,2), (18,12,3,2,1), (18,9,4,3,2), (12,9,6,4,3,2). - _Gus Wiseman_, Aug 01 2018
		

Crossrefs

Programs

  • PARI
    A074971(n) = { my(q=0); fordiv(n,i,my(p=1); fordiv(i,j,p *= (1 + 'x^j)); q += moebius(n/i)*p); polcoeff(q,n); }; \\ Antti Karttunen, Dec 19 2018

Formula

Coefficient of x^n in expansion of Sum_{i divides n} mu(n/i)*Product_{j divides i} (1+x^j).

A316433 Number of integer partitions of n whose length is equal to the LCM of all parts.

Original entry on oeis.org

1, 0, 1, 1, 1, 0, 2, 1, 4, 3, 4, 4, 8, 5, 7, 8, 10, 8, 13, 13, 20, 18, 25, 25, 36, 34, 48, 52, 64, 64, 85, 85, 108, 106, 129, 133, 160, 158, 189, 194, 229, 228, 276, 279, 332, 336, 394, 402, 476, 489, 572, 599, 699, 728, 845, 889, 1032, 1094, 1251, 1332, 1523
Offset: 1

Views

Author

Gus Wiseman, Jul 02 2018

Keywords

Examples

			The a(13) = 8 partitions are (4441), (55111), (322222), (332221), (333211), (622111), (631111), (7111111).
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],LCM@@#==Length[#]&]],{n,30}]
  • PARI
    a(n) = {my(nb = 0); forpart(p=n, if (lcm(Vec(p))==#p, nb++);); nb;} \\ Michel Marcus, Jul 03 2018

A319333 Heinz numbers of integer partitions whose sum is equal to their LCM.

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, 157, 163, 167, 173, 179, 181, 191, 193, 197, 198, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263
Offset: 1

Views

Author

Gus Wiseman, Sep 17 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 partitions whose Heinz numbers are in the sequence begins: (1), (2), (3), (4), (5), (6), (7), (8), (9), (10), (3,2,1), (11), (12).
		

Crossrefs

Programs

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

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

A305564 Number of finite sets of relatively prime positive integers with least common multiple n.

Original entry on oeis.org

1, 1, 1, 2, 1, 7, 1, 4, 2, 7, 1, 32, 1, 7, 7, 8, 1, 32, 1, 32, 7, 7, 1, 136, 2, 7, 4, 32, 1, 193, 1, 16, 7, 7, 7, 322, 1, 7, 7, 136, 1, 193, 1, 32, 32, 7, 1, 560, 2, 32, 7, 32, 1, 136, 7, 136, 7, 7, 1, 3464, 1, 7, 32, 32, 7, 193, 1, 32, 7, 193, 1, 2852, 1, 7
Offset: 1

Views

Author

Gus Wiseman, Jun 05 2018

Keywords

Examples

			The a(6) = 7 sets are {1,6}, {2,3}, {1,2,3}, {1,2,6}, {1,3,6}, {2,3,6}, {1,2,3,6}.
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Rest[Subsets[Divisors[n]]],And[GCD@@#==1,LCM@@#==n]&]],{n,100}]

A305565 Regular triangle where T(n,k) is the number of finite sets of positive integers with least common multiple n and greatest common divisor k.

Original entry on oeis.org

1, 1, 1, 1, 0, 1, 2, 1, 0, 1, 1, 0, 0, 0, 1, 7, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 4, 2, 0, 1, 0, 0, 0, 1, 2, 0, 1, 0, 0, 0, 0, 0, 1, 7, 1, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 32, 7, 2, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Jun 05 2018

Keywords

Examples

			The T(12,2) = 7 sets are {2,12}, {4,6}, {2,4,6}, {2,4,12}, {2,6,12}, {4,6,12}, {2,4,6,12}.
Triangle begins:
   1
   1  1
   1  0  1
   2  1  0  1
   1  0  0  0  1
   7  1  1  0  0  1
   1  0  0  0  0  0  1
   4  2  0  1  0  0  0  1
   2  0  1  0  0  0  0  0  1
   7  1  0  0  1  0  0  0  0  1
   1  0  0  0  0  0  0  0  0  0  1
  32  7  2  1  0  1  0  0  0  0  0  1
		

Crossrefs

Programs

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

Formula

If k divides n then T(n,k) = T(n/k,1) = A305564(n/k); otherwise T(n,k) = 0.

A305567 Irregular triangle where T(n,k) is the number of finite sets of positive integers with least common multiple n and greatest common divisor k, where k runs over all divisors of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 4, 2, 1, 1, 2, 1, 1, 7, 1, 1, 1, 1, 1, 32, 7, 2, 1, 1, 1, 1, 1, 7, 1, 1, 1, 7, 1, 1, 1, 8, 4, 2, 1, 1, 1, 1, 32, 2, 7, 1, 1, 1, 1, 1, 32, 7, 1, 2, 1, 1, 7, 1, 1, 1, 7, 1, 1, 1, 1, 1, 136, 32, 4, 7, 2, 1, 1, 1, 2
Offset: 1

Views

Author

Gus Wiseman, Jun 05 2018

Keywords

Examples

			Triangle begins:
   1
   1  1
   1  1
   2  1  1
   1  1
   7  1  1  1
   1  1
   4  2  1  1
   2  1  1
   7  1  1  1
   1  1
  32  7  2  1  1  1
   1  1
   7  1  1  1
   7  1  1  1
   8  4  2  1  1
   1  1
  32  2  7  1  1  1
   1  1
  32  7  1  2  1  1
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Subsets[Divisors[n]],And[GCD@@#==k,LCM@@#==n]&]],{n,30},{k,Divisors[n]}]
Showing 1-10 of 11 results. Next