已知数列{an},a1=1,a(n+1)=an^2+an+1,求通项公式an