HDU6214 2017青岛赛区网络赛

本弱并不会骚套路，当时并想了很久做不出。

将每个边的流量w变成w*(m+1)+1，因为边总数为m，所以后面那个1并不会影响最大流等于最小割的性质，再加上后面那个1体现了边数，跑出来的最大流为a*(m+1)+b，%掉m+1,就得到b，而每一种最大流等于最小割都可以得到相同的a，但b是不同的，由于最大流等于最小割，这个b最大就是最小割的边数。

#include<cstdio>#include<cstring>#define maxn 210#define maxl 1010#define inf 2000000001int n,m,S,T,cnt,ans;int cur[maxn],dis[maxn],a[maxn],ehead[maxn];struct ed{int to,nxt,w;} e[maxl<<1];void add(int u,int v,int w){e[++cnt].nxt=ehead[u];e[cnt].to=v;e[cnt].w=w;ehead[u]=cnt;}void prework(){scanf("%d%d",&n,&m);scanf("%d%d",&S,&T);for(int i=1;i<=n;i++)ehead[i]=0;int u,v,w;cnt=0;for(int i=1;i<=m;i++){scanf("%d%d%d",&u,&v,&w);add(u,v,w*(m+1)+1);add(v,u,0);}}bool bfs(){int u,v,head=0,tail=1;for(int i=1;i<=n;i++)dis[i]=-1;a[1]=S;dis[S]=0;while(head<tail){u=a[++head];for(int i=ehead[u];i;i=e[i].nxt){v=e[i].to;if(dis[v]<0 && e[i].w){dis[v]=dis[u]+1;a[++tail]=v;}}}if(dis[T]<0)return false;elsereturn true;}int min(int a,int b){if(a<b)return a;elsereturn b;}int find(int u,int low){int t,v;if(u==T)return low;for(int &i=cur[u];i;i=e[i].nxt){v=e[i].to;if(dis[v]==dis[u]+1 && e[i].w && (t=find(v,min(low,e[i].w)))){if(i&1)e[i+1].w+=t;elsee[i-1].w+=t;e[i].w-=t;return t;}}return 0;}void mainwork(){int t;ans=0;while(bfs()){for(int i=1;i<=n;i++)cur[i]=ehead[i];while((t=find(S,inf)))ans+=t;}}void print(){ans%=m+1;printf("%d\n",ans);}int main(){int t;scanf("%d",&t);for(int i=1;i<=t;i++){prework();mainwork();print();}return 0;}

据说以前有一道原题是跑完一遍以后吧满流设为1，否则为inf，再跑一遍，这样感觉没什么道理，最大流的时候最小割的数量是无法保证的，不同的跑法得到的满流边集都不一样，那道题是数据弱，这里就过不了了。然而又有人说反向边不改流量这样也可以过，这。。。就更没道理了，你不反悔那还叫网络流的算法？数据是会玩的。、、