zoukankan      html  css  js  c++  java
  • [bzoj1733][Usaco2005 feb]Secret Milking Machine 神秘的挤奶机_网络流

    [Usaco2005 feb]Secret Milking Machine 神秘的挤奶机

    题目大意翰正在制造一台新型的挤奶机,但他不希望别人知道.他希望尽可能久地隐藏这个秘密.他把挤奶机藏在他的农场里,使它不被发现.在挤奶机制造的过程中,他需要去挤奶机所在的地方T(1≤T≤200)次.他的农场里有秘密的地道,但约翰只在返回的时候用它.农场被划分成N(2≤N≤200)块区域,用1到200标号.这些区域被P(1≤P≤40000)条道路连接,每条路有一个小于10^6的长度L.两块区域之间可能有多条道路连接.为了减少被发现的可能,约翰不会两次经过农场上的任何一条道路.当然了,他希望走最短的路. 请帮助约翰寻找这T次从仓库走到挤奶机所在地的路线.仓库是区域1,挤奶机所在地是区域N.我们现在要求的是约翰经过的这些道路中最长的路的长度最小是多少,当然他不能重复走某一条路.请注意,我们要求的不是最短的总路程长度,而是所经过的直揍连接两个区域的道路中最长的道路的最小长度. 数据保证约翰可以找到T条没有重复的从仓库到挤奶机所在区域的路.

    数据范围:如题面。


    题解

    失了智........

    首先二分没问题,把最优化问题转换为判定性问题。

    接下来怎么办?

    其实只需要,原图怎么建图我们怎么建图,然后暴力跑最大流即可........

    真是有趣......

    代码:

    #include <bits/stdc++.h>
    
    #define N 210 
    
    #define M 40010 
    
    using namespace std;
    
    int head[N], to[M << 2], nxt[M << 2], val[M << 2], tot;
    
    queue<int >q;
    
    int dis[N], n, m, t, S, T;
    
    struct Edge {
    	int x, y, z;
    }e[M];
    
    char *p1, *p2, buf[100000];
    
    #define nc() (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2) ? EOF : *p1 ++ )
    
    int rd() {
    	int x = 0, f = 1;
    	char c = nc();
    	while (c < 48) {
    		if (c == '-')
    			f = -1;
    		c = nc();
    	}
    	while (c > 47) {
    		x = (((x << 2) + x) << 1) + (c ^ 48), c = nc();
    	}
    	return x * f;
    }
    
    inline void add(int x, int y, int z) {
    	to[ ++ tot] = y;
    	val[tot] = z;
    	nxt[tot] = head[x];
    	head[x] = tot;
    
    	to[ ++ tot] = x;
    	val[tot] = 0;
    	nxt[tot] = head[x];
    	head[x] = tot;
    }
    
    bool bfs() {
    	while (!q.empty())
    		q.pop();
    	memset(dis, -1, sizeof dis);
    	dis[S] = 0;
    	q.push(S);
    	while (!q.empty()) {
    	   int x = q.front();
    	   q.pop();
    	   for (int i = head[x]; i; i = nxt[i]) {
    		   if (dis[to[i]]  == -1 && val[i]) {
    			   dis[to[i]] = dis[x] + 1;
    			   if (to[i] == T)
    				   return true;
    			   q.push(to[i]);
    		   }
    	   }
       }
       return false;
    }
    
    int dinic(int x, int fl) {
    	int tmp = fl;
    	if (x == T)
    		return fl;
    	for (int i = head[x]; i; i = nxt[i]) {
    		if (dis[to[i]] == dis[x] + 1 && val[i]) {
    			int mdl = dinic(to[i], min(tmp, val[i]));
    			if (!mdl) dis[to[i]] = -1;
    			tmp -= mdl, val[i] -= mdl, val[i ^ 1] += mdl;
    			if(!tmp) break;
    		}
    	}
    	return fl - tmp;
    }
    
    void build(int x) {
    	for (int i = 1; i <= m; i ++ ) {
    		if (e[i].z <= x) {
    			add(e[i].x, e[i].y, 1);
    			add(e[i].y, e[i].x, 1);
    		}
    	}
    }
    
    bool check(int x) {
    	memset(head, 0, sizeof head);
    	tot = 1;
    	build(x);
    	int ans = 0;
    	while (bfs()) {
    		ans += dinic(1, 1 << 30);
    	}
    	return ans >= t;
    }
    
    int main() {
    	n = rd(), m = rd(), t = rd();
    	S = 1, T = n;
    	for (int i = 1; i <= m; i ++ ) {
    		e[i].x = rd(), e[i].y = rd(), e[i].z = rd();
    	}
    	int l = 1, r = 1000001;
    	while (l < r) {
    		int mid = (l + r) >> 1;
    		if (check(mid))
    			r = mid;
    		else
    			l = mid + 1;
    	}
    	printf("%d
    ", r);
    	return 0;
    }
    

    小结:我能总结些什么么......欧对,这个代码风格是真心好调真心好看,安利一波。

  • 相关阅读:
    洛谷P4382 [八省联考2018]劈配(网络流,二分答案)
    洛谷P3380 【模板】二逼平衡树(树套树,树状数组,线段树)
    C++实用整数快速输入输出模板(C++)
    洛谷P3348 [ZJOI2016]大森林(LCT,虚点,树上差分)
    洛谷P4338 [ZJOI2018]历史(LCT,树形DP,树链剖分)
    洛谷P3613 睡觉困难综合征(LCT,贪心)
    洛谷P3960 列队(NOIP2017)(Splay)
    洛谷P3275 [SCOI2011]糖果(差分约束,最长路,Tarjan,拓扑排序)
    博弈论总结(只会打表,永不证明)(博弈论)
    洛谷P1450 [HAOI2008]硬币购物(背包问题,容斥原理)
  • 原文地址:https://www.cnblogs.com/ShuraK/p/11235175.html
Copyright © 2011-2022 走看看