cp's OEIS Frontend

A228409 a(n) = 4*mu(n) + 5, where mu is the Moebius function (A008683).

%I A228409 #32 Oct 03 2017 03:20:02
%S A228409 9,1,1,5,1,9,1,5,5,9,1,5,1,9,9,5,1,5,1,5,9,9,1,5,5,9,5,5,1,1,1,5,9,9,
%T A228409 9,5,1,9,9,5,1,1,1,5,5,9,1,5,5,5,9,5,1,5,9,5,9,9,1,5,1,9,5,5,9,1,1,5,
%U A228409 9,1,1,5,1,9,5,5,9,1,1,5,5,9,1,5,9,9,9
%N A228409 a(n) = 4*mu(n) + 5, where mu is the Moebius function (A008683).
%C A228409 If n is prime (A000040), then a(n) = 1. The converse is not true: when n is the product of an odd number of distinct primes, mu(n) = -1 => a(n) = 1 (30 = 2*3*5, so a(30) = 1).
%C A228409 If n is semiprime (A001358), a(n) gives the number of divisors of n^2. In particular, if n = p^2 then n^2 = (p^2)^2 = p^4 has 5 divisors: p^4, p^3, p^2, p, 1. If n = pq (p,q distinct primes) then n^2 = (pq)^2 has 9 divisors: (pq)^2, qp^2, pq^2, p^2, q^2, pq, p, q, and 1.
%C A228409 a(n) = 1 if and only if n has an odd number of distinct prime factors, A030059. - _Jon Perry_, Nov 12 2013.
%H A228409 Antti Karttunen, <a href="/A228409/b228409.txt">Table of n, a(n) for n = 1..10000</a>
%H A228409 <a href="/index/Eu#epf">Index entries for sequences computed from exponents in factorization of n</a>
%e A228409 a(6) = 9; 4*mu(6) + 5 = 4*1 + 5 = 9.
%p A228409 with(numtheory); A228409:=n->4*mobius(n)+5; seq(A228409(n), n=1..100);
%t A228409 Table[4 MoebiusMu[n] + 5, {n, 100}]
%o A228409 (PARI) a(n)=4*moebius(n)+5 \\ _Charles R Greathouse IV_, Nov 12 2013
%o A228409 (Scheme) (define (A228409 n) (+ 5 (* 4 (A008683 n)))) ;; _Antti Karttunen_, Jul 26 2017
%Y A228409 Cf. A000040, A001358, A008683, A030059.
%K A228409 nonn,easy
%O A228409 1,1
%A A228409 _Wesley Ivan Hurt_, Nov 09 2013