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.

A332334 Let a(1) = a(2) = 1, and for n > 2 let a(n) = p where p is the largest prime such that p# divides phi(n), where phi is Euler's totient function and # is the primorial.

Original entry on oeis.org

1, 1, 2, 2, 2, 2, 3, 2, 3, 2, 2, 2, 3, 3, 2, 2, 2, 3, 3, 2, 3, 2, 2, 2, 2, 3, 3, 3, 2, 2, 5, 2, 2, 2, 3, 3, 3, 3, 3, 2, 2, 3, 3, 2, 3, 2, 2, 2, 3, 2, 2, 3, 2, 3, 2, 3, 3, 2, 2, 2, 5, 5, 3, 2, 3, 2, 3, 2, 2, 3, 2, 3, 3, 3, 2, 3, 5, 3, 3, 2, 3, 2, 2, 3, 2, 3, 2
Offset: 1

Views

Author

Keywords

Comments

Pollack and Pomerance show that the normal order of a(n) is log log n/log log log n. The maximal order is log n (for primorial primes A018239, by the prime number theorem) and the minimal order, for n > 2, is 2 (for products of Fermat primes A143512, apart from 1).

Crossrefs

Programs

  • PARI
    a(n)=my(ph=eulerphi(n)); my(p=1); forprime(q=2,, if(ph%q, return(p), p=q))