#include <bits/stdc++.h>
using namespace std;
struct Dbg {
  template <typename T>
  Dbg& operator,(const T& v) {
    cerr << v << ", ";
    return *this;
  }
} Dbg;
template <typename T>
ostream& operator<<(ostream& _o_, const vector<T>& _v_) {
  if (!_v_.empty()) {
    _o_ << '[';
    copy(_v_.begin(), _v_.end(), ostream_iterator<T>(_o_, ", "));
    _o_ << "\b\b]";
  }
  return _o_;
}
template <class C>
void mini(C& _a, C _b) {
  _a = min(_a, _b);
}
template <class C>
void maxi(C& _a, C _b) {
  _a = max(_a, _b);
}
int MOD;
map<long long, int> cnt;
long long pw(int a, int b) {
  long long res = 1;
  for (int i = 0; i < b; i++) {
    res *= a;
    res %= MOD;
  }
  return res;
}
int sub(int a, int b) {
  a -= b;
  if (a < 0) a += MOD;
  return a;
}
long long mul(int a, int b) {
  long long res = a;
  res *= b;
  res %= MOD;
  return res;
}
vector<int> scn(int n) {
  vector<int> a(n);
  for (int i = 0; i < n; i++) {
    cin >> a[i];
  }
  return a;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  cout.tie(0);
  int n, k;
  cin >> n >> MOD >> k;
  vector<int> a(n);
  for (int i = 0; i < n; i++) {
    cin >> a[i];
  }
  long long ans = 0;
  for (int i = 0; i < n; i++) {
    cnt[mul(a[i], (sub(pw(a[i], 3), k)))]++;
  }
  for (auto x : cnt) {
    ans += (x.second * (x.second - 1)) / 2;
  }
  cout << ans;
  return 0;
}
