# find all squares modulo a prime number

I have been using PARI-GP to find the squares of a prime number. For example

```
for(JJ=0,7,print1(JJ^2%8"\t"))
```

in PARI gives

```
0 1 4 1 0 1 4 1
```

from which I know that the squares modulo 8 are 1 and 4. Now that I need to find the squares modulo large prime, I was wondering is there a way SAGE can give me a complete list without repeating the values.