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

A002174 Values taken by reduced totient function psi(n).

Original entry on oeis.org

1, 2, 4, 6, 8, 10, 12, 16, 18, 20, 22, 24, 28, 30, 32, 36, 40, 42, 44, 46, 48, 52, 54, 56, 58, 60, 64, 66, 70, 72, 78, 80, 82, 84, 88, 90, 92, 96, 100, 102, 104, 106, 108, 110, 112, 116, 120, 126, 128, 130, 132, 136, 138, 140, 144, 148, 150, 156, 160, 162, 164, 166, 168
Offset: 1

Views

Author

Keywords

Comments

If p is a Sophie Germain prime (A005384), then 2p is here. - T. D. Noe, Aug 13 2008
Terms of A002322, sorted and multiple values taken just once. - Vladimir Joseph Stephan Orlovsky, Jul 21 2009
a(2445343) = 10^7, suggesting that Luca & Pomerance's lower bound may be closer to the truth than the upper bound. The fit exponent log a(n)/log n - 1 = 0.0957... in this case. - Charles R Greathouse IV, Jul 02 2017

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    lst={}; Do[AppendTo[lst, CarmichaelLambda[n]], {n, 6*7!}]; lst; Take[Union[lst], 123] (* Vladimir Joseph Stephan Orlovsky, Jul 21 2009 *)
    (* warning: there seems to be no guarantee that no terms near the end are omitted! - Joerg Arndt, Dec 23 2014 *)
    TakeWhile[Union@ Table[CarmichaelLambda@ n, {n, 10^6}], # <= 168 &] (* Michael De Vlieger, Mar 19 2016 *)
  • PARI
    list(lim)=my(v=List([1]),u,t); forprime(p=3,lim\3+1, u=List(); listput(u,p-1); while((t=u[#u]*p)<=lim, listput(u,t)); for(j=1,#v, for(i=1,#u, t=lcm(u[i],v[j]); if(t<=lim && t!=v[j], listput(v,t)))); v=List(Set(v))); forprime(p=lim\3+2,lim+1, listput(v,p-1)); v=List(Set(v)); for(i=1,#v, t=2*v[i]; if(t>lim, break); listput(v,t); while((t*=2)<=lim, listput(v,t))); Set(v) \\ Charles R Greathouse IV, Jun 23 2017
    
  • PARI
    is(n)=if(n%2, return(n==1)); my(f=factor(n),pe); for(i=1,#f~, if(n%(f[i,1]-1)==0, next); pe=f[i,1]^f[i,2]; forstep(q=2*pe+1,n+1,2*pe, if(n%(q-1)==0 && isprime(q), next(2))); return(0)); 1 \\ Charles R Greathouse IV, Jun 25 2017

Formula

n (log n)^0.086 << a(n) << n (log n)^0.36 where << is the Vinogradov symbol, see Luca & Pomerance. - Charles R Greathouse IV, Dec 28 2013

Extensions

More terms from T. D. Noe, Aug 13 2008

A079612 Largest number m such that a^n == 1 (mod m) whenever a is coprime to m.

Original entry on oeis.org

2, 24, 2, 240, 2, 504, 2, 480, 2, 264, 2, 65520, 2, 24, 2, 16320, 2, 28728, 2, 13200, 2, 552, 2, 131040, 2, 24, 2, 6960, 2, 171864, 2, 32640, 2, 24, 2, 138181680, 2, 24, 2, 1082400, 2, 151704, 2, 5520, 2, 1128, 2, 4455360, 2, 264, 2, 12720, 2, 86184, 2, 13920
Offset: 1

Views

Author

N. J. A. Sloane, Jan 29 2003

Keywords

Comments

a(m) divides the Jordan function J_m(n) for all n except when n is a prime dividing a(m) or m=2, n=4; it is the largest number dividing all but finitely many values of J_m(n). For m > 0, a(m) also divides Sum_{k=1}^n J_m(k) for n >= the largest exceptional value. - Franklin T. Adams-Watters, Dec 10 2005
The numbers m with this property are the divisors of a(n) that are not divisors of a(r) for r

References

  • R. C. Vaughan and T. D. Wooley, Waring's problem: a survey, pp. 285-324 of Surveys in Number Theory (Urbana, May 21, 2000), ed. M. A. Bennett et al., Peters, 2003. (The function K(n), see p. 303.)

Crossrefs

Cf. A006863 (bisection except for initial term); A059379 (Jordan function).

Programs

  • PARI
    a(n) = {if (n%2, 2, res = 1; forprime(p=2, n+1, if (!(n % (p-1)), t = valuation(n, p); if (p==2, if (t, res *= p^(t+2)), res *= p^(t+1)););); res;);} \\ Michel Marcus, May 12 2018

Formula

a(n) = 2 for n odd; for n even, a(n) = product of 2^(t+2) (where 2^t exactly divides n) and p^(t+1) (where p runs through all odd primes such that p-1 divides n and p^t exactly divides n).
From Antti Karttunen, Dec 19 2018: (Start)
a(n) = A185633(n)*(2-A000035(n)).
It also seems that for n > 1, a(n) = 2*A075180(n-1). (End)
We have 2*A075180(2n-1) = A006863(n) by definition, and A006863(n) = a(2n) by the comments in A006863. Hence a(n) = 2*A075180(n-1) for all even n. For all odd n > 1, we have a(n) = 2, which is also equal to 2*A075180(n-1). So the formula above is true. - Jianing Song, Apr 05 2021

Extensions

Edited by Franklin T. Adams-Watters, Dec 10 2005
Definition corrected by T. D. Noe, Aug 13 2008
Rather arbitrary term a(0) removed by Max Alekseyev, May 27 2010

A143407 Largest number k such that the reduced totient function psi(k) = A002174(n).

Original entry on oeis.org

2, 24, 240, 504, 480, 264, 65520, 16320, 28728, 13200, 552, 131040, 6960, 171864, 32640, 138181680, 1082400, 151704, 5520, 1128, 4455360, 12720, 86184, 13920, 1416, 6814407600, 65280, 776664, 18744, 20174525280, 39816, 36801600, 1992
Offset: 1

Author

T. D. Noe, Aug 13 2008

Keywords

Comments

For each of the values in A002174, there are only a finite number of numbers k such that psi(k)=A002174(n). This sequence gives the largest such k. Sequence A002396 gives the least k. The number of such k is given in A143408. When A002174(n) is twice a Sophie Germain prime, then a(n) is particularly small.

Crossrefs

Cf. A002322 (reduced totient function), A002174, A002396, A143408.

Formula

a(n) = A006863(A002174(n)/2) for n>1.

A141162 Smallest k such that lambda(k) = n, or 0 if there is no such k.

Original entry on oeis.org

1, 3, 0, 5, 0, 7, 0, 32, 0, 11, 0, 13, 0, 0, 0, 17, 0, 19, 0, 25, 0, 23, 0, 224, 0, 0, 0, 29, 0, 31, 0, 128, 0, 0, 0, 37, 0, 0, 0, 41, 0, 43, 0, 115, 0, 47, 0, 119, 0, 0, 0, 53, 0, 81, 0, 928, 0, 59, 0, 61, 0, 0, 0, 256, 0, 67, 0, 0, 0, 71, 0, 73, 0, 0, 0, 0, 0, 79, 0, 187, 0, 83, 0, 203, 0, 0, 0, 89, 0, 209, 0, 235, 0, 0, 0, 97, 0
Offset: 1

Author

Michel Lagneau, Mar 17 2011

Keywords

Comments

Sequence A002174 gives the n such that a(n) > 0. Removing the zeros from this sequence produces A002396. Note that some n appear only for large k. For example, 728 does not appear until k=49184. See A143407 for the largest k that produces a particular value of the lambda function. See A143408 for the number of times each value occurs. - T. D. Noe, Mar 17 2011

Examples

			a(8) = 32 because lambda(32) = 8.
		

Crossrefs

Cf. A002174, A002322 (Carmichael lambda function), A002396, A143407, A143408.

Programs

  • Maple
    with(numtheory):for k from 1 to 100 do:id:=0:for n from 1 to 1000 while(id=0)
      do: if lambda(n) = k then id:=1:printf(`%d, `,n):else fi:od:if id=0 then printf(`%d, `,0):else fi:od:
  • Mathematica
    nn = 100; t = Table[0, {nn}]; Do[c = CarmichaelLambda[k]; If[c <= nn && t[[c]] == 0, t[[c]] = k], {k, 1000}]; t

Formula

a(A002174(n)) = A002396(n).

A143417 Irregular triangle in which row n gives the k such that the reduced totient function psi(k) = A002174(n).

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 12, 24, 5, 10, 15, 16, 20, 30, 40, 48, 60, 80, 120, 240, 7, 9, 14, 18, 21, 28, 36, 42, 56, 63, 72, 84, 126, 168, 252, 504, 32, 96, 160, 480, 11, 22, 33, 44, 66, 88, 132, 264, 13, 26, 35, 39, 45, 52, 65, 70, 78, 90, 91, 104, 105, 112, 117, 130, 140, 144, 156
Offset: 1

Author

T. D. Noe, Aug 13 2008

Keywords

Comments

Row n consists of the divisors of A143407(n) that are not divisors of A143407(r) for rA143408(n).

Examples

			1,2; 3,4,6,8,12,24; 5,10,15,16,20,30,40,48,60,80,120,240
		

Crossrefs

Cf. A002322 (reduced totient function).
Showing 1-5 of 5 results.