%I #13 Dec 01 2021 11:07:13
%S 1,1,2,1,4,0,6,1,1,2,10,2,12,4,6,1,16,3,18,2,10,8,22,6,1,10,2,14,28,
%T 12,30,1,18,14,22,1,36,16,22,10,40,12,42,2,6,20,46,14,1,1,30,2,52,12,
%U 38,2,34,26,58,6,60,28,2,1,46,12,66,2,42,4,70,3,72,34,2,2,58,12,78,2,1,38,82,14,62,40,54,2,88,6,70,2
%N a(n) = gcd(A066503(n), A326143(n)) = gcd(n - A007947(n), sigma(n) - A007947(n) - n).
%H Antti Karttunen, <a href="/A326144/b326144.txt">Table of n, a(n) for n = 1..20000</a>
%F a(n) = gcd(A066503(n), A326143(n)) = gcd(n-A007947(n), A000203(n)-A007947(n)-n).
%t rad[n_] := Times @@ (First /@ FactorInteger@n);
%t a[n_] := GCD[n - rad[n], DivisorSigma[1, n] - rad[n] - n];
%t Array[a, 100] (* _Jean-François Alcover_, Dec 01 2021 *)
%o (PARI)
%o A007947(n) = factorback(factorint(n)[, 1]);
%o A066503(n) = (n - A007947(n));
%o A326143(n) = (sigma(n)-A007947(n)-n);
%o A326144(n) = gcd(A066503(n), A326143(n));
%Y Cf. A000203, A007947, A066503, A326142, A326143, A326145.
%Y Cf. also A009194, A325385, A325813, A325975, A326046, A326047, A326048, A326056, A326057, A326060, A326062, A326129, A326130, A326140.
%K nonn
%O 1,3
%A _Antti Karttunen_, Jun 09 2019