Note that the question is "how to determine" the sequence, rather than simply find the function which controls it. So below is a method I've used to try to find what controls the sequence.
If the sequence is 1, 21, 467, 5216, and we can take rough logs (base 10) of these numbers we get:
0, 1.32, 2.67, 3.72 which implies that the sequence could be logarithmic, because the ratios are roughly arithmetic with a common difference about 1.3 which is roughly the log of 20. This would generate the sequence 1, 20, 400, 8000. 8000 is much greater than 5216, so let's see what happens if we use 19 instead of 20: 1, 19, 361, 6859. The last term is still larger than we would like but the second term is too small.
Now try something else. Suppose the sequence is of the form anxnn+bn for each term, where integer n starts at 0, and x, an and bn are constants then:
b0=1, and the other terms are a1x+b1=21, a2x2+b2=467, a3x3+b3=5216.
Let's try some possible values for x (some figures are approximate) :
x=15:
15a1+b1=21; 225a2+b2=467; 3375a3+b3=5216, a1=(21-b1)/15=1.4-b1/15; a2=(467-b2)/225=2.076-b2/225; a3=(5216-b3)/3375=1.545-b3/3375;
x=16:
16a1+b1=21; 256a2+b2=467; 4096a3+b3=5216; a1=(21-b1)/16=1.313-b1/16; a2=(467-b2)/256=1.824-b2/256; a3=(5216-b3)/4096=1.545-b3/4096;
x=17:
17a1+b1=21; 289a2+b2=467; 4913a3+b3=5216; a1=(21-b1)/17=1.235-b1/17; a2=(467-b2)/289=1.616-b2/289; a3=(5216-b3)/4913=1.062-b3/4913;
x=18:
18a1+b1=21;324a2+b2=467; 4913a3+b3=5832; a1=(21-b1)/18=1.167-b1/18; a2=(467-b2)/324=1.441-b2/324; a3=(5216-b3)/5832=0.894-b3/5832.
From these examples we hope to get a common value for an, or a pattern in terms of n for an, and a pattern in terms of n for bn. Unfortunately no such patterns emerge.
It's always possible to find a formula for the nth term of a series when there are so few terms. The result is likely to be a formula with ungainly coefficients, but nevertheless it illustrates a method.
Let T(n)=a0+a1n+a2n2+a3n3, where an values are constants.
T(0)=a0=1;
T(1)=a0+a1+a2+a3=21, a1+a2+a3=20;
T(2)=a0+2a1+4a2+8a3=467, 2a1+4a2+8a3=466⇒a1+2a2+4a3=233;
T(3)= 3a1+9a2+27a3=5215.
Now we have a system of 3 equations and 3 unknowns.
(1) T(2)-T(1)=a2+3a3=213
(2) T(3)-3T(1)=6a2+24a3=5155
(3) (2)-6(1)=6a3=3877, a3=3877/6⇒a2=-3451/2⇒a1=3298/3.
So T(n)=1+3298n/3-3451n2/2+3877n3/6. This produces the series: 1, 21, 467, 5216, 18145, 43131, ...
However, it's unlikely to be the right solution.