博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
P4843 清理雪道
阅读量:5096 次
发布时间:2019-06-13

本文共 1946 字,大约阅读时间需要 6 分钟。

上下界网络流入门题...

每条边有一个下界流量 $1$,没有上界,求最小流

直接上下界最小流模板套进去就好了...

#include
#include
#include
#include
#include
#include
using namespace std;typedef long long ll;inline int read(){ int x=0,f=1; char ch=getchar(); while(ch<'0'||ch>'9') { if(ch=='-') f=-1; ch=getchar(); } while(ch>='0'&&ch<='9') { x=(x<<1)+(x<<3)+(ch^48); ch=getchar(); } return x*f;}const int N=1e5+7,M=4e6+7,INF=1e9+7;int fir[N],from[M],to[M],val[M],cntt=1;inline void add(int a,int b,int c){ from[++cntt]=fir[a]; fir[a]=cntt; to[cntt]=b; val[cntt]=c; from[++cntt]=fir[b]; fir[b]=cntt; to[cntt]=a; val[cntt]=0;}int n,SS,TT,sumflow[N],Ans;int dep[N],Fir[N],S,T;queue
q;bool BFS(){ for(int i=0;i<=n+3;i++) Fir[i]=fir[i],dep[i]=0; q.push(S); dep[S]=1; int x; while(!q.empty()) { x=q.front(); q.pop(); for(int i=fir[x];i;i=from[i]) { int &v=to[i]; if(dep[v]||!val[i]) continue; dep[v]=dep[x]+1; q.push(v); } } return dep[T]>0;}int DFS(int x,int mxf){ if(x==T||!mxf) return mxf; int fl=0,res; for(int &i=Fir[x];i;i=from[i]) { int &v=to[i]; if(dep[v]!=dep[x]+1||!val[i]) continue; if( res=DFS(v,min(mxf,val[i])) ) { mxf-=res; fl+=res; val[i]-=res; val[i^1]+=res; if(!mxf) break; } } return fl;}inline int Dinic() { int res=0; while(BFS()) res+=DFS(S,INF); return res; }inline void del(int x) { for(int i=fir[x];i;i=from[i]) val[i]=val[i^1]=0; }int main(){ n=read(); int a,b; S=0,T=n+1,SS=T+1,TT=SS+1; for(int i=1;i<=n;i++) { a=read(); for(int j=1;j<=a;j++) { b=read(),add(i,b,INF); sumflow[i]--; sumflow[b]++; } add(S,i,INF); add(i,T,INF); } for(int i=1;i<=n;i++) { if(!sumflow[i]) continue; if(sumflow[i]<0) add(i,TT,-sumflow[i]); else add(SS,i,sumflow[i]); } add(T,S,INF); int tmpS=S,tmpT=T; S=SS,T=TT; Dinic(); Ans=val[cntt]; val[cntt]=val[cntt^1]=0; del(SS),del(TT); S=tmpT,T=tmpS; printf("%d",Ans-Dinic()); return 0;}

 

转载于:https://www.cnblogs.com/LLTYYC/p/10805597.html

你可能感兴趣的文章
字符串处理
查看>>
ECharts(Enterprise Charts 商业产品图表库)初识
查看>>
LeetCode Factorial Trailing Zeroes (阶乘后缀零)
查看>>
hdu 5402 Travelling Salesman Problem (技巧,未写完)
查看>>
[AIR] 获取U盘,打开U盘
查看>>
HtmlUnitDriver 网页内容动态抓取
查看>>
ad logon hour
查看>>
获得进程可执行文件的路径: GetModuleFileNameEx, GetProcessImageFileName, QueryFullProcessImageName...
查看>>
证件照(1寸2寸)拍摄处理知识汇总
查看>>
罗马数字与阿拉伯数字转换
查看>>
Eclipse 反编译之 JadClipse
查看>>
asp.net 获取IP地理位置的几个主要接口
查看>>
Python入门-函数
查看>>
[HDU5727]Necklace(二分图最大匹配,枚举)
查看>>
距离公式汇总以及Python实现
查看>>
设计模式之装饰者模式
查看>>
【转】Linux内核调试方法总结
查看>>
一道不知道哪里来的容斥题
查看>>
Blender Python UV 学习
查看>>
window添加右键菜单
查看>>