A325975 a(n) = gcd(A325977(n), A325978(n)).

1, 1, 1, 1, 1, 6, 1, 1, 1, 2, 1, 4, 1, 2, 3, 1, 1, 3, 1, 2, 1, 2, 1, 12, 1, 2, 1, 4, 1, 6, 1, 1, 3, 2, 1, 1, 1, 2, 1, 2, 1, 6, 1, 4, 3, 2, 1, 4, 1, 1, 3, 2, 1, 6, 1, 8, 1, 2, 1, 12, 1, 2, 1, 1, 1, 6, 1, 2, 3, 2, 1, 3, 1, 2, 1, 4, 1, 6, 1, 2, 1, 2, 1, 4, 1, 2, 3, 4, 1, 18, 7, 4, 1, 2, 5, 12, 1, 1, 3, 1, 1, 6, 1, 2, 3

Offset: 1

Author

Antti Karttunen, Jun 02 2019

Keywords

• nonn

Comments

See comments in A325979 and A325981.

Links

• Antti Karttunen, Table of n, a(n) for n = 1..65537

Crossrefs

Cf. A034460, A048107, A325313, A325314, A325814, A325973, A325974, A325977, A325978, A325979, A325981.

Cf. also A009194, A325385, A325813, A326046, A326047, A326048, A326056, A326057, A326060, A326062.

Programs

PARI

A048250(n) = factorback(apply(p -> p+1,factor(n)[,1]));
A162296(n) = sumdiv(n, d, d*(1-issquarefree(d)));
A325313(n) = (A048250(n) - n);
A325314(n) = (n - A162296(n));
A034448(n) = { my(f=factorint(n)); prod(k=1, #f~, 1+(f[k, 1]^f[k, 2])); }; \\ After code in A034448
A034460(n) = (A034448(n) - n);
A048146(n) = (sigma(n)-A034448(n));
A325814(n) = (n-A048146(n));
A325977(n) = ((A034460(n)+A325313(n))/2);
A325978(n) = ((A325314(n)+A325814(n))/2);
A325975(n) = gcd(A325977(n), A325978(n));
\\ Or alternatively, as:
A325975(n) = (1/2)*gcd((A034460(n)+A325313(n)),(A325814(n)+A325314(n)));

Formula

a(n) = gcd(A325977(n), A325978(n)).

a(n) = (1/2)*gcd(A034460(n)+A325313(n), A325814(n)+A325314(n)).