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.

Showing 1-1 of 1 results.

A354866 Dirichlet inverse of A122111.

Original entry on oeis.org

1, -2, -4, 1, -8, 10, -16, -1, 7, 20, -32, -10, -64, 40, 46, 2, -128, -27, -256, -20, 92, 80, -512, 14, 37, 160, -17, -40, -1024, -150, -2048, -3, 184, 320, 202, 53, -4096, 640, 368, 28, -8192, -300, -16384, -80, -146, 1280, -32768, -26, 175, -129, 736, -160, -65536, 85, 404, 56, 1472, 2560, -131072, 242, -262144
Offset: 1

Views

Author

Antti Karttunen, Jun 09 2022

Keywords

Crossrefs

Cf. A122111, A354867, A354868 (parity), A354869 (positions of odd terms).

Programs

  • PARI
    A122111(n) = if(1==n,n,my(f=factor(n), es=Vecrev(f[,2]),is=concat(apply(primepi,Vecrev(f[,1])),[0]),pri=0,m=1); for(i=1, #es, pri += es[i]; m *= prime(pri)^(is[i]-is[1+i])); (m));
    memoA354866 = Map();
    A354866(n) = if(1==n,1,my(v); if(mapisdefined(memoA354866,n,&v), v, v = -sumdiv(n,d,if(dA122111(n/d)*A354866(d),0)); mapput(memoA354866,n,v); (v)));

Formula

a(1) = 1, and for n > 1, a(n) = -Sum_{d|n, dA122111(n/d) * a(d).
a(n) = A354867(n) - A122111(n).
Showing 1-1 of 1 results.