#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include "testlib.h" #define popcount __builtin_popcount using namespace std; using namespace atcoder; typedef long long ll; typedef pair P; struct unionfind{ int cnt; vector par, sz; unionfind() {} unionfind(int n):par(n), sz(n, 1), cnt(n){ for(int i=0; isz[y]) swap(x, y); par[x]=y; sz[y]+=sz[x]; } bool same(int x, int y){ return find(x)==find(y); } int size(int x){ return sz[find(x)]; } }; const int MAXN=100; const int MAXK=100; const ll INF=1e9; int n, k; vector gn[MAXN], gk[MAXK]; int vn[MAXN], vk[MAXK], ordn[MAXN], ordk[MAXK]; ll dp[MAXN][MAXK], dp2[MAXN][MAXK]; int a[MAXK], b[MAXK], c[MAXN], d[MAXN]; int main(int argc, char* argv[]){ registerValidation(argc, argv); k=inf.readInt(2, MAXK); inf.readEoln(); for(int i=0; i que; que.push(0); bool used[MAXK]={}; used[0]=1; int t=0; while(!que.empty()){ int x=que.front(); que.pop(); ordk[x]=t; vk[t++]=x; for(auto y:gk[x]){ if(used[y]) continue; que.push(y); used[y]=1; } } } { queue que; que.push(0); bool used[MAXN]={}; used[0]=1; int t=0; while(!que.empty()){ int x=que.front(); que.pop(); ordn[x]=t; vn[t++]=x; for(auto y:gn[x]){ if(used[y]) continue; que.push(y); used[y]=1; } } } for(int i=0; i=0; i--){ for(int j=k-1; j>=0; j--){ int x=vn[i], y=vk[j]; for(auto z:gn[x]){ if(ordn[z]0) n1--; if(j>0) k1--; if(k1>n1) continue; mcf_graph g(k1+n1+2); int s=k1+n1, t=s+1; for(int l=0; l=0){ g.add_edge(j1, i1+k1, 1, INF-dp[ordn[z]][ordk[w]]); } j1++; } i1++; } auto res=g.flow(s, t); if(res.first==k1){ dp[i][j]=max(dp[i][j], INF*k1-res.second); } } } ll ans=-INF; for(int i=0; i> st; st.insert({-INF, -1}); for(auto z:gn[x]){ if(ordn[z]first); st.insert({dp[ordn[z]][j]+1, z}); } int n1=gn[x].size(), k1=gk[y].size(); if(j>0) k1--; if(k1>n1) continue; if(k1==0){ for(auto z:gn[x]){ if(ordn[z] g0(k1+n1+2); int s=k1+n1, t=s+1; for(int l=0; l myon(n1); int i1=0; for(auto z:gn[x]){ int j1=0; for(auto w:gk[y]){ if(ordk[w]=0) g0.add_edge(j1, i1+k1, 1, INF-dp2[i][ordk[w]]); }else{ if(dp[ordn[z]][ordk[w]]>=0) g0.add_edge(j1, i1+k1, 1, INF-dp[ordn[z]][ordk[w]]); } j1++; } i1++; } auto res=g0.flow(s, t); ll mx=-INF; if(res.first==k1) mx=max(mx, INF*k1-res.second); if(j==0){ if(res.first==k1) ans=max(ans, mx); } auto edges=g0.edges(); for(auto e:edges){ if(e.from=k1 && e.to g(k1+n1+2); int s=k1+n1, t=s+1; for(int m=0; m=0) g.add_edge(j1, i1+k1, 1, INF-dp2[i][ordk[w]]); }else{ if(dp[ordn[z]][ordk[w]]>=0) g.add_edge(j1, i1+k1, 1, INF-dp[ordn[z]][ordk[w]]); } j1++; } i1++; } auto res=g.flow(s, t); if(res.first==k1){ dp2[ordn[z0]][j]=max(dp2[ordn[z0]][j], INF*k1-res.second); } } } } if(ans<0){ cout<<-1<