By way of a more targeted hint there is this MSE post which shows how to solve the general problem. You need to compute the cycle index of the permutation group of the grid acting on the squares.
Here is the (trivial) cycle index for $n=1$:
$$ a_{{1}}$$
and for $n=2$
$$ 1/8\,{a_{{1}}}^{4}+3/8\,{a_{{2}}}^{2}+1/4\,{a_{{1}}}^{2}a_{{2}}+1/4\,a_{{4}}$$
for $n=3$
$$1/8\,{a_{{1}}}^{9}+1/2\,{a_{{1}}}^{3}{a_{{2}}}^{3}+1/8\,a_{{1}}{a_{{2}}}^{4}+
1/4\,a_{{1}}{a_{{4}}}^{2} $$
for $n=4$
$$1/8\,{a_{{1}}}^{16}+3/8\,{a_{{2}}}^{8}+1/4\,{a_{{1}}}^{4}{a_{{2}}}^{6}+1/4\,{
a_{{4}}}^{4} $$
for $n=5$
$$1/8\,{a_{{1}}}^{25}+1/2\,{a_{{1}}}^{5}{a_{{2}}}^{10}+1/8\,a_{{1}}{a_{{2}}}^{
12}+1/4\,a_{{1}}{a_{{4}}}^{6} $$
and for $n=6$
$$1/8\,{a_{{1}}}^{36}+3/8\,{a_{{2}}}^{18}+1/4\,{a_{{1}}}^{6}{a_{{2}}}^{15}+1/4
\,{a_{{4}}}^{9}. $$
Shading exactly two of the squares you want to compute
$$[z^2] Z(G)(1+z)$$
where $Z(G)$ is the appropriate cycle index which gives the sequence
$$ 0, 2, 8, 21, 49, 93, 171, 278, 446, 660, \ldots$$
This is OEIS A014409 which seems to match your problem definition so this is probably correct.
Using pen and paper to verify the value $8$ for $n=3$ we get: colorings where a corner is colored and paired: 1 -- with an adjacent square, 2 -- with a square in the same row/column, 3 -- with the opposite corner, 4 -- with the midpoint of the opposite edge and 5 -- with the center square. Then there is 6 -- the center square paired with an adjacent square and 7 -- two opposite middle squares of two opposite sides and finally, 8 -- two side midpoint squares that share a corner.
Addendum. Inspired by MJDs work let me quote the general formula for the cycle index with $n$ even from the post I cited in the introduction, which is
$$1/8\,{a_{{1}}}^{{n}^{2}}+3/8\,{a_{{2}}}^{1/2\,{n}^{2}}+1/4\,{a_{{1}}}^{n}{a_{
{2}}}^{1/2\,{n}^{2}-1/2\,n}+1/4\,{a_{{4}}}^{1/4\,{n}^{2}}.$$
Substituting $1+z$ into this cycle index we get
$$1/8\,{(1+z)}^{{n}^{2}}+3/8\,{(1+z^2)}^{1/2\,{n}^{2}}+1/4\,{(1+z)}^{n}{(1+z^2)}^{1/2\,{n}^{2}-1/2\,n}+1/4\,{(1+z^4)}^{1/4\,{n}^{2}}.$$
Extracting the coefficient $[z^2]$ from this we see that the last term does not contribute (polynomial in $z^4$). The remaining contributions are (use the binomial theorem)
$$1/8 {n^2\choose 2} + 3/8 {1/2\times n^2\choose 1} + 1/4 {n\choose 2} +
1/4 {1/2 \times n^2 - 1/2\times n\choose 1}.$$
Expanding these binomial coefficients which are simple, fortunately, we obtain the formula
$$1/16\,{n}^{4}+3/8\,{n}^{2}-1/4\,n$$
which matches the result by MJD.