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 21-30 of 44 results. Next

A321588 Number of connected nonnegative integer matrices with sum of entries equal to n, no zero rows or columns, and distinct rows and columns.

Original entry on oeis.org

1, 1, 1, 9, 29, 181, 1285, 10635, 102355, 1118021, 13637175, 184238115, 2727293893, 43920009785, 764389610843, 14297306352937, 286014489487815, 6093615729757841, 137750602009548533, 3293082026520294529, 83006675263513350581, 2200216851785981586729, 61180266502369886181253
Offset: 0

Views

Author

Gus Wiseman, Nov 13 2018

Keywords

Comments

A matrix is connected if the positions in each row (or each column) of the nonzero entries form a connected hypergraph.

Examples

			The a(4) = 29 matrices:
4 31 13
.
3 21 21 20 12 12 11 110 11 110 101 101 1 10 10 02 011 011 01 01
1 10 01 11 10 01 20 101 02 011 110 011 3 21 12 11 110 101 21 12
.
11 11 10 10 01 01
10 01 11 01 11 10
01 10 01 11 10 11
		

Crossrefs

Programs

  • Mathematica
    prs2mat[prs_]:=Table[Count[prs,{i,j}],{i,Union[First/@prs]},{j,Union[Last/@prs]}];
    multsubs[set_,k_]:=If[k==0,{{}},Join@@Table[Prepend[#,set[[i]]]&/@multsubs[Drop[set,i-1],k-1],{i,Length[set]}]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[Length[Select[multsubs[Tuples[Range[n],2],n],And[Union[First/@#]==Range[Max@@First/@#],Union[Last/@#]==Range[Max@@Last/@#],UnsameQ@@prs2mat[#],UnsameQ@@Transpose[prs2mat[#]],Length[csm[Map[Last,GatherBy[#,First],{2}]]]==1]&]],{n,6}]
  • PARI
    permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
    K(q,t,wf)={prod(j=1, #q, wf(t*q[j]))-1}
    Q(m,n,wf=w->2)={my(s=0); forpart(p=m, s+=(-1)^#p*permcount(p)*exp(-sum(t=1, n, (-1)^t*x^t*K(p,t,wf)/t, O(x*x^n))) ); Vec((-1)^m*serchop(serlaplace(s),1), -n)}
    ConnectedMats(M)={my([m, n]=matsize(M), R=matrix(m, n)); for(m=1, m, for(n=1, n, R[m, n] = M[m, n] - sum(i=1, m-1, sum(j=1, n-1, binomial(m-1, i-1)*binomial(n, j)*R[i, j]*M[m-i, n-j])))); R}
    seq(n)={my(R=vectorv(n,m,Q(m,n,w->1/(1 - y^w) + O(y*y^n)))); for(i=2, #R, R[i] -= i*R[i-1]); Vec(1 + vecsum( vecsum( Vec( ConnectedMats( Mat(R))))))} \\ Andrew Howroyd, Jan 24 2024

Extensions

a(7) onwards from Andrew Howroyd, Jan 24 2024

A322136 Numbers whose number of prime factors counted with multiplicity exceeds half their sum of prime indices by at least 1.

Original entry on oeis.org

4, 8, 12, 16, 24, 32, 36, 40, 48, 64, 72, 80, 96, 108, 112, 120, 128, 144, 160, 192, 216, 224, 240, 256, 288, 320, 324, 336, 352, 360, 384, 400, 432, 448, 480, 512, 576, 640, 648, 672, 704, 720, 768, 800, 832, 864, 896, 960, 972
Offset: 1

Views

Author

Gus Wiseman, Nov 27 2018

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k). The sequence lists all Heinz numbers of integer partitions where the number of parts is at least 1 plus half the sum of parts.
Also Heinz numbers of integer partitions that are the vertex-degrees of some hypertree. We allow no singletons in a hypertree, so 2 is not included.

Examples

			The sequence of partitions with Heinz numbers in the sequence begins: (11), (111), (211), (1111), (2111), (11111), (2211), (3111), (21111), (111111), (22111), (31111), (211111), (22211), (41111), (32111), (1111111).
		

Crossrefs

Programs

  • Mathematica
    Select[Range[1000],PrimeOmega[#]>=(Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]+2)/2&]

A318403 Number of strict connected antichains of sets whose multiset union is an integer partition of n.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 4, 6, 8, 12, 13, 22, 31
Offset: 0

Views

Author

Gus Wiseman, Oct 12 2018

Keywords

Examples

			The a(1) = 1 through a(10) = 13 clutters:
  {{1}}  {{2}}  {{3}}    {{4}}    {{5}}    {{6}}      {{7}}
                {{1,2}}  {{1,3}}  {{1,4}}  {{1,5}}    {{1,6}}
                                  {{2,3}}  {{2,4}}    {{2,5}}
                                           {{1,2,3}}  {{3,4}}
                                                      {{1,2,4}}
                                                      {{1,2},{1,3}}
.
  {{8}}          {{9}}          {{10}}
  {{1,7}}        {{1,8}}        {{1,9}}
  {{2,6}}        {{2,7}}        {{2,8}}
  {{3,5}}        {{3,6}}        {{3,7}}
  {{1,2,5}}      {{4,5}}        {{4,6}}
  {{1,3,4}}      {{1,2,6}}      {{1,2,7}}
  {{1,2},{1,4}}  {{1,3,5}}      {{1,3,6}}
  {{1,2},{2,3}}  {{2,3,4}}      {{1,4,5}}
                 {{1,2},{1,5}}  {{2,3,5}}
                 {{1,2},{2,4}}  {{1,2,3,4}}
                 {{1,3},{1,4}}  {{1,2},{1,6}}
                 {{1,3},{2,3}}  {{1,2},{2,5}}
                                {{1,3},{1,5}}
		

Crossrefs

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    submultisetQ[M_,N_]:=Or[Length[M]==0,MatchQ[{Sort[List@@M],Sort[List@@N]},{{x_,Z___},{_,x_,W___}}/;submultisetQ[{Z},{W}]]];
    antiQ[s_]:=Select[Tuples[s,2],And[UnsameQ@@#,submultisetQ@@#]&]=={};
    Table[Length[Select[Join@@mps/@IntegerPartitions[n],And[UnsameQ@@#,And@@UnsameQ@@@#,Length[csm[#]]==1,antiQ[#]]&]],{n,8}]

A319079 Number of connected antichains of sets whose multiset union is an integer partition of n.

Original entry on oeis.org

1, 1, 2, 3, 4, 4, 8, 7, 12, 15, 19, 26, 43
Offset: 0

Views

Author

Gus Wiseman, Oct 12 2018

Keywords

Examples

			The a(10) = 19 clutters:
  {{10}}
  {{1,9}}
  {{2,8}}
  {{3,7}}
  {{4,6}}
  {{1,2,7}}
  {{1,3,6}}
  {{1,4,5}}
  {{2,3,5}}
  {{1,2,3,4}}
  {{5},{5}}
  {{1,2},{1,6}}
  {{1,2},{2,5}}
  {{1,3},{1,5}}
  {{1,4},{1,4}}
  {{2,3},{2,3}}
  {{1,2},{1,2},{1,3}}
  {{2},{2},{2},{2},{2}}
  {{1},{1},{1},{1},{1},{1},{1},{1},{1},{1}}
		

Crossrefs

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    submultisetQ[M_,N_]:=Or[Length[M]==0,MatchQ[{Sort[List@@M],Sort[List@@N]},{{x_,Z___},{_,x_,W___}}/;submultisetQ[{Z},{W}]]];
    antiQ[s_]:=Select[Tuples[s,2],And[UnsameQ@@#,submultisetQ@@#]&]=={};
    Table[Length[Select[Join@@mps/@IntegerPartitions[n],And[And@@UnsameQ@@@#,Length[csm[#]]==1,antiQ[#]]&]],{n,10}]

A321256 Regular triangle where T(n,k) is the number of non-isomorphic connected set systems of weight n with density -1 <= k <= n-2.

Original entry on oeis.org

1, 1, 0, 2, 0, 0, 4, 0, 0, 0, 6, 1, 0, 0, 0, 13, 5, 0, 0, 0, 0, 23, 12, 2, 0, 0, 0, 0, 49, 36, 11, 0, 0, 0, 0, 0, 100, 95, 39, 5, 0, 0, 0, 0, 0, 220, 262, 143, 32, 1, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Nov 01 2018

Keywords

Comments

A set system is a finite set of finite nonempty sets. The density of a set system is the sum of sizes of each part (weight) minus the number of parts minus the number of vertices.

Examples

			Triangle begins:
    1
    1   0
    2   0   0
    4   0   0   0
    6   1   0   0   0
   13   5   0   0   0   0
   23  12   2   0   0   0   0
   49  36  11   0   0   0   0   0
  100  95  39   5   0   0   0   0   0
  220 262 143  32   1   0   0   0   0   0
		

Crossrefs

First column is A321228. Row sums are A007718.

A321270 Number of connected multiset partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 5, 1, 5, 4, 7, 3, 11, 7, 10, 1, 15, 9, 22, 7, 19, 12, 30, 5, 22, 19, 28, 14, 42, 22, 56, 1, 33, 30, 42, 20, 77, 45
Offset: 1

Views

Author

Gus Wiseman, Nov 01 2018

Keywords

Comments

This multiset (row n of A305936) is generally not the same as the multiset of prime indices of n. For example, the prime indices of 12 are {1,1,2}, while a multiset whose multiplicities are {1,1,2} is {1,1,2,3}.

Examples

			The a(2) = 1 through a(12) = 3 connected multiset partitions:
  {{1}}  {{11}}    {{12}}  {{111}}      {{112}}    {{1111}}
         {{1}{1}}          {{1}{11}}    {{1}{12}}  {{1}{111}}
                           {{1}{1}{1}}             {{11}{11}}
                                                   {{1}{1}{11}}
                                                   {{1}{1}{1}{1}}
.
  {{123}}  {{1122}}      {{1112}}      {{11111}}          {{1123}}
           {{1}{122}}    {{1}{112}}    {{1}{1111}}        {{1}{123}}
           {{12}{12}}    {{11}{12}}    {{11}{111}}        {{12}{13}}
           {{2}{112}}    {{1}{1}{12}}  {{1}{1}{111}}
           {{1}{2}{12}}                {{1}{11}{11}}
                                       {{1}{1}{1}{11}}
                                       {{1}{1}{1}{1}{1}}
The a(18) = 9, a(27) = 28, and a(36) = 20 connected multiset partitions of {1,1,2,2,3}, {1,1,2,2,3,3}, and {1,1,2,2,3,4} respectively:
  {{1,1,2,2,3}}      {{1,1,2,2,3,3}}        {{1,1,2,2,3,4}}
  {{1},{1,2,2,3}}    {{1},{1,2,2,3,3}}      {{1},{1,2,2,3,4}}
  {{1,2},{1,2,3}}    {{1,1,2},{2,3,3}}      {{1,1,2},{2,3,4}}
  {{1,3},{1,2,2}}    {{1,1,3},{2,2,3}}      {{1,2},{1,2,3,4}}
  {{2},{1,1,2,3}}    {{1,2},{1,2,3,3}}      {{1,2,2},{1,3,4}}
  {{2,3},{1,1,2}}    {{1,2,2},{1,3,3}}      {{1,2,3},{1,2,4}}
  {{1},{1,2},{2,3}}  {{1,2,3},{1,2,3}}      {{1,3},{1,2,2,4}}
  {{1},{2},{1,2,3}}  {{1,3},{1,2,2,3}}      {{1,4},{1,2,2,3}}
  {{2},{1,2},{1,3}}  {{2},{1,1,2,3,3}}      {{2},{1,1,2,3,4}}
                     {{2,3},{1,1,2,3}}      {{2,3},{1,1,2,4}}
                     {{3},{1,1,2,2,3}}      {{2,4},{1,1,2,3}}
                     {{1},{1,2},{2,3,3}}    {{1},{1,2},{2,3,4}}
                     {{1},{1,3},{2,2,3}}    {{1},{2},{1,2,3,4}}
                     {{1},{2},{1,2,3,3}}    {{1,2},{1,3},{2,4}}
                     {{1,2},{1,3},{2,3}}    {{1,2},{1,4},{2,3}}
                     {{1},{2,3},{1,2,3}}    {{1},{2,3},{1,2,4}}
                     {{1},{3},{1,2,2,3}}    {{1},{2,4},{1,2,3}}
                     {{2},{1,2},{1,3,3}}    {{2},{1,2},{1,3,4}}
                     {{2},{1,3},{1,2,3}}    {{2},{1,3},{1,2,4}}
                     {{2},{2,3},{1,1,3}}    {{2},{1,4},{1,2,3}}
                     {{2},{3},{1,1,2,3}}
                     {{3},{1,2},{1,2,3}}
                     {{3},{1,3},{1,2,2}}
                     {{3},{2,3},{1,1,2}}
                     {{1},{2},{1,3},{2,3}}
                     {{1},{2},{3},{1,2,3}}
                     {{1},{3},{1,2},{2,3}}
                     {{2},{3},{1,2},{1,3}}
		

Crossrefs

A321584 Number of connected (0,1)-matrices with n ones and no zero rows or columns.

Original entry on oeis.org

1, 1, 2, 6, 27, 159, 1154, 9968, 99861, 1138234, 14544650, 205927012, 3199714508, 54131864317, 990455375968, 19488387266842, 410328328297512, 9205128127109576, 219191041679766542, 5521387415218119528, 146689867860276432637, 4099255234885039058842, 120199458455807733040338
Offset: 0

Views

Author

Gus Wiseman, Nov 13 2018

Keywords

Comments

A matrix is connected if the positions in each row (or each column) of the nonzero entries form a connected hypergraph.

Examples

			The a(4) = 27 matrices:
  [1111]
.
  [111][111][111][11][110][110][101][101][100][011][011][010][001]
  [100][010][001][11][101][011][110][011][111][110][101][111][111]
.
  [11][11][11][11][10][10][10][10][01][01][01][01]
  [10][10][01][01][11][11][10][01][11][11][10][01]
  [10][01][10][01][10][01][11][11][10][01][11][11]
.
  [1]
  [1]
  [1]
  [1]
		

Crossrefs

Programs

  • Mathematica
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[Length[Select[Subsets[Tuples[Range[n],2],{n}],And[Union[First/@#]==Range[Max@@First/@#],Union[Last/@#]==Range[Max@@Last/@#],Length[csm[Map[Last,GatherBy[#,First],{2}]]]==1]&]],{n,6}] (* Mathematica 7.0+ *)
  • PARI
    NonZeroCols(M)={my(C=Vec(M)); Mat(vector(#C, n, sum(k=1, n, (-1)^(n-k)*binomial(n,k)*C[k])))}
    ConnectedMats(M)={my([m,n]=matsize(M), R=matrix(m,n)); for(m=1, m, for(n=1, n, R[m,n] = M[m,n] - sum(i=1, m-1, sum(j=1, n-1, binomial(m-1,i-1)*binomial(n,j)*R[i,j]*M[m-i,n-j])))); R}
    seq(n)={my(M=matrix(n,n,i,j,sum(k=1, n, binomial(i*j,k)*x^k, O(x*x^n) ))); Vec(1 + vecsum(vecsum(Vec( ConnectedMats( NonZeroCols( NonZeroCols(M)~)) ))))} \\ Andrew Howroyd, Jan 17 2024

Extensions

a(7) onwards from Andrew Howroyd, Jan 17 2024

A321680 Number of non-isomorphic weight-n connected antichains (not necessarily strict) of multisets with multiset density -1.

Original entry on oeis.org

1, 1, 3, 4, 9, 14, 39, 80, 216, 538, 1460
Offset: 0

Views

Author

Gus Wiseman, Nov 16 2018

Keywords

Comments

The multiset density of a multiset partition is the sum of the numbers of distinct vertices in each part minus the number of parts minus the number of vertices.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(5) = 14 multiset trees:
  {{1}}  {{1,1}}    {{1,1,1}}      {{1,1,1,1}}        {{1,1,1,1,1}}
         {{1,2}}    {{1,2,2}}      {{1,1,2,2}}        {{1,1,2,2,2}}
         {{1},{1}}  {{1,2,3}}      {{1,2,2,2}}        {{1,2,2,2,2}}
                    {{1},{1},{1}}  {{1,2,3,3}}        {{1,2,2,3,3}}
                                   {{1,2,3,4}}        {{1,2,3,3,3}}
                                   {{1,1},{1,1}}      {{1,2,3,4,4}}
                                   {{1,2},{2,2}}      {{1,2,3,4,5}}
                                   {{1,3},{2,3}}      {{1,1},{1,2,2}}
                                   {{1},{1},{1},{1}}  {{1,2},{2,2,2}}
                                                      {{1,2},{2,3,3}}
                                                      {{1,3},{2,3,3}}
                                                      {{1,4},{2,3,4}}
                                                      {{3,3},{1,2,3}}
                                                      {{1},{1},{1},{1},{1}}
		

Crossrefs

A321681 Number of non-isomorphic weight-n connected strict antichains of multisets with multiset density -1.

Original entry on oeis.org

1, 1, 2, 3, 7, 13, 35, 77, 205, 517, 1399
Offset: 0

Views

Author

Gus Wiseman, Nov 16 2018

Keywords

Comments

The multiset density of a multiset partition is the sum of the numbers of distinct vertices in each part minus the number of parts minus the number of vertices.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(5) = 13 trees:
  {{1}}  {{1,1}}  {{1,1,1}}  {{1,1,1,1}}    {{1,1,1,1,1}}
         {{1,2}}  {{1,2,2}}  {{1,1,2,2}}    {{1,1,2,2,2}}
                  {{1,2,3}}  {{1,2,2,2}}    {{1,2,2,2,2}}
                             {{1,2,3,3}}    {{1,2,2,3,3}}
                             {{1,2,3,4}}    {{1,2,3,3,3}}
                             {{1,2},{2,2}}  {{1,2,3,4,4}}
                             {{1,3},{2,3}}  {{1,2,3,4,5}}
                                            {{1,1},{1,2,2}}
                                            {{1,2},{2,2,2}}
                                            {{1,2},{2,3,3}}
                                            {{1,3},{2,3,3}}
                                            {{1,4},{2,3,4}}
                                            {{3,3},{1,2,3}}
		

Crossrefs

A368096 Triangle read by rows where T(n,k) is the number of non-isomorphic set-systems of length k and weight n.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 4, 3, 1, 0, 1, 5, 8, 3, 1, 0, 1, 8, 18, 13, 3, 1, 0, 1, 9, 32, 37, 15, 3, 1, 0, 1, 13, 55, 96, 59, 16, 3, 1, 0, 1, 14, 91, 209, 196, 74, 16, 3, 1, 0, 1, 19, 138, 449, 573, 313, 82, 16, 3, 1, 0, 1, 20, 206, 863, 1529, 1147, 403, 84, 16, 3, 1
Offset: 0

Views

Author

Gus Wiseman, Dec 28 2023

Keywords

Comments

A set-system is a finite set of finite nonempty sets.
Conjecture: Column k = 2 is A101881.

Examples

			Triangle begins:
   1
   0   1
   0   1   1
   0   1   2   1
   0   1   4   3   1
   0   1   5   8   3   1
   0   1   8  18  13   3   1
   0   1   9  32  37  15   3   1
   0   1  13  55  96  59  16   3   1
   0   1  14  91 209 196  74  16   3   1
   0   1  19 138 449 573 313  82  16   3   1
   ...
Non-isomorphic representatives of the set-systems counted in row n = 5:
  .  {12345}  {1}{1234}  {1}{2}{123}  {1}{2}{3}{12}  {1}{2}{3}{4}{5}
              {1}{2345}  {1}{2}{134}  {1}{2}{3}{14}
              {12}{123}  {1}{2}{345}  {1}{2}{3}{45}
              {12}{134}  {1}{12}{13}
              {12}{345}  {1}{12}{23}
                         {1}{12}{34}
                         {1}{23}{24}
                         {1}{23}{45}
		

Crossrefs

Row sums are A283877, connected case A300913.
For multiset partitions we have A317533.
Counting connected components instead of edges gives A321194.
For set multipartitions we have A334550.
For strict multiset partitions we have A368099.
A000110 counts set-partitions, non-isomorphic A000041.
A003465 counts covering set-systems, unlabeled A055621.
A007716 counts non-isomorphic multiset partitions, connected A007718.
A049311 counts non-isomorphic set multipartitions, connected A056156.
A058891 counts set-systems, unlabeled A000612, connected A323818.
A316980 counts non-isomorphic strict multiset partitions, connected A319557.
A319559 counts non-isomorphic T_0 set-systems, connected A319566.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]] /@ Cases[Subsets[set],{i,_}];
    mpm[n_]:=Join@@Table[Union[Sort[Sort /@ (#/.x_Integer:>s[[x]])]& /@ sps[Range[n]]],{s,Flatten[MapIndexed[Table[#2,{#1}]&,#]]& /@ IntegerPartitions[n]}];
    brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{i,p[[i]]},{i,Length[p]}])], {p,Permutations[Union@@m]}]]];
    Table[Length[Union[brute /@ Select[mpm[n],UnsameQ@@#&&And@@UnsameQ@@@#&&Length[#]==k&]]], {n,0,5},{k,0,n}]
  • PARI
    WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, (-1)^(n-1)/n))))-1, -#v)}
    permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
    K(q, t, k)={WeighT(Vec(sum(j=1, #q, my(g=gcd(t, q[j])); g*x^(q[j]/g)) + O(x*x^k), -k))}
    G(n)={my(s=0); forpart(q=n, my(p=sum(t=1, n, y^t*subst(x*Ser(K(q, t, n\t))/t, x, x^t))); s+=permcount(q)*exp(p-subst(subst(p, x, x^2), y, y^2))); s/n!}
    T(n)={[Vecrev(p) | p <- Vec(G(n))]}
    { my(A=T(10)); for(n=1, #A, print(A[n])) } \\ Andrew Howroyd, Jan 11 2024

Extensions

Terms a(66) and beyond from Andrew Howroyd, Jan 11 2024
Previous Showing 21-30 of 44 results. Next