login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

a(n) = x/gcd(n,x), where x = 1+A060681(n).
1

%I #5 Jan 01 2021 11:44:08

%S 1,1,1,3,1,2,1,5,7,3,1,7,1,4,11,9,1,5,1,11,5,6,1,13,21,7,19,15,1,8,1,

%T 17,23,9,29,19,1,10,9,21,1,11,1,23,31,12,1,25,43,13,35,27,1,14,9,29,

%U 13,15,1,31,1,16,43,33,53,17,1,35,47,18,1,37,1,19,17,39,67,20,1,41,55,21,1,43,69,22,59,45,1,23

%N a(n) = x/gcd(n,x), where x = 1+A060681(n).

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

%F a(n) = (1+A060681(n)) / A323071(n).

%o (PARI)

%o A060681(n) = (n-if(1==n,n,n/vecmin(factor(n)[,1])));

%o A339913(n) = { my(x=1+A060681(n)); (x/gcd(n,x)); };

%Y Cf. A060681, A323071, A323072.

%K nonn

%O 1,4

%A _Antti Karttunen_, Jan 01 2021