Ich denke, man sollte im Falle p=q mit dem chinesischen Restsatz
arbeiten: Z/pqZ≅Z/pZ×Z/qZ. Für die invertierbaren Elemente gilt
dann (Z/pqZ)∗≅(Z/pZ)∗×(Z/qZ)∗.
Selbstinvers sind dann die Elemente, die den Paaren
(1,1),(1,−1),(−1,1) und (−1,−1) des Gruppenprodukts
entsprechen.
Beispiel: p=3 und q=7:
a≡1 mod 3, a≡1 mod 7 liefert a≡1 mod 21,
a≡1 mod 3, a≡−1 mod 7 liefert a≡13 mod 21,
a≡−1 mod 3, a≡1 mod 7 liefert a≡8 mod 21,
a≡−1 mod 3, a≡−1 mod 7 liefert a≡20≡−1 mod 21,
Also findet man als betragsmäßig kleinste selbstinverse Reste:
±1,±8.