mign := proc(a,n,m,k)
# n>>m>1, n>>k>1, |a|>1,  (n=20, m=3,k=3, a=100 I)
 (x^(n-k)+(a*x+1)^m)*(a*x+1)^k;
end: