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.

A317531 Expansion of Sum_{p prime, k>=1} x^(p^k)/(1 + x^(p^k)).

This page as a plain text file.
%I A317531 #11 Oct 05 2018 11:21:32
%S A317531 0,1,1,0,1,0,1,-1,2,0,1,-1,1,0,2,-2,1,-1,1,-1,2,0,1,-2,2,0,3,-1,1,-1,
%T A317531 1,-3,2,0,2,-2,1,0,2,-2,1,-1,1,-1,3,0,1,-3,2,-1,2,-1,1,-2,2,-2,2,0,1,
%U A317531 -2,1,0,3,-4,2,-1,1,-1,2,-1,1,-3,1,0,3,-1,2,-1,1,-3,4,0,1,-2,2,0,2,-2,1,-2,2,-1,2,0,2
%N A317531 Expansion of Sum_{p prime, k>=1} x^(p^k)/(1 + x^(p^k)).
%H A317531 Antti Karttunen, <a href="/A317531/b317531.txt">Table of n, a(n) for n = 1..65537</a>
%F A317531 G.f.: Sum_{k>=1} x^A246655(k)/(1 + x^A246655(k)).
%F A317531 L.g.f.: log(Product_{p prime, k>=1} (1 + x^(p^k))^(1/p^k)) = Sum_{n>=1} a(n)*x^n/n.
%F A317531 a(n) = Sum_{d|n} (-1)^(n/d+1)*A069513(d).
%F A317531 If n is odd, a(n) = A001222(n).
%t A317531 nmax = 95; Rest[CoefficientList[Series[Sum[Boole[PrimePowerQ[k]] x^k/(1 + x^k), {k, 1, nmax}], {x, 0, nmax}], x]]
%t A317531 nmax = 95; Rest[CoefficientList[Series[Log[Product[(1 + Boole[PrimePowerQ[k]] x^k)^(1/k), {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]]
%t A317531 Table[DivisorSum[n, (-1)^(n/# + 1) &, PrimePowerQ[#] &], {n, 95}]
%o A317531 (PARI) A317531(n) = sumdiv(n,d,((-1)^(n/d+1))*(1==omega(d))); \\ _Antti Karttunen_, Sep 30 2018
%Y A317531 Cf. A001222, A048272, A069513, A246655, A288571, A305614, A317528.
%K A317531 sign
%O A317531 1,9
%A A317531 _Ilya Gutkovskiy_, Jul 30 2018