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

A212180 Number of distinct second signatures (cf. A212172) represented among divisors of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 2, 1, 1, 1, 4, 1, 2, 1, 2, 1, 1, 1, 3, 2, 1, 3, 2, 1, 1, 1, 5, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 4, 2, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 2, 1, 1, 2, 6, 1, 1, 1, 2, 1, 1, 1, 5, 1, 1, 2, 2, 1, 1, 1, 4, 4, 1, 1, 2, 1, 1, 1, 3
Offset: 1

Views

Author

Matthew Vandermast, Jun 04 2012

Keywords

Comments

Completely determined by the exponents >=2 in the prime factorization of n (cf. A212172, A212173).
The fraction of the divisors of n which have a given second signature {S} is also a function of n's second signature. For example, if n has second signature {3,2}, it follows that 1/3 of n's divisors are squarefree. Squarefree numbers are represented with 0's in A212172, in accord with the usual OEIS custom of using 0 for nonexistent elements; in comments, their second signature is represented as { }.

Examples

			The divisors of 72 represent a total of 5 distinct second signatures (cf. A212172), as can be seen from the exponents >= 2, if any, in the canonical prime factorization of each divisor:
{ }: 1, 2 (prime), 3 (prime), 6 (2*3)
{2}: 4 (2^2), 9 (3^2), 12 (2^2*3), 18 (2*3^2)
{3}: 8 (2^3), 24 (2^3*3)
{2,2}: 36 (2^2*3^2)
{3,2}: 72 (2^3*3^2)
Hence, a(72) = 5.
		

Crossrefs

Programs

  • Mathematica
    Array[Length@ Union@ Map[Sort@ Select[FactorInteger[#][[All, -1]], # >= 2 &] &, Divisors@ #] &, 88] (* Michael De Vlieger, Jul 19 2017 *)
  • PARI
    A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); }; \\ This function from Charles R Greathouse IV, Aug 17 2011
    A057521(n) = { my(f=factor(n)); prod(i=1, #f~, if(f[i, 2]>1, f[i, 1]^f[i, 2], 1)); } \\ This function from Charles R Greathouse IV, Aug 13 2013
    A212173(n) = A046523(A057521(n));
    A212180(n) = { my(vals = Set()); fordiv(n, d, vals = Set(concat(vals, A212173(d)))); length(vals); }; \\ Antti Karttunen, Jul 19 2017
    
  • Python
    from sympy import factorint, divisors, prod
    def P(n): return sorted(factorint(n).values())
    def a046523(n):
        x=1
        while True:
            if P(n)==P(x): return x
            else: x+=1
    def a057521(n): return 1 if n==1 else prod(p**e for p, e in factorint(n).items() if e != 1)
    def a212173(n): return a046523(a057521(n))
    def a(n):
        l=[]
        for d in divisors(n):
            x=a212173(d)
            if not x in l:l+=[x, ]
        return len(l)
    print([a(n) for n in range(1, 51)]) # Indranil Ghosh, Jul 19 2017

A212644 If an integer's second signature (cf. A212172) is the n-th to appear among positive integers, a(n) = number of distinct second signatures represented among its divisors.

Original entry on oeis.org

1, 2, 3, 4, 5, 3, 6, 5, 7, 7, 6, 8, 9, 9, 9, 11, 12, 4, 10, 13, 10, 15, 7, 11, 15, 14, 18, 10, 12, 17, 18, 9, 21, 13, 15, 13, 19, 22, 14, 24, 16, 20, 14, 21, 26, 19, 10, 27, 19, 25, 16, 15, 23, 30, 24, 5, 21, 16, 30, 22, 30, 23, 16, 25, 34, 29, 9, 27, 22, 33
Offset: 1

Views

Author

Matthew Vandermast, Jun 07 2012

Keywords

Comments

Also, number of divisors of A181800(n) that are members of A181800.
Consider a member of A181800 with second signature {S} whose divisors represent a total of k distinct second signatures and a total of (j+k) distinct prime signatures (cf. A212642). Let m be any integer with second signature {S}. Then A212180(m) = k and A085082(m) is congruent to j modulo k. If {S} is the second signature of A181800(n), then A085082(m) is congruent to A212643(n) modulo a(n).

Examples

			The divisors of 72 represent 5 distinct second signatures (cf. A212172), as can be seen from the exponents >=2, if any, in the canonical prime factorization of each divisor:
{ }: 1, 2 (prime), 3 (prime), 6 (2*3)
{2}: 4 (2^2), 9 (3^2), 12 (2^2*3), 18 (2*3^2)
{3}: 8 (2^3), 24 (2^3*3)
{2,2}: 36 (2^2*3^2)
{3,2}: 72 (2^3*3^2)
Since 72 = A181800(8), a(8) = 5.
		

Crossrefs

Formula

a(n) = A212180(A181800(n)).

Extensions

Data corrected by Amiram Eldar, Jul 14 2019

A212643 Let b(n) and c(n) be the total numbers of distinct prime signatures and second signatures, respectively, represented among divisors of A181800(n) (first integers of each second signature; cf. A212172). b(n) mod c(n) = a(n).

Original entry on oeis.org

0, 1, 1, 1, 1, 0, 1, 4, 1, 5, 4, 1, 6, 5, 1, 7, 6, 2, 1, 8, 5, 7, 2, 1, 9, 6, 8, 2, 1, 10, 7, 1, 9, 2, 6, 1, 11, 8, 0, 10, 2, 7, 1, 12, 9, 18, 0, 11, 2, 8, 15, 1, 13, 10, 22, 0, 7, 14, 12, 2, 9, 20, 1, 14, 11, 26, 7, 8, 18, 13, 2, 10, 25, 1, 15, 15, 12, 30, 9
Offset: 1

Views

Author

Matthew Vandermast, Jun 05 2012

Keywords

Comments

Significance of the sequence: Consider a member of A181800 with second signature {S} whose divisors represent a total of k distinct second signatures and a total of (j+k) distinct prime signatures. For all integers n with second signature {S}, A212180(n) = k and A085082(n) is congruent to j modulo k; see examples.
Note: b(n) = A212642(n); c(n) = A212644(n).

Examples

			4 is the smallest integer with second signature {2}, and its divisors represent 3 distinct prime signatures and 2 distinct second signatures. 1 = 3 mod 2. Since 4 = A181800(2), a(2) = 1. For all integers m with second signature {2}, A085082(m) is congruent to 1 modulo 2.
10800 is the smallest integer with second signature {4,3,2}, and its divisors represent 28 distinct prime signatures and 14 distinct second signatures. 0 = 28 mod 14.  Since 10800 = A181800(39), a(39) = 0. For all integers m with second signature {4,3,2}, A085082(m) is congruent to 0 modulo 14.
		

Crossrefs

Formula

a(n) = A212642(n)-A212644(n), reduced modulo A212644(n).

A212173 First integer with same second signature as n (cf. A212172).

Original entry on oeis.org

1, 1, 1, 4, 1, 1, 1, 8, 4, 1, 1, 4, 1, 1, 1, 16, 1, 4, 1, 4, 1, 1, 1, 8, 4, 1, 8, 4, 1, 1, 1, 32, 1, 1, 1, 36, 1, 1, 1, 8, 1, 1, 1, 4, 4, 1, 1, 16, 4, 4, 1, 4, 1, 8, 1, 8, 1, 1, 1, 4, 1, 1, 4, 64, 1, 1, 1, 4, 1, 1, 1, 72, 1, 1, 4, 4, 1, 1, 1, 16, 16, 1, 1, 4
Offset: 1

Views

Author

Matthew Vandermast, Jun 03 2012

Keywords

Comments

Two integers have the same second signature iff the same exponents >= 2 occur in the canonical prime factorization of each integer, regardless of the order in which they occur in each factorization.

Examples

			12 = 2^2*3 has 1 exponent >= 2 in its prime factorization, namely, 2. Hence, its second signature is {2}.  The smallest number with second signature {2} is 4; hence, a(12) = 4.
		

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

Cf. A212172, A046523. All terms belong to A181800.

Programs

  • Magma
    A212173 := func; [A212173(n):n in[1..85]]; // Jason Kimberley, Jun 14 2012
    
  • Maple
    f:= proc(n) local E,i;
    E:= sort(select(`>`, map(t -> t[2], ifactors(n)[2]), 1),`>`);
    mul(ithprime(i)^E[i],i=1..nops(E))
    end proc:
    map(f, [$1..100]); # Robert Israel, Jul 19 2017
  • Mathematica
    Function[s, Sort[Apply[Join, Map[Function[k, Map[{#, First@ k} &, k]], Values@ s]]][[All, -1]]]@ KeySort@ PositionIndex@ Table[Sort@ DeleteCases[FactorInteger[n][[All, -1]], e_ /; e < 2] /. {} -> {1}, {n, 84}] (* Michael De Vlieger, Jul 19 2017 *)
  • PARI
    a(n) = {my(sn = vecsort(select(x->(x>=2), factor(n)[,2]))); for (i=1, n, if (vecsort(select(x->(x>=2),factor(i)[,2])) == sn, return(i)););} \\ Michel Marcus, Jul 19 2017
  • Python
    from functools import reduce
    from sympy import factorint
    from operator import mul
    def P(n): return sorted(factorint(n).values())
    def a046523(n):
        x=1
        while True:
            if P(n)==P(x): return x
            else: x+=1
    def a057521(n): return 1 if n==1 else reduce(mul, [1 if e==1 else p**e for p, e in factorint(n).items()])
    def a(n): return a046523(a057521(n))
    print([a(n) for n in range(1, 151)]) # Indranil Ghosh, Jul 19 2017
    

Formula

a(n) = A046523(A057521(n)) = A057521(A046523(n)).

A275812 Sum of exponents larger than one in the prime factorization of n: A001222(n) - A056169(n).

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Aug 11 2016

Keywords

Crossrefs

Differs from A212172 for the first time at n=36, where a(36)=4, while A212172(36)=2.

Programs

  • Mathematica
    Table[Total@ Map[Last, Select[FactorInteger@ n, Last@ # > 1 &] /. {} -> {{0, 0}}], {n, 120}] (* Michael De Vlieger, Aug 11 2016 *)
  • PARI
    a(n) = my(f = factor(n)); sum(k=1, #f~, if (f[k,2] > 1, f[k,2])); \\ Michel Marcus, Jul 19 2017
  • Perl
    sub a275812 { vecsum( grep {$> 1} map {$->[1]} factor_exp(shift) ); } # Dana Jacobsen, Aug 15 2016
    
  • Python
    from sympy import factorint, primefactors
    def a001222(n):
        return 0 if n==1 else a001222(n//primefactors(n)[0]) + 1
    def a056169(n):
        f=factorint(n)
        return 0 if n==1 else sum(1 for i in f if f[i]==1)
    def a(n):
        return a001222(n) - a056169(n)
    print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 19 2017
    

Formula

a(1) = 0, and for n > 1, if A067029(n)=1 [when n is one of the terms of A247180], a(n) = a(A028234(n)), otherwise a(n) = A067029(n)+a(A028234(n)).
a(n) = A001222(n) - A056169(n).
a(n) = A001222(A057521(n)). - Antti Karttunen, Jul 19 2017
From Amiram Eldar, Sep 28 2023: (Start)
Additive with a(p) = 0, and a(p^e) = e for e >= 2.
a(n) >= 0, with equality if and only if n is squarefree (A005117).
a(n) <= A001222(n), with equality if and only if n is powerful (A001694).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{p prime} (1/p^2 + 1/(p*(p-1))) = A085548 + A136141 = 1.22540408909086062637... . (End)
a(n) = A046660(n) + A056170(n). - Amiram Eldar, Jan 09 2024

A212171 Prime signature of n (nonincreasing version): row n of table lists positive exponents in canonical prime factorization of n, in nonincreasing order.

Original entry on oeis.org

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

Views

Author

Matthew Vandermast, Jun 03 2012

Keywords

Comments

Length of row n equals A001221(n).
The multiset of positive exponents in n's prime factorization completely determines a(n) for a host of OEIS sequences, including several "core" sequences. Of those not cross-referenced here or in A212172, many can be found by searching the database for A025487.
(Note: Differing opinions may exist about whether the prime signature of n should be defined as this multiset itself, or as a symbol or collection of symbols that identify or "signify" this multiset. The definition of this sequence is designed to be compatible with either view, as are the original comments. When n >= 2, the customary ways to signify the multiset of exponents in n's prime factorization are to list the constituent exponents in either nonincreasing or nondecreasing order; this table gives the nonincreasing version.)
Table lists exponents in the order in which they appear in the prime factorization of a member of A025487. This ordering is common in database comments (e.g., A008966).
Each possible multiset of an integer's positive prime factorization exponents corresponds to a unique partition that contains the same elements (cf. A000041). This includes the multiset of 1's positive exponents, { } (the empty multiset), which corresponds to the partition of 0.
Differs from A124010 from a(23) on, corresponding to the factorization of 18 = 2^1*3^2 which is here listed as row 18 = [2, 1], but as [1, 2] (in the order of the prime factors) in A124010 and also in A118914 which lists the prime signatures in nondecreasing order (so that row 12 = 2^2*3^1 is also [1, 2]). - M. F. Hasler, Apr 08 2022

Examples

			First rows of table read:
  1;
  1;
  2;
  1;
  1,1;
  1;
  3;
  2;
  1,1;
  1;
  2,1;
  ...
The multiset of positive exponents in the prime factorization of 6 = 2*3 is {1,1} (1s are often left implicit as exponents). The prime signature of 6 is therefore {1,1}.
12 = 2^2*3 has positive exponents 2 and 1 in its prime factorization, as does 18 = 2*3^2. Rows 12 and 18 of the table both read {2,1}.
		

Crossrefs

Cf. A025487, A001221 (row lengths), A001222 (row sums). A118914 gives the nondecreasing version. A124010 lists exponents in n's prime factorization in natural order, with A124010(1) = 0.
A212172 cross-references over 20 sequences that depend solely on n's prime exponents >= 2, including the "core" sequence A000688. Other sequences determined by the exponents in the prime factorization of n include:
Additive: A001221, A001222, A056169.
A highly incomplete selection of sequences, each definable by the set of prime signatures possessed by its members: A000040, A000290, A000578, A000583, A000961, A001248, A001358, A001597, A001694, A002808, A004709, A005117, A006881, A013929, A030059, A030229, A052486.

Programs

  • Magma
    &cat[Reverse(Sort([pe[2]:pe in Factorisation(n)])):n in[1..76]]; // Jason Kimberley, Jun 13 2012
    
  • PARI
    apply( {A212171_row(n)=vecsort(factor(n)[,2]~,,4)}, [1..40])\\ M. F. Hasler, Apr 19 2022

Formula

Row n of A118914, reversed.
Row n of A124010 for n > 1, with exponents sorted in nonincreasing order. Equivalently, row A046523(n) of A124010 for n > 1.

A212176 Row n of table lists exponents in canonical prime factorization of A181800(n) (n-th powerful number that is the first integer of its prime signature), in nonincreasing order.

Original entry on oeis.org

2, 3, 4, 5, 2, 2, 6, 3, 2, 7, 4, 2, 3, 3, 8, 5, 2, 4, 3, 9, 6, 2, 5, 3, 2, 2, 2, 10, 7, 2, 4, 4, 6, 3, 3, 2, 2, 11, 8, 2, 5, 4, 7, 3, 4, 2, 2, 12, 9, 2, 6, 4, 3, 3, 2, 8, 3, 5, 2, 2, 5, 5, 13, 10, 2, 7, 4, 4, 3, 2, 9, 3, 6, 2, 2, 6, 5, 14, 11, 2, 8, 4, 5, 3, 2, 3
Offset: 2

Views

Author

Matthew Vandermast, Jun 03 2012

Keywords

Comments

A212179(n) gives length of row n.
Table represents prime signature (cf. A212171) and second signature (cf. A212172) of A181800.

Examples

			Since 72 is a member of A181800, all positive exponents in its prime factorization (2^3*3^2) equal or exceed 2. Therefore, its second signature is the same as its prime signature, namely, {3,2} (nonincreasing version).  Since 72 = A181800 (8), row 8 represents the prime signature and second signature {3,2}.
		

Crossrefs

Formula

Row n is identical to row A181800(n) of tables A212171 and A212172.

A212181 Largest odd divisor of tau(n): a(n) = A000265(A000005(n)).

Original entry on oeis.org

1, 1, 1, 3, 1, 1, 1, 1, 3, 1, 1, 3, 1, 1, 1, 5, 1, 3, 1, 3, 1, 1, 1, 1, 3, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 9, 1, 1, 1, 1, 1, 1, 1, 3, 3, 1, 1, 5, 3, 3, 1, 3, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 3, 7, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 3, 3, 1, 1, 1, 5, 5, 1, 1, 3, 1, 1, 1, 1
Offset: 1

Views

Author

Matthew Vandermast, Jun 04 2012

Keywords

Comments

Completely determined by the exponents >=2 in the prime factorization of n (cf. A212172).
Not the same as the number of odd divisors of n (A001227(n)); see example.
Multiplicative because A000005 is multiplicative and A000265 is completely multiplicative. - Andrew Howroyd, Aug 01 2018
a(n) = 1 iff the number of divisors of n is a power of 2 (A036537). - Bernard Schott, Nov 04 2022

Examples

			48 has a total of 10 divisors (1, 2, 3, 4, 6, 8, 12, 16, 24 and 48). Since the largest odd divisor of 10 is 5, a(48) = 5.
		

Crossrefs

Cf. A000005, A000079, A000265, A036537, A108951, A212172, A295664, A331286 (applied to primorial inflation of n).

Programs

  • Mathematica
    Table[Block[{nd=DivisorSigma[0, n]}, nd/2^IntegerExponent[nd, 2]], {n, 100}] (* Indranil Ghosh, Jul 19 2017, after PARI code *)
  • PARI
    a(n) = my(nd = numdiv(n)); nd/2^valuation(nd, 2); \\ Michel Marcus, Jul 19 2017
    
  • Python
    from sympy import divisor_count, divisors
    def a(n): return [i for i in divisors(divisor_count(n)) if i%2][-1]
    print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jul 19 2017

Formula

a(n) = A000265(A000005(n)).
From Antti Karttunen, Jan 14 2020: (Start)
a(n) = A000005(n) / A000079(A295664(n)).
a(A108951(n)) = A331286(n).
(End)
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p odd prime} ((1 - 1/p)*(1 + Sum_{k>=1} a(k+1)/p^k)) = 2.076325817863586... . - Amiram Eldar, Oct 15 2022

A212642 a(n) = number of distinct prime signatures represented among divisors of A181800(n) (n-th powerful number that is the first integer of its prime signature).

Original entry on oeis.org

1, 3, 4, 5, 6, 6, 7, 9, 8, 12, 10, 9, 15, 14, 10, 18, 18, 10, 11, 21, 15, 22, 16, 12, 24, 20, 26, 22, 13, 27, 25, 19, 30, 28, 21, 14, 30, 30, 28, 34, 34, 27, 15, 33, 35, 37, 20, 38, 40, 33, 31, 16, 36, 40, 46, 15, 28, 30, 42, 46, 39, 43, 17, 39, 45, 55, 25, 35
Offset: 1

Views

Author

Matthew Vandermast, Jun 05 2012

Keywords

Comments

Also, number of divisors of A181800 that are members of A025487.
Consider a member of A181800 with second signature {S} whose divisors represent a total of k distinct second signatures and a total of (j+k) distinct prime signatures. Let n be any integer with second signature {S}. Then A212180(n) = k and A085082(n) is congruent to j modulo k. Cf. A212643, A212644.

Examples

			The divisors of 36 represent a total of 6 distinct prime signatures (cf. A085082), as can be seen from the positive exponents, if any, in the canonical prime factorization of each divisor:
{ }: 1 (multiset of positive exponents is the empty multiset)
{1}: 2 (2^1), 3 (3^1)
{1,1}: 6 (2^1*3^1)
{2}: 4 (2^2), 9 (3^2),
{2,1}: 12 (2^2*3^1), 18 (2^1*3^2)
{2,2}: 36 (2^2*3^2)
Since 36 = A181800(6), a(6) = 6.
		

Crossrefs

Formula

a(n) = A085082(A181800(n)).

A372603 The maximal exponent in the prime factorization of the powerful part of n.

Original entry on oeis.org

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

Views

Author

Amiram Eldar, May 07 2024

Keywords

Comments

First differs from A275812 at n = 36, and from A212172 at n = 37.

Crossrefs

Programs

  • Mathematica
    f[n_] := If[n == 1, 0, n]; a[n_] := f[Max[FactorInteger[n][[;; , 2]]]]; a[1] = 0; Array[a, 100]
  • PARI
    s(n) = if(n == 1, 0, n);
    a(n) = if(n>1, s(vecmax(factor(n)[,2])), 0);

Formula

a(n) = A051903(A057521(n)).
a(n) = A087156(A051903(n)).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 - 1/zeta(2) + Sum_{i>=2} (1 - 1/zeta(i)) = A033150 - A059956 = 1.09728403825134113562... .
Showing 1-10 of 25 results. Next