如何用mathematica有效地计算多项式展开系数？

Question

我想在Mathamatica中得到以下多项式展开的系数：

poly = 1 + Sum[a[i]*x^(i), {i, 1, imax}] ;
power = Expand[(poly)^(p), x];
coefs = CoefficientList[power,x];

上面的代码可以工作但是它会为大的imax和power p使用大量的内存。

现在我的问题是： 有没有一种有效的方法来获得任何功率系数x？

谢谢！