对于30%的数据,N ≤ 20,M ≤ 120。
对于100%的数据,N ≤ 200,M ≤ 20000。
/* 比较裸的一道费用流,就当打打板子。 需要注意的是一个点只能用一次,所以要拆点,在两点之间建一条容量为1的边。 */ #include<cstdio> #include<iostream> #define N 410 #define M 41000 #define inf 1000000000 using namespace std; int head[N],dis[N],q[N*10],vis[N],fa[N],n,m,S,T,cnt=1,timer,ans; struct node{ int v,pre,f,w; };node e[M]; void add(int u,int v,int f,int w){ e[++cnt].v=v;e[cnt].f=f;e[cnt].w=w;e[cnt].pre=head[u];head[u]=cnt; e[++cnt].v=u;e[cnt].f=0;e[cnt].w=-w;e[cnt].pre=head[v];head[v]=cnt; } bool spfa(){ for(int i=1;i<=T;i++)dis[i]=inf; int h=0,t=1; q[1]=S;dis[S]=0;vis[S]=1; while(h<t){ int x=q[++h];vis[x]=0; for(int i=head[x];i;i=e[i].pre){ int v=e[i].v; if(dis[v]>dis[x]+e[i].w&&e[i].f){ dis[v]=dis[x]+e[i].w; fa[v]=i; if(!vis[v]){ q[++t]=v; vis[v]=1; } } } } return dis[T]!=inf; } void updata(){ int tmp=T; while(tmp!=S){ int l=fa[tmp],v=e[l].v; e[l].f--;e[l^1].f++; tmp=e[l^1].v; } ans+=dis[T]; } int main(){ scanf("%d%d",&n,&m); S=0,T=n*2+1; add(S,1,inf,0);add(n*2,T,inf,0); for(int i=2;i<n;i++)add(i,i+n,1,0); add(1,1+n,inf,0);add(n,n+n,inf,0); for(int i=1;i<=m;i++){ int x,y,z;scanf("%d%d%d",&x,&y,&z); add(x+n,y,1,z); } while(spfa()){ timer++;updata(); } printf("%d %d",timer,ans); return 0; }
原文:http://www.cnblogs.com/harden/p/6260487.html