(PARI) a(n) = fromdigits(Vec(Pol(digits(11))^n)%10); \\
Seiichi Manyama, Mar 10 2023
(Python)
from math import comb, prod
from sympy.ntheory.modular import crt
from gmpy2 import digits
k, l = 0, len(s:=digits(n, 5))
for m in range(n+1):
t = digits(m, 5).zfill(l)
k = 10*k+crt([5, 2], [prod(comb(int(s[i]), int(t[i]))%5 for i in range(l))%5, int(not ~n & m)])[0]