Primality proof for n = 1621:

Take b = 2.

b^(n-1) mod n = 1.

5 is prime.
b^((n-1)/5)-1 mod n = 1488, which is a unit, inverse 195.

3 is prime.
b^((n-1)/3)-1 mod n = 1435, which is a unit, inverse 61.

(3^4 * 5) divides n-1.

(3^4 * 5)^2 > n.

n is prime by Pocklington's theorem.