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

A050326 Number of factorizations of n into distinct squarefree numbers > 1.

Original entry on oeis.org

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

Views

Author

Christian G. Bower, Oct 15 1999

Keywords

Comments

a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24=2^3*3 and 375=3*5^3 both have prime signature (3,1).
a(A212164(n)) = 0; a(A212166(n)) = 1; a(A006881(n)) = 2; a(A190107(n)) = 3; a(A085987(n)) = 4; a(A225228(n)) = 5; a(A179670(n)) = 7; a(A162143(n)) = 8; a(A190108(n)) = 11; a(A212167(n)) > 0; a(A212168(n)) > 1. - Reinhard Zumkeller, May 03 2013
The comment that a(A212164(n)) = 0 is incorrect. For example, 3600 belongs to A212164 but a(3600) = 1. The positions of zeros in this sequence are A293243. - Gus Wiseman, Oct 10 2017

Examples

			The a(30) = 5 factorizations are: 2*3*5, 2*15, 3*10, 5*6, 30. The a(180) = 5 factorizations are: 2*3*5*6, 2*3*30, 2*6*15, 3*6*10, 6*30. - _Gus Wiseman_, Oct 10 2017
		

Crossrefs

Cf. A001055, A005117, A045778, A046523, A050320, A050327, a(p^k)=0 (p>1), a(A002110)=A000110, a(n!)=A103775(n), A206778, A293243.

Programs

  • Haskell
    import Data.List (subsequences, genericIndex)
    a050326 n = genericIndex a050326_list (n-1)
    a050326_list = 1 : f 2 where
       f x = (if x /= s then a050326 s
                        else length $ filter (== x) $ map product $
                             subsequences $ tail $ a206778_row x) : f (x + 1)
             where s = a046523 x
    -- Reinhard Zumkeller, May 03 2013
  • Maple
    N:= 1000: # to get a(1)..a(N)
    A:= Vector(N):
    A[1]:= 1:
    for n from 2 to N do
      if numtheory:-issqrfree(n) then
         S:= [$1..N/n]; T:= n*S; A[T]:= A[T]+A[S]
        fi;
    od:
    convert(A,list); # Robert Israel, Oct 10 2017
  • Mathematica
    sqfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
    Table[Length[sqfacs[n]],{n,100}] (* Gus Wiseman, Oct 10 2017 *)

Formula

Dirichlet g.f.: prod{n is squarefree and > 1}(1+1/n^s).
a(n) = A050327(A101296(n)). - R. J. Mathar, May 26 2017

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).

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

A212166 Numbers k such that the maximum exponent in its prime factorization equals the number of positive exponents (A051903(k) = A001221(k)).

Original entry on oeis.org

1, 2, 3, 5, 7, 11, 12, 13, 17, 18, 19, 20, 23, 28, 29, 31, 36, 37, 41, 43, 44, 45, 47, 50, 52, 53, 59, 61, 63, 67, 68, 71, 73, 75, 76, 79, 83, 89, 92, 97, 98, 99, 100, 101, 103, 107, 109, 113, 116, 117, 120, 124, 127, 131, 137, 139, 147, 148, 149, 151, 153
Offset: 1

Views

Author

Matthew Vandermast, May 22 2012

Keywords

Examples

			36 = 2^2*3^2 has 2 positive exponents in its prime factorization. The maximal exponent in its prime factorization is also 2. Therefore, 36 belongs to this sequence.
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 844.

Crossrefs

Includes subsequences A000040, A006939, A138534, A181555, A181825.
Cf. A001221, A050326, A051903, A188654 (complement), A225230.

Programs

  • Haskell
    import Data.List (elemIndices)
    a212166 n = a212166_list !! (n-1)
    a212166_list = map (+ 1) $ elemIndices 0 a225230_list
    -- Reinhard Zumkeller, May 03 2013
    
  • Mathematica
    okQ[n_] := Module[{f = Transpose[FactorInteger[n]][[2]]}, Max[f] == Length[f]]; Select[Range[424], okQ] (* T. D. Noe, May 24 2012 *)
  • PARI
    is(k) = {my(e = factor(k)[, 2]); !(#e) || vecmax(e) == #e;} \\ Amiram Eldar, Sep 08 2024

Formula

A225230(a(n)) = 0; A050326(a(n)) = 1. - Reinhard Zumkeller, May 03 2013

A212164 Numbers k such that the maximum exponent in its prime factorization is greater than the number of positive exponents (A051903(k) > A001221(k)).

Original entry on oeis.org

4, 8, 9, 16, 24, 25, 27, 32, 40, 48, 49, 54, 56, 64, 72, 80, 81, 88, 96, 104, 108, 112, 121, 125, 128, 135, 136, 144, 152, 160, 162, 169, 176, 184, 189, 192, 200, 208, 216, 224, 232, 240, 243, 248, 250, 256, 272, 288, 289, 296, 297, 304, 320, 324, 328, 336
Offset: 1

Views

Author

Matthew Vandermast, May 22 2012

Keywords

Examples

			40 = 2^3*5^1 has 2 distinct prime factors, hence, 2 positive exponents in its prime factorization (namely, 3 and 1, although the 1 is often left implicit).   2 is less than the maximal exponent in 40's prime factorization, which is 3. Therefore, 40 belongs to the sequence.
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 844.

Crossrefs

Complement of A212167.
See also A212165, A212166, A212168.
Subsequence of A188654.

Programs

  • Haskell
    import Data.List (elemIndices)
    a212164 n = a212164_list !! (n-1)
    a212164_list = map (+ 1) $ findIndices (< 0) a225230_list
    -- Reinhard Zumkeller, May 03 2013
    
  • Mathematica
    okQ[n_] := Module[{f = Transpose[FactorInteger[n]][[2]]}, Max[f] > Length[f]]; Select[Range[1000], okQ] (* T. D. Noe, May 24 2012 *)
  • PARI
    is(k) = {my(e = factor(k)[, 2]); #e && vecmax(e) > #e;} \\ Amiram Eldar, Sep 08 2024

Formula

A225230(a(n)) < 0; A050326(a(n)) = 0. - Reinhard Zumkeller, May 03 2013

A212167 Numbers k such that the maximum exponent in its prime factorization is not greater than the number of positive exponents (A051903(k) <= A001221(k)).

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 26, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 50, 51, 52, 53, 55, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 82, 83
Offset: 1

Views

Author

Matthew Vandermast, May 22 2012

Keywords

Comments

Union of A212166 and A212168. Includes numerous subsequences that are subsequences of neither A212166 nor A212168.

Examples

			40 = 2^3*5^1 has 2 distinct prime factors, hence, 2 positive exponents in its prime factorization (although the 1 is often left implicit).  2 is less than the maximal exponent in 40's prime factorization, which is 3. Therefore, 40 does not belong to the sequence. But 10 = 2^1*5^1 and 20 = 2^2*5^1 belong, since the maximal exponents in their prime factorizations are 1 and 2 respectively.
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 844.

Crossrefs

Complement of A212164. See also A212165.
Subsequences (none of which are subsequences of A212166 or A212168) include A002110, A051451, A129912, A179983, A181826, A181827, A182862, A182863. Includes all members of A003418.

Programs

  • Haskell
    import Data.List (findIndices)
    a212167 n = a212167_list !! (n-1)
    a212167_list = map (+ 1) $ findIndices (>= 0) a225230_list
    -- Reinhard Zumkeller, May 03 2013
    
  • Maple
    isA212167 := proc(n)
        simplify(A051903(n) <= A001221(n)) ;
    end proc:
    for n from 1 to 1000 do
        if isA212167(n) then
            printf("%d,",n) ;
        end if;
    end do: # R. J. Mathar, Jan 06 2021
  • Mathematica
    okQ[n_] := Module[{f = Transpose[FactorInteger[n]][[2]]}, Max[f] <= Length[f]]; Select[Range[1000], okQ] (* T. D. Noe, May 24 2012 *)
  • PARI
    is(k) = {my(e = factor(k)[, 2]); !(#e) || vecmax(e) <= #e; } \\ Amiram Eldar, Sep 09 2024

Formula

A225230(a(n)) >= 0; A050326(a(n)) > 0. - Reinhard Zumkeller, May 03 2013

A212165 Numbers k such that the maximum exponent in its prime factorization is not less than the number of positive exponents (A051903(k) >= A001221(k)).

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 18, 19, 20, 23, 24, 25, 27, 28, 29, 31, 32, 36, 37, 40, 41, 43, 44, 45, 47, 48, 49, 50, 52, 53, 54, 56, 59, 61, 63, 64, 67, 68, 71, 72, 73, 75, 76, 79, 80, 81, 83, 88, 89, 92, 96, 97, 98, 99, 100, 101, 103, 104
Offset: 1

Views

Author

Matthew Vandermast, May 22 2012

Keywords

Comments

Union of A212164 and A212166. Includes numerous subsequences that are subsequences of neither A212164 nor A212166.
Includes all factorials except A000142(3) = 6.
Observation: all terms in DATA section are also the first 65 numbers n whose difference between the arithmetic derivative of n and the sum of the divisors of n is nonnegative. - Omar E. Pol, Dec 19 2012

Examples

			10 = 2^1*5^1 has 2 distinct prime factors, hence, 2 positive exponents in its prime factorization (although 1s are often left implicit).  2 is larger than the maximal exponent in 10's prime factorization, which is 1. Therefore, 10 does not belong to the sequence. But 20 = 2^2*5^1 and 40 = 2^3*5^1 belong, since the largest exponents in their prime factorizations are 2 and 3 respectively.
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 844.

Crossrefs

Complement of A212168.
See also A212167.
Subsequences (none of which are subsequences of A212164 or A212166) include A000079, A001021, A066120, A087980, A130091, A141586, A166475, A181818, A181823, A181824, A182763, A212169. Also includes all terms in A181813 and A181814.

Programs

  • Haskell
    import Data.List (findIndices)
    a212165 n = a212165_list !! (n-1)
    a212165_list = map (+ 1) $ findIndices (<= 0) a225230_list
    -- Reinhard Zumkeller, May 03 2013
    
  • Mathematica
    okQ[n_] := Module[{f = Transpose[FactorInteger[n]][[2]]}, Max[f] >= Length[f]]; Select[Range[1000], okQ] (* T. D. Noe, May 24 2012 *)
  • PARI
    is(k) = {my(e = factor(k)[, 2]); !(#e) || vecmax(e) >= #e;} \\ Amiram Eldar, Sep 08 2024

Formula

A225230(a(n)) <= 0. - Reinhard Zumkeller, May 03 2013

A211991 Difference between the arithmetic derivative of n and the sum of proper divisors of n.

Original entry on oeis.org

0, 0, 0, 1, 0, -1, 0, 5, 2, -1, 0, 0, 0, -1, -1, 17, 0, 0, 0, 2, -1, -1, 0, 8, 4, -1, 14, 4, 0, -11, 0, 49, -1, -1, -1, 5, 0, -1, -1, 18, 0, -13, 0, 8, 6, -1, 0, 36, 6, 2, -1, 10, 0, 15, -1, 28, -1, -1, 0, -16, 0, -1, 10, 129, -1, -17, 0, 14, -1, -15, 0, 33
Offset: 1

Views

Author

Omar E. Pol, Dec 18 2012

Keywords

Comments

Observations: at least the first 50 indices of nonnegative terms are also the first 50 terms of A212165. Also at least the first 28 indices of negative terms are also the first 28 terms of A212168, since A212168 is the complement of A212165.

Crossrefs

Programs

  • Mathematica
    dn[0] = 0; dn[1] = 0; dn[n_?Negative] := -dn[-n]; dn[n_] := Module[{f = Transpose[FactorInteger[n]]}, If[PrimeQ[n], 1, Plus @@ (n*f[[2]]/f[[1]])]]; Table[dn[n] - (DivisorSigma[1, n] - n), {n, 100}] (* T. D. Noe, Dec 27 2012 *)
  • PARI
    A003415(n) = {my(fac); if(n<1, 0, fac=factor(n); sum(i=1, matsize(fac)[1], n*fac[i, 2]/fac[i, 1]))}; \\ From A003415
    A211991(n) = (A003415(n) - (sigma(n)-n)); \\ Antti Karttunen, Mar 08 2018

Formula

a(n) = A003415(n) - A001065(n).
Sum_{k=1..n} a(k) ~ c * n^2, where c = (A136141 - A013661 + 1) / 2 = 0.0641113... . - Amiram Eldar, Mar 17 2024

A188654 Numbers k such that the maximum exponent in its prime factorization does not equal the number of positive exponents (A051903(k) <> A001221(k)).

Original entry on oeis.org

4, 6, 8, 9, 10, 14, 15, 16, 21, 22, 24, 25, 26, 27, 30, 32, 33, 34, 35, 38, 39, 40, 42, 46, 48, 49, 51, 54, 55, 56, 57, 58, 60, 62, 64, 65, 66, 69, 70, 72, 74, 77, 78, 80, 81, 82, 84, 85, 86, 87, 88, 90, 91, 93, 94, 95, 96, 102, 104, 105, 106, 108, 110, 111
Offset: 1

Views

Author

Reinhard Zumkeller, May 03 2013

Keywords

Crossrefs

Cf. A001221, A212166 (complement), A225230.
Union of A212164 and A212168.

Programs

  • Haskell
    import Data.List (findIndices)
    a188654 n = a188654_list !! (n-1)
    a188654_list = map (+ 1) $ findIndices (/= 0) a225230_list
    
  • Mathematica
    q[n_] := Module[{e = FactorInteger[n][[;;, 2]]}, Max[e] != Length[e]]; q[1] = False; Select[Range[120], q] (* Amiram Eldar, Sep 08 2024 *)
  • PARI
    is(k) = {my(e = factor(k)[, 2]); #e && vecmax(e) != #e;} \\ Amiram Eldar, Sep 08 2024

Formula

A051903(n) <> A001221(n);
A225230(a(n)) <> 0.

A225230 In the canonical prime factorization of n: (number of distinct primes) minus (largest prime exponent).

Original entry on oeis.org

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

Views

Author

Reinhard Zumkeller, May 03 2013

Keywords

Crossrefs

Programs

  • Haskell
    a225230 n = a001221 n - a051903 n
    
  • Mathematica
    a[n_] := Module[{e = FactorInteger[n][[;;, 2]]}, Length[e] - Max[e]]; Array[a, 100] (* Amiram Eldar, Sep 09 2024 *)
  • PARI
    a(n) = if (n>1, my(f=factor(n)); #f~ - vecmax(f[,2]), 0); \\ Michel Marcus, Jan 26 2022

Formula

a(n) = A001221(n) - A051903(n).
a(A212164(n)) < 0; a(A212165(n)) <= 0; a(A212166(n)) = 0; a(A188654(n)) <> 0; a(A212167(n)) >= 0; a(A212168(n)) > 0.
Showing 1-10 of 12 results. Next