ZOJ3496:Assignment——题解
http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3496
题目大意:A公司从S到T运货,每条路都有一个运货上限,而B公司则有p点费用点可以分配到个个路上,设分配到一条路上的费用点为x,该路走了y个货物,则需收费x*y,在以运货最大为要求下,B公司想要最大化收费,A公司则需要求出费用最大时的最小值(最小化费用)。
但是还没完,你还需要求出同条件下B公司想要最小化收费,A公司费用最小时的最大值(最大化费用)。
————————————————————————————
一道好题,顺便把我前面博客的bug找了出来。
首先需要想到,做到费用最大时,我们大可以把p点都加在流量最大的边上,这样一定是费用最大的(贪心策略)。
并且看到最大值最小,我们立刻想到了二分答案。
那么我们就可以二分所放p点的流量求解。
那么我们将所有边的容量都更新为min(容量,二分答案)重新跑一遍最大流,并且和之前的最大流相比,如果变小了的话就说明答案一定是小的,需要变大,否则答案就过大,需要变小。
此时我们已经做完了第一问,现在思考第二问。
我们用同样的思路二分,将所有路的容量增加一个我们二分答案为下界跑上下界网络流最大流,与原最大流比较即可,更新的方式和上面正好相反。
PS:一组数据如下:
1
5 11 2 1 2
3 4 12
4 1 11
0 1 8
2 3 11
4 0 12
3 1 9
2 4 11
4 3 8
0 3 6
2 0 8
2 1 11
ans:22 12
这组数据特别之处在于它有一条从u到v的边,同时还有从v到u的边,一旦有这种情况,只要在上下界网络流中处理不好就会出错。
显然我们只能删除超级源汇点(甚至你可以不删除,也不会对答案有什么影响),而不能将其他的边删除,否则就会造成残余网络中的流量无法留到汇点从而WA掉。
#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=;
const int M=;
const int INF=0x3f3f3f3f;
inline int read(){
int X=,w=;char ch=;
while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();
return w?-X:X;
}
struct node{
int nxt,to,w;
}edge[M<<];
struct line{
int u,v,w;
}e[M];
int head[N],du[N],cnt=-,ans,maxw;
int S,T,n,m,st,ed;
ll p;
inline void add(int u,int v,int w){
cnt++;
edge[cnt].to=v;
edge[cnt].w=w;
edge[cnt].nxt=head[u];
head[u]=cnt;
return;
}
int lev[N],cur[N],dui[N];
bool bfs(int k){
int r=;
for(int i=;i<=k;i++){
lev[i]=-;
cur[i]=head[i];
}
dui[]=S,lev[S]=;
int u,v;
for(int l=;l<=r;l++){
u=dui[l];
for(int e=head[u];e!=-;e=edge[e].nxt){
v=edge[e].to;
if(edge[e].w>&&lev[v]==-){
lev[v]=lev[u]+;
r++;
dui[r]=v;
if(v==T)return ;
}
}
}
return ;
}
int dinic(int u,int flow,int k){
if(u==k)return flow;
int res=,delta;
for(int &e=cur[u];e!=-;e=edge[e].nxt){
int v=edge[e].to;
if(edge[e].w>&&lev[u]<lev[v]){
delta=dinic(v,min(edge[e].w,flow-res),k);
if(delta>){
edge[e].w-=delta;
edge[e^].w+=delta;
res+=delta;
if(res==flow)break;
}
}
}
if(res!=flow)lev[u]=-;
return res;
}
//////////
bool check1(int k){
memset(head,-,sizeof(head));
cnt=-;
for(int i=;i<=m;i++){
add(e[i].u,e[i].v,min(k,e[i].w));
add(e[i].v,e[i].u,);
}
int awe=;S=st,T=ed;
while(bfs(n))awe+=dinic(S,INF,T);
return awe==ans;
}
ll solve1(){
int l=,r=maxw;
ll awe=;
while(l<=r){
int mid=(l+r)>>;
if(check1(mid)){
awe=mid;
r=mid-;
}else l=mid+;
}
return awe*p;
}
//////////////
bool check2(int k){
memset(head,-,sizeof(head));
memset(du,,sizeof(du));
cnt=-;
for(int i=;i<=m;i++){
if(e[i].w<k)return ;
add(e[i].u,e[i].v,e[i].w-k);
add(e[i].v,e[i].u,);
du[e[i].u]-=k;
du[e[i].v]+=k;
}
add(ed,st,INF);add(st,ed,);
S=n+;T=S+;
int awe=,full=;
for(int i=;i<=n;i++){
if(du[i]>){
add(S,i,du[i]);
add(i,S,);
full+=du[i];
}
else if(du[i]<){
add(i,T,-du[i]);
add(T,i,);
}
}
while(bfs(T))awe+=dinic(S,INF,T);
if(full!=awe)return ;
head[S]=head[T]=-;
S=st;T=ed;
awe=;
while(bfs(n))awe+=dinic(S,INF,T);
return awe==ans;
}
ll solve2(){
int l=,r=maxw;
ll awe=;
while(l<=r){
int mid=(l+r)>>;
if(check2(mid)){
awe=mid;
l=mid+;
}else r=mid-;
}
return awe*p;
}
int main(){
int t=read();
while(t--){
n=read(),m=read(),st=read()+,ed=read()+,p=read();
memset(head,-,sizeof(head));
maxw=;cnt=-;
for(int i=;i<=m;i++){
e[i].u=read()+,e[i].v=read()+,e[i].w=read();
add(e[i].u,e[i].v,e[i].w);
add(e[i].v,e[i].u,);
maxw=max(maxw,e[i].w);
}
S=st,T=ed;ans=;
while(bfs(n))ans+=dinic(S,INF,T);
printf("%lld %lld\n",solve1(),solve2());
}
return ;
}
ZOJ3496:Assignment——题解的更多相关文章
- ZOJ3496 Assignment
传送门 这题也是真恶心-- 题目大意是俩公司要运货,每条路有容量上限.然后B公司手里有p个--(技能点?)如果在一条路上放了x个技能点,这条路经过了y个货物,那么B公司就会收x*y的钱.现在要求的是, ...
- Lintcode208 Assignment Operator Overloading (C++ Only) solution 题解
[题目描述] Implement an assignment operator overloading method. Make sure that: The new data can be copi ...
- hdu4781 Assignment For Princess(构造)
题目链接:hdu4781 Assignment For Princess 题意:n个点m条边,每条有向边的权值分别是1,2,3…m,一个点能到达任意一个点,没有重边和自环,没有任何两条边的权值相同,任 ...
- 2017 google Round C APAC Test 题解
题解参考网上的答案,以及我自己的想法. 主要参考网站:http://codeforces.com/blog/entry/47181,http://codeforces.com/blog/entry/4 ...
- 2015 Multi-University Training Contest 1 题解&&总结
---------- HDU 5288 OO’s Sequence 题意 给定一个数列(长度<$10^5$),求有多少区间[l,r],且区间内有多少数,满足区间内其它数不是他的约数. 数的范围$ ...
- HDU 2853 Assignment(KM最大匹配好题)
HDU 2853 Assignment 题目链接 题意:如今有N个部队和M个任务(M>=N),每一个部队完毕每一个任务有一点的效率,效率越高越好.可是部队已经安排了一定的计划,这时须要我们尽量用 ...
- 题解——CodeForces 438D The Child and Sequence
题面 D. The Child and Sequence time limit per test 4 seconds memory limit per test 256 megabytes input ...
- HDU - 5289 Assignment (RMQ+二分)(单调队列)
题目链接: Assignment 题意: 给出一个数列,问其中存在多少连续子序列,使得子序列的最大值-最小值<k. 题解: RMQ先处理出每个区间的最大值和最小值(复杂度为:n×logn),相 ...
- HDU 5289 Assignment rmq
Assignment 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=5289 Description Tom owns a company and h ...
随机推荐
- node-redis使用记录
redis的高速存取性能让人印象深刻,虽然是分布式存储,但相比本地内存,性能毫不逊色. 之所以能做到这点,是由于redis的“单线程,多路复用IO”,同一时刻只有一个操作在进行. 而且多次建立从red ...
- hdu2147kiki's game(找规律)
kiki's game Time Limit: 5000/1000 MS (Java/Others) Memory Limit: 40000/10000 K (Java/Others)Total ...
- 提权基础-----mysql-udf提权
1.总结关于udf提权方法 通过弱口令,爆破,网站配置文件等方式得到mysql数据库帐号密码,---还要能外连 (1).将udf.dll代码的16进制数声明给my_udf_a变量 set @my_ud ...
- 使用unittest里面的discover()方法组织测试用例
import osimport unittest directory = os.getcwd()# 测试用例的目录organize = unittest.defaultTestLoader.disco ...
- Java开发工程师(Web方向) - 03.数据库开发 - 第5章.MyBatis
第5章--MyBatis MyBatis入门 Abstract: 数据库框架的工作原理和使用方法(以MyBatis为例) 面向对象的世界与关系型数据库的鸿沟: 面向对象世界中的数据是对象: 关系型数据 ...
- CsvHelper文档-2读
CsvHelper文档-2读 这个库默认不需要做任何设置就可以很容易的使用它.如果你的类属性名称直接匹配csv的标题名称,那么可以按照下面的实例来用: (以下所有的代码都需要引用using csvhe ...
- mysql数据库常用操作
目前最流行的数据库: oracle.mysql.sqlserver.db2.sqline --:单行注释 #:也是单行注释 /* 注释内容*/:多行注释 mysql -uroot -p密码:登录mys ...
- 概要梳理kafka知识点
主要是梳理一下kafka学习中的一些注意点,按照消息的流动方向进行梳理.详细的kafka介绍推荐看骑着龙的羊的系列博客,具体的某一块的知识点,可以参考我给出的一些参考文章. 1. kafka在系统中的 ...
- 【转载】2015Android 面试题 01
1.如何避免ANR? 答:ANR:Application Not Responding,五秒在Android中,活动管理器和窗口管理器这两个系统服务负责监视应用程序的响应. 当出现下列情况时,Andr ...
- sql STUFF 分组
---将sql想关的数据放到一个字段里 select r.Region_Name, )) [text()] from City c2 inner join Region as r2 on c2.Reg ...