|
%I #6 Dec 28 2020 09:51:22
%S 0,0,4,0,6,0,6,8,10,0,12,0,12,8,16,0,18,0,18,20,22,0,24,24,24,24,28,0,
%T 30,0,30,32,34,24,36,0,36,32,40,0,42,0,42,44,46,0,48,48,48,48,52,0,54,
%U 48,54,56,58,0,60,0,60,56,64,60,66,0,66,68,70,0,72,0,72,72,76,72,78,0,78,80,82,0,84,84,84,80
%N a(n) = 1 + n - A143771(n).
%H Antti Karttunen, <a href="/A339873/b339873.txt">Table of n, a(n) for n = 2..65539</a>
%o (PARI)
%o A143771(n) = my(d = divisors(n)); gcd(vector(#d, k, d[k]+n/d[k])); \\ From A143771.
%o A339873(n) = (1+n-A143771(n));
%Y Cf. A143771, A339914.
%K nonn
%O 2,3
%A _Antti Karttunen_, Dec 25 2020
|
