【最小割】【Dinic】【强联通分量缩点】bzoj1797 [Ahoi2009]Mincut 最小割
2024-08-26 22:25:07
结论:
满足条件一:当一条边的起点和终点不在 残量网络的 一个强联通分量中。且满流。
满足条件二:当一条边的起点和终点分别在 S 和 T 的强联通分量中。且满流。、
网上题解很多的。
#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
#include<queue>
using namespace std;
#define INF 2147483647
#define MAXN 4011
#define MAXM 120101
int v[MAXM],cap[MAXM],en,first[MAXN],next[MAXM];
int d[MAXN],cur[MAXN],cmp[MAXN],sum;
bool vis[MAXN];
queue<int>q;
vector<int>vs;
int n,m,S,T,A,B,C;
void Init_Dinic(){memset(first,-,sizeof(first)); en=;}
void AddEdge(const int &U,const int &V,const int &W)
{v[en]=V; cap[en]=W; next[en]=first[U]; first[U]=en++;
v[en]=U; next[en]=first[V]; first[V]=en++;}
bool bfs()
{
memset(d,-,sizeof(d)); q.push(S); d[S]=;
while(!q.empty())
{
int U=q.front(); q.pop();
for(int i=first[U];i!=-;i=next[i])
if(d[v[i]]==- && cap[i])
{
d[v[i]]=d[U]+;
q.push(v[i]);
}
}
return d[T]!=-;
}
int dfs(int U,int a)
{
if(U==T || !a) return a;
int Flow=,f;
for(int &i=cur[U];i!=-;i=next[i])
if(d[U]+==d[v[i]] && (f=dfs(v[i],min(a,cap[i]))))
{
cap[i]-=f; cap[i^]+=f;
Flow+=f; a-=f; if(!a) break;
}
if(!Flow) d[U]=-;
return Flow;
}
void max_flow()
{
while(bfs())
{
memcpy(cur,first,(n+)*sizeof(int));
while(dfs(S,INF));
}
}
void dfs(int U)
{
vis[U]=;
for(int i=first[U];i!=-;i=next[i]) if(cap[i]&&(!vis[v[i]])) dfs(v[i]);
vs.push_back(U);
}
void dfs2(int U)
{
cmp[U]=sum;
for(int i=first[U];i!=-;i=next[i]) if(cap[i^]&&(!cmp[v[i]])) dfs2(v[i]);
}
void scc()
{
for(int i=;i<=n;i++) if(!vis[i]) dfs(i);
vector<int>::iterator it=vs.end(); --it;
for(;;--it)
{
if(!cmp[*it]) {++sum; dfs2(*it);}
if(it==vs.begin()) break;
}
}
int main()
{
scanf("%d%d%d%d",&n,&m,&S,&T); Init_Dinic();
for(;m;--m)
{
scanf("%d%d%d",&A,&B,&C);
AddEdge(A,B,C);
}
max_flow(); scc();
for(int i=;i<en;i+=)
printf("%d %d\n",(!cap[i])&&cmp[v[i+]]!=cmp[v[i]],(!cap[i])&&cmp[v[i+]]==cmp[S]&&cmp[v[i]]==cmp[T]);
return ;
}
最新文章
- checkbox的全选与反选
- arcgis for flex展示GIS基本功能
- 设定自动获得DNS服务器地址
- Linux下多进程编程之exec函数语法及使用实例
- Poisson回归模型
- iOS学习之六种传值方式
- iOS:UIView的block函数实现转场动画---单视图
- spring读取prperties配置文件(1)
- novnc ignoring socket not reay
- java web实现 忘记密码(找回密码)功能及代码
- elike.python.function()
- Python2中while 1比while True更快
- (转)CentOS 7.0关闭默认防火墙启用iptables防火墙
- iOS上new Date异常解决办法
- MongoDB系列五(地理空间索引与查询).
- BZOJ 3817 Sum
- 优化算法系列-遗传算法(3)——NSGAII学习网址
- C#基础加强(6)之引用相等与运算符重载
- excel追加数据
- 20155311《网络对抗》Web安全基础实践