|
%I #9 Mar 30 2021 15:18:37
%S 1,1,1,5,1,7,1,3,5,11,1,13,1,5,2,17,1,19,1,7,11,23,1,25,13,9,7,29,1,
%T 31,1,11,17,35,3,37,1,13,5,41,1,43,1,15,23,47,1,49,25,17,13,53,1,55,7,
%U 19,29,59,1,61,1,21,8,65,11,67,1,23,35,71,1,73,1,25,19,77,13,79,1,27,41,83,1,85,43,29,11,89,1,91,23
%N a(n) = (1+n) / A143771(n).
%H Antti Karttunen, <a href="/A342918/b342918.txt">Table of n, a(n) for n = 1..16384</a>
%F a(n) = (1+n) / A143771(n).
%o (PARI)
%o A143771(n) = my(d = divisors(n)); gcd(vector(#d, k, d[k]+n/d[k])); \\ From A143771.
%o A342918(n) = ((1+n)/A143771(n));
%Y Cf. A143771, A339873.
%Y After n=1 differs from A342916 for the first time at n=44, where a(44) = 15, while A342916(44) = 5.
%K nonn
%O 1,4
%A _Antti Karttunen_, Mar 29 2021
|
