"orsoulx":
$ a_0=0; a_1=1; a_{k+1}=2a_k+a_{k-1} $
$n_k=2a_k^2 $ se $ k $ è pari; $ n_k=2a_k^2-1 $ se è dispari.

Ma non ho capito che razza di percorso logico hai fatto!
I primi termini della tua successione ${a_k}$ sono:
0, 1, 2, 5, 12, 29, 70, 169, 408, 985,...
Quelli della tua seconda successione vengono allora
0, 1, 8, 49, 288, 1681, 9800,57121,...
e sono effettivamente quelli della richiesta successione ${n_k}$.
"orsoulx":
Ricavati dalle terne pitagoriche 'quasi' isosceli (con i cateti che differiscono di 1)

Queste sono (iniziando dalla terna degenere $T_0$ = [0, 1, 1]):
[0, 1, 1]; [3, 4, 5]; [20, 21, 29]; [119, 120, 169]; [696, 697, 985]; ---
Come vettori tridimensionali sono la progressione geometrica
$T_r = T_0 · M^r$ (per ogni $r$ naturale).
dove $T_0 = [0, 1, 1]$ ed $M$ è la matrice (quadrata simmetrica):
!2, 1, 2|
|1, 2, 2|.
|2, 2, 3|
Vedo che nella tua prima successione ${a_k}$ i termini per k dispari– diciamo k = 2r + 1 – sono le "ipotenuse" delle dette terne pitagoriche $T_r$.
Ma i termini di indice pari (e quelli della tua seconda successione)... in che relazione stanno con le dette terne pitagoriche?
[La tua spiegazione mi risulta alquanto sibillina!]
E' però vero che, tanto le componenti delle terne pitagoriche con differenza costante tra i cateti quanto gli indici n richiesti dal quiz (e termini della tua seconda successione) (e pure le due successioni ${a_(2k)}$ e ${a_(2k+1)}$ dei termini rispettivamente di indice pari e dispari della tua ${a_k}$) verificano la medesima legge di ricorrenza:
$ y_(n+3) = 7·y_(n+2) – 7·y_(n+1) + y_n$.
Ossia: hanno polinomio caratteristico
$P_3(x) = x^3 - 7·x^2 + 7·x – 1 = (x -1)(x^2 -6x +1) =$
$= (x – 1)·[x – (3 + 2sqrt2)]·[x – (3 - 2sqrt2)]$.
Pertanto, il termine corrente delle succession delle componenti di ${T_r}$ e delle successioni ${a_(2k)}$, ${a_(2k+1)$ e ${n_k}$ espresso in funzione dell'indice è del tipo:
$y(k) = A + B·(3 + 2sqrt2)^k + C·(3 - 2sqrt(2)^k$.
Per la richiesta successione ${n_k}$ le costanti (A, B, C) valgono rispettivamente ($-1/2$, $1/4$, $1/4$).
In definitiva essa risulta:
$n_k =((3+2sqrt2)^k + (3 - 2sqrt2)^k)/4 - 1/2$.