kill(all)$

van=facebook.com/LearningCalculus$

f(k,p):=2^k-k=3*p^2$

P:primes(1,10000)$

length(P);

P[length(P)];

for k:1 thru 1000 do
    for i:1 thru length(P) do
        if(f(k,P[i])) then print("Answer: Solution: ", k,P [i])$

print("Answer: ", P[length(P)] * 1000, " combinations checked.")$

`;