login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Primes which are the fifth element of a generalized Wieferich sequence.
1

%I #16 Aug 04 2024 19:35:25

%S 359,6211,13477,19069,20431,22567,32933,39041

%N Primes which are the fifth element of a generalized Wieferich sequence.

%C A generalized Wieferich sequence is an increasing sequence of primes p[1],p[2],... such that a=p[n+1] is a Wieferich prime to base b=p[n], i.e., a^(b-1)=1 (mod b^2).

%C The present sequence is a subsequence of A179400. The numbers 728437, 1051139, 144 449603 and 5176 948723 found by M. Kamada and D. Broadhurst are also in this sequence. The last two numbers are even 7th elements of a generalized Wieferich sequence, cf. the link.

%H Kevin Acres, Mike Oakes, David Broadhurst, Makoto Kamada, <a href="/A179400/a179400.txt">1993/2011 puzzle</a>, digest of 15 messages in primenumbers Yahoo group, Jan 8 - Jan 9, 2011.

%H D. Broadhurst, <a href="http://groups.yahoo.com/group/primenumbers/message/22386">Re: 1993/2011 puzzle [and Puzzle 7]</a>, in primenumbers@yahoogroups.com, Jan 2011.

%e The smallest generalized Wieferich sequence of length 5 is (3, 11, 71, 331, 359): 3^10=1 (mod 11^2), 11^70=1 (mod 71^2), 71^330=1 (mod 331^2) and 331^358=1 (mod 359^2). Therefore, a(1)=359.

%e It happens that after a(2)=6211, there are four subsequent elements of A179400, namely 13477, 19069, 20431 and 22567, which are also member of a sequence of length 5 and thus in this sequence.

%o (PARI) fp(p)={ my(a=precprime(p)); while(Mod(a, p^2)^(p-1)-1 && a=precprime(a-1), ); a }

%o forprime(p=1, default(primelimit), my(a=p); for(c=1, 4, (a=fp(a))||next(2)); print1(p, ", "))

%K nonn

%O 1,1

%A _M. F. Hasler_, Jan 10 2011