|
%I #5 Oct 01 2013 21:35:30
%S 31,227,1051,3109,7151,15511,18127,30367,32143,32687,144719,151051,
%T 165311,186343,234191,302399,369997,371281,374239,407503,454303,
%U 509263,531263,537743,759359,1053007,1088081,1182287,1185601,1354321,1381441
%N Primes of the form x^5-y^4, where x,y >= 1.
%C If a prime has multiple representations of the format, it is entered only once.
%F If x^5-y^4 is prime for integers x,y list without duplicates.
%e 2^5 - 1^4 = 31.
%o (PARI) diffpowers(n,m) =
%o {
%o local(a,c=0,c2=0,j,k,y);
%o a=vector(floor(n^2/log(n^2)));
%o for(j=1,n,
%o for(k=1,n,
%o y=j^m-k^(m-1);
%o if(ispseudoprime(y),
%o c++;
%o \\ print(j","k","y);
%o a[c]=y;
%o );
%o );
%o );
%o a=vecsort(a);
%o for(j=2,length(a),
%o if(a[j]!=a[j-1]&&a[j]!=0,
%o c2++;
%o print1(a[j]",");
%o if(c2>100,break);
%o );
%o );
%o }
%K nonn
%O 1,1
%A _Cino Hilliard_, Jun 17 2009
|
