結果
| 問題 | No.3508 OR Mapping |
| コンテスト | |
| ユーザー |
 wsrtrt
|
| 提出日時 | 2026-04-18 12:37:04 |
| 言語 | C++23 (gcc 15.2.0 + boost 1.89.0) |
| 結果 |
AC
|
| 実行時間 | 1,867 ms / 2,000 ms |
| コード長 | 15,038 bytes |
| 記録 | |
| コンパイル時間 | 8,434 ms |
| コンパイル使用メモリ | 402,356 KB |
| 実行使用メモリ | 390,820 KB |
| 最終ジャッジ日時 | 2026-04-18 12:38:06 |
| 合計ジャッジ時間 | 52,435 ms |
|
ジャッジサーバーID （参考情報） |
judge2_0 / judge1_0 |
（要ログイン）
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 65 |
ソースコード
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#define INT(...) \
int __VA_ARGS__; \
IN(__VA_ARGS__)
#define LL(...) \
ll __VA_ARGS__; \
IN(__VA_ARGS__)
#define STR(...) \
string __VA_ARGS__; \
IN(__VA_ARGS__)
#define CHR(...) \
char __VA_ARGS__; \
IN(__VA_ARGS__)
#define DBL(...) \
double __VA_ARGS__; \
IN(__VA_ARGS__)
#define ll long long
#define yes cout<<"Yes"<<endl
#define no cout<<"No"<<endl
#define rep(i,a,b) for(int i=a;i<b;i++)
#define rrep(i,a,b) for(int i=a;i>=b;i--)
#define fore(i,a) for(auto &i:a)
#define all(x) (x).begin(),(x).end()
#define allr(x) (x).rbegin(),(x).rend()
#define SUM(v) accumulate(all(v), 0LL)
#define MIN(v) *min_element(all(v))
#define MAX(v) *max_element(all(v))
#define lb(c, x) distance((c).begin(), lower_bound(all(c), (x)))
#define ub(c, x) distance((c).begin(), upper_bound(all(c), (x)))
#define pii pair<int, int>
#define pll pair<long long,long long>
#define pb push_back
#define eb emplace_back
#define ff first
#define ss second
#define vi vector<int>
#define vll vector<long long>
#define vc vector<char>
#define vvi vector<vector<int>>
#define vec(type, name, ...) vector<type> name(__VA_ARGS__)
#define VEC(type, name, size) \
vector<type> name(size); \
IN(name)
int scan() { return getchar(); }
void scan(int &a) { cin >> a; }
void scan(long long &a) { cin >> a; }
void scan(char &a) { cin >> a; }
void scan(double &a) { cin >> a; }
void scan(string &a) { cin >> a; }
template <class T, class S> void scan(pair<T, S> &p) { scan(p.first), scan(p.second); }
template <class T> void scan(vector<T> &);
template <class T> void scan(vector<T> &a) {
for(auto &i : a) scan(i);
}
template <class T> void scan(T &a) { cin >> a; }
void IN() {}
template <class Head, class... Tail> void IN(Head &head, Tail &...tail) {
scan(head);
IN(tail...);
}
template <class T> void print(const T &a) { cout << a; }
void OUT() { cout << endl; }
template <class Head, class... Tail> void OUT(const Head &head, const Tail &...tail) {
print(head);
if(sizeof...(tail)) cout << ' ';
OUT(tail...);
}
#define vv(type, name, h, ...) vector<vector<type>> name(h, vector<type>(__VA_ARGS__))
#define VV(type, name, h, w) \
vector<vector<type>> name(h, vector<type>(w)); \
IN(name)
#define vvv(type, name, h, w, ...) vector<vector<vector<type>>> name(h, vector<vector<type>>(w, vector<type>(__VA_ARGS__)))
#define vvvv(type, name, a, b, c, ...) \
vector<vector<vector<vector<type>>>> name(a, vector<vector<vector<type>>>(b, vector<vector<type>>(c, vector<type>(__VA_ARGS__))))
template<typename T> using min_priority_queue = priority_queue<T, vector<T>, greater<T>>;
template <class T> pair<T, T> operator-(const pair<T, T> &x, const pair<T, T> &y) { return pair<T, T>(x.ff - y.ff, x.ss - y.ss); }
template <class T> pair<T, T> operator+(const pair<T, T> &x, const pair<T, T> &y) { return pair<T, T>(x.ff + y.ff, x.ss + y.ss); }
template <class T> pair<T, T> operator&(const pair<T, T> &l, const pair<T, T> &r) { return pair<T, T>(max(l.ff, r.ff), min(l.ss, r.ss)); }
template <class T> vector<T> &operator--(vector<T> &v) {
fore(e, v) e--;
return v;
}
template <class T> vector<T> operator--(vector<T> &v, int) {
auto res = v;
fore(e, v) e--;
return res;
}
template <class T> vector<T> &operator++(vector<T> &v) {
fore(e, v) e++;
return v;
}
template <class T> vector<T> operator++(vector<T> &v, int) {
auto res = v;
fore(e, v) e++;
return res;
}
template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; }
template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; }
#define UNIQUE(x) sort(all(x)), x.erase(unique(all(x)), x.end())
//座標圧縮
template <typename T> void zip(vector<T> &x) {
vector<T> y = x;
UNIQUE(y);
for(int i = 0; i < x.size(); ++i) { x[i] = lb(y, x[i]); }
}
template <class T> T ceil(T x, T y) {
assert(y >= 1);
return (x > 0 ? (x + y - 1) / y : x / y);
}
template <class T> T floor(T x, T y) {
assert(y >= 1);
return (x > 0 ? x / y : (x + y - 1) / y);
}
long long POW(long long x, int n) {
long long res = 1LL;
for(; n; n >>= 1, x *= x)
if(n & 1) res *= x;
return res;
}
//0^n=0
long long modpow(long long a, long long n, long long mod) {
a%=mod;
assert(a!=0||n!=0);
if(a==0)return 0;
long long res = 1;
while (n > 0) {
if (n & 1) res = res * a % mod;
a = a * a % mod;
n >>= 1;
}
return res;
}
//return 0<=a&&a<h&&0<=b&&b<w;
inline bool ingrid(ll a,ll b,ll h,ll w){return 0<=a&&a<h&&0<=b&&b<w;}
//return 0<=a&&a<n;
inline bool inl(int a,int n){return 0<=a&&a<n;}
// bit 演算系
ll pow2(int i) { return 1LL << i; }
int topbit(signed t) { return t == 0 ? -1 : 31 - __builtin_clz(t); }
int topbit(ll t) { return t == 0 ? -1 : 63 - __builtin_clzll(t); }
int lowbit(signed a) { return a == 0 ? 32 : __builtin_ctz(a); }
int lowbit(ll a) { return a == 0 ? 64 : __builtin_ctzll(a); }
// int allbit(int n) { return (1 << n) - 1; }
ll allbit(ll n) { return (1LL << n) - 1; }
int popcount(signed t) { return __builtin_popcount(t); }
int popcount(ll t) { return __builtin_popcountll(t); }
bool ispow2(int i) { return i && (i & -i) == i; }
int in() {
int x;
cin >> x;
return x;
}
ll lin() {
unsigned long long x;
cin >> x;
return x;
}
long long sqrtll(long long x) {
assert(x >= 0);
long long rev = sqrt(x);
while(rev * rev > x) --rev;
while((rev+1) * (rev+1)<=x) ++rev;
return rev;
}
int logN(long long n){
int ret=1;
while((1LL<<ret)<n)ret++;
return ret;
}
const double PI=3.1415926535897932384626433832795028841971;
const ll MOD998 = 998244353;
const int INFI = numeric_limits<int>::max() / 2; const long long INFL = numeric_limits<long long>::max() / 2;
#define inf INFINITY
template<class T>
void debug(vector<T> a){
rep(i,0,(int)a.size()){
cout<<a[i]<<' ';
}
cout<<endl;
return;
}
bool palindrome(const string& s){
return equal(all(s),s.rbegin());
}
template <std::uint_fast64_t Modulus> class modint {
using u64 = std::uint_fast64_t;
public:
u64 a;
constexpr modint(const u64 x = 0) noexcept : a(x % Modulus) {}
constexpr u64 &val() noexcept { return a; }
constexpr const u64 &val() const noexcept { return a; }
constexpr modint operator+(const modint rhs) const noexcept {
return modint(*this) += rhs;
}
constexpr modint operator-(const modint rhs) const noexcept {
return modint(*this) -= rhs;
}
constexpr modint operator*(const modint rhs) const noexcept {
return modint(*this) *= rhs;
}
constexpr modint operator/(const modint rhs) const noexcept {
return modint(*this) /= rhs;
}
constexpr modint &operator+=(const modint rhs) noexcept {
a += rhs.a;
if (a >= Modulus) {
a -= Modulus;
}
return *this;
}
constexpr modint &operator-=(const modint rhs) noexcept {
if (a < rhs.a) {
a += Modulus;
}
a -= rhs.a;
return *this;
}
constexpr modint &operator*=(const modint rhs) noexcept {
a = a * rhs.a % Modulus;
return *this;
}
constexpr modint &operator/=(modint rhs) noexcept {
u64 exp = Modulus - 2;
while (exp) {
if (exp % 2) {
*this *= rhs;
}
rhs *= rhs;
exp /= 2;
}
return *this;
}
friend bool operator==(const modint& a,const modint& b) { return a.val()==b.val(); }
friend bool operator!=(const modint& a,const modint& b) { return a.val()!=b.val(); }
};
using mint9=modint<998244353>;
using mint1=modint<1000000007>;
//costを指定しないと重みなし辺になります
struct Edge{
int from,to;
ll cost;
Edge()=default;
Edge(int from,int to,ll cost=1):from(from),to(to),cost(cost){}
operator int() const {return to;}
};
constexpr pii dx4[4] = {pii{-1, 0},pii{0, -1}, pii{0, 1}, pii{1, 0} };
constexpr pii dx8[8] = {{1, 0}, {1, 1}, {0, 1}, {-1, 1}, {-1, 0}, {-1, -1}, {0, -1}, {1, -1}};
constexpr pii dx[100]={{0,0},{-1,0},{0,-1},{-1,-1}};
template <typename G>
struct StronglyConnectedComponents {
const G &g;
vector<vector<int>> ng, rg; // normal-graphとreversed-graph
vector<int> compo, order, used; // compo…頂点がどの成分に属しているか、order…帰りがけ順で何番目か
StronglyConnectedComponents(G &g) : g(g), ng(g.size()), rg(g.size()), compo(g.size(), -1), used(g.size()) {
for (int i = 0; i < g.size(); i++)
for (auto t : g[i]) {
ng[i].emplace_back((int)t);
rg[(int)t].emplace_back(i); // 逆辺
}
}
// []で聞かれたときは成分番号を返す
int operator[](int k) { return compo[k]; }
void dfs(int now) {
if (used[now]) return;
used[now] = true;
for (auto to : ng[now]) dfs(to);
order.emplace_back(now);
}
void rdfs(int now, int count) {
// used代わりに結果を入れていく
if (compo[now] != -1) return;
compo[now] = count;
for (auto to : rg[now]) rdfs(to, count);
}
void build(vector<vector<int>> &ret) {
// 普通にDFS、辺の向きを変えてもう一度DFS
for (int i = 0; i < ng.size(); i++) dfs(i);
reverse(order.begin(), order.end());
int group = 0;
for (auto i : order)
if (compo[i] == -1) rdfs(i, group), group++;
// 縮めたグラフを構築する
ret.resize(group);
for (int i = 0; i < g.size(); i++) { // 全ての辺について
for (auto &to : g[i]) {
int s = compo[i], t = compo[to];
if (s != t) ret[s].emplace_back(t);
}
}
}
};
vector<int> topological_sort(vector<vector<int>> &g){
vector<int> ret;
queue<int> q;
int n=g.size();
vector<int> indeg(n);
rep(i,0,n){
fore(j,g[i]){
indeg[j]++;
}
}
rep(i,0,n){
if(indeg[i]==0){
q.push(i);
}
}
while(!q.empty()){
if(q.size()!=1){
cout<<"No"<<endl;
exit(0);
}
int tmp=q.front();
q.pop();
fore(i,g[tmp]){
indeg[i]--;
if(indeg[i]==0)q.push(i);
}
ret.push_back(tmp);
}
return ret;
}
int main(){
cin.tie(0)->sync_with_stdio(0);
int n,m,k;cin>>n>>m>>k;
vector<vector<int>> g(2*n);
vector<vector<int>> h(n);
vector<int> u(m),v(m);
rep(i,0,m){
cin>>u[i]>>v[i];u[i]--;v[i]--;
g[2*u[i]].pb(2*v[i]+1);
g[2*u[i]+1].pb(2*v[i]);
h[u[i]].pb(v[i]);
}
if(n==1){
cout<<"No"<<endl;
return 0;
}
StronglyConnectedComponents scc(g);
vector<vector<int>> tmp;
scc.build(tmp);
//頂点数1の連結成分があるかどうか
map<int,vector<int>> mp;
rep(i,0,2*n){
mp[scc[i]].push_back(i);
}
//全部の頂点に到達可能である必要がある
{
vector<int> seen(n);
seen[0]=1;
auto dfs=[&](auto&& self,int pos)->void {
fore(nxt,h[pos]){
if(seen[nxt])continue;
seen[nxt]=1;
self(self,nxt);
}
};
dfs(dfs,0);
rep(i,0,n){
if(seen[i]==0){
cout<<"No"<<endl;
return 0;
}
}
}
//奇数回で戻ってこれる or 例外
rep(i,0,n){
if(scc[2*i]!=scc[2*i+1]){
if(i!=0&&mp[scc[2*i]].size()==1){
continue;
}
cout<<"No"<<endl;
return 0;
}else continue;
}
{
//hをSCCして縮約する
map<int,int> cnt;
StronglyConnectedComponents scc_h(h);
vector<vector<int>> shr;
scc_h.build(shr);
rep(i,0,n){
cnt[scc_h[i]]++;
}
int flag{};
vector<int> ts=topological_sort(shr);
fore(pos,ts){
if(cnt[pos]==1){
if(flag==0){
cout<<"No"<<endl;
return 0;
}
flag=0;
}else{
flag=1;
}
}
}
cout<<"Yes"<<endl;
return 0;
}
/*
強連結成分は高々二つまで
奇数回で戻ってこれる
頂点2*iと2*i+1
頂点倍化して強連結成分分解
基本的には今日連結成分が2個以下かつ、奇数回で戻ってこれるならok
ただし、頂点1を含まない強連結成分の頂点数が1ならこれもok
普通に違う
「成分内の頂点は奇数回で自分に戻ってこれる」or「頂点1ではなく、成分の頂点数が1」
をすべての連結成分が満たし、とぽそが一意かつ、なんかいろいろならおｋ
7 8 3
1 2
2 3
3 1
1 4
4 5
5 6
6 7
7 5
*/