#include #pragma GCC optimize("Ofast") #define _GLIBCXX_DEBUG using namespace std; using std::cout; using std::cin; using std::endl; using ll=long long; using ld=long double; ll ILL=1167167167167167167; const int INF=2100000000; const ll mod=1e9+7; #define rep(i,a) for (ll i=0;i using _pq = priority_queue, greater>; template ll LB(vector &v,T a){return lower_bound(v.begin(),v.end(),a)-v.begin();} template ll UB(vector &v,T a){return upper_bound(v.begin(),v.end(),a)-v.begin();} template bool chmin(T &a,const T &b){if(a>b){a=b;return 1;}else return 0;} template bool chmax(T &a,const T &b){if(a void So(vector &v) {sort(v.begin(),v.end());} template void Sore(vector &v) {sort(v.begin(),v.end(),[](T x,T y){return x>y;});} void yneos(bool a){if(a) cout<<"Yes\n"; else cout<<"No\n";} namespace po167{ struct SCC{ std::vector> G; std::vector> G_rev; int size; SCC(int n):size(n),G(n),G_rev(n){} void add_edge(int from,int to){ assert(0<=to&&to &order,std::vector &seen,int ind,int &now){ if(seen[ind]!=0) return; seen[ind]=-1; for(auto x:G[ind]){ dfs1(order,seen,x,now); } order[now]=ind; now++; } void dfs2(std::vector &seen,int ind,int rank){ if(seen[ind]!=-1) return; seen[ind]=rank; for(auto x:G_rev[ind]){ dfs2(seen,x,rank); } } std::vector> tp_sort(){ std::vector seen(size,0),order(size); int now=0; for(int i=0;i=0;i--){ if(seen[order[i]]==-1){ dfs2(seen,order[i],now); now++; } } std::vector> ans(now); for(int i=0;i for_two_sat(){ std::vector seen(size,0),order(size); int now=0; for(int i=0;i=0;i--){ if(seen[order[i]]==-1){ dfs2(seen,order[i],now); now++; } } return seen; } }; struct two_SAT{ po167::SCC graph; std::vector ans; bool exists=false; int size; two_SAT(int n):size(n),ans(n),graph(2*n){} void add_clause(int i,bool f,int j,bool g){ graph.add_edge(i+size*(!f),j+size*(g)); graph.add_edge(j+size*(!g),i+size*(f)); } bool build(){ auto rank=graph.for_two_sat(); exists=true; for(int i=0;i>N>>M; int A=0,B=0; SCC G(N); rep(i,M){ int a,b; cin>>a>>b; a--,b--; G.add_edge(a,b); } auto p=G.tp_sort(); vector rev(N); rep(i,p.size()){ for(auto y:p[i]) rev[y]=i; } rep(i,p.size()){ int c=1,d=1; for(auto y:p[i]){ for(auto x:G.G[y]){ if(rev[y]