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.

A216157 Difference between the sum of the even divisors and the sum of the odd divisors of phi(n).

This page as a plain text file.
%I A216157 #7 Sep 03 2012 14:34:00
%S A216157 1,1,5,1,4,5,4,5,6,5,20,4,13,13,29,4,13,13,20,6,12,13,30,20,13,20,40,
%T A216157 13,24,29,30,29,52,20,65,13,52,29,78,20,32,30,52,12,24,29,32,30,61,52,
%U A216157 70,13,78,52,65,40,30,29,120,24,65,61,116,30,48,61,60,52
%N A216157 Difference between the sum of the even divisors and the sum of the odd divisors of phi(n).
%C A216157 phi(n) : A000010 is the Euler totient function, and even for n > 2.
%C A216157 If n prime, phi(n) = n-1 and a(n) = a((n-1)/2).
%H A216157 Michel Lagneau, <a href="/A216157/b216157.txt">Table of n, a(n) for n = 3..10000</a>
%e A216157 a(13) = 20 because the divisors of phi(13) = 12 are {1, 2, 3, 4, 6, 12} and (12 + 6 + 4 +2) - (3 + 1) = 20.
%p A216157 with(numtheory):for n from 3 to 100 do:x:=divisors(phi(n)):n1:=nops(x):s0:=0:s1:=0:for m from 1 to n1 do: if irem(x[m],2)=0 then s0:=s0+x[m]:else s1:=s1+x[m]:fi:od:if s0>s1  then printf(`%d, `,s0-s1):else fi:od:
%t A216157 Table[Total[Select[Divisors[EulerPhi[n]], EvenQ[#]&]]-Total[Select[Divisors[EulerPhi[n]], OddQ[#]&]], {n,3,80}]
%Y A216157 Cf. A000010, A215947.
%K A216157 nonn
%O A216157 3,3
%A A216157 _Michel Lagneau_, Sep 02 2012