%I #16 Aug 11 2024 10:58:13
%S 37,73,193,313,397,457,541,613,673,757,1153,1201,1321,1453,1621,1657,
%T 1753,1873,1993,2017,2137,2341,2473,2557,2593,2857,2917,3061,3217,
%U 3313,4057,4177,4273,4357,4441,4561,4933,5077,5101,5113,5233,5437,5581,5701
%N Primes of the form 2*p-1 where p is prime and p-1 is not squarefree.
%C Subsequence of A005383.
%H Robert Israel, <a href="/A153208/b153208.txt">Table of n, a(n) for n = 1..10000</a>
%e For p = 2 (the only case with p-1 odd), 2*p-1 = 3 is prime but p-1 = 1 is squarefree, so 3 is not in the sequence. For p = 19, 2*p-1 = 37 is prime and p-1 = 18 is not squarefree, so 37 is in the sequence.
%p R:= NULL; count:= 0: p:= 3:
%p while count < 100 do
%p p:= nextprime(p);
%p if isprime(2*p-1) and not numtheory:-issqrfree(p-1) then
%p R:= R, 2*p-1; count:= count+1;
%p fi
%p od:
%p R; # _Robert Israel_, Nov 22 2023
%t lst={}; Do[p = Prime[n]; If[ !SquareFreeQ[Floor[p/2]] && PrimeQ[Ceiling[p/2]], AppendTo[lst, p]], {n, 7!}]; lst
%t Select[2#-1&/@Select[Prime[Range[1000]],!SquareFreeQ[#-1]&],PrimeQ] (* _Harvey P. Dale_, Aug 11 2024 *)
%o (Magma) [ q: p in PrimesUpTo(2900) | not IsSquarefree(p-1) and IsPrime(q) where q is 2*p-1 ];
%Y Cf. A013929 (nonsquarefree numbers), A005383 (numbers n such that both n and (n+1)/2 are primes), A153207, A153209, A153210.
%K nonn
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Dec 20 2008
%E Edited by _Klaus Brockhaus_, Dec 24 2008
%E Mathematica updated by _Jean-François Alcover_, Jul 04 2013