Problem p03440 (C++)

Buggy submission (status: Time Limit Exceeded):

#include <algorithm>
#include <array>
#include <bitset>
#include <deque>
#include <iostream>
#include <map>
#include <math.h>
#include <queue>
#include <set>
#include <stack>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <unordered_map>
#include <vector>

long long mod = 1e9 + 7;
class Mint {
public:
  long long x;
  Mint(){}; // 引数なしでも定義できるように引数なしコンストラクタも用意しておく
  Mint(long long a) {
    x = a % mod;
    while (x < 0) {
      x += mod;
    }
  };
  Mint &operator+=(const Mint &a) {
    x += a.x;
    x %= mod;
    return *this;
  }
  Mint &operator-=(const Mint &a) {
    x += (mod - a.x);
    x %= mod;
    return *this;
  }
  Mint &operator*=(const Mint &a) {
    x *= a.x;
    x %= mod;
    return *this;
  }

  // a^n mod を計算する
  // Useaeg: Mint z = Mint(2).pow(n);
  Mint pow(long long n) const {
    if (n == 0)
      return Mint(1);
    Mint y = pow(n >> 1); // pow(n/2)を計算する
    y *= y;
    if (n % 2 == 1)
      y *= *this;
    return y;
  }

  // a^{-1} mod を計算する
  Mint modinv(const Mint &a) const { return a.pow(mod - 2); }

  Mint &operator/=(const Mint &a) {
    x *= modinv(a).x;
    x %= mod;
    return *this;
  }

  Mint operator+(Mint &a) const {
    Mint y(*this);
    y += a;
    return y;
  }

  Mint operator-(Mint &a) const {
    Mint y(*this);
    y -= a;
    return y;
  }

  Mint operator*(Mint &a) const {
    Mint y(*this);
    y *= a;
    return y;
  }

  Mint operator/(Mint &a) const {
    Mint y(*this);
    y /= a;
    return y.x;
  }

  // nCk @mod を計算する
  Mint nCk(Mint &n, const long long k) const {
    Mint y = Mint(1);
    Mint one = Mint(1);
    for (Mint i(0); (i.x) < k; i.x++) {
      y *= (n - i);
      y /= (i + one);
    }
    return y;
  }

  // nPk @mod を計算する
  Mint nPk(Mint &n, long long k) const {
    Mint y(1);
    for (Mint i(0); (i.x) < k; i.x++) {
      y *= (n - i);
    }
    return y;
  }
};

class UnionFind {
private:
  int n;
  std::vector<int> parent;

public:
  UnionFind(int num) : n(num) {
    parent = std::vector<int>(n);
    for (int i = 0; i < n; i++) {
      parent[i] = i;
    }
  }
  int GetRoot(int i) {
    if (parent[i] != i) {
      return UnionFind::GetRoot(parent[i]);
    } else {
      return i;
    }
  }

  bool SameGroup(int n, int m) {
    n = this->GetRoot(n);
    m = this->GetRoot(m);
    return n == m;
  }

  bool Unite(int n, int m) {
    n = this->GetRoot(n);
    m = this->GetRoot(m);
    if (n == m)
      return false;
    parent[n] = m;
    return true;
  }
};

int getgcd(int a, int b) {
  int max = std::max(a, b);
  int min = std::min(a, b);
  if (min == 0)
    return max;
  return getgcd(min, max % min);
}

int main() {
  int N, M;
  std::cin >> N >> M;
  std::vector<int> a(N);
  for (int i = 0; i < N; i++) {
    std::cin >> a[i];
  }
  UnionFind uf(N);
  int A, B;
  for (int i = 0; i < M; i++) {
    std::cin >> A >> B;
    uf.Unite(A, B);
  }
  std::vector<int> group[N];
  for (int i = 0; i < N; i++) {
    int root = uf.GetRoot(i);
    group[root].push_back(a[i]);
    // printf("root of %d = %d\n", i, root);
  }
  long long ans = 0;
  std::vector<int> amari;
  int groupnum = 0;
  for (int i = 0; i < N; i++) {
    if (group[i].size() > 0)
      groupnum++;
  }
  if (groupnum == 1) {
    printf("0\n");
    return 0;
  }

  for (int i = 0; i < N; i++) {
    if (group[i].size() == 0)
      continue;
    std::sort(group[i].begin(), group[i].end());
    ans += group[i][0];
    // printf("ans += %d\n", group[i][0]);
    for (int j = 1; j < group[i].size(); j++) {
      amari.push_back(group[i][j]);
      // printf("amari pushing %d\n", group[i][j]);
    }
  }
  if (amari.size() < groupnum - 2) {
    printf("Impossible\n");
  } else {
    std::sort(amari.begin(), amari.end());
    for (int j = 0; j < groupnum - 2; j++) {
      ans += amari[j];
    }
    std::cout << ans << std::endl;
  }

  return 0;
}

Fixed submission:

#include <algorithm>
#include <array>
#include <bitset>
#include <deque>
#include <iostream>
#include <map>
#include <math.h>
#include <queue>
#include <set>
#include <stack>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <unordered_map>
#include <vector>

long long mod = 1e9 + 7;
class Mint {
public:
  long long x;
  Mint(){}; // 引数なしでも定義できるように引数なしコンストラクタも用意しておく
  Mint(long long a) {
    x = a % mod;
    while (x < 0) {
      x += mod;
    }
  };
  Mint &operator+=(const Mint &a) {
    x += a.x;
    x %= mod;
    return *this;
  }
  Mint &operator-=(const Mint &a) {
    x += (mod - a.x);
    x %= mod;
    return *this;
  }
  Mint &operator*=(const Mint &a) {
    x *= a.x;
    x %= mod;
    return *this;
  }

  // a^n mod を計算する
  // Useaeg: Mint z = Mint(2).pow(n);
  Mint pow(long long n) const {
    if (n == 0)
      return Mint(1);
    Mint y = pow(n >> 1); // pow(n/2)を計算する
    y *= y;
    if (n % 2 == 1)
      y *= *this;
    return y;
  }

  // a^{-1} mod を計算する
  Mint modinv(const Mint &a) const { return a.pow(mod - 2); }

  Mint &operator/=(const Mint &a) {
    x *= modinv(a).x;
    x %= mod;
    return *this;
  }

  Mint operator+(Mint &a) const {
    Mint y(*this);
    y += a;
    return y;
  }

  Mint operator-(Mint &a) const {
    Mint y(*this);
    y -= a;
    return y;
  }

  Mint operator*(Mint &a) const {
    Mint y(*this);
    y *= a;
    return y;
  }

  Mint operator/(Mint &a) const {
    Mint y(*this);
    y /= a;
    return y.x;
  }

  // nCk @mod を計算する
  Mint nCk(Mint &n, const long long k) const {
    Mint y = Mint(1);
    Mint one = Mint(1);
    for (Mint i(0); (i.x) < k; i.x++) {
      y *= (n - i);
      y /= (i + one);
    }
    return y;
  }

  // nPk @mod を計算する
  Mint nPk(Mint &n, long long k) const {
    Mint y(1);
    for (Mint i(0); (i.x) < k; i.x++) {
      y *= (n - i);
    }
    return y;
  }
};

class UnionFind {
private:
  int n;
  std::vector<int> parent;

public:
  UnionFind(int num) : n(num) {
    parent = std::vector<int>(n);
    for (int i = 0; i < n; i++) {
      parent[i] = i;
    }
  }
  int GetRoot(int i) {
    int n = i;
    while (parent[n] != n) {
      n = parent[n];
    }
    parent[i] = n;
    return n;
  }

  bool SameGroup(int n, int m) {
    n = this->GetRoot(n);
    m = this->GetRoot(m);
    return n == m;
  }

  bool Unite(int n, int m) {
    n = this->GetRoot(n);
    m = this->GetRoot(m);
    if (n == m)
      return false;
    parent[n] = m;
    return true;
  }
};

int getgcd(int a, int b) {
  int max = std::max(a, b);
  int min = std::min(a, b);
  if (min == 0)
    return max;
  return getgcd(min, max % min);
}

int main() {
  int N, M;
  std::cin >> N >> M;
  std::vector<int> a(N);
  for (int i = 0; i < N; i++) {
    std::cin >> a[i];
  }
  UnionFind uf(N);
  int A, B;
  for (int i = 0; i < M; i++) {
    std::cin >> A >> B;
    uf.Unite(A, B);
  }
  std::vector<int> group[N];
  for (int i = 0; i < N; i++) {
    int root = uf.GetRoot(i);
    group[root].push_back(a[i]);
    // printf("root of %d = %d\n", i, root);
  }
  long long ans = 0;
  std::vector<int> amari;
  int groupnum = 0;
  for (int i = 0; i < N; i++) {
    if (group[i].size() > 0)
      groupnum++;
  }
  if (groupnum == 1) {
    printf("0\n");
    return 0;
  }

  for (int i = 0; i < N; i++) {
    if (group[i].size() == 0)
      continue;
    std::sort(group[i].begin(), group[i].end());
    ans += group[i][0];
    // printf("ans += %d\n", group[i][0]);
    for (int j = 1; j < group[i].size(); j++) {
      amari.push_back(group[i][j]);
      // printf("amari pushing %d\n", group[i][j]);
    }
  }
  if (amari.size() < groupnum - 2) {
    printf("Impossible\n");
  } else {
    std::sort(amari.begin(), amari.end());
    for (int j = 0; j < groupnum - 2; j++) {
      ans += amari[j];
    }
    std::cout << ans << std::endl;
  }

  return 0;
}