A348496 a(n) = gcd(A018804(n), A347130(n)) / A003557(n).

1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 4, 1, 3, 1, 2, 1, 21, 1, 36, 5, 1, 1, 1, 1, 15, 3, 4, 1, 1, 1, 1, 7, 1, 3, 1, 1, 3, 1, 1, 1, 1, 1, 12, 3, 5, 1, 10, 1, 3, 5, 4, 1, 9, 1, 1, 1, 1, 1, 12, 1, 3, 1, 1, 9, 1, 1, 12, 1, 1, 1, 1, 1, 3, 1, 4, 3, 1, 1, 2, 1, 1, 1, 4, 11, 15, 1, 35, 1, 3, 5, 36, 1, 1, 3, 1, 1, 3, 3, 1, 1

Offset

1,6

Links

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

Formula

a(n) = A348495(n) / A003557(n).

a(n) = gcd(A347128(n), A347129(n)).

Mathematica

Table[GCD[Total@ GCD[n, Range[n]], DivisorSum[n, #*(If[# < 2, 0, # Total[#2/#1 & @@@ FactorInteger[#]]] &[n/#]) &]]/Apply[Times, Map[#1^(#2 - 1) & @@ # &, FactorInteger[n]]], {n, 101}] (* Michael De Vlieger, Oct 21 2021 *)

Prog

(PARI)
\\ Needs also code from A348495:
A003557(n) = (n/factorback(factorint(n)[, 1]));
A348496(n) = (A348495(n)/A003557(n));

Crossrefs

Cf. A003557, A018804, A347128, A347129, A347130, A348495, A348501 [= a(A276086(n))].

Cf. also A348494, A348498.

Sequence in context: A165485 A179947 A336468 * A344757 A046623 A046602

Adjacent sequences: A348493 A348494 A348495 * A348497 A348498 A348499

Keyword

nonn

Author

Antti Karttunen, Oct 21 2021

Status

approved