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.

A356474 a(n) = phi(rad(prime(n)-1)), where phi = A000010 and rad = A007947.

This page as a plain text file.
%I A356474 #24 Aug 09 2022 09:02:16
%S A356474 1,1,1,2,4,2,1,2,10,6,8,2,4,12,22,12,28,8,20,24,2,24,40,10,2,4,32,52,
%T A356474 2,6,12,48,16,44,36,8,24,2,82,42,88,8,72,2,6,20,48,72,112,36,28,96,8,
%U A356474 4,1,130,66,8,44,24,92,72,32,120,24,78,80,12,172,56,10,178,120,60
%N A356474 a(n) = phi(rad(prime(n)-1)), where phi = A000010 and rad = A007947.
%H A356474 Michael De Vlieger, <a href="/A356474/b356474.txt">Table of n, a(n) for n = 1..10000</a>
%F A356474 a(n) = Product_{primes p dividing prime(n)-1} (p-1).
%F A356474 a(n) = phi(A077063(n)) = A173557(prime(n)-1).
%e A356474 prime(8) = 19, so a(8) = phi(rad(18)) = phi(6) = 2.
%t A356474 Array[EulerPhi[Times @@ FactorInteger[Prime[#] - 1][[All, 1]]] &, 74] (* _Michael De Vlieger_, Aug 09 2022 *)
%o A356474 (PARI) a(n) = my(f=factor(prime(n)-1)[, 1]); prod(k=1, #f, f[k]-1) \\ following _Michel Marcus_'s program for A173557
%Y A356474 Cf. A000010, A007947, A077063, A173557.
%K A356474 nonn,easy
%O A356474 1,4
%A A356474 _Jianing Song_, Aug 09 2022