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

A088529 Numerator of Bigomega(n)/Omega(n).

Original entry on oeis.org

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

Views

Author

Cino Hilliard, Nov 16 2003

Keywords

Examples

			bigomega(24) / omega(24) = 4/2 = 2, so a(24) = 2.
		

References

  • H. Z. Cao, On the average of exponents, Northeast. Math. J., Vol. 10 (1994), pp. 291-296.

Crossrefs

Cf. A001221, A001222, A070012, A070013, A070014, A088530 (gives the denominator).

Programs

  • Mathematica
    Table[Numerator[PrimeOmega[n]/PrimeNu[n]], {n, 2, 100}] (* Michael De Vlieger, Jul 12 2017 *)
  • PARI
    for(x=2,100,y=bigomega(x)/omega(x);print1(numerator(y)","))
    
  • Python
    from sympy import primefactors, Integer
    def bigomega(n): return 0 if n==1 else bigomega(Integer(n)/primefactors(n)[0]) + 1
    def omega(n): return Integer(len(primefactors(n)))
    def a(n): return (bigomega(n)/omega(n)).numerator
    print([a(n) for n in range(2, 51)]) # Indranil Ghosh, Jul 13 2017

Formula

Let B = number of prime divisors of n with multiplicity, O = number of distinct prime divisors of n. Then a(n) = numerator of B/O.
a(n) = A136565(n) = A181591(n) for n: 2 <= n < 24. - Reinhard Zumkeller, Nov 01 2010
Sum_{k=2..n} a(k)/A088530(k) ~ n + O(n/log(log(n))) (Duncan, 1970). - Amiram Eldar, Oct 14 2022
Sum_{k=2..n} a(k)/A088530(k) = n + c_1 * n/log(log(n)) + c_2 * n/log(log(n))^2 + O(n/log(log(n))^3), where c_1 = A136141 and c_2 = A272531 (Cao, 1994; Finch, 2020). - Amiram Eldar, Dec 15 2022

A088530 Denominator of bigomega(n)/omega(n).

Original entry on oeis.org

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

Views

Author

Cino Hilliard, Nov 16 2003

Keywords

Comments

a(n) is the denominator of A022559(n)/A000720(n). - Robert Israel, Jan 08 2024

Examples

			bigomega(24) / omega(24) = 4/2 = 2/1, so a(24) = 1.
		

Crossrefs

Cf. A001221, A001222, A000720, A022559, A070012, A070013, A070014, A088529 (gives the numerator).

Programs

  • Maple
    N:= 100:
    W:= ListTools:-PartialSums(map(numtheory:-bigomega,[$1..N])):
    seq(denom(W[i]/numtheory:-pi(i)),i=2..N); # Robert Israel, Jan 08 2024
  • Mathematica
    Table[Denominator[PrimeOmega[n]/PrimeNu[n]],{n,2,100}] (* Harvey P. Dale, Mar 22 2012 *)
  • PARI
    for(x=2,100,y=bigomega(x)/omega(x);print1(denominator(y)","))
    
  • Python
    from sympy import primefactors, Integer
    def bigomega(n): return 0 if n==1 else bigomega(Integer(n)/primefactors(n)[0]) + 1
    def omega(n): return Integer(len(primefactors(n)))
    def a(n): return (bigomega(n)/omega(n)).denominator
    print([a(n) for n in range(2, 51)]) # Indranil Ghosh, Jul 13 2017

Formula

Let B = number of prime divisors of n with multiplicity, O = number of distinct prime divisors of n. Then a(n) = denominator of B/O.

A070012 Floor of number of prime factors of n divided by the number of n's distinct prime factors.

Original entry on oeis.org

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

Views

Author

Rick L. Shepherd, Apr 11 2002

Keywords

Comments

a(n) is the integer part of the average of the exponents in the prime factorization of n.

Examples

			a(12)=1 because 12=2^2 * 3^1 and floor(bigomega(12)/omega(12)) = floor((2+1)/2) = 1.
a(36)=2 because 36=2^2 * 3^2 and floor(bigomega(36)/omega(36)) = floor((2+2)/2) = 2.
a(60)=1 because 60=2^2 * 3^1 * 5^1 and floor(bigomega(60)/omega(60)) = floor((2+1+1)/3) = 1.
36 is in A067340. 12 and 60 are in A070011.
		

Crossrefs

Cf. A001221 (omega(n)), A001222 (bigomega(n)), A067340 (ratio is an integer before floor applied), A070011 (ratio is not an integer), A070013 (ratio rounded), A070014 (ceiling of ratio), A046660 (bigomega(n)-omega(n)).

Programs

  • Mathematica
    A070012[n_]:=Floor[PrimeOmega[n]/PrimeNu[n]];Array[A070012,100]
  • PARI
    v=[]; for(n=2,150,v=concat(v,(bigomega(n)\omega(n)))); v

Formula

a(n) = floor(bigomega(n)/omega(n)) for n>=2.

A070014 Ceiling of number of prime factors of n divided by the number of n's distinct prime factors.

Original entry on oeis.org

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

Views

Author

Rick L. Shepherd, Apr 11 2002

Keywords

Comments

a(n) is the ceiling of the average of the exponents in the prime factorization of n.

Examples

			a(12) = 2 because 12 = 2^2 * 3^1 and ceiling(bigomega(12)/omega(12)) = ceiling((2+1)/2) = 2. a(36) = 2 because 36 = 2^2 * 3^2 and ceiling(bigomega(36)/omega(36)) = ceiling((2+2)/2) = 2. a(60) = 2 because 60 = 2^2 * 3^1 * 5^1 and ceiling(bigomega(60)/omega(60)) = ceiling((2+1+1)/3) = 2. 36 is in A067340. 12 and 60 are in A070011.
		

Crossrefs

Cf. A001221 (omega(n)), A001222 (bigomega(n)), A067340 (ratio is an integer before ceil is applied), A070011 (ratio is not an integer), A070012 (floor of ratio), A070013 (ratio rounded), A046660 (bigomega(n)-omega(n)), A088529, A088530.

Programs

  • Mathematica
    Table[Ceiling[PrimeOmega[n]/PrimeNu[n]], {n, 2, 106}] (* Michael De Vlieger, Jul 12 2017 *)
  • PARI
    v=[]; for(n=2,150,v=concat(v,ceil(bigomega(n)/omega(n)))); v
    
  • Python
    from sympy import primefactors, ceiling
    def bigomega(n): return 0 if n==1 else bigomega(n//primefactors(n)[0]) + 1
    def omega(n): return len(primefactors(n))
    def a(n): return ceiling(bigomega(n)/omega(n))
    print([a(n) for n in range(2, 51)]) # Indranil Ghosh, Jul 13 2017
  • Scheme
    (define (A070014 n) (let ((a (A001222 n)) (b (A001221 n))) (if (zero? (modulo a b)) (/ a b) (+ 1 (/ (- a (modulo a b)) b))))) ;; Antti Karttunen, Jul 12 2017
    

Formula

a(n) = ceiling(bigomega(n)/omega(n)) for n>=2.

A289621 Compound filter (omega & bigomega): a(1) = 0, for n > 1, a(n) = P(A001221(n), A001222(n)), where P(n,k) is sequence A000027 used as a pairing function.

Original entry on oeis.org

0, 1, 1, 2, 1, 5, 1, 4, 2, 5, 1, 8, 1, 5, 5, 7, 1, 8, 1, 8, 5, 5, 1, 12, 2, 5, 4, 8, 1, 13, 1, 11, 5, 5, 5, 12, 1, 5, 5, 12, 1, 13, 1, 8, 8, 5, 1, 17, 2, 8, 5, 8, 1, 12, 5, 12, 5, 5, 1, 18, 1, 5, 8, 16, 5, 13, 1, 8, 5, 13, 1, 17, 1, 5, 8, 8, 5, 13, 1, 17, 7, 5, 1, 18, 5, 5, 5, 12, 1, 18, 5, 8, 5, 5, 5, 23, 1, 8, 8
Offset: 1

Views

Author

Antti Karttunen, Jul 16 2017

Keywords

Crossrefs

Cf. A001221, A001222, A008966, A046660, A070012, A070013, A070014, A088529, A088530, A181591 (sequences with matching equivalence classes).

Programs

  • PARI
    A289621(n) = if(1==n,0,(1/2)*(2 + ((omega(n)+bigomega(n))^2) - omega(n) - 3*bigomega(n)));
    
  • Scheme
    (define (A289621 n) (if (= 1 n) 0 (* (/ 1 2) (+ (expt (+ (A001221 n) (A001222 n)) 2) (- (A001221 n)) (- (* 3 (A001222 n))) 2))))

Formula

a(1) = 0, for n > 1, a(n) = (1/2)*(2 + ((A001221(n)+A001222(n))^2) - A001221(n) - 3*A001222(n)).
Showing 1-5 of 5 results.