题目
思路
本题是最小割点模板,不知道大家有没有见过最小割边模板,那个是个网络瘤,作为网络瘤患者忠实爱好者,我们可以通过拆点把最小割点变成最小割边,相当于把x拆成x和x+n2个点,一个只有入边,一个只有出边,从只有出边的点向只有入边的点连一个容量为1的单向边,原来的边容量为inf。
code:
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
using namespace std;
int n,m,s,t,x,y,tot=2,head[2001],dep[2001],l,r,u[2001],r2;
long long w,ans;
long long mn(long long x,long long y)
{
return (x>y?y:x);
}
struct f{
int to,net;
long long w;
} a[100011];
void add(int x,int y,long long w)
{
a[tot].to=y,a[tot].w=w,a[tot].net=head[x],head[x]=tot++;
return;
}
bool bfs()
{
memset(dep,0,sizeof(dep));
l=0,r=0;
u[r++]=s;
dep[s]=1;
while (l<r)
{
r2=r;
for (int i=l;i<r;i++)
{
for (int j=head[u[i]];j;j=a[j].net)
{
if (a[j].w!=0&&dep[a[j].to]==0)
{
u[r2++]=a[j].to;
dep[a[j].to]=dep[u[i]]+1;
}
}
}
l=r,r=r2;
}
return dep[t];
}
long long dfs(int d,long long in)
{
if (d==t) return in;
long long out=0;
int uw=dep[d]+1;
for (int j=head[d];j&∈j=a[j].net)
{
if (a[j].w==0||dep[a[j].to]!=uw) continue;
long long s=dfs(a[j].to,mn(in,a[j].w));
out+=s,in-=s;
a[j].w-=s,a[j^1].w+=s;
}
if (out==0)
{
dep[d]=0;
}
return out;
}
int main()
{
scanf("%d%d%d%d",&n,&m,&s,&t);
s+=n;
w=1e30;
for (int i=1;i<=n;i++)
{
add(i,i+n,1);
add(i+n,i,0);
}
for (int i=0;i<m;i++)
{
scanf("%d%d",&x,&y);
add(x+n,y,w);
add(y,x+n,0);
add(y+n,x,w);
add(x,y+n,0);
}
while (bfs())
{
ans+=dfs(s,1e18);
}
cout<<ans;
return 0;
}