原贴链接：here。

现在问题已经得到解决（当然原贴里的人没帮上忙，我自己想的）：解决代码

但是很明显，在 $eps$ 足够小的时候这个代码的时间复杂度依然不可忍受（比如 $eps=10^{-16}$ 或更低），在这种情况下，我的代码对于 $a=3.1415926535897932$ 的情况跑了十多秒才出结果。

现在我的问题变成了：如何更快地出结果？

求助。