#include // #include //using namespace atcoder; // tabaicho see https://boostjp.github.io/tips/multiprec-int.html //#include // using namespace boost::multiprecision; // cpp_int // int128_t // int256_t // int512_t // int1024_t // uint128_t // uint256_t // uint512_t // uint1024_t #define int long long #define inf 1000000007 // #define inf 998244353 #define pa pair #define ppa pair #define ll long long #define PI 3.14159265358979323846 #define mp make_pair #define pb push_back #define EPS (1e-8) using namespace std; int dx[8]={0,1,0,-1,1,1,-1,-1}; int dy[8]={1,0,-1,0,-1,1,1,-1}; class pa3{ public: int x; int y,z; pa3(int x=0,int y=0,int z=0):x(x),y(y),z(z) {} bool operator < (const pa3 &p) const{ if(x!=p.x) return x (const pa3 &p) const{ if(x!=p.x) return x>p.x; if(y!=p.y) return y>p.y; return z>p.z; //return x != p.x ? x (const pa4 &p) const{ if(x!=p.x) return x>p.x; if(y!=p.y) return y>p.y; if(z!=p.z)return z>p.z; return w>p.w; //return x != p.x ? x (const pa2 &p) const{ return x != p.x ? xb) return Gcd(b,v); if(v==b) return b; if(b%v==0) return v; return Gcd(v,b%v); } int extgcd(int a, int b, int &x, int &y) { if (b == 0) { x = 1; y = 0; return a; } int d = extgcd(b, a%b, y, x); y -= a/b * x; return d; } pa operator+(const pa & l,const pa & r) { return {l.first+r.first,l.second+r.second}; } pa operator-(const pa & l,const pa & r) { return {l.first-r.first,l.second-r.second}; } ostream& operator<<(ostream& os, const vector& VEC){ for(auto v:VEC)os<& VEC){ for(auto v:VEC)os<& VEC){ for(auto v:VEC){ os<nn || nn<0) return 0; int r=pr[nn]*inv[rr]; r%=mod; r*=inv[nn-rr]; r%=mod; return r; } void gya(int ert){ pr[0]=1; for(int i=1;i<=ert;i++){ pr[i]=((ll)pr[i-1]*i)%mod; } inv[ert]=beki(pr[ert],mod-2,mod); for(int i=ert-1;i>=0;i--){ inv[i]=(ll)inv[i+1]*(i+1)%mod; } } // cin.tie(0); // ios::sync_with_stdio(false); //priority_queue,greater> pq; //sort(ve.begin(),ve.end(),greater()); // mt19937(clock_per_sec); // mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count()) ; struct DM_decomposition{ private: public: int V; const int MAX=1000000000; // < N vector> gr[2]; vectormm_use,group,vis; vectorssc_jun; /* void dfs_nibu(int r,int d){ if(nibu[r]>=0)return; nibu[r]=d; for(auto v:gr[0][r])dfs_nibu(v,1-d); } */ void dfs_toutatu(int r,int muki,int d){ if(group[r]>=0){ assert(group[r]==d); return; } //cout<,vector>> shori(int N,vectornibu,vectoreda,vectorMM){ V=N; gr[0].resize(N,{}); gr[1].resize(N,{}); mm_use.resize(N,0); group.resize(N,-1); vis.resize(N,0); /* for(auto v:eda){ gr[0][v.first].pb(v.second); gr[0][v.second].pb(v.first); } for(int i=0;i,vector>>ans(cnt+1); for(int i=0;i ZANYO[200020]; int label[200020]; queue qu_bfs; int dokomade[200020]; void add_edge(int s_point,int t_point, int capa){// (s!=t) s-->t ZANYO[s_point].pb((edge){t_point,capa,(int)ZANYO[t_point].size(),0}); ZANYO[t_point].pb((edge){s_point,0,(int)ZANYO[s_point].size()-1,1}); } void bfs_dinic(int s_point,int N){ for(int i=0;i=0)continue; label[z.first]=z.second; for(auto &v:ZANYO[z.first]){ if(v.cap==0) continue; if(v.cap<0){ cout<<"error"<=0) continue; qu_bfs.push(mp(v.to,z.second+1)); } } return; } int dfs_dinic(int s_point,int t_point, int ima_flow){ if(s_point==t_point) return ima_flow; while(1){ if(dokomade[s_point]>=(int)ZANYO[s_point].size()) break; int edanum=dokomade[s_point]; if(ZANYO[s_point][edanum].cap<0){ cout<<"minus"<0){ ZANYO[s_point][edanum].cap-=W; ZANYO[ZANYO[s_point][edanum].to][ZANYO[s_point][edanum].num].cap+=W; return W; } dokomade[s_point]++; } return 0; } int max_flow_dinic(int s_point,int t_point,int N){// Nは頂点数:存在する頂点番号が[0,N) (Nは多めにしてOK) int jougen=inf*1000000000ll; //容量の最大値 int ans_flow=0; while(1){ bfs_dinic(s_point,N); if(label[t_point]==-1) break; for(int i=0;i>n>>m>>l; vectoreda; for(int i=0;i>y>>yy; y--,yy--; eda.pb(mp(y,n+yy)); add_edge(y,n+yy,1); } for(int i=0;imm; for(int i=0;i=0)mm.pb(mp(i,bi.match[i])); //cout<tmp(n+m,0); for(int i=n;iiro(n+m,-1); for(int i=0;i<=G;i++){ for(auto w:V[i].first)iro[w]=i; for(auto w:V[i].second)iro[w]=i; } for(auto v:eda){ int p=iro[v.first]; int q=iro[v.second]; if(p!=q)cout<<"Yes"<>n; for(int i=0;i