The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A275812 Sum of exponents larger than one in the prime factorization of n: A001222(n) - A056169(n). 13
 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 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,4 LINKS Antti Karttunen, Table of n, a(n) for n = 1..10000 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 MATHEMATICA Table[Total@ Map[Last, Select[FactorInteger@ n, Last@ # > 1 &] /. {} -> {{0, 0}}], {n, 120}] (* Michael De Vlieger, Aug 11 2016 *) PROG (Scheme, two variants, the first one with memoizing definec-macro) (definec (A275812 n) (if (= 1 n) 0 (+ (if (> (A067029 n) 1) (A067029 n) 0) (A275812 (A028234 n))))) (define (A275812 n) (- (A001222 n) (A056169 n))) (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 (PARI) a(n) = my(f = factor(n)); sum(k=1, #f~, if (f[k, 2] > 1, f[k, 2])); \\ Michel Marcus, Jul 19 2017 CROSSREFS Cf. A001222, A028234, A056169, A056170, A057521, A067029, A247180. Differs from A212172 for the first time at n=36, where a(36)=4, while A212172(36)=2. Sequence in context: A216226 A123391 A212172 * A280683 A171871 A076260 Adjacent sequences:  A275809 A275810 A275811 * A275813 A275814 A275815 KEYWORD nonn AUTHOR Antti Karttunen, Aug 11 2016 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified June 7 01:26 EDT 2020. Contains 334836 sequences. (Running on oeis4.)