• 【BZOJ 3993】【SDOI 2015】星际战争


    http://www.lydsy.com/JudgeOnline/problem.php?id=3993
    调了好长时间啊
    这道题设时间为time,那么对于m个武器从S向这m个点连容量为time*Bi的边,代表能造成的总伤害。
    对于每个武器向每个能打到的机器人连容量为无穷的边。
    对每个机器人向T连自己的装甲量Ai为容量的边。
    又因为每个机器人的装甲必须被打完,所以设机器人到T的边的下界为Ai(上下界相等)。
    从T向S连容量无穷大的边,构成无源汇上下界网络流,所以构造附加网络然后二分时间time就可以了。
    注意附加网络有一些边是走不通的,所以把那些边删掉也可以
    (删掉后和原网络没什么区别啊Σ( ° △ °|||)︴那我还讲这么偏QAQ)。
    Dinic模板打残致命伤!!!以后改数组边界的时候注意要看程序中调用数组边界的常量是否也需要改。

    #include<cmath>
    #include<cstdio>
    #include<cstring>
    #include<algorithm>
    #define forU(i, a, b) for (int i = (a), UP = (b); i <= UP; ++i)
    #define forD(i, a, b) for (int i = (a), DOWN = (b); i >= DOWN; --i)
    using namespace std;
    
    inline double fabs(double x) {return x < 0 ? -x : x;}
    
    const int N = 55;
    const int M = N * N * 2;
    
    struct node {
    	int nxt, to; double c;
    	node(int _nxt = 0, int _to = 0, double _c = 0) : nxt(_nxt), to(_to), c(_c) {}
    } E[M];
    int cnt, point[N << 1], cur[N << 1];
    
    void ins(int u, int v, double c) {
    	E[++cnt] = node(point[u], v, c); point[u] = cnt;
    	E[++cnt] = node(point[v], u, 0); point[v] = cnt;
    }
    
    bool can[N][N];
    int q[N << 1], n, m, S, T, d[N << 1];
    double A[N], B[N];
    
    bool BFS() {
    	memset(d, 0, sizeof(int) * (T + 1));
    	int head = 0, tail = 1, u, v; q[1] = S; d[S] = 1;
    	while (head != tail) {
    		u = q[++head];
    		for (int i = point[u]; i; i = E[i].nxt)
    			if (fabs(E[i].c) > 1e-12 && !d[v = E[i].to]) {
    				d[v] = d[u] + 1;
    				q[++tail] = v;
    			}
    	}
    	return d[T];
    }
    
    double DFS(int u, double a) {
    	if (fabs(a) < 1e-12 || u == T) return a;
    	double f, flow = 0;
    	for (int &i = cur[u]; i; i = E[i].nxt)
    		if (d[E[i].to] == d[u] + 1 && fabs(f = DFS(E[i].to, min(a, E[i].c))) > 1e-12) {
    			a -= f; flow += f; E[i].c -= f; E[i ^ 1].c += f;
    			if (fabs(a) < 1e-12) break;
    		}
    	return flow;
    }
    
    double Dinic() {
    	double flow = 0;
    	while (BFS()) {
    		memcpy(cur, point, sizeof(int) * (T + 1));
    		flow += DFS(S, 100000000000000.0);
    	}
    	return flow;
    }
    
    double tot = 0;
    
    bool can2(double tim) {
    	memset(point, 0, sizeof(int) * (T + 1));
    	cnt = 1;
    	forU (i, 1, m) {
    		ins(S, i, tim * B[i]);
    		forU (j, 1, n)
    			if (can[i][j])
    				ins(i, j + m, 100000000000000.0);
    	}
    	forU (i, 1, n) ins(i + m, T, A[i]);
    	
    	return fabs(Dinic() - tot) < 1e-5;
    }
    
    int main() {
    	scanf("%d%d", &n, &m); S = n + m + 1; T = S + 1;
    	forU (i, 1, n) scanf("%lf", A + i), tot += A[i];
    	forU (i, 1, m) scanf("%lf", B + i);
    	int num;
    	double left = 0, right = 0, mid;
    	forU (i, 1, m)
    		forU (j, 1, n) {
    			scanf("%d", &num);
    			can[i][j] = num;
    		}
    	forU (i, 1, n) {
    		double fornow = 0;
    		forU (j, 1, m)
    			fornow = max(fornow, B[j]);
    		right += ceil(A[i] / fornow);
    	}
    	
    	while (right - left > 1e-5) {
    		mid = (left + right) / 2.0;
    		if (can2(mid)) right = mid;
    		else left = mid;
    	}
    	
    	printf("%lf
    ", left);
    	return 0;
    }
    
  • 相关阅读:
    P1135 奇怪的电梯
    pycharm设置快捷键在keymap下拉列表没有eclipse怎么办
    记录selenium简单实现自动点击操作
    selenium 批量下载文件,json,重命名
    python3.6+selenium使用chrome浏览器自动将文件下载到指定路径
    selenium + Java 设置文件默认下载路径
    详解介绍Selenium常用API的使用Java语言(完整版)
    Pycharm安装robot framework运行插件
    Python之robotframework+pycharm测试框架!
    基于Python3 Robot framework环境搭建
  • 原文地址:https://www.cnblogs.com/abclzr/p/6230778.html
Copyright © 2020-2023  润新知