ああ
a(n)=2^n/7 (4+3cos(n(Pi-Arccot(3/7^(1/2))))+7^(1/2)sin((n(Pi-Arccot(3/7^(1/2))))

です。
計算するとたしかに平方数になっていますが、チョットでかけるので