[go: up one dir, main page]

login
a(n) = gcd(A000120(n), A294898(n)) = gcd(A000120(n), sigma(n)-A005187(n)).
10

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

%S 1,1,2,1,2,2,3,1,1,2,1,2,3,1,2,1,2,1,3,2,1,1,2,2,1,1,2,3,4,4,5,1,2,2,

%T 1,1,3,1,2,2,1,3,2,1,4,4,1,2,1,1,2,3,4,4,1,1,2,2,1,4,5,1,2,1,2,2,3,2,

%U 3,1,2,1,3,1,2,3,2,4,1,2,1,1,2,1,2,4,1,1,4,2,5,4,1,1,2,2,3,1,2,1,4,4,1,1,2

%N a(n) = gcd(A000120(n), A294898(n)) = gcd(A000120(n), sigma(n)-A005187(n)).

%H Antti Karttunen, <a href="/A326130/b326130.txt">Table of n, a(n) for n = 1..65537</a>

%H <a href="/index/Bi#binary">Index entries for sequences related to binary expansion of n</a>

%H <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>

%F a(n) = gcd(A000120(n), A294898(n)) = gcd(A000120(n), A000203(n)-A005187(n)).

%o (PARI)

%o A005187(n) = { my(s=n); while(n>>=1, s+=n); s; };

%o A294898(n) = (A005187(n)-sigma(n));

%o A326130(n) = gcd(hammingweight(n), A294898(n));

%o (PARI) A326130(n) = gcd(hammingweight(n), sigma(n)-A005187(n));

%Y Cf. A000120, A000203, A005187, A033879, A326131.

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

%K nonn

%O 1,3

%A _Antti Karttunen_, Jun 09 2019