What is the remainder of dividing polynomial $$P(x)=x^{2016}-x^{2015}-1$$ with $x^2+1$?
So what I thought of doing is just dividing them the "school" way:
$(x^{2016}-x^{2015}-1)\div(x^2+1)=x^{2014}-x^{2013}-x^{2012}+x^{2011}\cdots$
But the problem is that this way it just goes on and on, and I can see that it should end up with $x$ as a remainder, but what's the way I can prove that, without knowing that it ends up in $x$?
The other way I thought of solving this is by maybe presenting that:
$P(x)\div(x^2+1)=Q(x)+y$,
so for $P(i)$
$i^{2016}-i^{2015}-1=y$,
where $i^2=-1$, and this way I get
$y=-i$
but that seems incorrect, since I can see that the regular division would end up to $x$.
Your approach works except that it is not fully general: because $x^2+1$ is quadratic, the quotient-remainder form should be $$ P(x) = Q(x)(x^2+1) + ax+b $$ for some unknown $a$ and $b$.
Now we can set $P(i) = Q(i)(0) + ai + b$, getting one equation for $a$ and $b$. We can also set $P(-i) = Q(i)(0) + a(-i) + b$, getting another equation, because $x^2+1$ is $0$ when $x=i$ or when $x=-i$.
Solving for $a$ and $b$, you should get $a=1$ and $b=0$.