login
Primes of the form 2*p-1 where p is prime and p-1 is not squarefree.
7

%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