CF786E ALT
题意
有一棵 \(n\) 个点的树和 \(m\) 个人,第 \(i\) 个人从 \(u_i\) 走到 \(v_i\)
现在要发宠物,要求一个人要么他自己发到宠物,要么他走的路径上的都有宠物。
求最小代价,并输出任意方案。
\(n,m \leq 20000\)
传送门
思路
对每个人和每条树边都建一个点。
源点向每个人连容量 \(1\) 的边,每条树边向汇点连容量 \(1\) 的边。每个人向他要走到的所有边连容量 \(+\infty\) 的边。
给人发就是割掉人与源的边,放边上就是割树边与汇的边,人与树边间的边不能割。当源汇不连通的时候,就是满足题意的,题目转化为最小割。
问题出在连边上,这是 \(nm\) 的,考虑如何优化连边。线段树?树上问题有点麻烦,想不到倍增。
对于一段相邻的\(2^i\)个点,都建一个点与它们所有相连,其实又有点像线段树,然后就会对于每条链得到一棵类似线段树的东西?大区间连向小区间,看起来就是。
然后对每个人,将他的路径在 LCA 处分成两条路径,这两条路径分别向对应区间覆盖。
那么我们得到的新图有 \(n\log n+m\) 个点和 \(n\log n+4m\) 条边,可以通过。
这里还要输出方案。
我们知道,最小割中的边一定满流。(全局最大流,分成两个集合,其间的边一定是满流的,而也就是最小割)因此我们从源点开始 dfs,只走没满流的边,并标记被 dfs 到的点。则图被分成两部分,一部分被访问过,一部分没被访问过。其中间那些边就是一个最小割。
然后,如果一个人代表的点没被访问过,则说明他所属那条边被割了。如果一条树边代表的点被访问过,则说明它所属那条边被割了。记录标号输出
#include <bits/stdc++.h>
using std::queue;
const int W=14,N=20005*(W+1),M=N+N*4;
int to[M<<1],w[M<<1],Next[M<<1],edge,n,x,y,f[20005][W+1],idn[N],deep[N],last[N],b[20005][W+1];
int s,t,m,cnt,ans,tag[N],cur[N];
queue <int> q;
void add(int x,int y,int z){
to[++edge]=y;
Next[edge]=last[x];
last[x]=edge;
w[edge]=z;
}
void dfs(int x,int fa){
f[x][0]=fa,deep[x]=deep[fa]+1;
for (int i=last[x];i;i=Next[i])
if (to[i]!=fa) {
idn[to[i]]=(i+1)/2;
dfs(to[i],x);
}
}
void Add(int x,int y,int w){
add(x,y,w);
add(y,x,0);
}
int lca(int x,int y){
if (deep[x]<deep[y]) std::swap(x,y);
for (int i=W;i>=0;i--)
if (deep[f[x][i]]>=deep[y]) x=f[x][i];
if (x==y) return x;
for (int i=W;i>=0;i--)
if (f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
return f[x][0];
}
void addedge(int x,int y,int id){
for (int i=W;i>=0;i--)
if (deep[f[x][i]]>=deep[y]){
Add(id,b[x][i],n+1);
x=f[x][i];
}
}
bool bfs(){
for (int i=0;i<=t;i++) cur[i]=last[i],deep[i]=0;
deep[s]=1;
q.push(s);
while (!q.empty()){
int x=q.front();
q.pop();
for (int i=last[x];i;i=Next[i])
if (w[i] && !deep[to[i]]){
deep[to[i]]=deep[x]+1;
q.push(to[i]);
}
}
return (deep[t]);
}
int c(int x) { return x&1?x+1:x-1; }
int dfs(int x){
if (x==t) return 1;
for (int i=cur[x];i;i=Next[i]){
cur[x]=i;
int u=to[i];
if (deep[u]>deep[x] && w[i]){
int di=dfs(to[i]);
if (di){
w[i]--;
w[c(i)]++;
return 1;
}
}
}
return 0;
}
void DFS(int x){
tag[x]=1;
for (int i=last[x];i;i=Next[i])
if (w[i] && !tag[to[i]]) DFS(to[i]);
}
int main(){
scanf("%d%d",&n,&m);
for (int i=1;i<n;i++){
scanf("%d%d",&x,&y);
add(x,y,0),add(y,x,0);
}
s=0;
dfs(1,0);
edge=0;memset(last,0,sizeof(last));
cnt=n;
for (int i=2;i<=n;i++) b[i][0]=i;
for (int i=1;i<=W;i++)
for (int j=1;j<=n;j++)
if (f[f[j][i-1]][i-1]){
f[j][i]=f[f[j][i-1]][i-1];
Add(++cnt,b[j][i-1],n+1);
Add(cnt,b[f[j][i-1]][i-1],n+1);
b[j][i]=cnt;
}
for (int i=1;i<=m;i++){
scanf("%d%d",&x,&y);
Add(s,++cnt,1);
idn[cnt]=i;
int l=lca(x,y);
addedge(x,l,cnt),addedge(y,l,cnt);
}
t=cnt+1;
for (int i=2;i<=n;i++) Add(i,t,1);
while (bfs())
while (dfs(s)==1) ans++;
printf("%d\n",ans);
DFS(s);
ans=0;
for (int i=last[s];i;i=Next[i])
if (tag[to[i]]!=tag[s]) ans++;
printf("%d ",ans);
for (int i=last[s];i;i=Next[i])
if (tag[to[i]]!=tag[s]) printf("%d ",idn[to[i]]);
puts("");
ans=0;
for (int i=last[t];i;i=Next[i])
if (tag[to[i]]!=tag[t]) ans++;
printf("%d ",ans);
for (int i=last[t];i;i=Next[i])
if (tag[to[i]]!=tag[t]) printf("%d ",idn[to[i]]);
return 0;
}
后记
追随神仙的脚步
抄袭于此Mrsrz
CF786E ALT的更多相关文章
- cf786E ALT (最小割+倍增优化建图)
如果把“我全都要”看作是我全不要的话,就可以用最小割解决啦 源点S,汇点T 我们试图让每个市民作为一个等待被割断的路径 把狗狗给市民:建边(S,i,1),其中i是市民 把狗狗给守卫:建边(j,T,1) ...
- 一句话题解&&总结
CF79D Password: 差分.两点取反,本质是匹配!最短路+状压DP 取反是套路,匹配是发现可以把操作进行目的化和阶段化,从而第二次转化问题. 且匹配不会影响别的位置答案 sequence 计 ...
- [No000093]按住Alt 再按数字键敲出任意汉字和字符!
1.在notepad里,(中文系统下) 按住Alt 然后按52946最后放开Alt 按住Alt 然后按45230最后放开Alt 按住Alt 然后按50403最后放开Alt 你会看到"我爱你& ...
- [No00008B]远程桌面发送“Ctrl+Alt+Delete”组合键调用任务管理器
向远程桌面发送"Ctrl+Alt+Delete"组合键的两种方法 1.在本地按下Ctrl+Alt+End,可以成功发送"Ctrl+Alt+Delete"组合键! ...
- 平常看到的Alt+xx 快捷键用法
1. 先按Alt, 哪一个菜单对应的字符是有划线的. 2. 输入对应的字符打开相应的菜单, 3 再输入相应的字符打开子菜单
- windows 中去除Ctrl+Alt+Del才能登录
安装windows 7后登录的时候有一样很麻烦的步骤是需要先按Ctrl+Alt+Del,才能输入用户密码进行登录.这里笔者介绍一下如何取消这个东西. 点击“开始菜单”,点击“控制面板”. [管理工具] ...
- TSql 巧用Alt 键
1,查看表的信息 在TSql 编辑器中,选中一个表,如图 点击Alt+F1,就可以查看表的属性定义 2,使用alt批量插入逗号 在Tsql中使用 in 子句,在(value_List)列表中,经常有很 ...
- title与alt的区别
html中的title属性和alt属性让人有些混淆. 以前不知道有title这个属性,第一次用到它时,就和alt产生了混淆.一位朋友告诉我说,alt是图片img标签里用的,title是超链接里用的, ...
- UML序列图总结(Loop、Opt、Par和Alt)
序列图主要用于展示对象之间交互的顺序. 序列图将交互关系表示为一个二维图.纵向是时间轴,时间沿竖线向下延伸.横向轴代表了在协作中各独立对象的类元角色.类元角色用生命线表示.当对象存在时,角色用一条虚线 ...
随机推荐
- MySQL存储引擎的介绍
数据库存储引擎是数据库底层软件组件,数据库管理系统使用数据引擎进行创建.查询.更新和删除数据操作.不同的存储引擎提供不同的存储机制.索引技巧.锁定水平等功能,使用不同的存储引擎还可以获得特定的功能. ...
- Fortify漏洞之Sql Injection(sql注入)
公司最近启用了Fortify扫描项目代码,报出较多的漏洞,安排了本人进行修复,近段时间将对修复的过程和一些修复的漏洞总结整理于此! 本篇先对Fortify做个简单的认识,同时总结一下sql注入的漏洞! ...
- Vue组件component创建及使用
组件化与模块化的区别 什么是组件:组件的出现,就是为了拆分Vue实例的代码量,能够让我们以不同的组件,来划分不同的功能模块 ,将来我们需要什么功能,就可以去调用对应的组件即可 组件化与模块化的不同: ...
- SQL SERVER-端口Port
Quick cheat sheet for port numbers used by SQL Server services or services that SQL Server may depen ...
- Linux之redis的安装,主从配置
一,redis安装 redis是一个key-value存储系统.和Memcached类似,它支持存储的value类型相对更多,包括string(字符串).list(链表).set(集合).zset(s ...
- 网络编程socketserver
一.网络编程回顾 tcp是流式传输,字节流,数据与数据之间没有边界 优点:不限定长度,可靠传输 缺点:慢,和一端的通信连接conn会一直占用通信资源 udp协议式面向数据包的传输 优点:快,由于不需要 ...
- sed初级教程
简介 sed是作为特殊目的的编辑器而创建的,用于专门执行脚本:与ed不同,它不能交互地使用.sed面向字符流.默认情况下,到sed的所有输入都会经过相应的处理,并转为标 准输出.输入文件本身不发生改变 ...
- Linux命令——finger
简介 查询并显示系统用户的相关信息. 最小化安装Linux可能没有改名了,需要单独安装. RHEL/CentOS yum install finger* -y Ubuntu apt-get inst ...
- JSON 的简介与使用
一.什么是JSON? JSON(javascript object notation)全称是 javascript 对象表示法,是一种数据交换的文本格式,用于读取结构化数据,提出目的是取代繁琐笨重的 ...
- 双端循环列表实现栈(python)
# -*- coding: utf-8 -*- from collections import deque class Node(object): def __init__(self, value=N ...