Sunday, 24 April 2022

Another sequence of note

In addition to the pentatrigesimal sequences I looked at last year — based upon using all numbers and latters except O — I have studied another sequence of numbers which I will tabulate below. Of the first 104 terms which I am tabulating, one is not yet known and is shaded in grey background with white text.

n

p

2

2

3

71

4

0

5

11

6

29

7

131

8

0

9

0

10

23

11

73

12

97

13

137

14

41

15

43

16

0

17

419

18

25667

19

59

20

1487

21

156217

22

79

23

3181

24

53

25

0

26

347

27

0

28

457

29

151

30

163

31

5581

32

0

33

197

34

1493

35

313

36

0

37

251

38

401

39

349

40

751

41

83

42

1319

43

6277

44

167

45

3319

46

67

47

18013

48

383

49

0

50

6521

51

4229

52

257

53

1571

54

389

55

839

56

157

57

16963

58

2333

59

479

60

173

61

37

62

757

63

3067

64

0

65

375017

66

19973

67

367

68

2767

69

2371

70

761

71

1583

72

227

73

110603

74

191

75

739

76

439

77

15361

78

1949

79

659

80

>399989

81

0

82

7607

83

2713

84

3917

85

2111

86

113

87

121487

88

577

89

571

90

5209

91

4421

92

13001

93

4903

94

170371

95

523

96

3343

97

1693

98

2801

99

5563

100

0

101

677

102

673

103

1549

104

263

105

4783

This sequence is — in essence — the inverse of OEIS sequence A066180. The nth term of this sequence is the smallest prime for which n is the smallest base yielding a generalised repunit prime. Alternatively, the nth term is defined as the first prime number yielding n in sequence A066180.

For bases that are perfect powers, generalised repunits can be factored algebraically and the sequence has the value 0. For base 65 — until the recent discovery of the probable prime 65375017-1/64 — and base 80, no known generalised repunit prime exists or existed without a smaller base yielding a generalised repunit prime, as can be seen from the table below:

base

primes

Smaller bases where Rp is prime

65

19

2, 10, 11, 12, 14, 19, 24, 40, 45, 46, 48

29

6, 40

631

39

80

3

2, 3, 5, 6, 8, 12, 14, 15, 17, 20, 21, 24, 27, 33, 38, 41, 50, 54, 57, 59, 62

7

2, 3, 5, 6, 13, 14, 26, 31, 38, 40, 46, 56, 60, 61, 66, 68, 72, 73

At present I am not sure over what range bases 65 and 80 have been checked, although data for adjacent bases from Henri and Renauld Lifchitz suggests they have probably been checked up to around four hundred thousand without additional generalised repunit primes being discovered.

As a last word, it might be noted that, of the first three hundred primes (up to 1987), eighty-three do not appear in this sequence at all, viz:

3

5

7

13

17

19

23

31

47

61

73

89

97

101

103

107

109

127

131

137

139

149

181

211

269

271

283

317

337

353

359

409

433

449

463

487

509

521

541

569

587

593

607

619

631

653

661

701

757

769

821

857

883

907

929

971

991

1013

1021

1031

1049

1061

1069

1087

1091

1151

1181

1193

1277

1279

1297

1303

1367

1409

1423

1487

1627

1699

1721

1759

1789

1861

1907

 

For these primes, the first base yielding a generalised repunit prime is also the first such base for a smaller prime. For instance:
  1. for 3, 5, and 7, base 2 is the first such base, but 22-1 is also prime
  2. for 103, 541, 1091 and 1367, base 3 is the first base yielding a prime but (371-1)/2 is also prime
  3. for 317 and 1031, base 10 is the first base yielding a prime, but (1023-1)/9 is also prime

4 comments:

123 said...

See http://www.fermatquotient.com/PrimSerien/GenRepu.txt, p=375017 is the smallest prime p such that 65 is the smallest such base b

Also, consider those generalised repunit primes:

(35^313-1)/34
(39^349-1)/38
(51^4229-1)/50
(91^4421-1)/90
(124^599-1)/123
(135^1171-1)/134
(142^1231-1)/141

For these generalised repunit primes, p is the smallest prime such that (b^p-1)/(b-1) is prime, also b is the smallest base such that (b^p-1)/(b-1) is prime, do there exist infinite many such generalised repunit primes?

jpbenney said...

Thanks!

80 is still missing and needing further search, but thanks for a few more terms beyond numbers I could calculate reasonably consistently. Actually, although a(106) and a(111) remain unknown, the next few numbers are mostly able to be calculated from your link:

a(107) = 24251
a(108) = 2477
a(109) = 13679
a(110) = 691
a(112) = 1697
a(113) = 6563
a(114) = 569
a(115) = 241
a(116) = 2503

before consecutive unknown values occur.

I wonder what the likely maximum of a(b) actually would be, so we might have an idea about the missing values for 80, 106, 111, 117 and 118?

123 said...

Well, a(80) should be > 500000 (I have not confirmed), and after a(105), by the data in http://www.fermatquotient.com/PrimSerien/GenRepu.txt, we have:

a(106) = ?
a(107) = 24251
a(108) = 2477
a(109) = 13679
a(110) = 691
a(111) = ?
a(112) = 1697
a(113) = 6563
a(114) = 709
a(115) = 241
a(116) = 2503
a(117) = ?
a(118) = 193
a(119) = 827
a(120) = 373
a(121) = 0
a(122) = 3803
a(123) = 563
a(124) = 599
a(125) = 0
a(126) = 20947
a(127) = 5281
a(128) = 0
a(129) = ?
a(130) = ?
a(131) = ?
a(132) = ?
a(133) = 3083
a(134) = 2843
a(135) = 1171
a(136) = 293
a(137) = 1009
a(138) = ?
a(139) = ?
a(140) = ?
a(141) = ?
a(142) = 1231
a(143) = ?
a(144) = 0
a(145) = ?
a(146) = 21961
a(147) = 983
a(148) = 1201
a(149) = 3251
a(150) = 3389
a(151) = 4831
a(152) = 270217
a(153) = 5099
a(154) = 8161
a(155) = 2087
a(156) = 199
a(157) = 2791
a(158) = 4003
a(159) = 1433
a(160) = 3989

Maybe you can consider to fill the gaps.

123 said...

New terms found:

a(129) = 51829
a(130) = 38329
a(132) = 53891

All other terms up to a(160) are > 100000