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.

A033949 Positive integers that do not have a primitive root.

Original entry on oeis.org

8, 12, 15, 16, 20, 21, 24, 28, 30, 32, 33, 35, 36, 39, 40, 42, 44, 45, 48, 51, 52, 55, 56, 57, 60, 63, 64, 65, 66, 68, 69, 70, 72, 75, 76, 77, 78, 80, 84, 85, 87, 88, 90, 91, 92, 93, 95, 96, 99, 100, 102, 104, 105, 108, 110, 111, 112, 114, 115, 116, 117, 119, 120, 123
Offset: 1

Views

Author

Calculated by Jud McCranie

Keywords

Comments

Numbers k such that the cyclotomic polynomial Phi(k,x) is reducible over Zp for all primes p. Harrison shows that this is equivalent to k > 2 and the discriminant of Phi(k,x), A004124(k), being a square. - T. D. Noe, Nov 06 2007
The multiplicative group modulo k is non-cyclic; the complement A033948. - Wolfdieter Lang, Mar 14 2012. See A281854 for the groups. - Wolfdieter Lang, Feb 04 2017
Numbers k with the property that there exists a positive integer m with 1 < m < k-1 and m^2 == 1 (mod k). - Reinhard Muehlfeld, May 27 2014
Also, numbers k for which A000010(k) > A002322(k), or equivalently A034380(k) > 1. - Ivan Neretin, Mar 28 2015
Numbers k of the form a + b + 2*sqrt(a*b + 1) for positive integers a,b such that a*b + 1 is a square. Proof: If 1 < m < k - 1 and m^2 == 1 (mod k), take a = (m^2 - 1)/k and b = ((k - m)^2 - 1)/k. Conversely, if k = a + b + 2*sqrt(a*b + 1), take m = a + sqrt(a*b + 1). - Tor Gunston, Apr 24 2021
Seems to be A050275 without the duplicates. - Charles R Greathouse IV, Feb 09 2025

References

  • I. Niven and H. S. Zuckerman, An Introduction to the Theory of Numbers, 4th edition, page 62, Theorem 2.25.

Crossrefs

Cf. A000010, A002322, A033948 (complement), A193305 (composites with primitive root).
Column k=1 of A277915, A281854.

Programs

  • Haskell
    a033949 n = a033949_list !! (n-1)
    a033949_list = filter
                   (\x -> any ((== 1) . (`mod` x) . (^ 2)) [2 .. x-2]) [1..]
    -- Reinhard Zumkeller, Dec 10 2014
    
  • Maple
    m := proc(n) local k, r; r := 1; if n = 2 then return false fi;
    for k from 1 to n do if igcd(n,k) = 1 then r := modp(r*k,n) fi od; r end:
    select(n -> m(n) = 1, [$1..123]); # Peter Luschny, May 25 2017
  • Mathematica
    Select[Range[2,130],!IntegerQ[PrimitiveRoot[#]]&] (* Harvey P. Dale, Oct 25 2011 *)
    a[n_] := Module[{j, l = {}}, While[Length[l] CarmichaelLambda[j], AppendTo[l, j]; Break[]]]]; l[[n]]]; Array[a, 100] (* Jean-François Alcover, May 29 2018, after Alois P. Heinz's Maple code for A277915 *)
  • PARI
    is(n)=n>7 && (!isprimepower(if(n%2,n,n/2)) || n>>valuation(n,2)==1) \\ Charles R Greathouse IV, Oct 08 2016
    
  • Python
    from itertools import count, islice
    from sympy.ntheory import sqrt_mod_iter
    def A033949_gen(): # generator of terms
        return filter(lambda n:max(filter(lambda k:k 1,count(3))
    A033949_list = list(islice(A033949_gen(),30)) # Chai Wah Wu, Oct 26 2022
    
  • Python
    from sympy import primepi, integer_nthroot
    def A033949(n):
        def f(x): return int(n+1+(x>=2)+(x>=4)+sum(primepi(integer_nthroot(x,k)[0])-1 for k in range(1,x.bit_length()))+sum(primepi(integer_nthroot(x>>1,k)[0])-1 for k in range(1,x.bit_length()-1)))
        m, k = n, f(n)
        while m != k: m, k = k, f(k)
        return m # Chai Wah Wu, Feb 24 2025
  • Sage
    [n for n in range(1,100) if not Integers(n).multiplicative_group_is_cyclic()]
    # Ralf Stephan, Mar 30 2014
    

Formula

Positive integers except 1, 2, 4 and numbers of the form p^i and 2p^i, where p is an odd prime and i >= 1.

A066498 Numbers k such that 3 divides phi(k).

Original entry on oeis.org

7, 9, 13, 14, 18, 19, 21, 26, 27, 28, 31, 35, 36, 37, 38, 39, 42, 43, 45, 49, 52, 54, 56, 57, 61, 62, 63, 65, 67, 70, 72, 73, 74, 76, 77, 78, 79, 81, 84, 86, 90, 91, 93, 95, 97, 98, 99, 103, 104, 105, 108, 109, 111, 112, 114, 117, 119, 122, 124, 126, 127, 129, 130, 133
Offset: 1

Views

Author

Benoit Cloitre, Jan 04 2002

Keywords

Comments

Numbers k such that x^3 == 1 (mod k) has solutions 1 < x < k.
Terms are multiple of 9 or of a prime of the form 6k+1.
If k is a term of this sequence, then G = is a non-abelian group of order 3k, where 1 < r < n and r^3 == 1 (mod k). For example, G can be the subgroup of GL(2, Z_k) generated by x = {{1, 1}, {0, 1}} and y = {{r, 0}, {0, 1}}. - Jianing Song, Sep 17 2019
The asymptotic density of this sequence is 1 (Dressler, 1975). - Amiram Eldar, Mar 21 2021

Examples

			If n < 7 then x^3 = 1 (mod n) has no solution 1 < x < n, but {2,4} are solutions to x^3 == 1 (mod 7), hence a(1) = 7.
		

Crossrefs

Complement of A088232.
A007645 gives the primes congruent to 1 mod 3.
Column k=2 of A277915.

Programs

  • Mathematica
    Select[Range[150], Divisible[EulerPhi[#], 3]&] (* Harvey P. Dale, Jan 12 2011 *)
  • PARI
    isok(k)={ eulerphi(k)%3 == 0 } \\ Harry J. Smith, Feb 18 2010

Extensions

Simpler definition from Yuval Dekel (dekelyuval(AT)hotmail.com), Oct 25 2003
Corrected and extended by Ray Chandler, Nov 05 2003

A066502 Numbers k such that 7 divides phi(k).

Original entry on oeis.org

29, 43, 49, 58, 71, 86, 87, 98, 113, 116, 127, 129, 142, 145, 147, 172, 174, 196, 197, 203, 211, 213, 215, 226, 232, 239, 245, 254, 258, 261, 281, 284, 290, 294, 301, 319, 337, 339, 343, 344, 348, 355, 377, 379, 381, 387, 392, 394, 406, 421, 422, 426, 430
Offset: 1

Views

Author

Benoit Cloitre, Jan 04 2002

Keywords

Comments

Related to the equation x^7 == 1 (mod k): sequence gives values of k such there are solutions 1 < x < k of x^7 == 1 (mod k).
If k is a term of this sequence, then G = is a non-abelian group of order 7k, where 1 < r < n and r^7 == 1 (mod k). For example, G can be the subgroup of GL(2, Z_k) generated by x = {{1, 1}, {0, 1}} and y = {{r, 0}, {0, 1}}. - Jianing Song, Sep 17 2019
The asymptotic density of this sequence is 1 (Dressler, 1975). - Amiram Eldar, May 23 2022

Examples

			x^7 == 1 (mod k) has solutions 1 < x < k for k = 29, 43, 49, ...
		

Crossrefs

Programs

  • Mathematica
    Select[Range[500],Divisible[EulerPhi[#],7]&] (* Harvey P. Dale, Apr 12 2012 *)
  • PARI
    isok(k) = { eulerphi(k)%7 == 0 } \\ Harry J. Smith, Feb 18 2010

Formula

a(n) are the numbers generated by 7^2 = 49 and all primes congruent to 1 mod 7 (A045465). Hence sequence gives all k such that k == 0 (mod A045465(n)) for some n > 1 or k == 0 (mod 49).

Extensions

Simpler definition from Yuval Dekel (dekelyuval(AT)hotmail.com), Oct 25 2003

A066500 Numbers k such that 5 divides phi(k).

Original entry on oeis.org

11, 22, 25, 31, 33, 41, 44, 50, 55, 61, 62, 66, 71, 75, 77, 82, 88, 93, 99, 100, 101, 110, 121, 122, 123, 124, 125, 131, 132, 142, 143, 150, 151, 154, 155, 164, 165, 175, 176, 181, 183, 186, 187, 191, 198, 200, 202, 205, 209, 211, 213, 217, 220, 225, 231, 241
Offset: 1

Views

Author

Benoit Cloitre, Jan 04 2002

Keywords

Comments

Related to the equation x^5 == 1 (mod k): sequence gives values of k such there are solutions 1 < x < k of x^5 == 1 (mod k).
If k is a term of this sequence, then G = is a non-abelian group of order 5k, where 1 < r < n and r^5 == 1 (mod k). For example, G can be the subgroup of GL(2, Z_k) generated by x = {{1, 1}, {0, 1}} and y = {{r, 0}, {0, 1}}. - Jianing Song, Sep 17 2019
The asymptotic density of this sequence is 1 (Dressler, 1975). - Amiram Eldar, May 23 2022

Examples

			x^5 == 1 (mod 11) has solutions 1 < x < 11, namely {3,4,5,9}.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[250], Divisible[EulerPhi[#], 5] &] (* Amiram Eldar, May 23 2022 *)
  • PARI
    isok(k) = { eulerphi(k)%5 == 0 } \\ Harry J. Smith, Feb 18 2010

Formula

a(n) are the numbers generated by 5^2 = 25 and all primes congruent to 1 mod 5 (A045453). Hence sequence gives all k such that k == 0 (mod A045453(n)) for some n > 1 or k == 0 (mod 25).

Extensions

Simpler definition from Yuval Dekel (dekelyuval(AT)hotmail.com), Oct 25 2003
Extended by Ray Chandler, Nov 06 2003

A305828 The n-th number m such that a nontrivial prime(n)-th root of unity modulo m exists.

Original entry on oeis.org

8, 9, 25, 58, 89, 158, 307, 571, 282, 466, 1303, 1043, 913, 1384, 1977, 1486, 3308, 2936, 3497, 3983, 2930, 4426, 2338, 3401, 4668, 6070, 6186, 6430, 10910, 4994, 9671, 4997, 9053, 10583, 13123, 10874, 22611, 14366, 30066, 8328, 10031, 15206, 11873, 17006
Offset: 1

Views

Author

Alois P. Heinz, Jun 10 2018

Keywords

Comments

For n > 1 the n-th number m such that phi(m) is divisible by the n-th prime.

Crossrefs

Main diagonal of A277915.

Formula

a(n) = A277915(n,n).
Showing 1-5 of 5 results.