%I #13 May 28 2024 11:45:15
%S 1,1,2,3,4,1,1,1,2,1,1,1,1,2,1,1,1,1,1,1,3,1,1,1,1,1,1,1,1,1,1,4,1,1,
%T 1,1,1,1,2,1,1,1,1,1,1,1,5,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,2,1,1,1,1,1,
%U 1,1,1,2,1,1,1,1,1,1,2,2,1,2,1,2,1,1,1,1,1,1,1,1,1,2,1,2,2
%N Greatest k >= 1 such that (p + 1 - 2^i) / 2^i is prime for i = 1..k and p is prime from A005385.
%C "k-safe primes" iff (p + 1 - 2^i) / 2^i is prime for i = 0..k, p a prime number. A000040 are 0-safe primes, A005385 are 1-safe primes, A066179 are 2-safe primes.
%e p = 11: (11 + 1 - 2^i) / 2^i is prime for i = 1..2, thus a(3) = 2.
%e p = 47: (47 + 1 - 2^i) / 2^i is prime for i = 1..4, thus a(5) = 4.
%o (PARI) isp(k) = (denominator(k) == 1) && isprime(k);
%o f(p) = my(i=1); while (isp((p+1-2^i)/2^i), i++); i-1;
%o apply(f, select(x->isp((x-1)/2), primes(1000))) \\ _Michel Marcus_, May 28 2024
%Y Cf. A000040, A005385, A066179.
%K nonn
%O 1,3
%A _Ctibor O. Zizka_, May 21 2024