luogu2402 奶牛隐藏

news/2024/7/5 20:29:45

题目描述

在一个农场里有n块田地。某天下午,有一群牛在田地里吃草,他们分散在农场的诸多田地上,农场由m条无向的路连接,每条路有不同的长度。
突然,天降大雨,奶牛们非常混乱,想要快点去躲雨。已知每个田地都建立有一个牛棚,但是每个牛棚只能容纳一定数量的牛躲雨,如果超过这个数量,那多出的牛只能去别的田地躲雨。奶牛们每移动1的距离花费1时间,奶牛们想知道它们全部都躲进牛棚,最少需要多少时间。(即最后一头奶牛最少要花多久才能躲进牛棚)。

输入输出格式

输入格式:

第一行输入两个整数N,M。N表示田地块数,M表示路径数。
接下来N行,每行两个整数S,P,分别表示该田地现在有几头牛以及该田地的牛棚最多可以容纳多少牛。
接下来M行,每行3个整数A,B,C,表示存在一条路径连接A,B,并且它的长度为C。

输出格式:

一个整数表示所有奶牛全都躲进牛棚所用的最少时间。如果无法使全部奶牛都躲进牛棚,输出-1。

输入输出样例

输入样例#1:

3 4 
7 2 
0 4 
2 6 
1 2 40 
3 2 70 
2 3 90 
1 3 120

输出样例#1:

110

说明

【样例解释】

1号点的两只牛直接躲进1号牛棚,剩下的5只中,4只跑去2号点,还有一只从1->2->3,3号点的2只牛也直接躲进去,这样最慢的牛花费的时间是110。

数据范围 : 对于100%的数据,N<=200 M<=1500


网络流,首先floyed跑出每两个农场间的最短路
考虑二分答案,每次只把路径长度小与二分的答案的边构图
网络流判断是否满流即可,注意开long long


# include <bits/stdc++.h>
# define IL inline
# define RG register
# define Fill(a, b) memset(a, b, sizeof(a))
# define Copy(a, b) memcpy(a, b, sizeof(a))
using namespace std;
typedef long long ll;
const int _(410), __(4e5 + 10);
const ll INF(1e18);

IL ll Read(){
    RG char c = getchar(); RG ll x = 0, z = 1;
    for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
    for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
    return x * z;
}
# define int ll
int n, m, num, w[__], fst[_], nxt[__], to[__], cnt, s[_], p[_];
int S, T, lev[_], cur[_], max_flow, ans, dis[_][_];
queue <int> Q;

IL void Add(RG int u, RG int v, RG int f){
    if(!f) return;
    w[cnt] = f; to[cnt] = v; nxt[cnt] = fst[u]; fst[u] = cnt++;
    w[cnt] = 0; to[cnt] = u; nxt[cnt] = fst[v]; fst[v] = cnt++;
}

IL int Dfs(RG int u, RG int maxf){
    if(u == T) return maxf;
    RG int ret = 0;
    for(RG int &e = cur[u]; e != -1; e = nxt[e]){
        if(lev[to[e]] != lev[u] + 1 || !w[e]) continue;
        RG int f = Dfs(to[e], min(w[e], maxf - ret));
        ret += f; w[e ^ 1] += f; w[e] -= f;
        if(ret == maxf) break;
    }
    return ret;
}

IL bool Bfs(){
    Fill(lev, 0); lev[S] = 1; Q.push(S);
    while(!Q.empty()){
        RG int u = Q.front(); Q.pop();
        for(RG int e = fst[u]; e != -1; e = nxt[e]){
            if(lev[to[e]] || !w[e]) continue;
            lev[to[e]] = lev[u] + 1;
            Q.push(to[e]);
        }
    }
    return lev[T];
}

IL bool Check(RG int lim){
    Fill(fst, -1); cnt = 0;
    for(RG int i = 1; i <= n; i++) Add(S, i, s[i]), Add(i + n, T, p[i]);
    for(RG int i = 1; i <= n; i++)
        for(RG int j = 1; j <= n; j++)
            if(dis[i][j] <= lim) Add(i, j + n, INF);
    for(max_flow = 0; Bfs(); ) Copy(cur, fst), max_flow += Dfs(S, INF);
    return max_flow == num;
}
# undef int
int main(RG int argc, RG char* argv[]){
# define int ll
    n = Read(); m = Read(); Fill(dis, 63); T = n + n + 1;
    for(RG int i = 1; i <= n; i++)
        dis[i][i] = 0, s[i] = Read(), p[i] = Read(), num += s[i];
    for(RG int i = 1, a, b, c; i <= m; i++){
        a = Read(); b = Read(); c = Read();
        if(c >= dis[a][b]) continue;
        dis[a][b] = dis[b][a] = c;
    }
    for(RG int k = 1; k <= n; k++)
        for(RG int i = 1; i <= n; i++)
            for(RG int j = 1; j <= n; j++)
                dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
    RG int l = 0, r = INF, ans = -1;
    while(l <= r){
        RG int mid = (l + r) >> 1;
        if(Check(mid)) ans = mid, r = mid - 1;
        else l = mid + 1;
    }
    printf("%lld\n", ans);
    return 0;
}

转载于:https://www.cnblogs.com/cjoieryl/p/8206360.html


http://www.niftyadmin.cn/n/4390454.html

相关文章

mysql注入内置函数,mysql注入内置函数

a. 原理由于load_file()函数相当于*nix中的cat函数,而在某些FreeBSD,Sunos系统中用cat / 命令可以列出根目录,而字符/ 的Ascii码是47,所以在MYSQL注入的时候......SQL注入常用_计算机软件及应用_IT/计算机_专业资料。MYSQL: MYSQL注释符: 上午插入记录的时候一直没有成功,郁闷不…

聊一聊JWT与session

前言 认证和授权&#xff0c;其实吧简单来说就是:认证就是让服务器知道你是谁&#xff0c;授权就是服务器让你知道你什么能干&#xff0c;什么不能干&#xff0c;认证授权俩种方式&#xff1a;Session-Cookie与JWT&#xff0c;下面我们就针对这两种方案就行阐述。 Session 工作…

Top-1 accuracy和Top-5 accuracy的概念及理解

官方解释&#xff0c;也是我所查到了的最多的解释 top-1 就是你预测的label取最后概率向量里面最大的那一个作为预测结果&#xff0c;如果你的预测结果中概率最大的那个分类正确&#xff0c;则预测正确。否则预测错误 top-5 就是最后概率向量最大的前五名中&#xff0c;只要…

手机麦克风声音太大_我的iphone手机话筒声音特别小,如何调大?

phone手机话筒声2113音特别小的原因有如下5261几种&#xff1a;一、手机网络问题。可能是运4102营商或者手机的1653网络有问题&#xff0c;可以更换时间、地点、联系人&#xff0c;重新进行通话操作&#xff0c;确认是否是因为网络原因而导致的通话质量不佳。二、声音被调小。这…

ubuntu系统下将chrome的语言设置为中文,并主动翻译英文界面

1. 打开命令行输入 $ cd /usr/share/applications $ sh -c "export LANGUAGEZH-CN.UTF-8 && /usr/bin/google-chrome-stable %U &" 会弹出chrome界面&#xff0c;接下来在弹出的界面进行操作 2. 按照其他博主在windows系统下操作那般&#xff0c;继续…

hector slam matlab,如何使用hector slam算法包

文/冷冬寒梅下面即为具体步骤step1&#xff0c;建立一个ros的hector_ws 的workspace&#xff0c;并在其中建立src文件夹建立ros的workspace的命令如下$mkdir -p ~/hector_ws/srcstep2&#xff0c;进入src文件夹&#xff0c;下载hector slam代码step3&#xff0c;编译&#xff0…

【力扣JavaScript】1047. 删除字符串中的所有相邻重复项

/*** param {string} s* return {string}*/ var removeDuplicates function(s) {let stack[];for(i of s){let prevstack.pop();if(prev!i){stack.push(prev);stack.push(i);}}return stack.join(); };

信息论与编码实验报告——MATLAB实现算术编码

一、实验内容 试用MATLAB编制算术编码算法实现程序。 二、实验过程 2.1 算术编码实现原理 算术编码的算法思想如下&#xff1a; &#xff08;1&#xff09;对一组信源符号按照符号的概率从大到小排序&#xff0c;将[0,1)设为当前分析区间。按信源符号的概率序列在当前分析区间划…