結果

問題 No.2948 move move rotti
コンテスト
ユーザー anago-pie
提出日時 2024-11-05 15:09:51
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 948 ms / 4,000 ms
コード長 4,953 bytes
コンパイル時間 3,758 ms
コンパイル使用メモリ 285,212 KB
実行使用メモリ 79,232 KB
最終ジャッジ日時 2024-11-05 15:10:04
合計ジャッジ時間 11,639 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define rep(i, n) for(int i=0; i<n; i++)
#define debug 0
#define YES cout << "Yes" << endl;
#define NO cout << "No" << endl;
using ll = long long;
using ld = long double;
const int mod = 998244353;
const int MOD = 1000000007;
const double pi = atan2(0, -1);
const int inf = 1 << 31 - 1;
const ll INF = 1LL << 63 - 1;
#include <time.h>
#include <chrono>

//vectorの中身を空白区切りで出力
template<typename T>
void printv(vector<T> v) {
	for (int i = 0; i < v.size(); i++) {
		cout << v[i];
		if (i < v.size() - 1) {
			cout << " ";
		}
	}
	cout << endl;
}

//vectorの中身を改行区切りで出力
template<typename T>
void print1(vector<T> v) {
	for (auto x : v) {
		cout << x << endl;
	}
}

//二次元配列を出力
template<typename T>
void printvv(vector<vector<T>> vv) {
	for (vector<T> v : vv) {
		printv(v);
	}
}

//vectorを降順にソート
template<typename T>
void rsort(vector<T>& v) {
	sort(v.begin(), v.end());
	reverse(v.begin(), v.end());
}

//昇順priority_queueを召喚
template<typename T>
struct rpriority_queue {
	priority_queue<T, vector<T>, greater<T>> pq;

	void push(T x) {
		pq.push(x);
	}

	void pop() {
		pq.pop();
	}

	T top() {
		return pq.top();
	}

	size_t size() {
		return pq.size();
	}

	bool empty() {
		return pq.empty();
	}
};

//mod mod下で逆元を算出する
//高速a^n計算(mod ver.)
ll power(ll a, ll n) {
	if (n == 0) {
		return 1;
	}
	else if (n % 2 == 0) {
		ll x = power(a, n / 2);
		x *= x;
		x %= mod;
		return x;
	}
	else {
		ll x = power(a, n - 1);
		x *= a;
		x %= mod;
		return x;
	}
}
//フェルマーの小定理を利用
ll modinv(ll p) {
	return power(p, mod - 2) % mod;
}

//Mexを求める
struct Mex {
	map<int, int> mp;
	set<int> s;
	Mex(int Max) {
		for (int i = 0; i <= Max; i++) {
			s.insert(i);
		}
	}

	int _mex = 0;
	void Input(int x) {
		mp[x]++;
		s.erase(x);
		if (_mex == x) {
			_mex = *begin(s);
		}
	}

	void Remove(int x) {
		if (mp[x] == 0) {
			cout << "Mex ERROR!: NO VALUE WILL BE REMOVED" << endl;
		}
		mp[x]--;
		if (mp[x] == 0) {
			s.insert(x);
			if (*begin(s) == x) {
				_mex = x;
			}
		}
	}

	int mex() {
		return _mex;
	}
};

//Union-Find
struct UnionFind {
	vector<int> par;
	UnionFind(int N) : par(N) {
		for (int i = 0; i < N; i++) {
			par[i] = -1;
		}
	}

	//root(x):xの根を求める関数
	int root(int x) {
		if (par[x] == -1) {
			return x;
		}
		else {
			return par[x] = root(par[x]);
		}
	}

	//isSame(x,y):xとyが同じグループならtrueを返す関数
	bool isSame(int x, int y) {
		if (root(x) == root(y)) {
			return true;
		}
		else {
			return false;
		}
	}

	//Union(x,y):xとyの根をつなげる関数
	void Union(int x, int y) {
		int X = root(x);
		int Y = root(y);
		if (X == Y) {
			return;
		}
		if (X < Y) {
			swap(X, Y);
		}
		par[X] = Y;
	}
};

//最大公約数(ユークリッドの互除法）
ll gcd(ll a, ll b) {
	if (b > a) {
		swap(a, b);
	}
	while (a % b != 0) {
		ll t = a;
		a = b;
		b = t % b;
	}
	return b;
}

//最小公倍数（gcdを定義しておく）
ll lcm(ll a, ll b) {
	ll g = gcd(a, b);
	ll x = (a / g) * b;
	return x;
}

struct SegTree {
	vector<ll> tree;
	ll n = 1;
	SegTree(int M) {
		while (n < M) {
			n *= 2;
		}
		vector<ll> t(n * 2, 0);
		swap(tree, t);
	}

	void update(ll i) {
		if (i > 0) {
			ll j = (i % 2 == 0 ? i - 1 : i + 1);
			tree[(i - 1) / 2] = tree[i] + tree[j];
			update((i - 1) / 2);
		}
	}

	void add(ll r) {
		tree[n - 1 + r]++;
		update(n - 1 + r);
	}

	ll sum(ll l, ll r, ll d, ll u, ll now) {
		if (l >= u || r <= d) {
			return 0;
		}
		else if (l <= d && r >= u) {
			return tree[now];
		}
		else {
			return sum(l, r, d, (d + u) / 2, now * 2 + 1) + sum(l, r, (d + u) / 2, u, now * 2 + 2);
		}
	}

	ll query(ll l, ll r) {
		return sum(l, r + 1, 0, n, 0);
	}
};

int main() {
	int N, M, K;
	cin >> N >> M >> K;
	vector<vector<int>> path(N);
	vector<int> start(K);
	rep(i, K) {
		int a;
		cin >> a;
		a--;
		start[i] = a;
	}
	rep(i, M) {
		int u, v;
		cin >> u >> v;
		u--; v--;
		path[u].push_back(v);
		path[v].push_back(u);
	}

	vector<vector<vector<unordered_set<bitset<15>>>>> memo(N, vector<vector<unordered_set<bitset<15>>>>(N, vector<unordered_set<bitset<15>>>(N)));
	vector<vector<unordered_set<int>>> rec(N, vector<unordered_set<int>>(N));
	vector<bitset<15>> bin(N + 1);
	rep(i, N+1) {
		bin[i] = 1<<i;
	}
	rep(k, K) {
		int& s = start[k];
		bitset<15> bs = bin[s];
		memo[s][s][0].insert(bs);
		rep(m, N) {
			rep(g, N) {
				if (!memo[s][g][m].empty()) {
					rec[g][m].insert(k);
					for (bitset<15> b : memo[s][g][m]) {
						for (int x : path[g]) {
							if (!b[x]) {
								//cout << "k=" << k << " s=" << s << " x=" << x << " m=" << m << endl;
								memo[s][x][m + 1].insert(b | bin[x]);
							}
						}
					}
				}
			}
		}
	}

	bool ok = false;
	rep(g, N) {
		rep(m, N) {
			if (rec[g][m].size() == K) {
				ok = true;
			}
		}
	}

	cout << (ok ? "Yes" : "No") << endl;
}
0