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.

A136566 a(n) = sum of the exponents occurring only once each in the prime-factorization of n.

This page as a plain text file.
%I A136566 #18 Sep 22 2017 03:51:23
%S A136566 0,1,1,2,1,0,1,3,2,0,1,3,1,0,0,4,1,3,1,3,0,0,1,4,2,0,3,3,1,0,1,5,0,0,
%T A136566 0,0,1,0,0,4,1,0,1,3,3,0,1,5,2,3,0,3,1,4,0,4,0,0,1,2,1,0,3,6,0,0,1,3,
%U A136566 0,0,1,5,1,0,3,3,0,0,1,5,4,0,1,2,0,0,0,4,1,2,0,3,0,0,0,6,1,3,3,0,1,0,1,4,0
%N A136566 a(n) = sum of the exponents occurring only once each in the prime-factorization of n.
%H A136566 Diana Mecum and Michael De Vlieger, <a href="/A136566/b136566.txt">Table of n, a(n) for n = 1..10000</a> (first 1000 terms from Diana Mecum)
%e A136566 4200 = 2^3 * 3^1 * 5^2 * 7^1. The exponents of the prime factorization are therefore 3,1,2,1. The exponents occurring exactly once are 2 and 3. So a(4200) = 2+3 = 5.
%t A136566 Table[Total@ Flatten@ Select[Split[Sort[FactorInteger[n][[All, -1]]]], Length@ # == 1 &] - Boole[n == 1], {n, 105}] (* _Michael De Vlieger_, Sep 21 2017 *)
%o A136566 (PARI) a(n) = my(f=factor(n)[,2]); sum(k=1, #f, f[k]*(#select(x->(x==f[k]), f) == 1)); \\ _Michel Marcus_, Sep 22 2017
%Y A136566 Cf. A136565, A136567.
%K A136566 nonn
%O A136566 1,4
%A A136566 _Leroy Quet_, Jan 07 2008
%E A136566 More terms from _Diana L. Mecum_, Jul 17 2008