EulerPhi[n^k] = EulerPhi[n] n^{k - 1}
EulerPhi[n^k]/n^{k - 1} = EulerPhi[n] =1210=2*5*11²
Prime(x)^{k-1}*(Prime(x)-1)=1210 -> k=3, x=5 und Prime(5)=11
Lösungen:
±1331=11³
±2662=2*11³
bis 100000 keine weiteren Lösungen...
Anderer Weg:
EulerPhi[n] EulerPhi[m] = (EulerPhi[m n] EulerPhi[ggT[n, m]])/ggT[n, m]
... reicht für heute...