WC模拟(1.2) T2 变量

标签: 省选 网络流 最大流 最小割
123人阅读 评论(0) 收藏 举报
分类:

变量

题目背景:

1.2 WC模拟T2

分析:网络流,最大流最小割

 

其实这个题····一看就是网络流最小割······结果我考场上脑子里全部是拆点,然后开心的没有想出建图······尴尬脸。

考虑,显然,W是来搞笑的,最后乘上去就好了,对于每一个变量建一个点i,显然每个等式后面的三项是可以直接提出来一起考虑的,对于每一个答案w[i]统计其贡献cnt[i],对于cnt[i]的正负性,讨论后选择在à i,和à T中的一条边上加上2 * abs(cnt[i]),然后对于一个等式a * |w[x] - w[y]|,我们在x, y之间连2 * a的无向边,表示,如果选择两个不同的方向则一定也要割掉中间的一条边,否则选择相同两个。然后解决变量间的限制,对于x < y直接强制xTINF的边,SyINF的边就可以了,如果对于x <= y,等价于不能存在x = W, y = -W的情况,那么所以直接xyINF的有向边,对于x == y,直接xyINF的边,yxINF的边就好了。然后原图最小割加上一开始默认值的答案就可以了。

Source:

 

/*
	created by scarlyw
*/
#include <cstdio>
#include <string>
#include <algorithm>
#include <cstring>
#include <iostream>
#include <cmath>
#include <cctype>
#include <vector>
#include <set>
#include <queue>
#include <ctime>
#include <bitset>

inline char read() {
	static const int IN_LEN = 1024 * 1024;
	static char buf[IN_LEN], *s, *t;
	if (s == t) {
		t = (s = buf) + fread(buf, 1, IN_LEN, stdin);
		if (s == t) return -1;
	}
	return *s++;
}

///*
template<class T>
inline void R(T &x) {
	static char c;
	static bool iosig;
	for (c = read(), iosig = false; !isdigit(c); c = read()) {
		if (c == -1) return ;
		if (c == '-') iosig = true;	
	}
	for (x = 0; isdigit(c); c = read()) 
		x = ((x << 2) + x << 1) + (c ^ '0');
	if (iosig) x = -x;
}
//*/

const int OUT_LEN = 1024 * 1024;
char obuf[OUT_LEN], *oh = obuf;
inline void write_char(char c) {
	if (oh == obuf + OUT_LEN) fwrite(obuf, 1, OUT_LEN, stdout), oh = obuf;
	*oh++ = c;
}

template<class T>
inline void W(T x) {
	static int buf[30], cnt;
	if (x == 0) write_char('0');
	else {
		if (x < 0) write_char('-'), x = -x;
		for (cnt = 0; x; x /= 10) buf[++cnt] = x % 10 + 48;
		while (cnt) write_char(buf[cnt--]);
	}
}

inline void flush() {
	fwrite(obuf, 1, oh - obuf, stdout);
}

/*
template<class T>
inline void R(T &x) {
	static char c;
	static bool iosig;
	for (c = getchar(), iosig = false; !isdigit(c); c = getchar())
		if (c == '-') iosig = true;	
	for (x = 0; isdigit(c); c = getchar()) 
		x = ((x << 2) + x << 1) + (c ^ '0');
	if (iosig) x = -x;
}
//*/

const int MAXN = 500 + 10;
const int INF = 1000000000;

struct node {
	int to, w, rev;
	node(int to = 0, int w = 0, int rev = 0)
		 : to(to), w(w), rev(rev) {}
} ;

int x, y, z, a, b, c, d, e, f, r, n, val, p, q, s, t, ans, data_case;
int temp[MAXN], dis[MAXN], cnt[MAXN];
std::vector<node> edge[MAXN];

inline void add_edge(int x, int y, int w) {
	edge[x].push_back(node(y, w, edge[y].size()));
	edge[y].push_back(node(x, 0, edge[x].size() - 1));
}

inline void build_graph() {
	R(n), R(val), R(p), R(q), s = 0, t = n + 1, ans = 0;
	for (int i = s; i <= t; ++i) edge[i].clear(), cnt[i] = 1;
	for (int i = 1; i <= p; ++i) {
		R(x), R(y), R(z), R(a), R(b), R(c), R(d), R(e), R(f);
		add_edge(x, y, 2 * a), add_edge(y, x, 2 * a);
		add_edge(y, z, 2 * b), add_edge(z, y, 2 * b);
		add_edge(z, x, 2 * c), add_edge(x, z, 2 * c);
		cnt[x] += (d - f), cnt[y] += (e - d), cnt[z] += (f - e);
	}
	for (int i = 1; i <= n; ++i)
		(cnt[i] > 0) ? (ans -= cnt[i], add_edge(i, t, 2 * cnt[i]))
			 : (ans += cnt[i], add_edge(s, i, -2 * cnt[i]));
	for (int i = 1; i <= q; ++i) {
		R(x), R(y), R(r);
		switch (r) {
			case 0: add_edge(x, y, INF);
					break ;
			case 1: add_edge(x, y, INF), add_edge(y, x, INF);
					break ;
			case 2: add_edge(s, y, INF), add_edge(x, t, INF);
		}
	}
	std::cerr << ans << '\n';
}

inline bool bfs(int s, int t) {
	memset(dis, -1, sizeof(int) * (n + 5));
	std::queue<int> q;
	dis[s] = 0, q.push(s);
	while (!q.empty()) {
		int cur = q.front();
		q.pop();
		for (int p = 0; p < edge[cur].size(); ++p) {
			node *e = &edge[cur][p];
			if (dis[e->to] == -1 && e->w > 0) {
				dis[e->to] = dis[cur] + 1;
				if (e->to == t) return true ;
				q.push(e->to);
			}
		}
	}
	return false;
}

inline int dfs(int cur, int low, int t) {
	if (cur == t) return low;
	int delta = 0;
	for (int &p = temp[cur]; p < edge[cur].size(); ++p) {
		node *e = &edge[cur][p];
		if (dis[e->to] == dis[cur] + 1 && e->w > 0) {
			int ret = dfs(e->to, std::min(low - delta, e->w), t);
			e->w -= ret, edge[e->to][e->rev].w += ret;
			if ((delta += ret) == low) break ;
		}
	}
	return delta;
}

inline int max_flow(int s, int t) {
	int ans = 0;
	while (bfs(s, t)) {
		int ret;
		while (memset(temp, 0, sizeof(int) * (n + 5)), 
			ret = dfs(s, INF, t)) 
				ans += ret;
	}
	return ans;
}

inline void solve() {
	build_graph();
	W((long long)(ans + max_flow(s, t)) * val), write_char('\n');
}

int main() {
	freopen("variable.in", "r", stdin);
	freopen("variable.out", "w", stdout);
	R(data_case);
	while (data_case--) std::cerr << data_case << '\n', solve();
	flush();
	return 0;
}
查看评论

WC模拟(1.10) T2 Permutation

结论推理
  • scar_lyw
  • scar_lyw
  • 2018年01月17日 20:23
  • 105

WC模拟(1.4) T2 轰炸

tarjan缩点 + DP最长路
  • scar_lyw
  • scar_lyw
  • 2018年01月11日 15:09
  • 94

2017.1.17WC模拟赛总结

boom zero了,很不开心,总结一下问题。 T1boom zero了,原因是常数太大自己又没有去卡。 T2boom zero了,原因是没看好题,不符合它输出要求。 T3boom zero了,...
  • WerKeyTom_FTD
  • WerKeyTom_FTD
  • 2017年01月17日 14:47
  • 512

WC模拟(1.2) T1 串

结论分析
  • scar_lyw
  • scar_lyw
  • 2018年01月03日 14:08
  • 113

WC2018口胡记

作为一个noip考挂选手,我并不能去WC2018。 当天下午一点半左右问@tututu要来题目,一开始只有T1的复述版本,感觉是个数据结构题,但并没有什么思路,而且部分分表也没有,似乎只会O(n2)...
  • zxin__
  • zxin__
  • 2018年02月10日 23:13
  • 305

WC2017总结

~~~~~~吾之于WC,如蛙之于井底,仰而望天之高,跃而不及地。Part 1突然开始的WC模拟      ~~~~~~noip2016在11月20日结束,经过一段时间的乱搞,我们于12月开始,为期3个...
  • rzO_KQP_Orz
  • rzO_KQP_Orz
  • 2017年02月09日 20:17
  • 4092

WC模拟(12.20) T2 随机二分图(bzoj5006)

记忆化搜索 + hash状态剪枝
  • scar_lyw
  • scar_lyw
  • 2017年12月21日 11:06
  • 314

WC模拟(1.8) T2 送你一个集合

二进制
  • scar_lyw
  • scar_lyw
  • 2018年01月11日 15:59
  • 186

WC模拟(1.6) T2 送你一棵圣诞树

树状数组 + 线段树 + LCA
  • scar_lyw
  • scar_lyw
  • 2018年01月11日 15:48
  • 115

WC模拟(1.2) T3 取石子

结论分析
  • scar_lyw
  • scar_lyw
  • 2018年01月03日 14:13
  • 151
    个人资料
    专栏达人 持之以恒
    等级:
    访问量: 2万+
    积分: 2169
    排名: 2万+
    博客专栏
    最新评论