# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a339181
Showing 1-1 of 1

%I A339181 #10 Nov 30 2020 01:05:23
%S A339181 17,31,151,241,577,3001,3571,4801,12097,21121,23761,28513,61441,65521,
%T A339181 77761,113023,126001,171697,174721,178753,193441,244901,287281,364801,
%U A339181 582427,616897,677321,976501,1016401,1425601,1431847,2015441,2080801,2483713,2672671,3089371,4321931,4667921,5177761
%N A339181 Primes p such that p == 1 (mod A001414(p-1)) and p == 1 (mod A001414(p+1)).
%C A339181 Members p of A339180 such that p == 1 (mod A001414(p-1)).
%H A339181 Robert Israel, <a href="/A339181/b339181.txt">Table of n, a(n) for n = 1..198</a>
%e A339181 a(3) = 151 is in the sequence because 151 is prime, A001414(150)=2+3+5+5=15, A001414(152)=2+2+2+19=25, 151 == 1 (mod 15) and 151 == 1 (mod 25).
%p A339181 spf:= n -> add(t[1]*t[2], t=ifactors(n)[2]):
%p A339181 select(p -> isprime(p) and p mod spf(p-1) = 1 and p mod spf(p+1) = 1, [seq(i, i=3..6*10^6, 2)]);
%Y A339181 Cf. A001414, A339180.
%K A339181 nonn
%O A339181 1,1
%A A339181 _J. M. Bergot_ and _Robert Israel_, Nov 26 2020

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE