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.

Previous Showing 11-20 of 4244 results. Next

A007755 Smallest number m such that the trajectory of m under iteration of Euler's totient function phi(n) [A000010] contains exactly n distinct numbers, including m and the fixed point.

Original entry on oeis.org

1, 2, 3, 5, 11, 17, 41, 83, 137, 257, 641, 1097, 2329, 4369, 10537, 17477, 35209, 65537, 140417, 281929, 557057, 1114129, 2384897, 4227137, 8978569, 16843009, 35946497, 71304257, 143163649, 286331153, 541073537, 1086374209, 2281701377, 4295098369
Offset: 1

Views

Author

Pepijn van Erp [ vanerp(AT)sci.kun.nl ]

Keywords

Comments

Least integer k such that the number of iterations of Euler phi function needed to reach 1 starting at k (k is counted) is n.
a(n) is smallest number in the class k(n) which groups families of integers which take the same number of iterations of the totient function to evolve to 1. The maximum is 2*3^(n-1).
Shapiro shows that the smallest number is greater than 2^(n-1). Catlin shows that if a(n) is odd and composite, then its factors are among the a(k), k < n. For example a(12) = a(5) a(8). There is a conjecture that all terms of this sequence are odd. - T. D. Noe, Mar 08 2004
The indices of odd prime terms are given by n=A136040(k)+2 for k=1,2,3,.... - T. D. Noe, Dec 14 2007
Shapiro mentions on page 30 of his paper the conjecture that a(n) is prime for each n > 1, but a(13) is composite and so the conjecture fails. - Charles R Greathouse IV, Oct 28 2011

Examples

			a(3) = 3 because trajectory={3,2,1}. n=1: a(1)=1 because trajectory={1}
		

References

  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 83, p. 29, Ellipses, Paris 2008. Also Entry 137, p. 47.
  • R. K. Guy, Unsolved Problems in Number Theory, 2nd Ed. New York: Springer-Verlag, p. 97, 1994, Section B41.

Crossrefs

Cf. A000010, A003434, A049108, A092873 (prime factors of a(n)), A060611, A098196, A227946.
A060611 has the same initial terms but is a different sequence.

Programs

  • Haskell
    a007755 = (+ 1) . fromJust . (`elemIndex` a003434_list) . (subtract 1)
    -- Reinhard Zumkeller, Feb 08 2013, Jul 03 2011
  • Mathematica
    f[n_] := Length[ NestWhileList[ EulerPhi, n, Unequal, 2]] - 1; a = Table[0, {30}]; Do[b = f[n]; If[a[[b]] == 0, a[[b]] = n; Print[n, " = ", b]], {n, 1, 22500000}] (* Robert G. Wilson v *)

Formula

a(n) = smallest m such that A049108(m)=n.
Alternatively, a(n) = smallest m such that A003434(m)=n-1.
a(n+2) ~ 2^n.

Extensions

More terms from David W. Wilson, May 15 1997
Additional comments from James S. Cronen (cronej(AT)rpi.edu)

A057635 a(n) is the largest m such that phi(m) = n, where phi is Euler's totient function = A000010, or a(n) = 0 if no such m exists.

Original entry on oeis.org

2, 6, 0, 12, 0, 18, 0, 30, 0, 22, 0, 42, 0, 0, 0, 60, 0, 54, 0, 66, 0, 46, 0, 90, 0, 0, 0, 58, 0, 62, 0, 120, 0, 0, 0, 126, 0, 0, 0, 150, 0, 98, 0, 138, 0, 94, 0, 210, 0, 0, 0, 106, 0, 162, 0, 174, 0, 118, 0, 198, 0, 0, 0, 240, 0, 134, 0, 0, 0, 142, 0, 270, 0, 0, 0, 0, 0, 158, 0, 330, 0
Offset: 1

Views

Author

Jud McCranie, Oct 10 2000

Keywords

Comments

To check that a property P holds for all EulerPhi(x) not exceeding n, for n with a(n) > 0, it suffices to check P for all EulerPhi(x) with x not exceeding a(n). - Joseph L. Pe, Jan 10 2002
The Alekseyev link in A131883 establishes the following explicit relationship between A131883, A036912 and A057635: for t belonging to A036912, we have t = A131883(A057635(t)-1). In other words, A036912(n) = A131883(A057635(A036912(n))-1) for all n.
From Jianing Song, Feb 16 2019: (Start)
Let f(n) = exp(gamma)*log(log(n)) + 2.5/log(log(n)), then a(n) < n*f(n^2) for all n > 1, where gamma = A001620.
Proof. Without loss of generality we suppose log(log(n)) > n_0 = sqrt(2.5/exp(gamma)) = 1.18475..., then f(n), n/f(n) and N(n) = ceiling(n*f(n^2)) are all monotonically increasing functions of n, and we have f(n) < 2*exp(gamma)*log(log(n)).
By the formula (3.41) in Theorem 15 by J. Barkley Rosser and Lowell Schoenfeld we have phi(k) > k/f(k) for k != 1, 2, 223092870. N(31802157) = 223092869 < 223092870, N(31802158) = 223092877 > 223092870, so N(n) != 223092870 (N(n) is increasing). So phi(N(n)) > N(n)/f(N(n)) > (n*f(n^2))/f(n*f(n^2)) (n/f(n) is increasing and log(log(n*f(n^2))) > n_0).
Note that f(n^2) < 2*exp(gamma)*log(log(n^2)) < 2*exp(gamma)*(log(n^2)/e) = 4*exp(gamma-1)*log(n) < 4*exp(gamma-2)*n < n, so n*f(n^2) < n^2, f(n*f(n^2)) < f(n^2) (f(n) is increasing and log(log(n*f(n^2))) > n_0), so phi(N(n)) > n. As a result, a(n) <= N(n) - 1 < n*f(n^2).
Conjecturally a(n) < n*f(n) for all n > 2. (End)

Examples

			m = 12 is the largest value of m such that phi(m) = 4, so a(4) = 12.
		

Crossrefs

Cf. A006511 (largest k for which A000010(k) = A002202(n)).

Programs

  • Mathematica
    a = Table[0, {100}]; Do[ t = EulerPhi[n]; If[t < 101, a[[t]] = n], {n, 1, 10^6}]; a
  • PARI
    a(n) = if(n%2, 2*(n==1), forstep(k=floor(exp(Euler)*n*log(log(n^2))+2.5*n/log(log(n^2))), n, -1, if(eulerphi(k)==n, return(k)); if(k==n, return(0)))) \\ Jianing Song, Feb 15 2019
    
  • PARI
    apply( {A057635(n,m=istotient(n))=if(!m, 0, n>1, m=log(log(n)*2); m=bitand(n*(exp(Euler)*m+2.5/m)\1,-2); while(eulerphi(m)!=n, m-=2); m, 2)}, [1..99]) \\ If n is known to be a totient, a nonzero 2nd arg can be given to avoid the check. - M. F. Hasler, Aug 13 2021
    
  • PARI
    a(n) = invphiMax(n); \\ Amiram Eldar, Nov 14 2024 using Max Alekseyev's invphi.gp

Formula

a(2n+1) = 0 for n > 0, and a(2n) = 0 iff 2n is in A005277.

Extensions

Edited and escape clause added to definition by M. F. Hasler, Aug 13 2021

A109395 Denominator of phi(n)/n = Product_{p|n} (1 - 1/p); phi(n)=A000010(n), the Euler totient function.

Original entry on oeis.org

1, 2, 3, 2, 5, 3, 7, 2, 3, 5, 11, 3, 13, 7, 15, 2, 17, 3, 19, 5, 7, 11, 23, 3, 5, 13, 3, 7, 29, 15, 31, 2, 33, 17, 35, 3, 37, 19, 13, 5, 41, 7, 43, 11, 15, 23, 47, 3, 7, 5, 51, 13, 53, 3, 11, 7, 19, 29, 59, 15, 61, 31, 7, 2, 65, 33, 67, 17, 69, 35, 71, 3, 73, 37, 15, 19, 77, 13, 79, 5, 3
Offset: 1

Views

Author

Franz Vrabec, Aug 26 2005

Keywords

Comments

a(n)=2 iff n=2^k (k>0); otherwise a(n) is odd. If p is prime, a(p)=p; the converse is false, e.g.: a(15)=15. It is remarkable that this sequence often coincides with A006530, the largest prime P dividing n. Theorem: a(n)=P if and only if for every prime p < P in n there is some prime q in n with p|(q-1). - Franz Vrabec, Aug 30 2005

Examples

			a(10) = 10/gcd(10,phi(10)) = 10/gcd(10,4) = 10/2 = 5.
		

Crossrefs

Cf. A076512 for the numerator.
Phi(m)/m = k: A000079 \ {1} (k=1/2), A033845 (k=1/3), A000244 \ {1} (k=2/3), A033846 (k=2/5), A000351 \ {1} (k=4/5), A033847 (k=3/7), A033850 (k=4/7), A000420 \ {1} (k=6/7), A033848 (k=5/11), A001020 \ {1} (k=10/11), A288162 (k=6/13), A001022 \ {1} (12/13), A143207 (k=4/15), A033849 (k=8/15), A033851 (k=24/35).

Programs

Formula

a(n) = n/gcd(n, phi(n)) = n/A009195(n).
From Antti Karttunen, Feb 09 2019: (Start)
a(n) = denominator of A173557(n)/A007947(n).
a(2^n) = 2 for all n >= 1.
(End)
From Amiram Eldar, Jul 31 2020: (Start)
Asymptotic mean of phi(n)/n: lim_{m->oo} (1/m) * Sum_{n=1..m} A076512(n)/a(n) = 6/Pi^2 (A059956).
Asymptotic mean of n/phi(n): lim_{m->oo} (1/m) * Sum_{n=1..m} a(n)/A076512(n) = zeta(2)*zeta(3)/zeta(6) (A082695). (End)

A077816 Wieferich numbers (1): n > 1 such that 2^A000010(n) == 1 (mod n^2).

Original entry on oeis.org

1093, 3279, 3511, 7651, 10533, 14209, 17555, 22953, 31599, 42627, 45643, 52665, 68859, 94797, 99463, 127881, 136929, 157995, 228215, 298389, 410787, 473985, 684645, 895167, 1232361, 2053935, 2685501, 3697083, 3837523, 6161805, 11512569
Offset: 1

Views

Author

Reinhard Zumkeller, Nov 17 2002

Keywords

Comments

A077815(a(n)) = 1.
The only known primes are a(1)=A001220(1)=1093 and a(3)=A001220(2)=3511, the Wieferich primes.
If there are finitely many Wieferich primes (A001220), this sequence is finite. In particular, unless there are other Wieferich primes besides 1093 and 3511, this sequence consists of 104 terms with the largest being 16547533489305 (Agoh et al., 1997).
a(105)=A001220(3) in the sense that either both numbers are well-defined and equal, or else neither number exists. - Jeppe Stig Nielsen, Oct 16 2016

Examples

			A077815(3279) = 2^A000010(3279) mod 3279^2 = 2^2184 mod 10751841 = 1, therefore 3279 is a term.
		

Crossrefs

For another definition of Wieferich numbers, see A182297.
Cf. A001220.

Programs

  • Magma
    [n: n in [1..8*10^5] | 2^EulerPhi(n) mod n^2 eq 1]; // Vincenzo Librandi, Dec 05 2015
  • Mathematica
    Reap[For[k = 1, k <= 10^8, k++, If[PowerMod[2, EulerPhi[k], k^2] == 1, Print[k]; Sow[k]]]][[2, 1]] (* Jean-François Alcover, Nov 17 2021 *)
  • PARI
    for(n=2, 10^9, if(Mod(2, n^2)^(eulerphi(n))==1, print1(n, ", "))); \\ Felix Fröhlich, May 27 2014
    

Extensions

More terms from Emeric Deutsch, Mar 05 2005
More terms from Sam Handler (sam_5_5_5_0(AT)yahoo.com), Jun 18 2005

A097945 a(n) = mu(n)*phi(n) where mu(n) is the Mobius function (A008683) and phi(n) is the Euler totient function (A000010).

Original entry on oeis.org

1, -1, -2, 0, -4, 2, -6, 0, 0, 4, -10, 0, -12, 6, 8, 0, -16, 0, -18, 0, 12, 10, -22, 0, 0, 12, 0, 0, -28, -8, -30, 0, 20, 16, 24, 0, -36, 18, 24, 0, -40, -12, -42, 0, 0, 22, -46, 0, 0, 0, 32, 0, -52, 0, 40, 0, 36, 28, -58, 0, -60, 30, 0, 0, 48, -20, -66, 0, 44, -24, -70, 0, -72, 36, 0, 0, 60, -24, -78, 0, 0, 40, -82, 0
Offset: 1

Views

Author

Gerald McGarvey, Sep 04 2004

Keywords

Comments

Also, a(n) = mu(n)*uphi(n) where mu(n) is the Mobius function (A008683) and uphi(n) is the unitary totient function (A047994), since phi(n) = uphi(n) when n is squarefree, while mu(n) = 0 when n is not squarefree. - Franklin T. Adams-Watters, May 14 2006
Conjecture: Sum_{n>=1} mu(n)/phi(n) = Sum_{n>=1} a(n)/phi(n)^2 = 0. It is true that Sum_{n>=1} mu(n)/phi(n)^s = 0 at least for s > 1 since: phi(2)=1, phi is multiplicative, so for n's that are squarefree, the phi(n) values can be partitioned in pairs where phi(m)=phi(2m) and mu(m) = -mu(2m). So Sum_{i=1..n} mu(i)/phi(i)^s < Sum_{j=floor(n/2)..n} 1/phi(j)^s, which approaches 0 as n increases since (1) n^(1-e) < phi(n) < n for any e > 0 and n > N(e) and (2) Sum_{i..n} 1/n^s converges for s > 1. Conjecture: Sum_{n>=1} mu(n)/phi(n)^z = 0 for Re(z) > 1.
Multiplicative with a(p^1) = 1-p, a(p^e) = 0, e > 1. - Mitch Harris, May 24 2005
Row sums of triangle A143153 = a signed version of the sequence such that parity = (-) iff A008683(n) = (+); 0 or (+): (1, 1, 2, 0, 4, -2, 6, 0, 0, -4, 10, 0, 12, -6, 0, 0, 0, ...). - Gary W. Adamson, Jul 27 2008
Dirichlet inverse of A003958. - R. J. Mathar, Jul 08 2011

Crossrefs

Programs

  • Maple
    with(numtheory):
    a:= n-> mobius(n)*phi(n):
    seq(a(n), n=1..100);  # Alois P. Heinz, Aug 06 2012
  • Mathematica
    Table[ MoebiusMu[n]EulerPhi[n], {n, 85}] (* Robert G. Wilson v, Sep 06 2004 *)
  • PARI
    a(n)=moebius(n)*eulerphi(n) \\ Charles R Greathouse IV, Feb 21 2013
    
  • PARI
    for(n=1, 100, print1(direuler(p=2, n, (1 - p*X + X))[n], ", ")) \\ Vaclav Kotesovec, Jun 14 2020

Formula

Dirichlet g.f.: Product_{primes p} (1-p^(1-s)+p^(-s)). - R. J. Mathar, Aug 29 2011
Sum_{d|n} abs(a(d)) = rad(n) = A007947(n). - Rémy Sigrist, Nov 05 2017
Sum_{k=1..n} abs(a(k)) ~ c * n^2, where c = A065464/2 = (1/2) * Product_{primes p} (1 - 2/p^2 + 1/p^3) = 0.21412475283854722... Equivalently, c = A065463 * 3 / Pi^2. - Vaclav Kotesovec, Jun 14 2020
From Antti Karttunen, Aug 20 2021: (Start)
a(n) = mu(n)*A000010(n) = mu(n)*A003958(n) = mu(n)*A047994(n) = mu(n)*A173557(n), where mu is Möbius mu function (A008683).
a(n) = A008966(n) * A023900(n) = abs(mu(n)) * A023900(n).
a(n) = A322581(n) - A003958(n).
(End)

Extensions

More terms from Robert G. Wilson v, Sep 06 2004
Edited by N. J. A. Sloane, May 20 2006

A053287 Euler totient function (A000010) of 2^n - 1.

Original entry on oeis.org

1, 2, 6, 8, 30, 36, 126, 128, 432, 600, 1936, 1728, 8190, 10584, 27000, 32768, 131070, 139968, 524286, 480000, 1778112, 2640704, 8210080, 6635520, 32400000, 44717400, 113467392, 132765696, 533826432, 534600000, 2147483646, 2147483648, 6963536448, 11452896600
Offset: 1

Views

Author

Labos Elemer, Mar 03 2000

Keywords

Comments

Number of elements of multiplicative order 2^n - 1 in GF(2^n).
n divides a(n) because 2^a(n) mod 2^n - 1 is 1, 2^n mod 2^n - 1 is 1, so n | a(n). A011260(n) = a(n)/n. - Jinyuan Wang, Oct 31 2018
The set {a(n)/(2^n-1)} is dense in [0, 1] (Luca, 2003). - Amiram Eldar, Mar 04 2021

Crossrefs

Programs

Formula

a(n) = A000010(A000225(n)).
a(A000079(n-1)) = A058891(n).
a(n) = A000010(2^n-1) or also a(n) = A062401(2^(n-1)) = phi(sigma(2^(n-1))). - Labos Elemer, Jul 19 2004

A055653 Sum of phi(d) [A000010] over all unitary divisors d of n (that is, gcd(d,n/d) = 1).

Original entry on oeis.org

1, 2, 3, 3, 5, 6, 7, 5, 7, 10, 11, 9, 13, 14, 15, 9, 17, 14, 19, 15, 21, 22, 23, 15, 21, 26, 19, 21, 29, 30, 31, 17, 33, 34, 35, 21, 37, 38, 39, 25, 41, 42, 43, 33, 35, 46, 47, 27, 43, 42, 51, 39, 53, 38, 55, 35, 57, 58, 59, 45, 61, 62, 49, 33, 65, 66, 67, 51, 69, 70, 71, 35, 73
Offset: 1

Views

Author

Labos Elemer, Jun 07 2000

Keywords

Comments

Phi-summation over d-s if runs over all divisors is n, so these values do not exceed n. Compare also other "Phi-summations" like A053570, A053571, or distinct primes dividing n, etc.
a(n) is also the number of solutions of x^(k+1)=x mod n for some k>=1. - Steven Finch, Apr 11 2006
An integer a is called regular (mod n) if there is an integer x such that a^2 x == a (mod n). Then a(n) is also the number of regular integers a (mod n) such that 1 <= a <= n. - Laszlo Toth, Sep 04 2008
Equals row sums of triangle A157361 and inverse Mobius transform of A114810. - Gary W. Adamson, Feb 28 2009
a(m) = m iff m is squarefree, a(A005117(n)) = A005117(n). - Reinhard Zumkeller, Mar 11 2012
Apostol & Tóth call this ϱ(n), i.e., varrho(n). - Charles R Greathouse IV, Apr 23 2013

Examples

			n=1260 has 36 divisors of which 16 are unitary ones: {1,4,5,7,9,20,28,35,36,45,63,140,180,252,315,1260}.
EulerPhi values of these divisors are: {1,2,4,6,6,8,12,24,12,24,36,48,48,72,144,288}.
The sum is 735, thus a(1260)=735.
Or, 1260=2^2*3^2*5*7, thus a(1260) = (1 + 2^2 - 2)*(1 + 3^2 - 3)*(1 + 5 - 5^0)*(1 + 7 - 7^0) = 735.
		

References

  • J. Morgado, Inteiros regulares módulo n, Gazeta de Matematica (Lisboa), 33 (1972), no. 125-128, 1-5. [From Laszlo Toth, Sep 04 2008]
  • J. Morgado, A property of the Euler phi-function concerning the integers which are regular modulo n, Portugal. Math., 33 (1974), 185-191.

Crossrefs

Programs

  • Haskell
    a055653 = sum . map a000010 . a077610_row
    -- Reinhard Zumkeller, Mar 11 2012
    
  • Maple
    A055653 := proc(n) local ans, i:ans := 1: for i from 1 to nops(ifactors(n)[ 2 ]) do ans := ans*(1+ifactors(n)[ 2 ][ i ] [ 1 ]^ifactors(n)[ 2 ] [ i ] [ 2 ]-ifactors(n)[ 2 ][ i ] [ 1 ]^(ifactors(n)[ 2 ] [ i ] [ 2 ]-1)): od: RETURN(ans) end:
  • Mathematica
    a[n_] := Total[EulerPhi[Select[Divisors[n], GCD[#, n/#] == 1 &]]]; Array[a, 73] (* Jean-François Alcover, May 03 2011 *)
    f[p_, e_] := p^e - p^(e-1) + 1; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 10 2020 *)
  • PARI
    a(n) = sumdiv(n, d, if(gcd(n/d, d)==1, eulerphi(d))); \\ Charles R Greathouse IV, Feb 19 2013, corrected by Antti Karttunen, Sep 03 2018
    
  • PARI
    a(n)=my(f=factor(n));prod(i=1,#f[,1],f[i,1]^f[i,2]-f[i,1]^(f[i,2]-1)+1) \\ Charles R Greathouse IV, Feb 19 2013

Formula

If n = product p_i^e_i, a(n) = product (1+p_i^e_i-p_i^(e_i-1)). - Vladeta Jovovic, Apr 19 2001
Dirichlet g.f.: zeta(s)*zeta(s-1)*product_{primes p} (1+p^(-2s)-p^(1-2s)-p^(-s)). - R. J. Mathar, Oct 24 2011
Dirichlet convolution square of A318661(n)/A318662(n). - Antti Karttunen, Sep 03 2018
Sum_{k=1..n} a(k) ~ c * Pi^2 * n^2 / 12, where c = Product_{primes p} (1 - 1/p^2 - 1/p^3 + 1/p^4) = A330523 = 0.535896... - Vaclav Kotesovec, Dec 17 2019

A006511 Largest inverse of totient function (A000010): a(n) is the largest x such that phi(x) = m, where m = A002202(n) is the n-th number in the range of phi.

Original entry on oeis.org

2, 6, 12, 18, 30, 22, 42, 60, 54, 66, 46, 90, 58, 62, 120, 126, 150, 98, 138, 94, 210, 106, 162, 174, 118, 198, 240, 134, 142, 270, 158, 330, 166, 294, 276, 282, 420, 250, 206, 318, 214, 378, 242, 348, 354, 462, 254, 510, 262, 414, 274, 278, 426, 630, 298, 302
Offset: 1

Views

Author

Keywords

Comments

Always even, as phi(2n) = phi(n) when n is odd. - Alain Jacques (thegentleway(AT)bigpond.com), Jun 15 2006

References

  • J. W. L. Glaisher, Number-Divisor Tables. British Assoc. Math. Tables, Vol. 8, Camb. Univ. Press, 1940, p. 64.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

For records see A036913, A132154, A036912.

Programs

  • Mathematica
    phiinv[n_, pl_] := Module[{i, p, e, pe, val}, If[pl=={}, Return[If[n==1, {1}, {}]]]; val={}; p=Last[pl]; For[e=0; pe=1, e==0||Mod[n, (p-1)pe/p]==0, e++; pe*=p, val=Join[val, pe*phiinv[If[e==0, n, n*p/pe/(p-1)], Drop[pl, -1]]]]; Sort[val]]; phiinv[n_] := phiinv[n, Select[1+Divisors[n], PrimeQ]]; Last/@Select[phiinv/@Range[1, 200], #!={}&] (* phiinv[n, pl] = list of x with phi(x)=n and all prime divisors of x in list pl. phiinv[n] = list of x with phi(x)=n *)
  • PARI
    g(n) = if(n%2, 2*(n==1), forstep(k = floor(exp(Euler)*n*log(log(n^2))+2.5*n/log(log(n^2))), n, -1, if(eulerphi(k)==n, return(k)); if(k==n, return(0)))); \\ A057635
    lista(nn) = for(m = 1, nn, if(istotient(m), print1(g(m), ", "))); \\ Jinyuan Wang, Aug 29 2019
    
  • PARI
    lista(nmax) = my(s); for(n = 1, nmax, s = invphiMax(n); if(s > 0, print1(s, ", "))); \\ Amiram Eldar, Nov 14 2024, using Max Alekseyev's invphi.gp
  • Perl
    use ntheory ":all"; my $k=1; for my $i (1..100) { my @v; do{@v=inverse_totient($k++)} until @v; print "$i $v[-1]\n"; } # Dana Jacobsen, Mar 04 2019
    

Formula

a(n) = A057635(A002202(n)). - T. D. Noe

A135303 a(n) = phi(2*n+1)/(2*A003558(n)), where phi = A000010.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 3, 2, 1, 1, 1, 2, 3, 1, 1, 1, 2, 1, 1, 2, 1, 1, 3, 4, 1, 1, 1, 4, 1, 1, 1, 1, 1, 4, 1, 4, 3, 3, 1, 2, 2, 1, 1, 2, 1, 3, 1, 4, 1, 3, 2, 1, 2, 1, 9, 6, 1, 3, 1, 2, 1, 1, 1, 4, 1, 1, 5, 2, 3, 3, 1, 2, 1, 2, 1, 1, 6, 1, 1, 2
Offset: 1

Views

Author

N. J. A. Sloane, Dec 05 2007

Keywords

Comments

The Froemke-Grossman 1988 reference is the earliest I have seen. a(n) is the charm bracelet function b(2,2*n+1) in their notation. - N. J. A. Sloane, Feb 28 2023
This sequence is called the "coach numbers" ("c(2*n+1)"), and was studied by J. Pedersen, Byron Walden, Victor Quintanar-Zilinskas and Linda Velarde of Santa Clara University. Coach Theorem: Let b = 2*n+1 > 1 and let phi(b) be the Euler totient function. Let Sigma(b) be the complete symbol of b, let c be the number of coaches in Sigma(b), and let k = Sum_{i=1..r} k(i). Then phi(b) = 2 * c * k [Hilton & Pedersen, p. 262]. The complete symbols for b = 17 and 43 are shown in the examples. - Gary W. Adamson, Aug 15 2012
Conjecture relating to primes with more than one coach: The combined set of integers in the top rows of all coaches of these primes is composed of a permutation of the first q odd integers, where prime p = (4q-1) or (4q+1), (q > 0). Example: As shown for 17, this prime has two coaches with the top rows [1] and [3, 7, 5]. 43 has three coaches with q = 11. The top rows are [1, 21, 11], [3, 5, 19], [7, 9, 17, 13, 15]. The comment of Sep 08 2012 in A216371 applies to primes with one coach, in which case "all coaches" is reduced to one and the set of q odd integers is in the top row of the coach. - Gary W. Adamson, Sep 10 2012
Conjecture [Carl Schick]: If 2*n+1 is prime, then these are the number of distinct cycles of f(k) = |(2*n+1) - 2*k| beginning at an odd number 0 < k < 2*n. - Jonathan Skowera, Aug 03 2013 [See also the Brändli and Beyne link, eq. (2). - Wolfdieter Lang, Feb 08 2020]
From Gary W. Adamson, Oct 04 2019: (Start)
Conjecture of Aug 03 2013 proved by Jean Pedersen. By way of example, take A003558(5) = 11, such that
2^5 == -1 (mod 11). Then Pedersen on p. 98 has:
11 - 1 = 2^1 * 5 (pick "1", odd, the putative seed number)
11 - 5 = 2^1 * 3 (then subtract 3 in the next row)
11 - 3 = 2^3 * 1 (cycle ends). Then Pedersen constructs the "coach" (p. 98) for N= 11: [1, 5, 3]
.......... [1, 1, 3]. The top row represents the angles on the tape used to construct an 11-gon at the operative crease lines beginning with Pi/11. (extract the (1,5,3) column). Then extract the exponents of 2: (1,1,3); which are the bottom row terms. The final result is that at successive creases on the tape are at angles of j*Pi/11, j = (1,5,3); alternatively at the top of the tape, then the bottom. The code U(1), D(1), U(3) is understood to be those numbers of bisections at each vertex. The total numbers of bisections = 5 = (1 + 1 + 3), shown to be the entry for N = 11 in A003558. (End)

Examples

			Refer to A003558 for the J. Pedersen definition of a Coach. a(8) for b = 17 = 2 since 17 has two possible Coaches:
17: [1] and [3, 7, 5]
    [4]     [1, 1, 2];
where sum of the bottom row terms = k = 4 = A003558(8). For b = 43, a(21) = 3 since there are three possible coaches for 43:
43: [1, 21, 11]  [3, 5, 19]  [7, 9, 17, 13, 15]
    [1,  1,  5], [3, 1,  3], [2, 1,  1,  1,  2],
where k = sum of terms in bottom rows of all possible coaches = 7 = A003558(21). For the coach with a "1" in the top row, the numbers of terms in the rows ("j" in A003558), = A179480(22) = 3. Note that the parity of numbers of terms in the bottom coach rows is the same.
From _Gary W. Adamson_, Aug 24 2019: (Start)
An alternative to the coach method of Pedersen and Hilton involves the doubling sequence, mod n; (43 in this case). The top row begins (1, 2, 4, 8, 16, ...) but the next number is 11, not 32. 32 == -11 (mod 43). We pick the least (in absolute value) of the two candidates (11 and 32): 11. The top row ends when the rightmost term is (n-1)/2 = 21. In subsequent rows the leftmost term is the least odd number not previously used, in this case 3. Continue with the doubling sequence and stop when the next row produces a term already used.
"20" ends row 2 since (2 * 20) = 40 == -3 (mod 43). 3 has been used so that row ends and our next row begins with the next unused odd term, a 7. That row ends with 18 since 2 * 18 = 36 == -7 (mod 43).
The entire set is complete when every term (1 through (n-1)/2) is present without duplication. In this method, k is likewise 7 but is represented by the numbers of terms in the top row. Pedersen's [1, 21, 11] appears as the only odd terms of the top row. [3,  5, 19] appears as the odd terms of the middle row, and [7, 9, 17, 13, 15] are the only odd terms of the bottom row. The three completed rows are:
  [1,  2,  4,  8, 16, 11, 21;
   3,  6, 12, 19,  5, 10, 20;
   7, 14, 15, 13, 17,  9, 18]
  It appears that the numbers of rows is equal to Pedersen's
  number of coaches. Another example is the complete system of coaches shown on p. 261 of (Hilton and Pedersen):
  31: [1, 15], [3, 7], [5, 13, 9, 11]
      [1,  4], [2, 3], [1,  1, 1,  2]
  The alternative system, called an r-t table in A065941, is
     [1,  2,  4, 8, 15;
      3,  6, 12, 7, 14;
      5, 10, 11, 9, 13]
  The odd terms of the top row (1, 15) appear in the leftmost coach. The odd terms (3, 7) appear in the middle coach, and (5, 11, 9, 13) are shown in the rightmost coach. (End)
Pedersen's coaches can be derived from the alternative system, doubling (mod N) since her coaches are simply another version: (repeated bisections (mod N)). First write out the doubling terms (mod N). Say N = 23: 1, 2, 4, 8, 7, 9, 5, 10, 3, 6, 11, representing the trajectory of terms 2*cos(j*Pi/23), using (x^2-2); j = 1, 2, 4, .... Begin with ("1"), then jump to the next odd term (11), then to each odd term in succession going left, getting: 23: (1, 11, 3, 5, 9, 7); the top row in Pedersen's coach. ....(1, 2, 2, 1, 1, 4) is the bottom row for 23 as shown on p. 105. Those terms are the numbers of term jumps in the previous operation; for example (1 to 11) = 1, (11 to 3) = 2, (3 to 5) = 2; and so on.  Note that the number of terms in the doubling trajectory (11) matches the sum of terms in the bottom row of the coach, satisfying 2^11 == 1 (mod 23). - _Gary W. Adamson_, Oct 23 2019
		

References

  • Froemke, Jon, and Jerrold W. Grossman. "An algebraic approach to some number-theoretic problems arising from paper-folding regular polygons." The American mathematical monthly 95.4 (1988): 289-307. See Appendix.
  • Peter Hilton & Jean Pedersen, A Mathematical Tapestry: Demonstrating the Beautiful Unity of Mathematics; Cambridge University Press, 2010, pages 260-264.
  • Carl Schick, Trigonometrie und unterhaltsame Zahlentheorie, Bokos Druck, Zürich, 2003 (ISBN 3-9522917-0-6). Tables 3.1 to 3.10, for odd p = 3..113, pp. 158-166.

Crossrefs

Cf. A216371 (odd primes with one coach).

Programs

  • Maple
    A135303 := proc(n)
        numtheory[phi](2*n+1)/2/A003558(n) ;
    end proc:
    seq(A135303(n),n=1..40) ; # R. J. Mathar, Dec 01 2014
  • Mathematica
    Array[EulerPhi[#2]/(2 If[#2 > 1 && GCD[#1, #2] == 1, Min[MultiplicativeOrder[#1, #2, {-1, 1}]], 0]) & @@ {2, 2 # + 1} &, 105] (* Michael De Vlieger, Oct 25 2019 *)
  • PARI
    isok8(m, n) = my(md = Mod(2, 2*n+1)^m); (md==1) || (md==-1);
    A003558(n) = my(m=1); while(!isok8(m, n) , m++); m;
    a(n) = eulerphi(2*n+1)/(2*A003558(n)); \\ Michel Marcus, Jun 11 2020

Formula

a(n) = "c", a Coach number; = A000010(n)/(2*A003558(n-1)/2)); or phi(2*n+1) = 2 * c * k, with c = Coach numbers, k = A003558.

Extensions

Title changed by Wolfdieter Lang and M. F. Hasler, Feb 20 2020

A307868 Decimal expansion of the asymptotic mean of phi(k)/psi(k), where phi(k) is Euler totient function (A000010) and psi(k) is Dedekind psi function (A001615).

Original entry on oeis.org

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

Views

Author

Amiram Eldar, May 02 2019

Keywords

Comments

Also, the asymptotic mean of A162511. - Amiram Eldar, Sep 18 2022

Examples

			0.47168061361299786807523563308048208742592638200698...
		

Crossrefs

Programs

  • Mathematica
    $MaxExtraPrecision = 1000; m = 1000; c = LinearRecurrence[{-2, 1, 2}, {0, -4, 6}, m]; RealDigits[(2/3) * Exp[NSum[Indexed[c, n]*(PrimeZetaP[n] - 1/2^n)/n, {n, 2, m}, NSumTerms -> m, WorkingPrecision -> m]], 10, 100][[1]]
  • PARI
    prodeulerrat(1 - 2/(p*(p+1))) \\ Vaclav Kotesovec, Sep 19 2020

Formula

Equals lim_{m->oo} (1/m)*Sum_{k=1..m} phi(k)/psi(k).
Equals Product_{p prime} (1 - 2/(p * (p+1))).
Equals A065472 / zeta(2). - Amiram Eldar, Sep 18 2022

Extensions

More digits from Vaclav Kotesovec, Sep 19 2020
Previous Showing 11-20 of 4244 results. Next