結果

問題 No.1779 Magical Swap
ユーザー tnakao0123tnakao0123
提出日時 2021-12-09 14:25:29
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,534 bytes
コンパイル時間 501 ms
コンパイル使用メモリ 54,444 KB
実行使用メモリ 23,720 KB
最終ジャッジ日時 2023-09-24 07:27:53
合計ジャッジ時間 2,472 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 5 ms
12,308 KB
testcase_01 AC 4 ms
12,348 KB
testcase_02 AC 7 ms
12,404 KB
testcase_03 WA -
testcase_04 AC 92 ms
22,008 KB
testcase_05 AC 36 ms
16,580 KB
testcase_06 AC 42 ms
17,088 KB
testcase_07 AC 75 ms
20,616 KB
testcase_08 AC 7 ms
12,608 KB
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 AC 37 ms
12,400 KB
testcase_16 AC 96 ms
23,680 KB
testcase_17 WA -
testcase_18 AC 5 ms
12,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

/* -*- coding: utf-8 -*-
 *
 * 1779.cc:  No.1779 Magical Swap - yukicoder
 */

#include<cstdio>
#include<set>
#include<algorithm>
 
using namespace std;

/* constant */

const int MAX_N = 100000;

/* typedef */

typedef multiset<int> msi;

template <const int MAX_N>
struct UFT {
  int links[MAX_N], ranks[MAX_N], sizes[MAX_N];
  UFT() {}

  void init(int n) {
    for (int i = 0; i < n; i++) links[i] = i, ranks[i] = sizes[i] = 1;
  }

  int root(int i) {
    int i0 = i;
    while (links[i0] != i0) i0 = links[i0];
    return (links[i] = i0);
  }

  int rank(int i) { return ranks[root(i)]; }
  int size(int i) { return sizes[root(i)]; }
  bool same(int i, int j) { return root(i) == root(j); }

  int merge(int i0, int i1) {
    int r0 = root(i0), r1 = root(i1), mr;
    if (r0 == r1) return r0;
    if (ranks[r0] == ranks[r1]) {
      links[r1] = r0;
      sizes[r0] += sizes[r1];
      ranks[r0]++;
      mr = r0;
    }
    else if (ranks[r0] > ranks[r1]) {
      links[r1] = r0;
      sizes[r0] += sizes[r1];
      mr = r0;
    }
    else {
      links[r0] = r1;
      sizes[r1] += sizes[r0];
      mr = r1;
    }
    return mr;
  }
};

/* global variables */

int as[MAX_N + 1], bs[MAX_N + 1];
msi sas[MAX_N + 1], sbs[MAX_N + 1];
UFT<MAX_N+1> uft;

/* subroutines */

void mergemsi(msi &s0, msi &s1) {
  if (s0.size() < s1.size()) swap(s0, s1);
  for (auto v: s1) s0.insert(v);
  s1.clear();
}

void printmsi(msi &s) {
  for (auto v: s) printf("%d ", v);
}

bool eqmsi(msi &s0, msi &s1) {
  //printmsi(s0), printf("<=> "), printmsi(s1), putchar('\n');
  if (s0.size() != s1.size()) return false;
  msi::iterator sit0 = s0.begin();
  msi::iterator sit1 = s1.begin();
  for (; sit0 != s0.end(); sit0++, sit1++)
    if (*sit0 != *sit1) return false;
  return true;
}

/* main */

int main() {
  int tn;
  scanf("%d", &tn);

  while (tn--) {
    int n;
    scanf("%d", &n);
    for (int i = 1; i <= n; i++) scanf("%d", as + i);
    for (int i = 1; i <= n; i++) scanf("%d", bs + i);

    for (int i = 1; i <= n; i++) {
      sas[i].clear(), sas[i].insert(as[i]);
      sbs[i].clear(), sbs[i].insert(bs[i]);
    }

    uft.init(n + 1);
    for (int u = 2; u * 2 <= n; u++)
      for (int v = u * 2; v <= n; v += u)
	if (! uft.same(u, v)) {
	  int p = uft.merge(u, v), q = (p != u) ? u : v;
	  mergemsi(sas[p], sas[q]);
	  mergemsi(sbs[p], sbs[q]);
	}

    bool ok = true;
    for (int i = 1; ok && i <= n; i++)
      if (uft.root(i) == i)
	ok = eqmsi(sas[i], sbs[i]);

    if (ok) puts("Yes");
    else puts("No");
  }
  return 0;
}
0