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-3 of 3 results.

A014197 Number of numbers m with Euler phi(m) = n.

Original entry on oeis.org

2, 3, 0, 4, 0, 4, 0, 5, 0, 2, 0, 6, 0, 0, 0, 6, 0, 4, 0, 5, 0, 2, 0, 10, 0, 0, 0, 2, 0, 2, 0, 7, 0, 0, 0, 8, 0, 0, 0, 9, 0, 4, 0, 3, 0, 2, 0, 11, 0, 0, 0, 2, 0, 2, 0, 3, 0, 2, 0, 9, 0, 0, 0, 8, 0, 2, 0, 0, 0, 2, 0, 17, 0, 0, 0, 0, 0, 2, 0, 10, 0, 2, 0, 6, 0, 0, 0, 6, 0, 0, 0, 3
Offset: 1

Views

Author

Keywords

Comments

Carmichael conjectured that there are no 1's in this sequence. - Jud McCranie, Oct 10 2000
Number of cyclotomic polynomials of degree n. - T. D. Noe, Aug 15 2003
Let v == 0 (mod 24), w = v + 24, and v < k < q < w, where k and q are integer. It seems that, for most values of v, there is no b such that b = a(k) + a(q) and b > a(v) + a(w). The first case where b > a(v) + a(w) occurs at v = 888: b = a(896) + a(900) = 15 + 4, b > a(888) + a(912), or 19 > 8 + 7. The first case where v < n < w and a(n) > a(v) + a(w) occurs at v = 2232: a(2240) > a(2232) + a(2256), or 27 > 7 + 8. - Sergey Pavlov, Feb 05 2017
One elementary result relating to phi(m) is that if m is odd, then phi(m)=phi(2m) because 1 and 2 both have phi value 1 and phi is multiplicative. - Roderick MacPhee, Jun 03 2017

References

  • Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004, Section B39, pp. 144-146.
  • Joe Roberts, Lure of The Integers, The Mathematical Association of America, 1992, entry 32, page 182.

Crossrefs

Cf. A000010, A002202, A032446 (bisection), A049283, A051894, A055506, A057635, A057826, A058277 (nonzero terms), A058341, A063439, A066412, A070243 (partial sums), A070633, A071386 (positions of odd terms), A071387, A071388 (positions of primes), A071389 (where prime(n) occurs for the first time), A082695, A097942 (positions of records), A097946, A120963, A134269, A219930, A280611, A280709, A280712, A296655 (positions of positive even terms), A305353, A305656, A319048, A322019.
For records see A131934.
Column 1 of array A320000.

Programs

  • GAP
    a := function(n)
    local S, T, R, max, i, k, r;
    S:=[];
    for i in DivisorsInt(n)+1 do
        if IsPrime(i)=true then
            S:=Concatenation(S,[i]);
        fi;
    od;
    T:=[];
    for k in [1..Size(S)] do
        T:=Concatenation(T,[S[k]/(S[k]-1)]);
    od;
    max := n*Product(T);
    R:=[];
    for r in [1..Int(max)] do
        if Phi(r)=n then
            R:=Concatenation(R,[r]);
        fi;
    od;
    return Size(R);
    end; # Miles Englezou, Oct 22 2024
  • Magma
    [#EulerPhiInverse(n): n in [1..100]]; // Marius A. Burtea, Sep 08 2019
    
  • Maple
    with(numtheory): A014197:=n-> nops(invphi(n)): seq(A014197(n), n=1..200);
  • Mathematica
    a[1] = 2; a[m_?OddQ] = 0; a[m_] := Module[{p, nmax, n, k}, p = Select[ Divisors[m]+1, PrimeQ]; nmax = m*Times @@ (p/(p - 1)); n = m; k = 0; While[n <= nmax, If[EulerPhi[n] == m, k++]; n++]; k]; Array[a, 92] (* Jean-François Alcover, Dec 09 2011, updated Apr 25 2016 *)
    With[{nn = 116}, Function[s, Function[t, Take[#, nn] &@ ReplacePart[t, Map[# -> Length@ Lookup[s, #] &, Keys@ s]]]@ ConstantArray[0, Max@ Keys@ s]]@ KeySort@ PositionIndex@ Array[EulerPhi, Floor[nn^(3/2)] + 10]] (* Michael De Vlieger, Jul 19 2017 *)
  • PARI
    A014197(n,m=1) = { n==1 && return(1+(m<2)); my(p,q); sumdiv(n, d, if( d>=m && isprime(d+1), sum( i=0,valuation(q=n\d,p=d+1), A014197(q\p^i,p))))} \\ M. F. Hasler, Oct 05 2009
    
  • PARI
    a(n) = invphiNum(n); \\ Amiram Eldar, Nov 15 2024 using Max Alekseyev's invphi.gp
    
  • Python
    from sympy import totient, divisors, isprime, prod
    def a(m):
        if m == 1: return 2
        if m % 2: return 0
        X = (x + 1 for x in divisors(m))
        nmax=m*prod(i/(i - 1) for i in X if isprime(i))
        n=m
        k=0
        while n<=nmax:
            if totient(n)==m:k+=1
            n+=1
        return k
    print([a(n) for n in range(1, 51)]) # Indranil Ghosh, Jul 18 2017, after Mathematica code
    

Formula

Dirichlet g.f.: Sum_{n>=1} a(n)*n^-s = zeta(s)*Product_(1+1/(p-1)^s-1/p^s). - Benoit Cloitre, Apr 12 2003
Limit_{n->infinity} (1/n) * Sum_{k=1..n} a(k) = zeta(2)*zeta(3)/zeta(6) = 1.94359643682075920505707036... (see A082695). - Benoit Cloitre, Apr 12 2003
From Christopher J. Smyth, Jan 08 2017: (Start)
Euler transform = Product_{n>=1} (1-x^n)^(-a(n)) = g.f. of A120963.
Product_{n>=1} (1+x^n)^a(n)
= Product_{n>=1} ((1-x^(2n))/(1-x^n))^a(n)
= Product_{n>=1} (1-x^n)^(-A280712(n))
= Euler transform of A280712 = g.f. of A280611.
(End)
a(A000010(n)) = A066412(n). - Antti Karttunen, Jul 18 2017
From Antti Karttunen, Dec 04 2018: (Start)
a(A000079(n)) = A058321(n).
a(A000142(n)) = A055506(n).
a(A017545(n)) = A063667(n).
a(n) = Sum_{d|n} A008683(n/d)*A070633(d).
a(n) = A056239(A322310(n)).
(End)

A082695 Decimal expansion of zeta(2)*zeta(3)/zeta(6).

Original entry on oeis.org

1, 9, 4, 3, 5, 9, 6, 4, 3, 6, 8, 2, 0, 7, 5, 9, 2, 0, 5, 0, 5, 7, 0, 7, 0, 3, 6, 2, 5, 7, 4, 7, 6, 3, 4, 3, 7, 1, 8, 7, 8, 5, 8, 5, 0, 1, 7, 6, 7, 8, 0, 5, 7, 1, 6, 0, 2, 6, 6, 3, 5, 6, 8, 8, 9, 0, 0, 5, 3, 4, 9, 5, 0, 6, 9, 3, 5, 5, 4, 0, 5, 3, 9, 4, 8, 1, 7, 9, 1, 0, 0, 8, 2, 1, 1, 1, 1, 3, 0, 1, 0, 6, 9, 0, 5
Offset: 1

Views

Author

Benoit Cloitre, Apr 12 2003

Keywords

Comments

Equals the Dirichlet zeta-function Sum_{n>=1} A001615(n)/n^s at s=3. - R. J. Mathar, Apr 02 2011
Dressler shows that this is the average value of A014197, that is, the number of values m such that phi(m) <= n is asymptotically n times this constant. Erdős had shown earlier that this limit exists. - Charles R Greathouse IV, Nov 26 2013
From Stanislav Sykora, Nov 14 2014: (Start)
Equals lim_{n->infinity} (Sum_{k=1..n} k/phi(k))/n, i.e., the limit mean value of k/phi(k), where phi(k) is Euler's totient function.
Also equals lim_{n->infinity} (Sum_{k=1..n} 1/phi(k))/log(n).
Proofs are trivial using the formulas for Sum_{k=1..n} k/phi(k) and Sum_{k=1..n} 1/phi(k) listed in the Wikipedia link.
For the limit mean value of phi(k)/k, see A059956. (End)
The asymptotic mean of A005361. - Amiram Eldar, Apr 13 2020

Examples

			1.94359643682075920505707036257476343718785850176780571602663568890 ...
		

References

  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Section 2.7, p. 116.
  • Joe Roberts, Lure of the Integers, Mathematical Association of America, 1992. See p. 74.

Crossrefs

Programs

  • Mathematica
    First@RealDigits[ Zeta[2]*Zeta[3]/Zeta[6], 10, 100]
    RealDigits[ 315 Zeta[3]/(2 Pi^4), 10, 111][[1]] (* Robert G. Wilson v, Aug 11 2014 *)
  • PARI
    zeta(3)*315/2/Pi^4

Formula

Decimal expansion of Product_{p prime} (1+1/p/(p-1)) = zeta(2)*zeta(3)/zeta(6) = 1.94359643682075920505707...
The sum of the reciprocals of the powerful numbers, A001694. - T. D. Noe, May 03 2006
Equals A013661 * A002117 / A013664 = 1 / A068468 = zeta(3) * 315/(2*Pi^4) = zeta(3) * A157292.
Equals Sum_{k>=1} mu(k)^2/(k*phi(k)) (the sum of reciprocals of the squarefree numbers multiplied by their Euler totient function values, A000010). - Amiram Eldar, Aug 18 2020

Extensions

New definition from Eric W. Weisstein, May 04 2006

A082696 Continued fraction expansion of Product_{p prime} (1+1/(p*(p-1))).

Original entry on oeis.org

1, 1, 16, 1, 2, 1, 2, 3, 1, 1, 3, 2, 1, 8, 1, 1, 1, 1, 1, 1, 1, 32, 33, 2, 14, 1, 1, 1, 1, 1, 6, 2, 1, 3, 1, 3, 5, 3, 2, 13, 21, 1, 1, 1, 1, 8, 3, 27, 1, 2, 5, 1, 6, 2, 4, 3, 1, 1, 3, 10, 3, 5, 12, 48, 1, 2, 2, 1, 1, 1, 2, 2, 1, 11, 1, 1, 3, 7, 1, 30, 42, 1, 1, 1, 12, 1, 7, 1, 41, 1, 1, 1, 3, 5, 2, 1, 2, 7
Offset: 0

Views

Author

Benoit Cloitre, Apr 12 2003

Keywords

Crossrefs

Cf. A014197, A070243, A082695 (decimal expansion).

Programs

  • PARI
    contfrac(prodeulerrat(1+1/(p*(p-1)))) \\ Amiram Eldar, Jun 13 2021

Extensions

Offset changed by Andrew Howroyd, Jul 06 2024
Showing 1-3 of 3 results.