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

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A326140 a(n) = gcd(A318878(n), A318879(n)). 9

%I #8 Jun 10 2019 17:56:38

%S 1,1,2,1,4,2,6,1,5,2,10,2,12,2,6,1,16,1,18,2,10,2,22,2,19,2,14,6,28,6,

%T 30,1,18,2,22,1,36,2,22,2,40,2,42,2,12,2,46,2,41,1,30,6,52,2,38,2,34,

%U 2,58,6,60,2,22,1,46,6,66,2,42,2,70,1,72,2,26,6,58,2,78,2,41,2,82,2,62,2,54,2,88,6,70,2,58,2

%N a(n) = gcd(A318878(n), A318879(n)).

%H Antti Karttunen, <a href="/A326140/b326140.txt">Table of n, a(n) for n = 1..20000</a>

%o (PARI) A326140(n) = { my(t=0, u=0); fordiv(n,d, d -= 2*eulerphi(d); if(d<0, t -= d, u += d)); gcd(t,u); };

%o (PARI)

%o A318878(n) = sumdiv(n,d,d=(2*eulerphi(d))-d; (d>0)*d);

%o A318879(n) = sumdiv(n,d,d=d-(2*eulerphi(d)); (d>0)*d);

%o A326140(n) = gcd(A318878(n), A318879(n));

%Y Cf. A033879, A083254, A318878, A318879, A326141.

%Y Cf. also A009194, A325385, A325813, A325975, A326046, A326047, A326048, A326056, A326057, A326060, A326062, A326129, A326130, A326144.

%K nonn

%O 1,3

%A _Antti Karttunen_, Jun 09 2019

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.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 08:59 EDT 2024. Contains 371268 sequences. (Running on oeis4.)