%I #13 Apr 09 2019 03:20:10
%S 2,11,1093,3511,20771,534851,1006003,3152573
%N Primes p such that Omega(p + 1)^(p - 1) == 1 (mod p^2), where Omega is A001222.
%C a(9) > 807795277 if it exists.
%C a(9) > 3.5*10^12 if it exists. - _Giovanni Resta_, Apr 09 2019
%e A001222(20772) = 5 and 5^(20771-1) == 1 (mod 20771^2), so 20771 is a term of the sequence.
%t Select[Prime@ Range@ 230000, PowerMod[ PrimeOmega[# + 1], #-1, #^2] == 1 &] (* _Giovanni Resta_, Apr 09 2019 *)
%o (PARI) forprime(p=1, , if(Mod(bigomega(p+1), p^2)^(p-1)==1, print1(p, ", ")))
%Y Cf. A001220, A001222, A260377, A267487.
%K nonn,hard,more
%O 1,1
%A _Felix Fröhlich_, Mar 16 2019
|