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

 

Logo

Please make a donation to keep the OEIS running. We are now in our 56th year. In the past year we added 10000 new sequences and reached almost 9000 citations (which often say "discovered thanks to the OEIS").
Other ways to donate

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A297168 Difference between A156552 and its Moebius transform: a(n) = A156552(n) - A297112(n). 10
0, 0, 0, 1, 0, 3, 0, 3, 2, 5, 0, 7, 0, 9, 6, 7, 0, 9, 0, 11, 10, 17, 0, 15, 4, 33, 6, 19, 0, 17, 0, 15, 18, 65, 12, 19, 0, 129, 34, 23, 0, 29, 0, 35, 14, 257, 0, 31, 8, 17, 66, 67, 0, 21, 20, 39, 130, 513, 0, 35, 0, 1025, 22, 31, 36, 53, 0, 131, 258, 33, 0, 39, 0, 2049, 18, 259, 24, 101, 0, 47, 14, 4097, 0, 59, 68, 8193, 514, 71, 0, 37, 40 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,6

LINKS

Antti Karttunen, Table of n, a(n) for n = 1..8192

Index entries for sequences computed from indices in prime factorization

FORMULA

a(n) = -Sum_{d|n, d<n} A008683(n/d)*A156552(d).

a(n) = Sum_{d|n, d<n} A297112(d).

For n > 1, a(n) = Sum_{d|n, 1<d<n} 2^A033265(A156552(d)).

a(n) = A156552(n) - A297112(n).

a(1) = 0, for n > 1, a(n) = A156552(n) - 2^A297167(n).

MATHEMATICA

With[{s = Array[Total@ MapIndexed[#1 2^(First@ #2 - 1) &, Flatten@ Map[ConstantArray[2^(PrimePi@ #1 - 1), #2] & @@ # &, FactorInteger@ #]] - Boole[# == 1]/2 &, 91]}, Table[-DivisorSum[n, MoebiusMu[n/#] s[[#]] &, # < n &], {n, Length@ s}]] (* Michael De Vlieger, Mar 13 2018 *)

PROG

(PARI)

A064989(n) = {my(f); f = factor(n); if((n>1 && f[1, 1]==2), f[1, 2] = 0); for (i=1, #f~, f[i, 1] = precprime(f[i, 1]-1)); factorback(f)};

A156552(n) = if(1==n, 0, if(!(n%2), 1+(2*A156552(n/2)), 2*A156552(A064989(n))));

A297112(n) = sumdiv(n, d, moebius(n/d)*A156552(d));

A297168(n) = (A156552(n)-A297112(n));

\\ Or alternatively as:

A297168(n) = -sumdiv(n, d, (d<n)*moebius(n/d)*A156552(d));

(PARI)

A061395(n) = if(1==n, 0, primepi(vecmax(factor(n)[, 1])));

A297167(n) = if(1==n, 0, (A061395(n) + (bigomega(n)-omega(n)) - 1));

A297112(n) = if(1==n, 0, 2^A297167(n));

A297168(n) = sumdiv(n, d, (d<n)*A297112(d)); \\ Antti Karttunen, Mar 13 2018

(Scheme)

(define (A297168 n) (- (A156552 n) (A297112 n)))

(define (A297168 n) (if (= 1 n) 0 (- (A156552 n) (A000079 (A297167 n)))))

CROSSREFS

Cf. A008683, A033265, A156552, A297112, A297112, A297113, A297167, A297169, A300827.

Sequence in context: A271860 A219428 A016035 * A112470 A331924 A201582

Adjacent sequences:  A297165 A297166 A297167 * A297169 A297170 A297171

KEYWORD

nonn

AUTHOR

Antti Karttunen, Feb 27 2018

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.

License Agreements, Terms of Use, Privacy Policy. .

Last modified November 23 21:51 EST 2020. Contains 338603 sequences. (Running on oeis4.)