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.

A130421 Numbers k such that 4^k == 2 (mod k).

Original entry on oeis.org

1, 2, 14, 1022, 20066, 80519, 107663, 485918, 1284113, 1510313, 2531678, 3677198, 3933023, 4557713, 8277458, 8893262, 21122318, 24849833, 26358638, 39852014, 42448478, 71871113, 76712318, 80646143, 98058097, 104832833, 106694033, 131492498, 144322478, 146987033, 164360606, 168204191, 175126478, 176647378, 188997463, 196705598
Offset: 1

Views

Author

Jon E. Schoenfield, May 26 2007

Keywords

Comments

Some terms above 10^15: 26435035805519327, 158975398896178078, 64044049390757098943, 1063423126446412987081943, 1091220919655042176978844783, 81074850280355100090334498663, 6317483763950169936179578094903, 5672799393875320397186007124651847, 170923900137537174138295268515194974, 195746953975871672436191077726091399305155458, 325665752547333314939363628501536024940097079718953, 953533053776414279913696071891872697927468471633033, 85791212788381063775490416118630897060666265030605503, 334519297382630382793758729321508383611586565722054114034741260213364710519401967713. - Max Alekseyev, Jun 18 2014

Crossrefs

Cf. A006935 (odd terms times 2), A130422, A347906 (odd terms), A347908 (even terms).

Programs

  • Mathematica
    Join[{1,2},Select[Range[107000000],PowerMod[4,#,#]==2&]] (* Harvey P. Dale, Jun 13 2013 *)

Extensions

Terms a(28) onward from Max Alekseyev, Jun 18 2014
b-file corrected by Max Alekseyev, Oct 09 2016

A347906 Numbers k such that 2^(2*k-1) == 1 (mod k).

Original entry on oeis.org

1, 80519, 107663, 1284113, 1510313, 3933023, 4557713, 24849833, 71871113, 80646143, 98058097, 104832833, 106694033, 146987033, 168204191, 188997463, 205428713, 332693873, 333681761, 336327863, 380284847, 533039513, 552913169, 711999113, 725943719, 805031663, 1000519033, 1069441313, 1476327353, 1610020913
Offset: 1

Views

Author

Jianing Song, Sep 18 2021

Keywords

Comments

Odd numbers k such that ord(2,k) divides 2*k-1, where ord(2,k) is the multiplicative order of 2 modulo k.
Numbers k such that 2*k is in A006935. For k > 1, k is a term if and only if 2*k is an even pseudoprime to base 2.
Odd terms in A130421. Complement of A347908 in A130421.
Terms > 1 must be composite, since for odd primes p we have 2^(2*p-1) == 2 (mod p). If k > 1 is a term, then 2*k-1 must also be composite, since ord(2,k) | (2*k-1) and ord(2,k) <= eulerphi(k) <= k-1 < 2*k-1.
If k > 1 is a term, then (2^(2*k-1) - 1)/k is composite. Proof: since 2*k-1 is composite, write 2*k-1 = u*v, u >= v > 1, then (2^(2*k-1) - 1)/k = (2^u - 1)*(2^(u*(v-1)) + ... + 2^u + 1)/k. Since k | 2^(2*k-1) - 1, there exist positive integers a,b such that a*b = k and that a | 2^u - 1 and b | 2^(u*(v-1)) + ... + 2^u + 1. Note that (2^u - 1)/a, (2^(u*(v-1)) + ... + 2^u + 1)/b >= (2^u - 1)/k >= (2^sqrt(2*k-1) - 1)/k > 1, so (2^(2*k-1) - 1)/k is the product of two integers > 1, so it is composite.
2^t - 1 is a term if and only if 2^(t+1) == 3 (mod t) (t = 1, 111481, 465793, ... in A296370).

Examples

			80519 is a term since 80519 divides 2^161037 - 1 (the multiplicative order of 2 modulo 80519 is 261, which is a divisor of 161037). Note that 2 * 80519 = 161038 = A006935(2) is the smallest even pseudoprime to base 2.
		

Crossrefs

Cf. A347907 (a similar sequence).

Programs

  • PARI
    isA347906(k) = if(k%2 && !isprime(k), Mod(2, k)^(2*k-1)==1, 0)

Formula

a(n) = A006935(n)/2.

A347907 Numbers k such that 2^(4*k-1) == 1 (mod k).

Original entry on oeis.org

1, 7, 511, 10033, 242959, 1265839, 1838599, 4138729, 4446631, 10561159, 13179319, 19926007, 21224239, 38356159, 65746249, 72161239, 82180303, 87563239, 88323689, 98352799, 124563313, 153394537, 158525689, 219011569, 248520769, 348485359, 498260329, 636381799, 638395369, 685333399, 689463889
Offset: 1

Views

Author

Jianing Song, Sep 18 2021

Keywords

Comments

Odd numbers k such that ord(2,k) divides 4*k-1, where ord(2,k) is the multiplicative order of 2 modulo k.
Numbers k such that 2*k is in A130421.
Terms > 7 must be composite, since for odd primes p we have 2^(4*p-1) == 8 (mod p). If k > 1 is a term, then 4*k-1 must also be composite, since ord(2,k) | (4*k-1) and ord(2,k) <= eulerphi(k) <= k-1 < 4*k-1.
If k > 1 is a term, then (2^(4*k-1) - 1)/k is composite. Proof: since 4*k-1 is composite, write 4*k-1 = u*v, u >= v > 1. Proof: since 4*k-1 is composite, write 4*k-1 = u*v, u >= v > 1, then (2^(4*k-1) - 1)/k = (2^u - 1)*(2^(u*(v-1)) + ... + 2^u + 1)/k. Since k | 2^(4*k-1) - 1, there exist positive integers a,b such that a*b = k and that a | 2^u - 1 and b | 2^(u*(v-1)) + ... + 2^u + 1. Note that (2^u - 1)/a, (2^(u*(v-1)) + ... + 2^u + 1)/b >= (2^u - 1)/k >= (2^sqrt(4*k-1) - 1)/k > 1, so (2^(4*k-1) - 1)/k is the product of two integers > 1, so it is composite.
2^t - 1 is a term if and only if 2^(t+2) == 5 (mod t) (t = 1, 3, 9, 871, 2043, 2119, 8769, ...).

Examples

			7 is a term since 7 divides 2^27 - 1.
		

Crossrefs

Cf. A347906 (a similar sequence).

Programs

  • Mathematica
    Join[{1},Parallelize[Select[Range[69*10^7],PowerMod[2,4#-1,#]==1&]]] (* Harvey P. Dale, Apr 16 2023 *)
  • PARI
    isA347907(k) = if(k%2 && (!isprime(k) || k==7), Mod(2, k)^(4*k-1)==1, 0)

Formula

a(n) = A347908(n)/2.
Showing 1-3 of 3 results.