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 31-40 of 51 results. Next

A322436 Number of pairs of factorizations of n into factors > 1 where no factor of the second properly divides any factor of the first.

Original entry on oeis.org

1, 1, 1, 3, 1, 3, 1, 5, 3, 3, 1, 8, 1, 3, 3, 11, 1, 8, 1, 8, 3, 3, 1, 18, 3, 3, 5, 8, 1, 12, 1, 15, 3, 3, 3, 31, 1, 3, 3, 18, 1, 12, 1, 8, 8, 3, 1, 39, 3, 8, 3, 8, 1, 18, 3, 18, 3, 3, 1, 42, 1, 3, 8, 33, 3, 12, 1, 8, 3, 12, 1, 67, 1, 3, 8, 8, 3, 12, 1, 39, 11
Offset: 1

Views

Author

Gus Wiseman, Dec 08 2018

Keywords

Examples

			The a(12) = 8 pairs of factorizations:
  (2*2*3)|(2*2*3)
  (2*2*3)|(2*6)
  (2*2*3)|(3*4)
  (2*2*3)|(12)
    (2*6)|(12)
    (3*4)|(3*4)
    (3*4)|(12)
     (12)|(12)
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    divpropQ[x_,y_]:=And[x!=y,Divisible[x,y]];
    Table[Length[Select[Tuples[facs[n],2],!Or@@divpropQ@@@Tuples[#]&]],{n,100}]

A322438 Number of unordered pairs of factorizations of n into factors > 1 where no factor of one properly divides any factor of the other.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 4
Offset: 1

Views

Author

Gus Wiseman, Dec 08 2018

Keywords

Comments

First differs from A322437 at a(144) = 4, A322437(144) = 3.
First differs from A379958 at a(120) = 2, A379958(120) = 1.

Examples

			The a(240) = 5 pairs of factorizations::
  (4*4*15)|(4*6*10)
    (6*40)|(15*16)
    (8*30)|(12*20)
   (10*24)|(15*16)
   (12*20)|(15*16)
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    divpropQ[x_,y_]:=And[x!=y,Divisible[x,y]];
    Table[Length[Select[Subsets[facs[n],{2}],And[!Or@@divpropQ@@@Tuples[#],!Or@@divpropQ@@@Reverse/@Tuples[#]]&]],{n,100}]
  • PARI
    factorizations(n, m=n, f=List([]), z=List([])) = if(1==n, listput(z,Vec(f)); z, my(newf); fordiv(n, d, if((d>1)&&(d<=m), newf = List(f); listput(newf,d); z = factorizations(n/d, d, newf, z))); (z));
    is_proper_ndf_pair(fac1,fac2) = { for(i=1,#fac1,for(j=1,#fac2,if((fac1[i]!=fac2[j]) && (!(fac1[i]%fac2[j]) || !(fac2[j]%fac1[i])),return(0)))); (1); };
    number_of_proper_ndfpairs(z) = sum(i=1,#z,sum(j=i+1,#z,is_proper_ndf_pair(z[i],z[j])));
    A322438(n) = number_of_proper_ndfpairs(Vec(factorizations(n))); \\ Antti Karttunen, Jan 24 2025

Extensions

Data section extended up to a(144) by Antti Karttunen, Jan 24 2025

A320802 Number of non-isomorphic aperiodic multiset partitions of weight n whose dual is also an aperiodic multiset partition.

Original entry on oeis.org

1, 1, 2, 8, 26, 89, 274, 908, 2955, 9926, 34021, 119367, 428612, 1574222, 5914324, 22699632, 88997058, 356058538, 1453059643, 6044132792, 25612530061, 110503625785, 485161109305, 2166488899640, 9835209048655, 45370059225137, 212582814591083, 1011306624492831
Offset: 0

Views

Author

Gus Wiseman, Nov 06 2018

Keywords

Comments

Also the number of nonnegative integer matrices with sum of entries equal to n and no zero rows or columns where the multiset of rows and the multiset of columns are both aperiodic, up to row and column permutations.
A multiset is aperiodic if its multiplicities are relatively prime.
The dual of a multiset partition has, for each vertex, one part consisting of the indices (or positions) of the parts containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,2}} is {{1},{1,2,2}}.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.
Also the number of non-isomorphic aperiodic multiset partitions of weight n whose parts have relatively prime periods, where the period of a multiset is the GCD of its multiplicities.

Examples

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

Crossrefs

Formula

Second Moebius transform of A007716, or Moebius transform of A303546, where the Moebius transform of a sequence b is a(n) = Sum_{d|n} mu(d) * b(n/d).

Extensions

a(26)-a(27) from Jinyuan Wang, Jun 27 2020

A320803 Number of non-isomorphic multiset partitions of weight n in which all parts are aperiodic multisets.

Original entry on oeis.org

1, 1, 3, 7, 21, 56, 174, 517, 1664, 5383, 18199, 62745, 223390, 813425, 3040181, 11620969, 45446484, 181537904, 740369798, 3079779662, 13059203150, 56406416004, 248027678362, 1109626606188, 5048119061134, 23342088591797, 109648937760252, 523036690273237
Offset: 0

Views

Author

Gus Wiseman, Nov 06 2018

Keywords

Comments

A multiset is aperiodic if its multiplicities are relatively prime.
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(4) = 21 multiset partitions with aperiodic parts:
  {{1}}  {{1,2}}    {{1,2,2}}      {{1,2,2,2}}
         {{1},{1}}  {{1,2,3}}      {{1,2,3,3}}
         {{1},{2}}  {{1},{2,3}}    {{1,2,3,4}}
                    {{2},{1,2}}    {{1},{1,2,2}}
                    {{1},{1},{1}}  {{1,2},{1,2}}
                    {{1},{2},{2}}  {{1},{2,3,3}}
                    {{1},{2},{3}}  {{1},{2,3,4}}
                                   {{1,2},{3,4}}
                                   {{1,3},{2,3}}
                                   {{2},{1,2,2}}
                                   {{3},{1,2,3}}
                                   {{1},{1},{2,3}}
                                   {{1},{2},{1,2}}
                                   {{1},{2},{3,4}}
                                   {{1},{3},{2,3}}
                                   {{2},{2},{1,2}}
                                   {{1},{1},{1},{1}}
                                   {{1},{1},{2},{2}}
                                   {{1},{2},{2},{2}}
                                   {{1},{2},{3},{3}}
                                   {{1},{2},{3},{4}}
		

Crossrefs

Programs

  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, 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)={EulerT(Vec(sum(j=1, #q, gcd(t, q[j])*x^lcm(t, q[j])) + O(x*x^k), -k))}
    a(n)={if(n==0, 1, my(mbt=vector(n, d, moebius(d)), s=0); forpart(q=n, s+=permcount(q)*polcoef(exp(x*Ser(dirmul(mbt, sum(t=1, n, K(q, t, n)/t)))), n)); s/n!)} \\ Andrew Howroyd, Jan 16 2023

Extensions

Terms a(11) and beyond from Andrew Howroyd, Jan 16 2023

A320809 Number of non-isomorphic multiset partitions of weight n in which each part and each part of the dual, as well as the multiset union of the parts, is an aperiodic multiset.

Original entry on oeis.org

1, 1, 2, 5, 13, 40, 99, 344, 985, 3302, 10583
Offset: 0

Views

Author

Gus Wiseman, Nov 07 2018

Keywords

Comments

Also the number of nonnegative integer matrices up to row and column permutations with sum of elements equal to n and no zero rows or columns, in which (1) the positive entries in each row and column are relatively prime and (2) the column sums are relatively prime.
The dual of a multiset partition has, for each vertex, one part consisting of the indices (or positions) of the parts containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,2}} is {{1},{1,2,2}}.
A multiset is aperiodic if its multiplicities are relatively prime.
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(4) = 13 multiset partitions:
  {{1}}  {{1,2}}    {{1,2,3}}      {{1,2,3,4}}
         {{1},{2}}  {{1},{2,3}}    {{1},{2,3,4}}
                    {{2},{1,2}}    {{1,2},{3,4}}
                    {{1},{2},{2}}  {{1,3},{2,3}}
                    {{1},{2},{3}}  {{2},{1,2,2}}
                                   {{3},{1,2,3}}
                                   {{1},{1},{2,3}}
                                   {{1},{2},{3,4}}
                                   {{1},{3},{2,3}}
                                   {{2},{2},{1,2}}
                                   {{1},{2},{2},{2}}
                                   {{1},{2},{3},{3}}
                                   {{1},{2},{3},{4}}
		

Crossrefs

A321390 Third Moebius transform of A007716. Number of non-isomorphic aperiodic multiset partitions of weight n whose parts have relatively prime periods and whose dual is also an aperiodic multiset partition.

Original entry on oeis.org

1, 1, 1, 7, 24, 88, 265, 907, 2929, 9918, 33931, 119366, 428314, 1574221, 5913415, 22699536, 88994103, 356058537, 1453049451, 6044132791, 25612496016, 110503624870, 485160989937, 2166488899639, 9835208617114, 45370059225048
Offset: 0

Views

Author

Gus Wiseman, Nov 08 2018

Keywords

Comments

The Moebius transform c of a sequence b is c(n) = Sum_{d|n} mu(d) * b(n/d).
Also the number of nonnegative integer matrices with sum of entries equal to n and no zero rows or columns where the multiset of rows and the multiset of columns are both aperiodic and the nonzero entries are relatively prime, up to row and column permutations.
A multiset is aperiodic if its multiplicities are relatively prime. The period of a multiset is the GCD of its multiplicities.
The dual of a multiset partition has, for each vertex, one part consisting of the indices (or positions) of the parts containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,2}} is {{1},{1,2,2}}.
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(4) = 24 multiset partitions:
  {{1}}  {{1},{2}}  {{1,2,2}}      {{1,2,2,2}}
                    {{1},{1,1}}    {{1,2,3,3}}
                    {{1},{2,2}}    {{1},{1,1,1}}
                    {{1},{2,3}}    {{1},{1,2,2}}
                    {{2},{1,2}}    {{1},{2,2,2}}
                    {{1},{2},{2}}  {{1,2},{2,2}}
                    {{1},{2},{3}}  {{1},{2,3,3}}
                                   {{1,2},{3,3}}
                                   {{1},{2,3,4}}
                                   {{1,3},{2,3}}
                                   {{2},{1,2,2}}
                                   {{3},{1,2,3}}
                                   {{1},{1},{1,1}}
                                   {{1},{1},{2,2}}
                                   {{1},{1},{2,3}}
                                   {{1},{2},{1,2}}
                                   {{1},{2},{2,2}}
                                   {{1},{2},{3,3}}
                                   {{1},{2},{3,4}}
                                   {{1},{3},{2,3}}
                                   {{2},{2},{1,2}}
                                   {{1},{2},{2},{2}}
                                   {{1},{2},{3},{3}}
                                   {{1},{2},{3},{4}}
		

Crossrefs

A108572 Number of partitions of n which, as multisets, are nontrivial repetitions of a multiset.

Original entry on oeis.org

0, 0, 0, 1, 0, 3, 0, 4, 2, 7, 0, 13, 0, 15, 8, 21, 0, 37, 0, 44, 16, 56, 0, 93, 6, 101, 29, 137, 0, 217, 0, 230, 57, 297, 20, 450, 0, 490, 102, 643, 0, 918, 0, 1004, 202, 1255, 0, 1783, 14, 1992, 298, 2438, 0, 3364, 61, 3734, 491, 4565, 0, 6251, 0, 6842, 818
Offset: 1

Views

Author

Len Smiley, Jul 25 2005

Keywords

Comments

The singleton and the all-ones partitions are ignored, so that a(n)=0 if n is prime. If a partition is listed as m_1^am_2^bm_3^c..., then it is counted exactly when gcd(a,b,c,...)>1. These are equinumerous (conjugate) with those partitions for which gcd(m_1,m_2,...)>1 (less 1, the singleton), hence the formula.

Examples

			a(25) = 6: 1^(15)2^5 = 5{1, 1, 1, 2}, 1^52^(10) = 5{1, 2, 2}, 1^(10)3^5 = 5{3, 1, 1}, 2^53^5 = 5{3, 2}, 1^44^4 = 5{4, 1}, 5^5 = 5{5}.
Note that A000041(25)=P(25)=1958, only 6 of which satisfy the criterion.
		

Crossrefs

Programs

  • Maple
    with(combinat):PartMulti:=proc(n::nonnegint) local count,a,i,j,b,m,k,part_vec;
    bigcount:=0; if isprime(n) then return(bigcount) else ps:=partition(n); b:=nops(ps);
    for m from 2 to b-1 do p:=ps[m]; a:=nops(p); part_vec:=array(1..n);
    for k from 1 to n do part_vec[k]:=0 od;
    for i from 1 to a do j:=p[i]; part_vec[j]:=part_vec[j]+1 od;
    g:=0; for j from 1 to n do g:=igcd(g,part_vec[j]) od;
    if g>1 then bigcount:=bigcount+1 fi od; return(bigcount) end if end proc;
    seq(PartMulti(q),q=1..49);
  • Mathematica
    Table[Length[Select[IntegerPartitions[n],And[Length[#]1]&]],{n,20}] (* Gus Wiseman, Dec 06 2018 *)

Formula

a(n) = A018783(n)-1, n>1. - Vladeta Jovovic, Jul 28 2005

Extensions

More terms from Gus Wiseman, Dec 06 2018

A305253 Number of connected factorizations of n into factors greater than 1 whose distinct factors are pairwise indivisible.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, May 28 2018

Keywords

Comments

Given a finite multiset S of positive integers greater than one, let G(S) be the simple labeled graph with vertex set S and edges between any two vertices with a common divisor greater than 1. For example, G({6,14,15,35}) is a 4-cycle. This sequence counts factorizations S whose distinct factors are pairwise indivisible and such that G(S) is a connected graph.

Examples

			The a(360) = 8 factorizations: (360), (4*90), (10*36), (12*30), (15*24), (18*20), (4*6*15), (6*6*10).
		

Crossrefs

Programs

  • Mathematica
    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]]]]]]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    sacs[n_]:=Select[facs[n],Function[f,Length[zsm[f]]==1&&Select[Tuples[Union[f],2],UnsameQ@@#&&Divisible@@#&]=={}]]
    Table[Length[sacs[n]],{n,500}]
  • PARI
    is_connected(facs) = { my(siz=length(facs)); if(1==siz,1,my(m=matrix(siz,siz,i,j,(gcd(facs[i],facs[j])!=1))^siz); for(n=1,siz,if(0==vecmin(m[n,]),return(0))); (1)); };
    A305253aux(n, m, facs) = if(1==n, is_connected(Vec(facs)), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m)&&factorback(apply(x -> (x==d)||(x%d),Vec(facs))), newfacs = List(facs); listput(newfacs,d); s += A305253aux(n/d, d, newfacs))); (s));
    A305253(n) = if(1==n,0,A305253aux(n, n, List([]))); \\ Antti Karttunen, Dec 06 2018

Formula

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

Extensions

Definition clarified by Gus Wiseman, more terms from Antti Karttunen, Dec 06 2018

A317852 Number of plane trees with n nodes where the sequence of branches directly under any given node is aperiodic, meaning its cyclic permutations are all different.

Original entry on oeis.org

1, 1, 1, 3, 8, 26, 76, 247, 783, 2565, 8447, 28256, 95168, 323720, 1108415, 3821144, 13246307, 46158480, 161574043, 567925140, 2003653016, 7092953340, 25186731980, 89690452750, 320221033370, 1146028762599, 4110596336036, 14774346783745, 53203889807764, 191934931634880
Offset: 1

Views

Author

Gus Wiseman, Sep 05 2018

Keywords

Comments

Also the number of plane trees with n nodes where the sequence of branches directly under any given node has relatively prime run-lengths.

Examples

			The a(5) = 8 locally aperiodic plane trees:
  ((((o)))),
  (((o)o)), ((o(o))), (((o))o), (o((o))),
  ((o)oo), (o(o)o), (oo(o)).
The a(6) = 26 locally aperiodic plane trees:
  (((((o)))))  ((((o)o)))  (((o)oo))  ((o)ooo)
               (((o(o))))  ((o(o)o))  (o(o)oo)
               ((((o))o))  ((oo(o)))  (oo(o)o)
               ((o((o))))  (((o)o)o)  (ooo(o))
               ((((o)))o)  ((o(o))o)
               (o(((o))))  (o((o)o))
               (((o))(o))  (o(o(o)))
               ((o)((o)))  (((o))oo)
                           (o((o))o)
                           (oo((o)))
                           ((o)(o)o)
                           ((o)o(o))
                           (o(o)(o))
		

Crossrefs

Programs

  • Mathematica
    aperQ[q_]:=Array[RotateRight[q,#]&,Length[q],1,UnsameQ];
    aperplane[n_]:=If[n==1,{{}},Join@@Table[Select[Tuples[aperplane/@c],aperQ],{c,Join@@Permutations/@IntegerPartitions[n-1]}]];
    Table[Length[aperplane[n]],{n,10}]
  • PARI
    Tfm(p, n)={sum(d=1, n, moebius(d)*(subst(1/(1+O(x*x^(n\d))-p), x, x^d)-1))}
    seq(n)={my(p=O(1)); for(i=1, n, p=1+Tfm(x*p, i)); Vec(p)} \\ Andrew Howroyd, Feb 08 2020

Extensions

a(16)-a(17) from Robert Price, Sep 15 2018
Terms a(18) and beyond from Andrew Howroyd, Feb 08 2020

A327501 Maximum divisor of n that is 1 or not a perfect power.

Original entry on oeis.org

1, 2, 3, 2, 5, 6, 7, 2, 3, 10, 11, 12, 13, 14, 15, 2, 17, 18, 19, 20, 21, 22, 23, 24, 5, 26, 3, 28, 29, 30, 31, 2, 33, 34, 35, 18, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 7, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 2, 65, 66, 67, 68, 69
Offset: 1

Views

Author

Gus Wiseman, Sep 16 2019

Keywords

Comments

First differs from A052410 at a(36) = 18, A052410(36) = 6.
The number of divisors that are 1 or not a perfect power is given by A327502.
A multiset is aperiodic if its multiplicities are relatively prime. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). Heinz numbers of aperiodic multisets are numbers that are not perfect powers (A007916).
a(n) = n iff n is in A175082. - Bernard Schott, Sep 20 2019

Examples

			The divisors of 36 that are not perfect powers are {1, 2, 3, 6, 12, 18}, so a(36) = 18.
		

Crossrefs

See link for additional cross-references.

Programs

  • Magma
    [1] cat [Max([d:d in Divisors(n)| d gt 1 and not IsPower(d)]):n in [2..70]]; // Marius A. Burtea, Sep 20 2019
  • Mathematica
    Table[Max[Select[Divisors[n],GCD@@Last/@FactorInteger[#]==1&]],{n,100}]
  • PARI
    isp(n) = !ispower(n) && (n>1); \\ A007916
    a(n) = if (n==1, 1, vecmax(select(x->isp(x), divisors(n)))); \\ Michel Marcus, Sep 18 2019
    
Previous Showing 31-40 of 51 results. Next