A300566 Numbers z such that there is a solution to x^4 + y^5 = z^6 with x, y, z >= 1.
8748, 10368, 342732
Offset: 1
Examples
a(1) = 8748 = 2^2*3^7 is in the sequence because 8748^6 = (2^3*3^8)^5 + (2^3*3^10)^4, using 2^3 + 1 = 3^2. Similarly, all z = 4*3^(10k-3) are in the sequence for k >= 1, with x = 8*3^(15k-5) and y = 8*3^(12k-4). a(2) = 10368 = 2^7*3^4 is in the sequence because 10368^6 = (2^8*3^5)^5 + (2^10*3^6)^4, using 3 + 1 = 2^2. Similarly, any z = 2^7*3^(10k+4) is in the sequence for k >= 0, with x = 2^10*3^(15k+6) and y = 2^8*3^(12k+5). z = 342732 = 2^2*3*13^4 is in the sequence because (2^2*3*13^4)^6 = (2^3*13^5)^5 + (2^3*5*13^6)^4, using 2^3*13 + 5^4 = 3^6. z = 810000 = 2^4*3^4*5^4 is in the sequence because z^6 = x^4 + y^5 with x = 2^5*3^6*5^6 and y = 2^4*3^5*5^5 (using 1 + 3*5 = 2^4). z = 1361367 = 3^4*7^5 is in the sequence because z^6 = x^4 + y^5 with x = 3^5*7^8 and y = 2*3^4*7^6.
Crossrefs
Programs
-
PARI
is(z)=for(y=1,sqrtnint(-1+z=z^6,5),ispower(z-y^5,4)&&return(y)) /* Code below for illustration only, not guaranteed to give a complete list. Half-integral values give the additional term 31505922 for b = 63/2. Third-integral values give the additional solution z = 342732 for b = 26/3. */ S=[]; N=1e5; forstep(b=1,9,1/3, forstep(a=1,N,1/3, issquare(b^10+a^4<<2,&r)&& !frac(z=b^5/2+r/2)&& !print1(z",")&&S=setunion(S,[z])); print1([b])); S
Comments