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-8 of 8 results.

A181796 a(n) = number of divisors of n whose canonical prime factorizations contain no repeated positive exponents (cf. A130091).

Original entry on oeis.org

1, 2, 2, 3, 2, 3, 2, 4, 3, 3, 2, 5, 2, 3, 3, 5, 2, 5, 2, 5, 3, 3, 2, 7, 3, 3, 4, 5, 2, 4, 2, 6, 3, 3, 3, 7, 2, 3, 3, 7, 2, 4, 2, 5, 5, 3, 2, 9, 3, 5, 3, 5, 2, 7, 3, 7, 3, 3, 2, 7, 2, 3, 5, 7, 3, 4, 2, 5, 3, 4, 2, 10, 2, 3, 5, 5, 3, 4, 2, 9, 5, 3, 2, 7, 3, 3, 3, 7, 2, 7, 3, 5, 3, 3, 3, 11, 2, 5, 5, 7, 2, 4, 2, 7, 4
Offset: 1

Views

Author

Matthew Vandermast, Nov 22 2010

Keywords

Comments

The canonical factorization of n into prime powers can be written as Product p(i)^e(i), for example. A host of equivalent notations can also be used (for another example, see Weisstein link). a(n) depends only on prime signature of n (cf. A025487).
a(n) >= A085082(n). (A085082(n) equals the number of members of A025487 that divide A046523(n), and each member of A025487 is divisible by at least one member of A130091 that divides no smaller member of A025487.) a(n) > A085082(n) iff n has in its canonical prime factorization at least two exponents greater than 1.
a(n) = number of such divisors of n that in their prime factorization all exponents are unique. - Antti Karttunen, May 27 2017
First differs from A335549 at a(90) = 7, A335549(90) = 8. First differs from A335516 at a(180) = 9, A335516(180) = 10. - Gus Wiseman, Jun 28 2020

Examples

			12 has a total of six divisors (1, 2, 3, 4, 6 and 12). Of those divisors, the number 1 has no prime factors, hence, no positive exponents at all (and no repeated positive exponents) in its canonical prime factorization. The lists of positive exponents for 2, 3, 4, 6 and 12 are (1), (1), (2), (1,1) and (2,1) respectively (cf. A124010). Of all six divisors, only the number 6 (2^1*3^1) has at least one positive exponent repeated (namely, 1). The other five do not; hence, a(12) = 5.
For n = 90 = 2 * 3^2 * 5, the divisors that satisfy the condition are: 1, 2, 3, 3^2, 5, 2 * 3^2, 3^2 * 5, altogether 7, (but for example 90 itself is not included), thus a(90) = 7.
		

Crossrefs

Diverges from A088873 at n=24 and from A085082 at n=36. a(36) = 7, while A085082(36) = 6.
Partitions with distinct multiplicities are A098859.
Sorted prime signature is A118914.
Unsorted prime signature is A124010.
a(n) is the number of divisors of n in A130091.
Factorizations with distinct multiplicities are A255231.
The largest of the counted divisors is A327498.
Factorizations using the counted divisors are A327523.

Programs

  • Mathematica
    Table[DivisorSum[n, 1 &, Length@ Union@ # == Length@ # &@ FactorInteger[#][[All, -1]] &], {n, 105}] (* Michael De Vlieger, May 28 2017 *)
  • PARI
    no_repeated_exponents(n) = { my(es = factor(n)[, 2]); if(length(Set(es)) == length(es),1,0); }
    A181796(n) = sumdiv(n,d,no_repeated_exponents(d)); \\ Antti Karttunen, May 27 2017
    
  • Python
    from sympy import factorint, divisors
    def ok(n):
        f=factorint(n)
        ex=[f[i] for i in f]
        for i in ex:
            if ex.count(i)>1: return 0
        return 1
    def a(n): return sum([1 for i in divisors(n) if ok(i)]) # Indranil Ghosh, May 27 2017

Formula

a(A000079(n)) = a(A002110(n)) = n+1.
a(A006939(n)) = A000110(n+1).
a(A181555(n)) = A002720(n).

A336866 Number of integer partitions of n without all distinct multiplicities.

Original entry on oeis.org

0, 0, 0, 1, 1, 2, 4, 5, 9, 15, 21, 28, 46, 56, 80, 114, 149, 192, 269, 337, 455, 584, 751, 943, 1234, 1527, 1944, 2422, 3042, 3739, 4699, 5722, 7100, 8668, 10634, 12880, 15790, 19012, 23093, 27776, 33528, 40102, 48264, 57469, 68793, 81727, 97372, 115227
Offset: 0

Views

Author

Gus Wiseman, Aug 09 2020

Keywords

Examples

			The a(0) = 0 through a(9) = 15 partitions (empty columns shown as dots):
  .  .  .  (21)  (31)  (32)  (42)    (43)    (53)     (54)
                       (41)  (51)    (52)    (62)     (63)
                             (321)   (61)    (71)     (72)
                             (2211)  (421)   (431)    (81)
                                     (3211)  (521)    (432)
                                             (3221)   (531)
                                             (3311)   (621)
                                             (4211)   (3321)
                                             (32111)  (4221)
                                                      (4311)
                                                      (5211)
                                                      (32211)
                                                      (42111)
                                                      (222111)
                                                      (321111)
		

Crossrefs

A098859 counts the complement.
A130092 gives the Heinz numbers of these partitions.
A001222 counts prime factors with multiplicity.
A013929 lists nonsquarefree numbers.
A047966 counts uniform partitions.
A047967 counts non-strict partitions.
A071625 counts distinct prime multiplicities.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A327498 gives the maximum divisor with distinct prime multiplicities.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],!UnsameQ@@Length/@Split[#]&]],{n,0,30}]

Formula

a(n) = A000041(n) - A098859(n).

A327498 Maximum divisor of n whose prime multiplicities are distinct (A130091).

Original entry on oeis.org

1, 2, 3, 4, 5, 3, 7, 8, 9, 5, 11, 12, 13, 7, 5, 16, 17, 18, 19, 20, 7, 11, 23, 24, 25, 13, 27, 28, 29, 5, 31, 32, 11, 17, 7, 18, 37, 19, 13, 40, 41, 7, 43, 44, 45, 23, 47, 48, 49, 50, 17, 52, 53, 54, 11, 56, 19, 29, 59, 20, 61, 31, 63, 64, 13, 11, 67, 68, 23
Offset: 1

Views

Author

Gus Wiseman, Sep 16 2019

Keywords

Comments

A number's prime multiplicities are also called its (unsorted) prime signature.
Every positive integer appears a finite number of times in the sequence; a prime p occurs 2^(PrimePi(p) - 1) times. - David A. Corneth, Sep 17 2019

Examples

			The divisors of 60 whose prime multiplicities are distinct are {1, 2, 3, 4, 5, 12, 20}, so a(60) = 20, the largest of these divisors.
		

Crossrefs

See link for additional cross-references.

Programs

  • Mathematica
    Table[Max[Select[Divisors[n],UnsameQ@@Last/@FactorInteger[#]&]],{n,100}]
  • PARI
    a(n) = {my(m = Map(), f = factor(n), res = 1); forstep(i = #f~, 1, -1, forstep(j = f[i, 2], 1, -1, if(!mapisdefined(m, j), mapput(m, j, j); res*=f[i, 1]^j; next(2)))); res} \\ David A. Corneth, Sep 17 2019
    
  • PARI
    A351564(n) = issquarefree(factorback(apply(e->prime(e),(factor(n)[,2]))));
    A327498(n) = fordiv(n,d,if(A351564(n/d), return(n/d))); \\ Antti Karttunen, Apr 02 2022

Formula

a(A130091(n)) = n and a(A130092(n)) < n. - Ivan N. Ianakiev, Sep 17 2019
a(n) = n / A327499(n). - Antti Karttunen, Apr 02 2022

A259936 Number of ways to express the integer n as a product of its unitary divisors (A034444).

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, 2, 1, 2, 1, 2, 1, 5, 1, 1, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 5, 1, 2, 2, 1, 2, 5, 1, 2, 2, 5, 1, 2, 1, 2, 2, 2, 2, 5, 1, 2, 1, 2, 1, 5, 2, 2, 2, 2, 1, 5, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 5
Offset: 1

Views

Author

Geoffrey Critzer, Jul 09 2015

Keywords

Comments

Equivalently, a(n) is the number of ways to express the cyclic group Z_n as a direct sum of its Hall subgroups. A Hall subgroup of a finite group G is a subgroup whose order is coprime to its index.
a(n) is the number of ways to partition the set of distinct prime factors of n.
Also the number of singleton or pairwise coprime factorizations of n. - Gus Wiseman, Sep 24 2019

Examples

			a(60) = 5 because we have: 60 = 4*3*5 = 4*15 = 3*20 = 5*12.
For n = 36, its unitary divisors are 1, 4, 9, 36. From these we obtain 36 either as 1*36 or 4*9, thus a(36) = 2. - _Antti Karttunen_, Oct 21 2017
		

Crossrefs

Differs from A050320 for the first time at n=36.
Differs from A354870 for the first time at n=210, where a(210) = 15, while A354870(210) = 12.
Related classes of factorizations:
- No conditions: A001055
- Strict: A045778
- Constant: A089723
- Distinct multiplicities: A255231
- Singleton or coprime: A259936
- Relatively prime: A281116
- Aperiodic: A303386
- Stable (indivisible): A305149
- Connected: A305193
- Strict relatively prime: A318721
- Uniform: A319269
- Intersecting: A319786
- Constant or distinct factors coprime: A327399
- Constant or relatively prime: A327400
- Coprime: A327517
- Not relatively prime: A327658
- Distinct factors coprime: A327695

Programs

  • Maple
    map(combinat:-bell @ nops @ numtheory:-factorset, [$1..100]); # Robert Israel, Jul 09 2015
  • Mathematica
    Table[BellB[PrimeNu[n]], {n, 1, 75}]
    (* second program *)
    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],Length[#]==1||CoprimeQ@@#&]],{n,100}] (* Gus Wiseman, Sep 24 2019 *)
  • PARI
    a(n) = my(t=omega(n), x='x, m=contfracpnqn(matrix(2, t\2, y, z, if( y==1, -z*x^2, 1 - (z+1)*x)))); polcoeff(1/(1 - x + m[2, 1]/m[1, 1]) + O(x^(t+1)), t) \\ Charles R Greathouse IV, Jun 30 2017

Formula

a(n) = A000110(A001221(n)).
a(n > 1) = A327517(n) + 1. - Gus Wiseman, Sep 24 2019

Extensions

Incorrect comment removed by Antti Karttunen, Jun 11 2022

A327499 Quotient of n over the maximum divisor of n whose prime multiplicities are distinct.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 3, 1, 1, 1, 1, 1, 3, 2, 1, 1, 1, 2, 1, 1, 1, 6, 1, 1, 3, 2, 5, 2, 1, 2, 3, 1, 1, 6, 1, 1, 1, 2, 1, 1, 1, 1, 3, 1, 1, 1, 5, 1, 3, 2, 1, 3, 1, 2, 1, 1, 5, 6, 1, 1, 3, 10, 1, 1, 1, 2, 1, 1, 7, 6, 1, 1, 1, 2, 1, 3, 5, 2, 3, 1, 1, 2, 7, 1, 3, 2, 5, 1, 1, 1, 1, 2, 1, 6, 1, 1, 15
Offset: 1

Views

Author

Gus Wiseman, Sep 16 2019

Keywords

Comments

A number's prime multiplicities are also called its (unsorted) prime signature.

Examples

			The maximum such divisor of 60 is 20, so a(60) = 3.
		

Crossrefs

See link for additional cross-references.

Programs

  • Mathematica
    Table[n/Max[Select[Divisors[n],UnsameQ@@Last/@FactorInteger[#]&]],{n,100}]
  • PARI
    A351564(n) = issquarefree(factorback(apply(e->prime(e),(factor(n)[,2]))));
    A327499(n) = fordiv(n,d,if(A351564(n/d), return(d))); \\ Antti Karttunen, Apr 02 2022

Formula

a(n) = n/A327498(n).

Extensions

Data section extended up to 105 terms by Antti Karttunen, Apr 02 2022

A327500 Number of steps to reach a fixed point starting with n and repeatedly taking the quotient by the maximum divisor whose prime multiplicities are distinct (A327498, A327499).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Sep 16 2019

Keywords

Comments

A number's prime multiplicities are also called its (unsorted) prime signature. Numbers whose prime multiplicities are distinct are A130091.

Examples

			We have 9282 -> 546 -> 42 -> 6 -> 2 -> 1, so a(9282) = 5.
		

Crossrefs

See link for additional cross-references.
Position of first appearance of n is A002110(n).
Cf. also A327503.

Programs

  • Mathematica
    Table[Length[FixedPointList[#/Max[Select[Divisors[#],UnsameQ@@Last/@FactorInteger[#]&]]&,n]]-2,{n,100}]
  • PARI
    A351564(n) = issquarefree(factorback(apply(e->prime(e),(factor(n)[,2]))));
    A327499(n) = fordiv(n,d,if(A351564(n/d), return(d)));
    A327500(n) = { my(u=A327499(n)); if(u==n, 0, 1+A327500(u)); }; \\ Antti Karttunen, Apr 02 2022

Extensions

Data section extended up to 105 terms by Antti Karttunen, Apr 02 2022

A327523 Number of factorizations of the n-th number with distinct prime multiplicities A130091(n) into numbers > 1 with distinct prime multiplicities.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 3, 2, 1, 3, 1, 5, 1, 3, 1, 3, 1, 5, 2, 3, 3, 1, 1, 7, 1, 5, 1, 1, 3, 3, 1, 9, 2, 3, 3, 1, 5, 5, 1, 1, 3, 11, 1, 3, 1, 11, 1, 3, 3, 1, 9, 5, 1, 5, 1, 3, 14, 1, 3, 3, 1, 1, 5, 1, 11, 1, 9, 1, 3, 3, 2, 3, 3, 1, 15, 1, 5, 5, 1, 1, 20, 3, 3, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Sep 16 2019

Keywords

Comments

A number's prime multiplicities are also called its (unsorted) prime signature.

Examples

			The a(57) = 14 factorizations of 96 together with the corresponding multiset partitions of {1,1,1,1,1,2}:
  (2*2*2*2*2*3)  {{1}{1}{1}{1}{1}{2}}
  (2*2*2*3*4)    {{1}{1}{1}{2}{11}}
  (2*2*2*12)     {{1}{1}{1}{112}}
  (2*2*3*8)      {{1}{1}{2}{111}}
  (2*2*24)       {{1}{1}{1112}}
  (2*3*4*4)      {{1}{2}{11}{11}}
  (2*3*16)       {{1}{2}{1111}}
  (2*4*12)       {{1}{11}{112}}
  (2*48)         {{1}{11112}}
  (3*4*8)        {{2}{11}{111}}
  (3*32)         {{2}{11111}}
  (4*24)         {{11}{1112}}
  (8*12)         {{111}{112}}
  (96)           {{111112}}
		

Crossrefs

See link for additional cross-references.

Programs

  • Mathematica
    nn=100;
    facsusing[s_,n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facsusing[Select[s,Divisible[n/d,#]&],n/d],Min@@#>=d&]],{d,Select[s,Divisible[n,#]&]}]];
    y=Select[Range[nn],UnsameQ@@Last/@FactorInteger[#]&];
    Table[Length[facsusing[Rest[y],n]],{n,y}]

A319272 Numbers whose prime multiplicities are distinct and whose prime indices are term of the sequence.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 12, 16, 17, 18, 19, 20, 23, 24, 25, 27, 28, 31, 32, 37, 40, 44, 45, 48, 49, 50, 53, 54, 56, 59, 61, 63, 64, 67, 68, 71, 72, 75, 76, 80, 81, 83, 88, 89, 92, 96, 97, 98, 99, 103, 107, 108, 112, 121, 124, 125, 127, 128, 131, 135, 136, 144, 147, 148
Offset: 1

Views

Author

Gus Wiseman, Sep 16 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n.
Also Matula-Goebel numbers of rooted trees in which the multiplicities in the multiset of branches directly under any given node are distinct.

Examples

			36 is not in the sequence because 36 = 2^2 * 3^2 does not have distinct prime multiplicities.
The sequence of terms of the sequence followed by their Matula-Goebel trees begins:
   1: o
   2: (o)
   3: ((o))
   4: (oo)
   5: (((o)))
   7: ((oo))
   8: (ooo)
   9: ((o)(o))
  11: ((((o))))
  12: (oo(o))
  16: (oooo)
  17: (((oo)))
  18: (o(o)(o))
  19: ((ooo))
  20: (oo((o)))
  23: (((o)(o)))
  24: (ooo(o))
  25: (((o))((o)))
  27: ((o)(o)(o))
  28: (oo(oo))
  31: (((((o)))))
		

Crossrefs

Programs

  • Mathematica
    mgsiQ[n_]:=Or[n==1,And[UnsameQ@@Last/@FactorInteger[n],And@@Cases[FactorInteger[n],{p_,_}:>mgsiQ[PrimePi[p]]]]];
    Select[Range[100],mgsiQ]
  • PARI
    is(n)={my(f=factor(n)); if(#Set(f[,2])<#f~, 0, for(i=1, #f~, if(!is(primepi(f[i,1])), return(0))); 1)}
    { select(is, [1..200]) } \\ Andrew Howroyd, Mar 01 2020

Extensions

Terms a(53) and beyond from Andrew Howroyd, Mar 01 2020
Showing 1-8 of 8 results.