cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A348501 a(n) = A348496(A276086(n)), where A348496(n) = A348495(n) / A003557(n).

This page as a plain text file.
%I A348501 #16 Nov 13 2021 22:33:34
%S A348501 1,1,1,5,1,21,1,1,1,1,3,3,1,3,1,3,1,3,1,1,1,1,1,3,1,1,5,35,3,3,1,3,5,
%T A348501 1,1,3,3,1,1,13,3,3,1,3,1,3,1,3,13,3,1,1,1,3,3,1,5,5,39,3,1,3,1,3,1,3,
%U A348501 3,9,3,3,9,9,1,3,1,3,1,3,1,3,1,57,1,3,57,9,15,15,3,9,1,1,1,1,1,3,1,1,1,1,3,3
%N A348501 a(n) = A348496(A276086(n)), where A348496(n) = A348495(n) / A003557(n).
%H A348501 Antti Karttunen, <a href="/A348501/b348501.txt">Table of n, a(n) for n = 0..11550</a>
%H A348501 <a href="/index/Pri#primorialbase">Index entries for sequences related to primorial base</a>
%F A348501 a(n) = A348496(A276086(n)).
%t A348501 Block[{b = MixedRadix[Reverse@ Prime@ Range@ 12]}, Map[Function[n, GCD[Total@ GCD[n, Range[n]], DivisorSum[n, #*(If[# < 2, 0, # Total[#2/#1 & @@@ FactorInteger[#]]] &[n/#]) &]]/Apply[Times, Map[#1^(#2 - 1) & @@ # &, FactorInteger[n]]]], Table[Function[k, Times @@ Power @@@ # &@ Transpose@ {Prime@ Range@ Length@ k, Reverse@ k}]@ IntegerDigits[n, b], {n, 0, 101}] ] ] (* _Michael De Vlieger_, Oct 21 2021 *)
%o A348501 (PARI)
%o A348501 \\ Needs also code from A348496:
%o A348501 A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
%o A348501 A348501(n) = A348496(A276086(n));
%Y A348501 Cf. A003557, A276086, A348495, A348496.
%Y A348501 Cf. also A348500, A348502.
%K A348501 nonn,base
%O A348501 0,4
%A A348501 _Antti Karttunen_, Oct 21 2021