#include <bits/stdc++.h>
#pragma GCC optimize(2)
using namespace std;
const int INF = 0x3f3f3f3f;
const int MOD = 1e9 + 7;
const double eps = 1e-6;
const double pi = 3.141593;
int gcd(int a, int b) { return b == 0 ? a : gcd(b, a % b); }
long long qpow(long long x, long long y) {
  long long ans = 1;
  for (; y > 0; y >>= 1) {
    if (y & 1) ans *= x;
    x *= x;
  }
  return ans;
}
long long qpow(long long x, long long y, int MOD) {
  long long ans = 1;
  for (; y > 0; y >>= 1) {
    if (y & 1) ans = ans * x % MOD;
    x = x * x % MOD;
  }
  return ans;
}
void exgcd(long long a, long long b, long long &x, long long &y) {
  if (!b) {
    x = 1;
    y = 0;
    return;
  }
  exgcd(b, a % b, x, y);
  long long X = x;
  x = y;
  y = X - a / b * y;
}
long long inv(long long x, int mod = MOD) { return qpow(x, mod - 2) % mod; }
long long m_m(long long a) { return (a % MOD + MOD) % MOD; }
int main() {
  ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
  int n, p, k;
  cin >> n >> p >> k;
  map<int, int> mp;
  for (int i = (0); i < (n); ++i) {
    int x;
    cin >> x;
    x = qpow(x, 4, p) - 1ll * k * x % p;
    if (x < 0) x += p;
    ++mp[x];
  }
  long long ans = 0;
  for (auto &it : mp) {
    ans += 1ll * it.second * (it.second - 1) / 2;
  }
  cout << ans;
  return 0;
}
