Submission #994488

---

Source Code Expand

#include <bits/stdc++.h>

using namespace std;

#define REP(i,a,b) for(int i=a;i<(int)b;i++)
#define rep(i,n) REP(i,0,n)
#define all(c) (c).begin(), (c).end()
#define zero(a) memset(a, 0, sizeof a)
#define minus(a) memset(a, -1, sizeof a)
#define watch(a) { cout << #a << " = " << a << endl; }
template<class T1, class T2> inline bool minimize(T1 &a, T2 b) { return b < a && (a = b, 1); }
template<class T1, class T2> inline bool maximize(T1 &a, T2 b) { return a < b && (a = b, 1); }

typedef long long ll;
int const inf = 1<<29;

namespace tree {
struct union_find {
  vector<int> par, rank, size;
  int compnum;

  union_find(int N) {
    compnum = N;
    par.resize(N), rank.resize(N), size.resize(N);
    for(int i=0; i<N; i++) {
      par[i] = i;
      rank[i] = 0;
      size[i] = 1;
    }
  }

  int root(int x) {
    return par[x] == x ? x : par[x] = root(par[x]);
  }

  void unite(int x, int y) {
    x = root(x), y = root(y);
    if(x == y) return;
    if(rank[x] < rank[y]) {
      par[x] = y, size[y] += size[x];
    } else {
      par[y] = x, size[x] += size[y];
      if(rank[x] == rank[y]) rank[x]++;
    }
    compnum--;
  }

  int operator[](int x) { return root(x); }
  void operator()(int x, int y) { return unite(x, y); }

  bool same(int x, int y) { return root(x) == root(y); }
  int size_of(int x) { return size[root(x)]; }
  int num_of_comps() { return compnum; }
};
}

int main() {

  int N, M; cin >> N >> M;
  vector<int> X(N);
  rep(i, N) {
    cin >> X[i];
  }

  sort(all(X));

  tree::union_find uf(N);

  // 元のindexが同じUFをする
  rep(i, N - 1) {
    if(X[i] == X[i + 1]) {
      uf(i, i + 1);
    }
  }

  map<int, int> st;

  // MODが同じ数で、元の値ごとに数え上げる
  unordered_map<int, unordered_map<int, bool>> odd;
  rep(i, N) {
    odd[X[i] % M][X[i]] ^= 1;
    // Mでmodをとり
    X[i] %= M;
    // 端数の一致する対称を詰める
    st[(M - X[i]) % M] ++;
  }

  auto tmpST = st;

  unordered_map<int, int> oddcnt;
  for(auto const& e: odd) {
    int ocnt = 0;
    for(auto const& u: e.second) {
      if(u.second) ocnt++;
    }
    oddcnt[e.first] = ocnt;
  }

  ll ans = 0;

  rep(i, N) {
    auto tmp = (M - X[i]) % M;

    if(st[tmp] == 0) continue;

    // 自分は選んだので消す
    st[tmp]--;

    // 自分の端数対称が一致するものを見つけては消す
    if(st[X[i]] > 0) {
      st[X[i]]--;
      ans ++;
    }
    else {
      // 復活
      st[tmp]++;
    }
  }

  // 残った (M - X[i]) % M を用いて、元の数が同じものを出来る限り生成する

  // MODが同じ値の個数のうち、既に使ったMODの数をできるだけ各値の奇数の端数から削る。

  // MODが同じ値の個数 / 2を単純に足した場合、解以上になることを確かめる

  for(auto const& e: st) {
    // 残った (M - X[i]) % M のmapについて、
    // まず、X[i] % M に直す。この個数は未使用分の個数。そのうち使えるものだけを取り除きたい。
    // 上で、奇数の端数から削っていた、と考える
    // つまり、元のMOD一致ペア数から、奇数の数/2を引いて、0とmaxを取る。これをmpnumとする。
    // mpnumをe.second/2から引いて足し合わせる
    auto x = (M - e.first) % M;
    auto mpnum = max(0, (tmpST[x] - oddcnt[x] / 2));
    ans += e.second / 2 - mpnum;
  }

  cout << abs(ans) << endl;
  
  return 0;
}

Submission Info

Judge Result