(PARI) list(lim)=my(upper=log(lim*sqrt(5))\log((1+sqrt(5))/2)+1, t, tt, v=List([1])); if(fibonacci(t)>lim, t--); for(i=3, upper, t=fibonacci(i); for(j=2, i-1, tt=t+fibonacci(j); if(tt>lim, break); if(ispower(tt), listput(v, tt)))); Set(v) \\
Charles R Greathouse IV, May 03 2016