经典的最小割建图,用 ∞ ∞ 边体现限制条件。
inline char gc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int getint(){
char ch=gc(); int res=0,ff=1;
while(!isdigit(ch)) ch=='-'?ff=-1:0, ch=gc();
while(isdigit(ch)) res=(res<<3)+(res<<1)+ch-'0', ch=gc();
return res*ff;
}
const int maxn=100005,maxe=200005;
int N,m,S,T;
struct Edge{
int from,to,cap,flow;
Edge(int t1=0,int t2=0,int t3=0,int t4=0){ from=t1;to=t2;cap=t3;flow=t4; }
} Es[maxe];
int fir[maxn],nxt[maxe],tot=1,pos_fir[maxn],d[maxn];
void add(int x,int y,int z){
Es[++tot]=Edge(x,y,z,0);
nxt[tot]=fir[x]; fir[x]=tot;
Es[++tot]=Edge(y,x,0,0);
nxt[tot]=fir[y]; fir[y]=tot;
}
queue<int> que;
bool Bfs(){
memset(d,63,sizeof(d)); int INF=d[0];
d[S]=0; que.push(S);
while(!que.empty()){
int x=que.front(); que.pop();
for(int j=fir[x];j;j=nxt[j])
if(Es[j].cap>Es[j].flow&&d[Es[j].to]==INF) d[Es[j].to]=d[x]+1, que.push(Es[j].to);
}
return d[T]!=INF;
}
int Dfs(int x,int flw){
if(x==T||flw==0) return flw;
int t,res=0;
for(int &j=pos_fir[x];j;j=nxt[j])
if(d[x]+1==d[Es[j].to]&&(t=Dfs(Es[j].to,min(Es[j].cap-Es[j].flow,flw)))>0){
res+=t; flw-=t; Es[j].flow+=t; Es