polynomial.<< (k) is at fault, fix it and then go with expan_powers2(p)