bzoj
luogu
sol
很显然的建边方式:普通的桥建双向流量\(inf\)的边,危桥建双向流量\(2\)的边。
然后分别从\(a_1\),\(b_1\)向\(a_2\),\(b_2\)跑最大流,保证流量分别要大于等于\(a_n\),\(b_n\)。
多源汇且源汇一一对应的网络流是一类NPC问题,但在这里只有两源两汇,可以采用两次\(Dinic\)的方式,第一次\(S\)连\(a_1\),\(b_1\),\(T\)连\(a_2\),\(b_2\)跑最大流,第二次\(S\)连\(a_1\),\(b_2\),\(T\)连\(a_2\),\(b_1\)再跑一遍,只要两遍都能保证满流就行了。
code
#include
#include
#include
#include
using namespace std;
const int N = 60;
const int inf = 1e9;
struct edge{int to,nxt,w;}a[N*N<<2];
int n,a1,a2,an,b1,b2,bn,head[N],cnt,dep[N],cur[N];
char map[N][N];queueQ;
void link(int u,int v,int w)
{
a[++cnt]=(edge){v,head[u],w};
head[u]=cnt;
a[++cnt]=(edge){u,head[v],w};
head[v]=cnt;
}
bool bfs(int s,int t)
{
memset(dep,0,sizeof(dep));
dep[s]=1;Q.push(s);
while (!Q.empty())
{
int u=Q.front();Q.pop();
for (int e=head[u];e;e=a[e].nxt)
if (a[e].w&&!dep[a[e].to])
dep[a[e].to]=dep[u]+1,Q.push(a[e].to);
}
return dep[t];
}
int dfs(int u,int f,int t)
{
if (u==t) return f;
for (int &e=cur[u];e;e=a[e].nxt)
if (a[e].w&&dep[a[e].to]==dep[u]+1)
{
int tmp=dfs(a[e].to,min(a[e].w,f),t);
if (tmp) {a[e].w-=tmp;a[e^1].w+=tmp;return tmp;}
}
return 0;
}
int Dinic(int s,int t)
{
int res=0;
while (bfs(s,t))
{
for (int i=1;i<=t;++i) cur[i]=head[i];
while (int tmp=dfs(s,inf,t)) res+=tmp;
}
return res;
}
int solve(int s1,int s2,int t1,int t2)
{
int s=n+1,t=n+2;
memset(head,0,sizeof(head));cnt=1;
for (int i=1;i