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

A252594 Records in A072994.

Original entry on oeis.org

1, 2, 8, 16, 27, 32, 54, 64, 100, 128, 200, 243, 256, 400, 486, 500, 512, 800, 972, 1000, 1024, 1600, 1944, 2000, 2048, 3200, 3888, 4000, 4096, 4624, 6400, 7776, 8000, 8192, 9248, 12100, 12500, 12800, 15552, 16000, 16384, 18496, 24200, 25000, 25600, 31104, 32000, 32768, 36992, 48400, 50000
Offset: 1

Views

Author

Robert G. Wilson v, Dec 18 2014

Keywords

Comments

Consider the function f(n) = the number of positive integers k < n such that k^n (mod n) == 1. This sequence lists the values of n at which f(n) reaches a new maximum.
All powers of two are present except its square. f(2^n) (with exception noted) = 2^(n-1) = 2^n/2.
All powers of two multiplied by 100, 1000 and 100000, but not 10000, are also present.
Terms other than the above are 27, 54, 243, 486, 500, 972, 1944, 3888, 4624, 7776, 9248, 12100, 12500, 15552, 18496, 24200, 25000, 31104, 36992, 48400, 50000, ..., .
Conjecture: f(x)/x -> 5/12.

Crossrefs

Cf. A072994.

Programs

  • Mathematica
    f[n_] := (d = If[ OddQ@ n, 1, 2];   d*Length@ Select[ Range[ n/d], PowerMod[#, n, n] == 1 &]); f[1] = 1; k = 1; mx = 0; lst = {}; While[k < 10000001, a = f@ k; If[a > mx, mx = a; AppendTo[lst, k]; Print[{a, k}]]; k++]; lst

A251642 Records values in A072994.

Original entry on oeis.org

1, 2, 4, 8, 9, 16, 18, 32, 40, 64, 80, 81, 128, 160, 162, 200, 256, 320, 324, 400, 512, 640, 648, 800, 1024, 1280, 1296, 1600, 2048, 2176, 2560, 2592, 3200, 4096, 4352, 4400, 5000, 5120, 5184, 6400, 8192, 8704, 8800, 10000, 10240, 10368, 12800, 16384, 17408, 17600, 20000, 20480, 22000, 25000, 25600, 32768, 34816, 36992, 40000, 40960, 43520, 44000, 50000, 51200, 65536, 69632, 73984, 80000, 81920, 87040, 88000, 100000, 102400, 131072, 139264, 147968
Offset: 1

Views

Author

Robert G. Wilson v, Dec 06 2014

Keywords

Crossrefs

Cf. A072994.

A050399 Least k such that n = A009195(k) (= gcd(phi(k), k)).

Original entry on oeis.org

1, 4, 9, 8, 25, 18, 49, 16, 27, 50, 121, 36, 169, 98, 225, 32, 289, 54, 361, 100, 147, 242, 529, 72, 125, 338, 81, 196, 841, 450, 961, 64, 1089, 578, 1225, 108, 1369, 722, 507, 200, 1681, 294, 1849, 484, 675, 1058, 2209, 144, 343, 250, 2601, 676, 2809
Offset: 1

Views

Author

Christian G. Bower, Nov 15 1999

Keywords

Comments

Coincides with A072995 for many terms, but differs, e.g., in n = 20, 40, 52, ... in addition to the zeros in A072995. See also the comments in A072994. - M. F. Hasler, Feb 23 2014
a(n) <= n^2. - Robert G. Wilson v, Feb 27 2014

Programs

  • Mathematica
    t = Table[0, {10000}]; k = 1; While[k < 100000001, a = GCD[k, EulerPhi@ k]; If[a < 10001 && t[[a]] == 0, t[[a]] = k; Print[{a, k}]]; k++]; t (* Robert G. Wilson v, Feb 27 2014 *)
  • PARI
    A050399=n->for(k=1,oo,gcd(eulerphi(k),k)==n&&return(k)) \\ M. F. Hasler, Feb 23 2014

A072995 Least k > 0 such that the number of solutions to x^k == 1 (mod k) 1 <= x <= k is equal to n, or 0 if no such k exists.

Original entry on oeis.org

1, 4, 9, 8, 25, 18, 49, 16, 27, 50, 121, 36, 169, 98, 225, 32, 289, 54, 361, 110, 147, 242, 529, 72, 125, 338, 81, 196, 841, 0, 961, 64, 1089, 578, 1225, 108, 1369, 722, 507, 100, 1681, 0, 1849, 484, 675, 1058, 2209, 144, 343, 250, 2601, 1378, 2809
Offset: 1

Views

Author

Benoit Cloitre, Aug 21 2002

Keywords

Comments

A072989 lists the indices for which a(n) differs from A050399(n), e.g., in n = 20, 40, 52, ... in addition to the zeros in this sequence (n = 30, 42, 66, 70, 78, 90, ...). See also A009195 vs. A072994. [Corrected and extended by M. F. Hasler, Feb 23 2014]
The sequence seems difficult to extend, as the next term a(30) is larger than 5100. However, a(32)=64, a(64)=128 and a(128)=256 can be easily calculated. It thus appears that a(2^k)=2^(k+1), for k=1,2,3,.... Is this known to be true? - John W. Layman, Aug 05 2003 -- Answer: It's true. One could have defined the sequence so that a(1)=2: then it would be true for 2^0 also. - Don Reble, Feb 23 2014
a(30), if it exists, is greater than 400000. - Ryan Propper, Sep 10 2005
a(30) doesn't exist: If N is even, and divisible by D different odd primes, but not divisible by 2^D, then a(N) doesn't exist. - Don Reble, Feb 23 2014 [This and the preceding comment refer to the former definition lacking the clause "0 if no such k exists". - Ed.]
Conjecture: a(n)=0 iff n/2 is in A061346. - Robert G. Wilson v, Feb 23 2014
[n=420 seems to be a counterexample to the above conjecture. - M. F. Hasler, Feb 24 2014]
From Robert G. Wilson v, Mar 05 2014: (Start)
Observation:
If n = 1 then a(n) = 1 by definition;
If, but not iff, n (an even number) is a member of A238367 then a(n) = 0;
If n (an even number not in A238367) is {684, 954, ...}, then a(n) = 0;
If n (an odd number) is {273, 399, 651, 741, 777, 903, ...}, then a(n) = 0;
If p is a prime [A000040] and e is its exponent, then a(p^e) = p^(e+1);
If p is a prime then a(2p^e) = 2p^(e+1);
If p is a prime then a(n) # p since the f(p)=1.
(End)
Often A072995(n) equals A050399(n). They differ at n: 20, 30, 40, 42, 52, 60, 66, 68, 70, 78, 80, 84, 90, 100, 102, 104, 110, 114, 116, 120, 126, 130, 132, ... - Robert G. Wilson v, Dec 06 2014
When A072995(n)>0 and does not equal A050399(n): 20, 40, 52, 60, 68, 80, 84, 100, 104, 116, 120, 132, 136, 140, 148, 156, 160, 164, 168, 171, 180, 200, ... - Robert G. Wilson v, Dec 06 2014
When a(n) > 1, then 2n <= a(n) <= n^2. - Robert G. Wilson v, Dec 10 2014

Crossrefs

Cf. A072994.

Programs

  • Mathematica
    t = Table[0, {1000}]; f[n_] := (d = If[EvenQ@ n, 2, 1]; d*Length@ Select[ Range[ n/d], PowerMod[#, n, n] == 1 &]); f[1] = 1; k = 1; While[k < 520001, If[ PrimeQ@ k, k++]; a = f@ k; If[a < 1001 && t[[a]] == 0, t[[a]] = k; Print[{a, k}]]; k++]; t (* Robert G. Wilson v, Dec 12 2014 *)
  • PARI
    A072995(n)=(n%2||n%2^(omega(n)-1)==0)&&for(k=1,9e9,A072994(k)==n&&return(k)) \\ M. F. Hasler, Feb 23 2014

Formula

First occurrence of n in A072994.

Extensions

More terms from Don Reble, Feb 23 2014
Edited, at the suggestion of Don Reble, by M. F. Hasler, Feb 23 2014

A196082 Greatest residue of x^n (mod n), x=0..n-1.

Original entry on oeis.org

0, 1, 2, 1, 4, 4, 6, 1, 8, 9, 10, 9, 12, 11, 14, 1, 16, 10, 18, 16, 20, 20, 22, 16, 24, 25, 26, 25, 28, 25, 30, 1, 32, 33, 34, 28, 36, 36, 38, 25, 40, 36, 42, 37, 44, 41, 46, 33, 48, 49, 50, 48, 52, 28, 54, 49, 56, 57, 58, 45, 60, 59, 62, 1, 64, 64, 66, 64, 68
Offset: 1

Views

Author

Michel Lagneau, Sep 27 2011

Keywords

Comments

a(n) = 1 if n is of the form 2^p and a(n) = n-1 if n prime.

Examples

			a(18) = 10 because x^18 == 0, 1, 9, 10  (mod 18) => 10 is the greatest residue.
		

Crossrefs

Programs

  • Mathematica
    Table[Max[Union[PowerMod[Range[0, n - 1], n, n]]], {n, 100}]

A072989 Numbers m>0 such that the number of solutions to x^m==1 (mod m), 1<=x<=m, is not equal to gcd(m, phi(m)).

Original entry on oeis.org

20, 30, 40, 42, 52, 60, 66, 68, 70, 78, 80, 84, 90, 100, 102, 104, 110, 114, 116, 120, 126, 130, 132, 136, 138, 140, 148, 150, 154, 156, 160, 164, 168, 170, 171, 174, 180, 182, 186, 190, 198, 200, 204, 208, 210, 212, 220, 222, 228, 230, 232, 234, 238, 240
Offset: 1

Views

Author

Benoit Cloitre, Aug 21 2002

Keywords

Comments

Conjecture: limit of a(n)/n is zero.
This conjecture is certainly wrong as stated, because sequences "Numbers such that..." have lim a(n)/n >= 1 and a(n) > n for all indices following the first one for which this holds, as here: a(1) > 1. - M. F. Hasler, Feb 24 2014

Crossrefs

Programs

  • PARI
    isok(m) = sum(x=1, m, Mod(x, m)^m==1) != gcd(m, eulerphi(m)); \\ Michel Marcus, Feb 18 2021

Formula

Equals { m>0 | A009195(m) != A072994(m) }. - M. F. Hasler, Feb 23 2014
Showing 1-6 of 6 results.