码迷,mamicode.com
首页 > 其他好文 > 详细

【COGS】894. 追查坏牛奶

时间:2015-03-01 17:03:14      阅读:177      评论:0      收藏:0      [点我收藏+]

标签:

http://cojs.tk/cogs/problem/problem.php?pid=894

题意:n个点m条边的加权网络,求最少边数的按编号字典序最小的最小割。(n<=32, m<=1000)

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
struct Gr {
	static const int N=33, M=1005, oo=~0u>>1;
	struct E { int next, from, to, cap; }e[M<<1];
	int ihead[N], cnt, d[N], p[N], gap[N], cur[N];
	Gr() { memset(ihead, 0, sizeof ihead); cnt=1; }
	Gr & operator=(const Gr &g) {
		memcpy(ihead, g.ihead, sizeof g.ihead);
		memcpy(e, g.e, sizeof g.e);
		cnt=g.cnt;
		return *this;
	}
	void add(int u, int v, int cap) {
		e[++cnt]=(E){ihead[u], u, v, cap}; ihead[u]=cnt;
		e[++cnt]=(E){ihead[v], v, u, 0}; ihead[v]=cnt;
	}
	ll isap(int s, int t, int n) {
		for(int i=0; i<=n; ++i) d[i]=0, gap[i]=0, cur[i]=ihead[i];
		gap[0]=n; int u=s, i, f; ll ret=0;
		while(d[s]<n) {
			for(i=cur[u]; i; i=e[i].next) if(e[i].cap && d[e[i].to]+1==d[u]) break;
			if(i) {
				cur[u]=i; p[e[i].to]=i; u=e[i].to;
				if(u==t) {
					for(f=oo; u!=s; u=e[p[u]].from) f=min(f, e[p[u]].cap);
					for(u=t; u!=s; u=e[p[u]].from) e[p[u]].cap-=f, e[p[u]^1].cap+=f;
					ret+=f;
				}
			}
			else {
				if(!(--gap[d[u]])) break;
				d[u]=n; cur[u]=ihead[u];
				for(i=ihead[u]; i; i=e[i].next) if(e[i].cap && d[e[i].to]+1<d[u]) d[u]=d[e[i].to]+1;
				++gap[d[u]];
				if(u!=s) u=e[p[u]].from;
			}
		}
		return ret;
	}
}g, G;
int n, m, arr[Gr::M], tot;
int main() {
	freopen("milk6.in", "r", stdin);
	freopen("milk6.out", "w", stdout);
	scanf("%d%d", &n, &m);
	for(int i=1; i<=m; ++i) {
		int x, y, cap;
		scanf("%d%d%d", &x, &y, &cap);
		g.add(x, y, cap*(m+1)+1);
	}
	ll ans=0, f;
	G=g;
	f=G.isap(1, n, n);
	for(int i=1; i<=m; ++i) {
		int id=i<<1, cap=g.e[id].cap;
		g.e[id].cap=g.e[id^1].cap=0;
		G=g;
		ll mn=G.isap(1, n, n);
		if(mn+cap==f) arr[++tot]=i, ans+=(cap-1)/(m+1), f=mn;
		else g.e[id].cap=cap, g.e[id^1].cap=0;
	}
	printf("%lld %d ", ans, tot);
	for(int i=1; i<=tot; ++i) printf("%d ", arr[i]);
	return 0;
}

  

神题...在uoj群被神犇们吊打&裱

= =这么水的usaco training我都不会做= =

首先求最少边的话可以直接将权值变为$w*(m+1)+1$,然后就是最小割。理由很简单,最小割是$\sum_{i为一条割边} (w[i]*(m+1)+1) = (m+1)\sum_{i为一条割边} w[i] + c$首先我们分离出了原来没有改变权值的割,显然这样做$\sum_{i为一条割边} w[i]$是原图的最小割,因为乘上$(m+1)$后始终大于边数,和边数$c$都是常数= =。其次我们在求改变权值后的最小割后,比较完原图最小割后还比较了右边常数$c$,而$c$正是割边数。而$(m+1)$这个乘数也是因为割边数<=m最大可能到了m,可能会影响到原来图上的边权值,所以乘上$(m+1)$来消除影响。

那么现在问题变为求普通的字典序最少的最小割。我们考虑从小到大枚举边然后删边。

如果删掉$i$这条边时,最小割变小的值恰好为$i$改变后的权值,那么显然$i$是割边,此时维护的最小割减小,删掉$i$不恢复。

否则把删掉的边恢复。

【COGS】894. 追查坏牛奶

标签:

原文地址:http://www.cnblogs.com/iwtwiioi/p/4307241.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!