|
%I #8 May 26 2017 16:37:37
%S 0,1,7,22,23,67,29,122,79,173,67,408,107,277,328,531,214,742,191,949,
%T 530,631,277,1894,498,905,781,1598,467,2704,497,2149,1178,1600,1228,
%U 4188,743,1771,1656,4282,949,4658,947,3572,3163,2557,1129,8005,1597,4373,2855,4953,1487,7141,2704,7384,3242,4097,1771,14539,1955,4561,5462,8520,3745
%N Compound filter (2-adic valuation of phi(n) & sigma(n)): a(n) = P(A053574(n), A000203(n)), where P(n,k) is sequence A000027 used as a pairing function.
%H Antti Karttunen, <a href="/A286572/b286572.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = (1/2)*(2 + ((A053574(n)+A000203(n))^2) - A053574(n) - 3*A000203(n)).
%o (PARI)
%o A000203(n) = sigma(n);
%o A053574(n) = valuation(eulerphi(n), 2);
%o A286572(n) = (1/2)*(2 + ((A053574(n)+A000203(n))^2) - A053574(n) - 3*A000203(n));
%o (Scheme) (define (A286572 n) (* (/ 1 2) (+ (expt (+ (A053574 n) (A000203 n)) 2) (- (A053574 n)) (- (* 3 (A000203 n))) 2)))
%o (Python)
%o from sympy import totient, divisor_sigma
%o def T(n, m): return ((n + m)**2 - n - 3*m + 2)/2
%o def a007814(n): return 1 + bin(n - 1)[2:].count("1") - bin(n)[2:].count("1")
%o def a(n): return T(a007814(totient(n)), divisor_sigma(n)) # _Indranil Ghosh_, May 26 2017
%Y Cf. A000010, A000027, A000203, A053574, A286360, A286460, A286568.
%K nonn
%O 1,3
%A _Antti Karttunen_, May 26 2017
|
