以知多项式3x^3+ax^2+bx+1能被x^2+1整除,且商式为3x+1,求(-a)^b的值

3x^3+ax^2+bx+1能被x^2+1整除,且商式为3x+1
(x^2+1)(3x+1)
=3x^3+x^2+3x+1
=3x^3+ax^2+bx+1
a=1
b=3
(-a)^b=(-1)^3=-1

Matlab的做法：
p1=[1 0 1]
p2=[3 1]
p=conv（p1，p2) %多项式的乘法用卷积函数conv（）;；; ;
ans=
3 1 3 1
所以
a=1,b=3;
(-a)^b=-1^3=-1 。

因为原式能被x^2+1整除，且商式为3x+1
所以原式为3x^3+x^2+3x+1
所以a=1 b=3
所以(-a)^b的值为-1