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.

Showing 1-2 of 2 results.

A061026 Smallest number m such that phi(m) is divisible by n, where phi = Euler totient function A000010.

Original entry on oeis.org

1, 3, 7, 5, 11, 7, 29, 15, 19, 11, 23, 13, 53, 29, 31, 17, 103, 19, 191, 25, 43, 23, 47, 35, 101, 53, 81, 29, 59, 31, 311, 51, 67, 103, 71, 37, 149, 191, 79, 41, 83, 43, 173, 69, 181, 47, 283, 65, 197, 101, 103, 53, 107, 81, 121, 87, 229, 59, 709, 61, 367, 311, 127, 85
Offset: 1

Views

Author

Melvin J. Knight (knightmj(AT)juno.com), May 25 2001

Keywords

Comments

Conjecture: a(n) is odd for all n. Verified up to n <= 3*10^5. - Jianing Song, Feb 21 2021
The conjecture above is false because a(16842752) = 33817088; see A002181 and A143510. - Flávio V. Fernandes, Oct 08 2023

Examples

			a(48) = 65 because phi(65) = phi(5)*phi(13) = 4*12 = 48 and no smaller integer m has phi(m) divisible by 48.
		

Crossrefs

Cf. A233516, A233517 (records).
Cf. A005179 (analog for number of divisors), A070982 (analog for sum of divisors).

Programs

  • Mathematica
    a = ConstantArray[1, 64]; k = 1; While[Length[vac = Rest[Flatten[Position[a, 1]]]] > 0, k++; a[[Intersection[Divisors[EulerPhi[k]], vac]]] *= k]; a  (* Ivan Neretin, May 15 2015 *)
  • PARI
    a(n) = my(s=1); while(eulerphi(s)%n, s++); s;
    vector(100, n, a(n))
    
  • Python
    from sympy import totient as phi
    def a(n):
      k = 1
      while phi(k)%n != 0: k += 1
      return k
    print([a(n) for n in range(1, 65)]) # Michael S. Branicky, Feb 21 2021

Formula

Sequence is unbounded; a(n) <= n^2 since phi(n^2) is always divisible by n.
If n+1 is prime then a(n) = n+1.
a(n) = min{ k : phi(k) == 0 (mod n) }.
a(n) = a(2n) for odd n > 1. - Jianing Song, Feb 21 2021

A233517 Record values in A061026, the smallest number m such that n divides phi(m), where phi is Euler's totient function.

Original entry on oeis.org

1, 3, 7, 11, 29, 53, 103, 191, 311, 709, 1021, 1091, 1597, 2339, 3547, 5449, 8243, 9337, 13711, 16673, 17579, 18899, 25367, 37217, 62207, 74441, 87869, 94439, 94789, 96353, 114013, 171167, 229981, 397253, 424769, 432781, 496747, 542599, 583397, 673451, 741677
Offset: 1

Views

Author

T. D. Noe, Feb 12 2014

Keywords

Comments

See A233516 for the n that produce these values. After the initial 1, these numbers are prime.

Crossrefs

Programs

  • Mathematica
    t2 = {{1, 1}}; Do[k = 1; While[Mod[EulerPhi[k], n] > 0, k++]; If[k > t2[[-1,2]], AppendTo[t2, {n, k}]; Print[{n, k}]], {n, 2, 10^3}]; Transpose[t2][[2]]
  • PARI
    lista(cmax) = {my(v = vector(cmax), c = 0, k = 1, d, vm = 0); while(c < cmax, d = divisors(eulerphi(k)); for(i = 1, #d, if(d[i] <= cmax && v[d[i]] == 0, c++; v[d[i]] = k)); k++); for(i = 1, cmax, if(v[i] > vm, vm = v[i]; print1(v[i], ", ")));} \\ Amiram Eldar, May 26 2024
Showing 1-2 of 2 results.