#pragma GCC optimize("Ofast") #pragma GCC optimize("unroll-loops") #pragma GCC optimize("inline") #include using namespace std; template struct cLtraits_identity{ using type = T; } ; template using cLtraits_try_make_signed = typename conditional< is_integral::value, make_signed, cLtraits_identity >::type; template struct cLtraits_common_type{ using tS = typename cLtraits_try_make_signed::type; using tT = typename cLtraits_try_make_signed::type; using type = typename common_type::type; } ; void*wmem; char memarr[96000000]; template inline auto min_L(S a, T b) -> typename cLtraits_common_type::type{ return (typename cLtraits_common_type::type) a <= (typename cLtraits_common_type::type) b ? a : b; } template inline auto max_L(S a, T b) -> typename cLtraits_common_type::type{ return (typename cLtraits_common_type::type) a >= (typename cLtraits_common_type::type) b ? a : b; } template inline void walloc1d(T **arr, int x, void **mem = &wmem){ static int skip[16] = {0, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1}; (*mem) = (void*)( ((char*)(*mem)) + skip[((unsigned long long)(*mem)) & 15] ); (*arr)=(T*)(*mem); (*mem)=((*arr)+x); } template inline void walloc1d(T **arr, int x1, int x2, void **mem = &wmem){ walloc1d(arr, x2-x1, mem); (*arr) -= x1; } inline int my_getchar_unlocked(){ static char buf[1048576]; static int s = 1048576; static int e = 1048576; if(s == e && e == 1048576){ e = fread_unlocked(buf, 1, 1048576, stdin); s = 0; } if(s == e){ return EOF; } return buf[s++]; } inline void rd(int &x){ int k; int m=0; x=0; for(;;){ k = my_getchar_unlocked(); if(k=='-'){ m=1; break; } if('0'<=k&&k<='9'){ x=k-'0'; break; } } for(;;){ k = my_getchar_unlocked(); if(k<'0'||k>'9'){ break; } x=x*10+k-'0'; } if(m){ x=-x; } } struct MY_WRITER{ char buf[1048576]; int s; int e; MY_WRITER(){ s = 0; e = 1048576; } ~MY_WRITER(){ if(s){ fwrite_unlocked(buf, 1, s, stdout); } } } ; MY_WRITER MY_WRITER_VAR; void my_putchar_unlocked(int a){ if(MY_WRITER_VAR.s == MY_WRITER_VAR.e){ fwrite_unlocked(MY_WRITER_VAR.buf, 1, MY_WRITER_VAR.s, stdout); MY_WRITER_VAR.s = 0; } MY_WRITER_VAR.buf[MY_WRITER_VAR.s++] = a; } inline void wt_L(char a){ my_putchar_unlocked(a); } inline void wt_L(int x){ int s=0; int m=0; char f[10]; if(x<0){ m=1; x=-x; } while(x){ f[s++]=x%10; x/=10; } if(!s){ f[s++]=0; } if(m){ my_putchar_unlocked('-'); } while(s--){ my_putchar_unlocked(f[s]+'0'); } } inline void wt_L(const char c[]){ int i=0; for(i=0;c[i]!='\0';i++){ my_putchar_unlocked(c[i]); } } template struct maxflow{ int node; int st; int ed; int*es; int*emem; int**edge; int**rev; int*level; int*qq; T**flow; T eps; void malloc(int N){ int i; es = (int*)std::malloc(N*sizeof(int)); emem = (int*)std::malloc(N*sizeof(int)); level = (int*)std::malloc(N*sizeof(int)); qq = (int*)std::malloc(N*sizeof(int)); edge = (int**)std::malloc(N*sizeof(int*)); rev = (int**)std::malloc(N*sizeof(int*)); flow = (T**)std::malloc(N*sizeof(T*)); for(i=(0);i<(N);i++){ emem[i] = 0; edge[i] = rev[i] = NULL; flow[i] = NULL; } } void malloc(int N, int init_flag){ int i; es = (int*)std::malloc(N*sizeof(int)); emem = (int*)std::malloc(N*sizeof(int)); level = (int*)std::malloc(N*sizeof(int)); qq = (int*)std::malloc(N*sizeof(int)); edge = (int**)std::malloc(N*sizeof(int*)); rev = (int**)std::malloc(N*sizeof(int*)); flow = (T**)std::malloc(N*sizeof(T*)); for(i=(0);i<(N);i++){ emem[i] = 0; edge[i] = rev[i] = NULL; flow[i] = NULL; } if(init_flag){ init(N); } } void walloc(int N, void**mem = &wmem){ int i; walloc1d(&es, N, mem); walloc1d(&emem, N, mem); walloc1d(&level, N, mem); walloc1d(&qq, N, mem); walloc1d(&edge, N, mem); walloc1d(&rev, N, mem); walloc1d(&flow, N, mem); (*mem) = (flow + N); } void walloc(int N, int init_flag, void**mem = &wmem){ int i; walloc1d(&es, N, mem); walloc1d(&emem, N, mem); walloc1d(&level, N, mem); walloc1d(&qq, N, mem); walloc1d(&edge, N, mem); walloc1d(&rev, N, mem); walloc1d(&flow, N, mem); (*mem) = (flow + N); if(init_flag){ init(N); } } void levelize(void){ int i; int j; int k; int t; int q_st = 0; int q_ed = 1; for(i=(0);i<(node);i++){ level[i] = -1; } level[st] = 0; qq[0] = st; while(q_st != q_ed){ i = qq[q_st++]; t = level[i] + 1; for(j=(0);j<(es[i]);j++){ if(flow[i][j] > eps){ k = edge[i][j]; if(level[k]!=-1){ continue; } level[k] = t; qq[q_ed++] = k; if(k==ed){ return; } } } } } S pushflow(int i, S lim){ int j; int k; int ji; S s; S t; S res = 0; if(i==ed){ return lim; } for(j=(0);j<(es[i]);j++){ if(flow[i][j] > eps){ k = edge[i][j]; if(level[k] != level[i]+1){ continue; } s =min_L(lim, (S)flow[i][j]); t = pushflow(k, s); if(!t){ continue; } res += t; lim -= t; ji = rev[i][j]; flow[i][j] -= t; flow[k][ji] += t; if(!lim){ break; } } } if(lim){ level[i] = -1; } return res; } S solve(int st_, int ed_){ S res = 0; st = st_; ed = ed_; for(;;){ levelize(); if(level[ed] == -1){ break; } res += pushflow(st, numeric_limits::max()); } return res; } void init(int N){ int i; node = N; for(i=(0);i<(N);i++){ es[i] = 0; } eps = (T)1e-9; } void memoryExpand(int i, int sz){ if(sz <= emem[i]){ return; } sz =max_L(sz,max_L(3, emem[i]*2)); emem[i]=sz; edge[i] = (int*)realloc(edge[i], sz*sizeof(int)); rev[i] = (int*)realloc(rev[i], sz*sizeof(int)); flow[i] = (T*)realloc(flow[i], sz*sizeof(T)); } void addEdge(int n1, int n2, T f1, T f2 = 0){ int s1 = es[n1]++; int s2 = es[n2]++; if(s1 >= emem[n1]){ memoryExpand(n1, es[n1]); } if(s2 >= emem[n2]){ memoryExpand(n2, es[n2]); } edge[n1][s1]=n2; edge[n2][s2]=n1; flow[n1][s1]=f1; flow[n2][s2]=f2; rev[n1][s1]=s2; rev[n2][s2]=s1; } void addEdgeAdv(int n1, int n2, T f1, T f2 = 0){ int s1 = es[n1]++; int s2 = es[n2]++; edge[n1][s1]=n2; edge[n2][s2]=n1; flow[n1][s1]=f1; flow[n2][s2]=f2; rev[n1][s1]=s2; rev[n2][s2]=s1; } void setGraph(int N, int M, int n1[], int n2[], T f1[], T f2[]){ int i; node = N; for(i=(0);i<(N);i++){ es[i] = 0; } for(i=(0);i<(M);i++){ es[n1[i]]++; es[n2[i]]++; } for(i=(0);i<(N);i++){ memoryExpand(i, es[i]); } for(i=(0);i<(N);i++){ es[i] = 0; } for(i=(0);i<(M);i++){ addEdgeAdv(n1[i], n2[i], f1[i], f2[i]); } eps = (T)1e-9; } void setGraph_w(int N, int M, int n1[], int n2[], T f1[], T f2[], void **mem = wmem){ int i; int j; int k; node = N; for(i=(0);i<(N);i++){ es[i] = emem[i] = 0; } for(i=(0);i<(M);i++){ es[n1[i]]++; es[n2[i]]++; } edge[0] = (int*)(*mem); int Oq6TK1Zh = N; for(i=(1);i<(Oq6TK1Zh);i++){ edge[i] = edge[i-1] + es[i-1]; } rev[0] = edge[N-1] + es[N-1]; int jG1yfsum = N; for(i=(1);i<(jG1yfsum);i++){ rev[i] = rev[i-1] + es[i-1]; } flow[0] = (T*)(rev[N-1] + es[N-1]); int n23oAcQn = N; for(i=(1);i<(n23oAcQn);i++){ flow[i] = flow[i-1] + es[i-1]; } *mem = (void*)(flow[N-1] + es[N-1]); for(i=(0);i<(N);i++){ es[i] = 0; } for(i=(0);i<(M);i++){ addEdgeAdv(n1[i], n2[i], f1[i], f2[i]); } eps = (T)1e-9; } } ; int N; int A[100000]; int B[100000]; int main(){ wmem = memarr; int i; int j; int k; int node; int st; int ed; int f; maxflow flow; rd(N); { int Lj4PdHRW; for(Lj4PdHRW=(0);Lj4PdHRW<(N);Lj4PdHRW++){ rd(A[Lj4PdHRW]);A[Lj4PdHRW] += (-1); rd(B[Lj4PdHRW]);B[Lj4PdHRW] += (-1); } } node = 2*N; st = node++; ed = node++; flow.walloc(node,1); for(i=(0);i<(N);i++){ flow.addEdge(st, i, 1); } for(i=(0);i<(N);i++){ flow.addEdge(N+i, ed, 1); } for(i=(0);i<(N);i++){ flow.addEdge(i, N+A[i], 1); flow.addEdge(i, N+B[i], 1); } f = flow.solve(st,ed); if(f==N){ wt_L("Yes"); wt_L('\n'); for(i=(0);i<(N);i++){ for(j=(0);j<(flow.es[i]);j++){ k = flow.edge[i][j] - N; if(flow.flow[i][j] == 0 && 0 <= k && k < N){ wt_L(k+1); wt_L('\n'); } } } } else{ wt_L("No"); wt_L('\n'); } return 0; } // cLay version 20210717-1 [beta] // --- original code --- // int N, A[1d5], B[]; // { // int i, j, k; // int node, st, ed, f; // maxflow flow; // rd(N,(A--,B--)(N)); // node = 2*N; // st = node++; // ed = node++; // flow.walloc(node,1); // rep(i,N) flow.addEdge(st, i, 1); // rep(i,N) flow.addEdge(N+i, ed, 1); // rep(i,N) flow.addEdge(i, N+A[i], 1), flow.addEdge(i, N+B[i], 1); // f = flow.solve(st,ed); // if(f==N){ // wt("Yes"); // rep(i,N) rep(j,flow.es[i]){ // k = flow.edge[i][j] - N; // if(flow.flow[i][j] == 0 && 0 <= k < N) wt(k+1); // } // } else { // wt("No"); // } // }