The OEIS is supported by the many generous donors to the OEIS Foundation.

 Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 59th year, we have over 358,000 sequences, and we’ve crossed 10,300 citations (which often say “discovered thanks to the OEIS”). Other ways to Give
 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 | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified December 5 18:02 EST 2022. Contains 358588 sequences. (Running on oeis4.)