PellP = DifferenceRoot[ Function[{y, n}, {-y[n] - 2 y[1 + n] + y[2 + n] == 0, y[0] == 0, y[1] == 1}]];
DifferenceRootReduce[ PellP[n] == (-(1 - Sqrt[2])^n + (1 + Sqrt[2])^n)/(2 Sqrt[2]), n]
DifferenceRootReduce[PellP[n + 1] PellP[n - 1] - PellP[n]^2 == (-1)^n, n]
DifferenceRootReduce[\!\( \*UnderoverscriptBox[\(\[Sum]\), \(i = 0\), \(4\ n + 1\)]\(PellP[ i]\)\) == (PellP[2 n + 1] + PellP[2 n])^2, n]