A348501 -id:A348501 - cp's OEIS frontend

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

Antti Karttunen, Oct 21 2021

nonn

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

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

Cf. also A348494, A348498.

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 *)

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

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

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

A348500 a(n) = A348494(A276086(n)).

Original entry on oeis.org

1, 1, 1, 5, 1, 7, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 3, 5, 1, 1, 1, 3, 1, 1, 13, 1, 1, 1, 1, 1, 1, 1, 1, 13, 3, 1, 1, 1, 1, 3, 1, 1, 1, 13, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 19, 1, 1, 19, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1

Offset: 0

Antti Karttunen, Oct 21 2021

Antti Karttunen, Table of n, a(n) for n = 0..11550
Index entries for sequences related to primorial base

Cf. A003557, A276086, A348492, A348494.

Cf. also A348501.

Block[{b = MixedRadix[Reverse@ Prime@ Range@ 12]}, Map[GCD[Total@GCD[#1, Range[#1]], #1 Total[#2/#1 & @@@ #2]]/Apply[Times, Map[#1^(#2 - 1) & @@ # &, #2]] & @@ {#, FactorInteger[#]} &, Table[Function[k, Times @@ Power @@@ # &@ Transpose@ {Prime@ Range@ Length@ k, Reverse@ k}]@ IntegerDigits[n, b], {n, 0, 105}] ] ] (* Michael De Vlieger, Oct 21 2021 *)

\\ Needs also code from A348494:
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A348500(n) = A348494(A276086(n));

a(n) = A348494(A276086(n)), where A348494(n) = A348492(n) / A003557(n).

A348502 a(n) = A348498(A276086(n)).

Original entry on oeis.org

0, 1, 1, 5, 1, 7, 1, 7, 8, 31, 13, 41, 1, 9, 11, 37, 8, 47, 3, 11, 14, 43, 19, 53, 2, 13, 17, 49, 11, 59, 1, 9, 10, 41, 17, 55, 12, 59, 71, 247, 106, 317, 19, 73, 92, 289, 127, 359, 26, 87, 113, 331, 148, 401, 33, 101, 134, 373, 169, 443, 1, 11, 13, 47, 10, 61, 17, 69, 86, 277, 121, 347, 12, 83, 107, 319, 71, 389, 31

Offset: 0

Antti Karttunen, Oct 23 2021

Antti Karttunen, Table of n, a(n) for n = 0..11550
Index entries for sequences related to primorial base

Cf. A276086, A348498, A348500, A348501.

Cf. also A342002 for a similar scatter plot.

b = MixedRadix[Reverse@ Prime@ Range@ 12]; f[n_] := If[n < 2, 0, n Total[#2/#1 & @@@ FactorInteger[n]]]; Map[Function[n, GCD[f[n], DivisorSum[n, # f[n/#] &]]*Apply[Times, FactorInteger[n][[All, 1]]]/n], Array[Times @@ Power @@@ # &@ Transpose@ {Prime@ Range@ Length@ #, Reverse@ #} &@ IntegerDigits[#, b] &, 79, 0]] (* Michael De Vlieger, Oct 25 2021 *)

\\ Needs also code from A348498:
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A348502(n) = A348498(A276086(n));

cp's OEIS Frontend

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

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A348500 a(n) = A348494(A276086(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A348502 a(n) = A348498(A276086(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI