洛谷P4382 [八省联考2018]劈配(网络流,二分答案)
说不定比官方sol里的某理论最优算法还优秀一点?
所以\(n,m\)说不定可以出到\(1000\)?
无所谓啦,反正是个得分题。Orz良心出题人,暴力有70分2333
思路分析
正解的思路很巧妙,其实我并不觉得这是个正儿八经的网络流或者二分图匹配的题目,主要还是个思维+建图模型+乱搞。。。。。。
\(C=1\)时我们就可以对于每个人直接匹配而不会影响到后面的选择了。但是\(C>1\)的话,可能某一个人可以选多个导师,当他随便选了其中一个以后,可能影响到后面某个人使其选不到本来的最优解。而此时后面那个人就要让前面那个人改变决策,做出反悔。
这时候有没有想到网络流算法的反悔边的应用呢?(因为蒟蒻只会网络流所以就用网络流来跑二分图匹配)建一个二分图,左边\(n\)个点代表选手,右边\(m\)个点代表导师(战队),一开始所有导师向汇点\(T\)连流量为战队人数上限的边。每次从前往后枚举选手,找到能够选择的最高志愿,动态地向该档志愿能选择的所有导师连容量为\(1\)的边,然后增广一次(这时会随意选择一个可行的导师)。因为增广后会留下反悔边,所以这样做就很方便啦。因为每次是动态加边直接从每个选手开始增广 ,所以貌似连源点都不需要(这就是蒟蒻一直觉得这题像一个乱搞题的原因)
这时候留下一个问题,如何方便地判断某时刻某战队是否已经不可选了呢?暴力存图实在是太麻烦啦!而且我们并不需要知道整个图是什么样子,只要知道这一些bool型状态,暴力存图实在是浪费时空。大概算法复杂度优化的瓶颈就在这里吧。一开始蒟蒻yy各种极其错误的方法来乱搞。后来发现,将问题带入到我们建的网络流模型里面,不就是要知道每个代表导师的点到汇点是否存在增广路吗?具体方法,从汇点开始dfs搜索,如果某条边的反向边有流量就搜过去,把搜到的点标记起来就好啦!
对于第二问,显然可以想到二分。先把每次每个战队能否选择的状态保存下来,对于每个选手,二分他之前的排名,根据是否存在一个可以选择而且令他满意的战队来决定下面往哪一半区间接着找。
复杂度分析
有三步。首先要暴力找到能满足的最高志愿,复杂度\(O(nm)\)
接着,增广(蒟蒻直接dfs去找而不是dinic)和更新当前战队是否可选状态,因为是在图上搞,边数不超过\(O(nC)\),所以这部分的复杂度都是\(O(n^2C)\)的(据说在二分图上使用dinic可以使增广的复杂度降到\(O(n\sqrt nC)\)?但是这是动态加边,每次都要跑,那应该是不能优化的吧)
最后,二分最少上升的排名,复杂度\(O(nm\log n)\)
综上,因为n与m同阶,所以总复杂度应该是\(O(n^2(C+\log n))\)级别的。由于是网络流,所以复杂度非常玄学,也不至于删边那么麻烦,蒟蒻随意加了一点点剪枝,还有register,inline,快读快写,就只要8ms,开O2 4ms。居然没有0msTAT,欢迎超越
最后提一个惨痛的经历:蒟蒻调这题一天多了,所有的算法性错误都更正了,还是WA,最后才发现重建图的时候反向边的流量忘记赋成0来覆盖掉原来的信息了TAT
#include<cstdio>
#include<cstdlib>
#include<cstring>
#define R register int
#define I inline void
#define pc(C) *po++=C
#define add(X,Y,F)\
to[++p]=Y,ne[p]=he[X],f[he[X]=p]=F,\
to[++p]=X,ne[p]=he[Y],f[he[Y]=p]=0//反边赋0!
const int SZ=1<<20,N=209,M=N<<1,L=N*22;//边数最多2(nC+m)
char ibuf[SZ],obuf[SZ],*pi=ibuf-1,*po=obuf;
I in(R&x){
while(*++pi<'-');
x=*pi&15;
while(*++pi>'-')x*=10,x+=*pi&15;
}
I out(R x){
if(x>9)out(x/10);
pc(x%10|'0');
}
int T,he[M],ne[L],to[L],f[L],he1[N],ne1[N];
int a[N][N],b[N],up[N];//up记录最少上升排名
bool c[N][N],vis[M];//c存下战队是否可选
void remain(R x){//更新状态,利用好反向边!
vis[x]=1;
for(R i=he[x];i;i=ne[i])
if(f[i^1]&&!vis[to[i]])remain(to[i]);
}
bool flow(R x){//dfs增广
if(vis[x])return x==T;
vis[x]=1;
for(R i=he[x];i;i=ne[i])
if(f[i]&&flow(to[i])){
--f[i];++f[i^1];
return 1;
}
return 0;
}
int main(){
fread(ibuf,1,SZ,stdin);
R TT,CC,n,m,p,i,j,x,s,l,r,mi;
in(TT);in(CC);
while(TT--){
in(n);in(m);
T=n+m+1;p=vis[T]=1;
for(j=1;j<=m;++j)
in(b[j]),add(j+n,T,b[j]);
for(i=1;i<=n;++i)
for(j=1;j<=m;++j)in(a[i][j]);
for(j=1;j<=m;++j)
c[1][j]=(bool)b[j];//初始状态就是收不收人(不会来个导师不收人吧233)
for(i=1;i<=n;++i){
memset(he1,0,(n+1)<<2);
for(j=1;j<=m;++j)//链表把同等级志愿挂一起,像邻接表一样
ne1[j]=he1[a[i][j]],he1[a[i][j]]=j;
for(x=1;x<=m;++x){
for(j=he1[x];j;j=ne1[j])
if(c[i][j])break;
if(j){
for(j=he1[x];j;j=ne1[j])
if(c[i][j])add(i,j+n,1);//贪心选择,动态连边
memset(vis,0,T);flow(i);
break;
}
}
out(x);pc(' ');
if(x<=m){
memset(vis,0,T);remain(T);
memcpy(c[i+1]+1,vis+n+1,m);
}//没选到肯定状态不变啦,直接copy上一次的
else memcpy(c[i+1]+1,c[i]+1,N);
in(s);if(x<=s){up[i]=0;continue;}//判掉一开始就满意的
l=0;r=i-1;
F:while(l!=r){
mi=(l+r+1)>>1;//注意二分区间是左闭右开的,所以要写成这样
for(x=s;x;--x)
for(j=he1[x];j;j=ne1[j])
if(c[mi][j]){l=mi;goto F;}
r=mi-1;
}
up[i]=i-l;
}
pc('\n');for(i=1;i<=n;++i){out(up[i]);pc(' ');}pc('\n');
memset(he,0,(T+1)<<2);//应该只有he需要清空
}
fwrite(obuf,1,po-obuf,stdout);
return 0;
}
洛谷P4382 [八省联考2018]劈配(网络流,二分答案)的更多相关文章
- P4382 [八省联考2018]劈配
题目链接 题意分析 受到了\(olinr\ \ julao\)的影响 写了匈牙利算法 首先 我们对于每一个人 从高到低枚举志愿 如果当前志愿的老师有剩余的话 那么我们就选 否则的话 我们看看谁的那个志 ...
- BZOJ.5251.[八省联考2018]劈配mentor(最大流)
BZOJ 洛谷 对于每个人,每次枚举一个志愿看是否能增广即可. 对于第二问,可以保留第一问中\(n\)次增广前后的\(n\)张图,二分,在对应图上看是否能增广即可. 貌似匈牙利的某种写法比网络流优多了 ...
- [BZOJ5251][九省联考2018]劈配(网络流)
5251: [2018多省省队联测]劈配 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 33 Solved: 22[Submit][Status][ ...
- 洛谷P4383 [八省联考2018]林克卡特树lct(DP凸优化/wqs二分)
题目描述 小L 最近沉迷于塞尔达传说:荒野之息(The Legend of Zelda: Breath of The Wild)无法自拔,他尤其喜欢游戏中的迷你挑战. 游戏中有一个叫做“LCT” 的挑 ...
- 洛谷 4383 [八省联考2018]林克卡特树lct——树形DP+带权二分
题目:https://www.luogu.org/problemnew/show/P4383 关于带权二分:https://www.cnblogs.com/flashhu/p/9480669.html ...
- 洛谷.4383.[八省联考2018]林克卡特树lct(树形DP 带权二分)
题目链接 \(Description\) 给定一棵边带权的树.求删掉K条边.再连上K条权为0的边后,新树的最大直径. \(n,K\leq3\times10^5\). \(Solution\) 题目可以 ...
- dp凸优化/wqs二分学习笔记(洛谷4383 [八省联考2018]林克卡特树lct)
qwq 安利一个凸优化讲的比较好的博客 https://www.cnblogs.com/Gloid/p/9433783.html 但是他的暴力部分略微有点问题 qwq 我还是详细的讲一下这个题+这个知 ...
- [八省联考2018] 劈配 mentor
Description 一年一度的综艺节目<中国新代码>又开始了.Zayid 从小就梦想成为一名程序员,他觉得这是一个展示自己的舞台,于是他毫不犹豫地报名了. Input 轻车熟路的Zay ...
- BZOJ5251 八省联考2018劈配(网络流)
劈配,匹配,网络流.那么考虑怎么跑网络流. 先看第一问.首先套路的建出超源超汇.不用想也知道导师向汇连容量为战队人数上限的边.特别地,给出局也建一个点,向汇连容量inf的边(似乎没有必要).对于一个新 ...
随机推荐
- array_map,array_filter,array_walk区别
这几个方法都有遍历数组,操作的过程. 很容易搞混.尤其很多人养成的编程习惯都是使用foreach来循环遍历.. 就举个简单例子. $a = array( 0,1,2,3,4,5,6,7,8,9 ); ...
- [SCOI2007]修车 BZOJ1070
很久之前写的题了,今天翻出来写一篇博客复习一下... 分析: 考虑,T <= 1000,并不能针对这一维处理,所以考虑将,每个人拆点,之后,拆完之后表示,这个人第n-j+1个修k这辆车,也就是, ...
- Django Rest Framework源码剖析(一)-----认证
一.简介 Django REST Framework(简称DRF),是一个用于构建Web API的强大且灵活的工具包. 先说说REST:REST是一种Web API设计标准,是目前比较成熟的一套互联网 ...
- 20155223 实验5 MSF基础应用
20155223 实验5 MSF基础应用 基础问题回答 用自己的话解释什么是exploit,payload,encode? exploit:漏洞攻击.一个exploit程序肯定会触发系统的一个或多个漏 ...
- 20155328 《网络对抗》 实验八:Web基础
20155328 <网络对抗> 实验八:Web基础 实验内容及过程记录 一.Web前端HTML 我们的kali是默认安装好了apache的.首先输入netstat -tupln |grep ...
- SVD(奇异值分解)Python实现
注:在<SVD(奇异值分解)小结 >中分享了SVD原理,但其中只是利用了numpy.linalg.svd函数应用了它,并没有提到如何自己编写代码实现它,在这里,我再分享一下如何自已写一个S ...
- dubbo见解
调用关系说明 服务容器负责启动,加载,运行服务提供者. 服务提供者在启动时,向注册中心注册自己提供的服务. 服务消费者在启动时,向注册中心订阅自己所需的服务. 注册中心返回服务提供者地址列表给消费者, ...
- Git管理多个SSH密钥,Git多帐号配置
版权声明:转载须标明本文转自严振杰的博客 https://blog.csdn.net/yanzhenjie1003/article/details/69487932版权声明:转载必须注明本文转自严振杰 ...
- Golang Context 详细介绍
Golang context 本文包含对context实现上的分析和使用方式,分析部分源码讲解比价多,可能会比较枯燥,读者可以直接跳过去阅读使用部分. ps: 作者本着开源分享的精神撰写本篇文章,如果 ...
- LeetCode-765.情侣牵手
N 对情侣坐在连续排列的 2N 个座位上,想要牵到对方的手. 计算最少交换座位的次数,以便每对情侣可以并肩坐在一起. 一次交换可选择任意两人,让他们站起来交换座位. 人和座位用 0 到 2N-1 的整 ...