Problem p03354 (C++)

Buggy submission (status: Runtime Error):
/*
                  友利奈緒ぱわーでACしてくれ！！！！！！！！！！！！
                  Charlotteまだ見てない人は見ような！！！！！！！

                  ／　 /:／　　　　　　　　　　　 　 ＼:ヽ　＼
  　　　　　　　 /　　／　　　　　　 　 　 　 　 　 　 　 ＼　　ヽ
  　　　　　＼__L 彡　　　　　　　　　　　　　　　　　　　 ヽ　　',
  　　　　　　　ﾌ /　　　/　　/　 /　　　,　　 　 　 　 ヽ　　', 　ｊ
  　　　　　　　ﾚ　　　 /　　/ 　/　 　 　 　 　 　 　 　 ', 　ヽ　j
  　　　 　 　 // 　 　 !　/_l＿/__ノ 　　{ ＼}　　　　 　 ',　　',/
  　　　　　　i ﾊ 　 　 {　从j　ﾊ　i　　　 ',　 ト-､}　 i i　 }　　jj
  　　　　　　( ﾊ　　　(　 川　ハ ﾊ　 !　 '　ﾊ　 ﾊ　j j　,' 　 八
  　　　 ＜＿,ｨ∧　　　斗芹ﾃミxハ ﾊ　 ﾚ　} /__ﾚﾚ　/ 　 ∧ ＼　ﾉ
  　　　 　 　 ∨　＼ 　ゝ　辷:ソ　　 　)　芹ﾚ心ヾレ′　/ ト--　´
  　　　＼ｰ‐'　ﾉ　　 ＼ゝ　　　　　　 　 　 ゞ:_ソ "/ 　／ 　ヽ
  　　　　 　 フ　　　　i八 " "　　　,　　　　　　　ム彡　　　　＼
  　　　　 　/　　/　　j　 ト　　　　　　 　 　 " "イト＜　 　 　 ＼＼
  　　　　 /　　/　　　j　/／ヽ､ 　 ∩　　 　 イ　{ {　　￣ フフへ ＼＼
  　　　　(　／　　　 ,/ / 　　i　＞――＜ニニニﾆ┐　 〃／: : ヽ　ヽヽ
  　　　　／　　ィT´/ /┌―　￣￣ /::: , ,)　 　 　 〃／／: : : : : :}　ヽヽ)
  　　　/ /　 ∧ヾi┌― { {￣　　　 ﾉ:::トﾟ<　 　 　 ∥//: : : : : : : : i　 ﾉ ﾉ
  . 　　( 人　 {:ヽヾi∨　∧V　　　 /:::/　､ヽ　　　 ∥/: : : : : : : : : :}　/
  　　　V　( ∧: :＼'∨　∧V　　 ﾉ:::/　 ∧ 〉＿__∥: : : : : :／: : : :ﾚ
  　　 　 　 /∧: : : : ∨　∧V┬ｲ:::ﾉ　 〈 TT　|　|{{: : : : ／: : : : : ∧　　 ﾉ
  　　　 　 (　ハ: : : : :∨'T∧Vi　i (　 　 V　!　!　{{: : ／ : : : : : :
  /　ｰ　ノ 　　　　　V　ハ: : : : :∨ﾍ∧V　i ゝゝ　 i／⌒＼{／ : : : : : : : /ｰ― ´
  　　　　　 ∨　ヽ: : : : :∨ﾍ∧∨i　〉 〉 /　/二　 ): : : : : : : : :/

          (https://seesaawiki.jp/asciiart/d/Charlotte より)
*/
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using ull = unsigned long long;
// #define int ll
// #define DEBUG 42
inline void nprint() {}
template <class Head, class... Tail>
inline void nprint(Head &&head, Tail &&...tail) {
  cout << head << endl;
  nprint(move(tail)...);
}
#ifdef DEBUG
#define eprint(...) nprint(__VA_ARGS__)
#else
#define eprint(...)                                                            \
  if (0 == 1)                                                                  \
    cout << 1 << endl;
#endif
#define Yes(a) cout << (a ? "Yes" : "No") << endl
#define YES(a) cout << (a ? "YES" : "NO") << endl
#define POSSIBLE(a) cout << (a ? "POSSIBLE" : "IMPOSSIBLE") << endl
using cmp = complex<double>;
using vb = vector<bool>;
using vvb = vector<vb>;
using vi = vector<int>;
using vvi = vector<vi>;
using vl = vector<ll>;
using vvl = vector<vl>;
template <class T> using V = vector<T>;
template <class T> using VV = vector<V<T>>;
#define fi first
#define se second
#define maxs(x, y) (x = max(x, y))
#define mins(x, y) (x = min(x, y))
using pii = pair<int, int>;
using pll = pair<ll, ll>;
#define FOR(i, a, b) for (ll i = (a); i < (ll)(b); ++i)
#define REP(i, n) FOR(i, 0, n)
#define FORS(i, a, b) FOR(i, a, b + 1)
#define REPS(i, n) REP(i, n + 1)
#define RFOR(i, a, b) for (ll i = (ll)(b)-1; i >= a; --i)
#define RREP(i, n) RFOR(i, 0, n)
#define RREPS(i, n) RREP(i, n + 1)
#define RFORS(i, a, b) RFOR(i, a, b + 1)
#define ALL(obj) (obj).begin(), (obj).end()
#define RALL(obj) (obj).rbegin(), (obj).rend()
#define PERM(c)                                                                \
  sort(ALL(c));                                                                \
  for (bool cp = true; cp; cp = next_permutation(ALL(c)))
#define eb(val) emplace_back(val)
const double PI = acos(-1), EPS = 1e-10;
constexpr ll MOD = 1E9 + 7;
constexpr int dx[] = {1, 0, -1, 0};
constexpr int dy[] = {0, 1, 0, -1};
template <typename T1, typename T2>
ostream &operator<<(ostream &s, const pair<T1, T2> &p) {
  return s << "(" << p.first << ", " << p.second << ")";
}
template <class T> istream &operator>>(istream &is, vector<T> &st) {
  for (size_t i = 0; i < st.size(); ++i)
    is >> st[i];
  return is;
}
template <class T> istream &operator>>(istream &is, vector<vector<T>> &st) {
  for (size_t i = 0; i < st.size(); ++i)
    is >> st[i];
  return is;
}
template <class T> ostream &operator<<(ostream &os, const vector<T> &st) {
  for (size_t i = 0; i < st.size(); ++i) {
    if (i == st.size() - 1)
      os << st[i];
    else
      os << st[i] << " ";
  }
  return os;
}
template <class T>
ostream &operator<<(ostream &os, const vector<vector<T>> &st) {
  for (size_t i = 0; i < st.size(); ++i) {
    os << st[i];
    if (i != st.size() - 1)
      os << endl;
  }
  return os;
}
struct union_find {
  vector<long long>
      data; // そのノードがrootであれば-（集合のサイズ）を表し、そうでなければ親ノードidを示す。
  vector<long long> diff_weight;
  union_find(long long size, long long abel_sum_unity = 0)
      : data(size, -1), diff_weight(size, abel_sum_unity) {
  } // コンストラクタ,size個のノードを作成し、-1で初期化
  // 根を探す,経路圧縮
  inline long long root(long long a) {
    if (data[a] < 0)
      return a;
    long long r = root(data[a]);
    diff_weight[a] += diff_weight[data[a]];
    return data[a] = r;
    // return data[a] < 0 ? a : data[a]=root(data[a]),diff_weight[a] +=
    // diff_weight[data[a]];
  }
  // a,bが同じ集合に属するか調べる
  inline bool same(long long a, long long b) { return root(a) == root(b); }
  // 小さい方の木の根に大きい方の根を上書きする
  inline bool unite(long long a, long long b, long long w = 1) {
    w += weight(a) - weight(b);
    a = root(a), b = root(b);
    if (a == b)
      return false; // 根が同じなら終了,uniteが失敗したことを出力
    if (-data[a] < -data[b])
      swap(a, b), w = -w; //-data[a]≥-data[b]となるようにswap
    // 集合サイズの大きい方に小さい方をつける。
    data[a] += data[b]; // 集合サイズを統合
    data[b] = a;        // bの親をaにする。
    diff_weight[b] = w;
    return true; // unite成功を出力
  }
  // aが属する集合のサイズを出力
  inline long long size(long long a) { return -data[root(a)]; }
  // root(a)からaまでの距離を返す。
  inline long long weight(long long a) {
    root(a);
    return diff_weight[a];
  }
  inline long long diff(long long a, long long b) {
    return weight(b) - weight(a);
  }
};
signed main() {
  cin.tie(0);
  ios::sync_with_stdio(false);
  cout << fixed << setprecision(10);
  int n, m;
  cin >> n >> m;
  vi p(n);
  cin >> p;
  REP(i, n) p[i]--;
  union_find uf(n);
  REP(i, m) {
    int x, y;
    cin >> x >> y;
    x--, y--;
    uf.unite(x, y);
  }
  int ans = 0;
  REP(i, n) {
    if (uf.same(i, p[i] - 1))
      ++ans;
  }
  // vvi gr;
  // vector<set<int>> sts;
  // map<int,int> memo;
  // int id = 0;
  // REP(i,n){
  //   if(memo.find(uf.root(i))==memo.end()){
  //     memo[uf.root(i)] = id++;
  //     gr.eb(vi());
  //     sts.eb(set<int>());
  //     gr.back().eb(i);
  //     sts.back().insert(p[i]);
  //   }else{
  //     gr[memo[uf.root(i)]].eb(i);
  //     sts[memo[uf.root(i)]].insert(p[i]);
  //   }
  // }
  // int ans = 0;
  // REP(i,gr.size())REP(j,gr[i].size()){
  //   if(sts[i].find(gr[i][j])!=sts[i].end()) ++ans;
  // }
  cout << ans << endl;
}


Fixed submission:
/*
                  友利奈緒ぱわーでACしてくれ！！！！！！！！！！！！
                  Charlotteまだ見てない人は見ような！！！！！！！

                  ／　 /:／　　　　　　　　　　　 　 ＼:ヽ　＼
  　　　　　　　 /　　／　　　　　　 　 　 　 　 　 　 　 ＼　　ヽ
  　　　　　＼__L 彡　　　　　　　　　　　　　　　　　　　 ヽ　　',
  　　　　　　　ﾌ /　　　/　　/　 /　　　,　　 　 　 　 ヽ　　', 　ｊ
  　　　　　　　ﾚ　　　 /　　/ 　/　 　 　 　 　 　 　 　 ', 　ヽ　j
  　　　 　 　 // 　 　 !　/_l＿/__ノ 　　{ ＼}　　　　 　 ',　　',/
  　　　　　　i ﾊ 　 　 {　从j　ﾊ　i　　　 ',　 ト-､}　 i i　 }　　jj
  　　　　　　( ﾊ　　　(　 川　ハ ﾊ　 !　 '　ﾊ　 ﾊ　j j　,' 　 八
  　　　 ＜＿,ｨ∧　　　斗芹ﾃミxハ ﾊ　 ﾚ　} /__ﾚﾚ　/ 　 ∧ ＼　ﾉ
  　　　 　 　 ∨　＼ 　ゝ　辷:ソ　　 　)　芹ﾚ心ヾレ′　/ ト--　´
  　　　＼ｰ‐'　ﾉ　　 ＼ゝ　　　　　　 　 　 ゞ:_ソ "/ 　／ 　ヽ
  　　　　 　 フ　　　　i八 " "　　　,　　　　　　　ム彡　　　　＼
  　　　　 　/　　/　　j　 ト　　　　　　 　 　 " "イト＜　 　 　 ＼＼
  　　　　 /　　/　　　j　/／ヽ､ 　 ∩　　 　 イ　{ {　　￣ フフへ ＼＼
  　　　　(　／　　　 ,/ / 　　i　＞――＜ニニニﾆ┐　 〃／: : ヽ　ヽヽ
  　　　　／　　ィT´/ /┌―　￣￣ /::: , ,)　 　 　 〃／／: : : : : :}　ヽヽ)
  　　　/ /　 ∧ヾi┌― { {￣　　　 ﾉ:::トﾟ<　 　 　 ∥//: : : : : : : : i　 ﾉ ﾉ
  . 　　( 人　 {:ヽヾi∨　∧V　　　 /:::/　､ヽ　　　 ∥/: : : : : : : : : :}　/
  　　　V　( ∧: :＼'∨　∧V　　 ﾉ:::/　 ∧ 〉＿__∥: : : : : :／: : : :ﾚ
  　　 　 　 /∧: : : : ∨　∧V┬ｲ:::ﾉ　 〈 TT　|　|{{: : : : ／: : : : : ∧　　 ﾉ
  　　　 　 (　ハ: : : : :∨'T∧Vi　i (　 　 V　!　!　{{: : ／ : : : : : :
  /　ｰ　ノ 　　　　　V　ハ: : : : :∨ﾍ∧V　i ゝゝ　 i／⌒＼{／ : : : : : : : /ｰ― ´
  　　　　　 ∨　ヽ: : : : :∨ﾍ∧∨i　〉 〉 /　/二　 ): : : : : : : : :/

          (https://seesaawiki.jp/asciiart/d/Charlotte より)
*/
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using ull = unsigned long long;
// #define int ll
// #define DEBUG 42
inline void nprint() {}
template <class Head, class... Tail>
inline void nprint(Head &&head, Tail &&...tail) {
  cout << head << endl;
  nprint(move(tail)...);
}
#ifdef DEBUG
#define eprint(...) nprint(__VA_ARGS__)
#else
#define eprint(...)                                                            \
  if (0 == 1)                                                                  \
    cout << 1 << endl;
#endif
#define Yes(a) cout << (a ? "Yes" : "No") << endl
#define YES(a) cout << (a ? "YES" : "NO") << endl
#define POSSIBLE(a) cout << (a ? "POSSIBLE" : "IMPOSSIBLE") << endl
using cmp = complex<double>;
using vb = vector<bool>;
using vvb = vector<vb>;
using vi = vector<int>;
using vvi = vector<vi>;
using vl = vector<ll>;
using vvl = vector<vl>;
template <class T> using V = vector<T>;
template <class T> using VV = vector<V<T>>;
#define fi first
#define se second
#define maxs(x, y) (x = max(x, y))
#define mins(x, y) (x = min(x, y))
using pii = pair<int, int>;
using pll = pair<ll, ll>;
#define FOR(i, a, b) for (ll i = (a); i < (ll)(b); ++i)
#define REP(i, n) FOR(i, 0, n)
#define FORS(i, a, b) FOR(i, a, b + 1)
#define REPS(i, n) REP(i, n + 1)
#define RFOR(i, a, b) for (ll i = (ll)(b)-1; i >= a; --i)
#define RREP(i, n) RFOR(i, 0, n)
#define RREPS(i, n) RREP(i, n + 1)
#define RFORS(i, a, b) RFOR(i, a, b + 1)
#define ALL(obj) (obj).begin(), (obj).end()
#define RALL(obj) (obj).rbegin(), (obj).rend()
#define PERM(c)                                                                \
  sort(ALL(c));                                                                \
  for (bool cp = true; cp; cp = next_permutation(ALL(c)))
#define eb(val) emplace_back(val)
const double PI = acos(-1), EPS = 1e-10;
constexpr ll MOD = 1E9 + 7;
constexpr int dx[] = {1, 0, -1, 0};
constexpr int dy[] = {0, 1, 0, -1};
template <typename T1, typename T2>
ostream &operator<<(ostream &s, const pair<T1, T2> &p) {
  return s << "(" << p.first << ", " << p.second << ")";
}
template <class T> istream &operator>>(istream &is, vector<T> &st) {
  for (size_t i = 0; i < st.size(); ++i)
    is >> st[i];
  return is;
}
template <class T> istream &operator>>(istream &is, vector<vector<T>> &st) {
  for (size_t i = 0; i < st.size(); ++i)
    is >> st[i];
  return is;
}
template <class T> ostream &operator<<(ostream &os, const vector<T> &st) {
  for (size_t i = 0; i < st.size(); ++i) {
    if (i == st.size() - 1)
      os << st[i];
    else
      os << st[i] << " ";
  }
  return os;
}
template <class T>
ostream &operator<<(ostream &os, const vector<vector<T>> &st) {
  for (size_t i = 0; i < st.size(); ++i) {
    os << st[i];
    if (i != st.size() - 1)
      os << endl;
  }
  return os;
}
struct union_find {
  vector<long long>
      data; // そのノードがrootであれば-（集合のサイズ）を表し、そうでなければ親ノードidを示す。
  vector<long long> diff_weight;
  union_find(long long size, long long abel_sum_unity = 0)
      : data(size, -1), diff_weight(size, abel_sum_unity) {
  } // コンストラクタ,size個のノードを作成し、-1で初期化
  // 根を探す,経路圧縮
  inline long long root(long long a) {
    if (data[a] < 0)
      return a;
    long long r = root(data[a]);
    diff_weight[a] += diff_weight[data[a]];
    return data[a] = r;
    // return data[a] < 0 ? a : data[a]=root(data[a]),diff_weight[a] +=
    // diff_weight[data[a]];
  }
  // a,bが同じ集合に属するか調べる
  inline bool same(long long a, long long b) { return root(a) == root(b); }
  // 小さい方の木の根に大きい方の根を上書きする
  inline bool unite(long long a, long long b, long long w = 1) {
    w += weight(a) - weight(b);
    a = root(a), b = root(b);
    if (a == b)
      return false; // 根が同じなら終了,uniteが失敗したことを出力
    if (-data[a] < -data[b])
      swap(a, b), w = -w; //-data[a]≥-data[b]となるようにswap
    // 集合サイズの大きい方に小さい方をつける。
    data[a] += data[b]; // 集合サイズを統合
    data[b] = a;        // bの親をaにする。
    diff_weight[b] = w;
    return true; // unite成功を出力
  }
  // aが属する集合のサイズを出力
  inline long long size(long long a) { return -data[root(a)]; }
  // root(a)からaまでの距離を返す。
  inline long long weight(long long a) {
    root(a);
    return diff_weight[a];
  }
  inline long long diff(long long a, long long b) {
    return weight(b) - weight(a);
  }
};
signed main() {
  cin.tie(0);
  ios::sync_with_stdio(false);
  cout << fixed << setprecision(10);
  int n, m;
  cin >> n >> m;
  vi p(n);
  cin >> p;
  REP(i, n) p[i]--;
  union_find uf(n);
  REP(i, m) {
    int x, y;
    cin >> x >> y;
    x--, y--;
    uf.unite(x, y);
  }
  int ans = 0;
  REP(i, n) {
    if (uf.same(i, p[i]))
      ++ans;
  }
  // vvi gr;
  // vector<set<int>> sts;
  // map<int,int> memo;
  // int id = 0;
  // REP(i,n){
  //   if(memo.find(uf.root(i))==memo.end()){
  //     memo[uf.root(i)] = id++;
  //     gr.eb(vi());
  //     sts.eb(set<int>());
  //     gr.back().eb(i);
  //     sts.back().insert(p[i]);
  //   }else{
  //     gr[memo[uf.root(i)]].eb(i);
  //     sts[memo[uf.root(i)]].insert(p[i]);
  //   }
  // }
  // int ans = 0;
  // REP(i,gr.size())REP(j,gr[i].size()){
  //   if(sts[i].find(gr[i][j])!=sts[i].end()) ++ans;
  // }
  cout << ans << endl;
}