BZOJ 1758 / Luogu P4292 [WC2010]重建计划 (分数规划(二分/迭代) + 长链剖分/点分治)
题意
分析
求这个平均值的最大值就是分数规划,二分一下就变成了求一条长度在[L,R]内路径的权值和最大.有淀粉质的做法但是我没写,感觉常数会很大.这道题可以用长链剖分做.
先对树长链剖分. 我们像做dsu on tree一样先做重儿子,用线段树继承重儿子的全部信息,然后做其他轻儿子
查询的时候枚举一下路径的长度len,一边单点O(1)O(1)O(1)查询长度为len的最大权值,一边线段树O(logn)O(logn)O(logn)查询长度为[L-len,R-len]的区间即可
时间复杂度我不会证(后面有)…反正枚举轻儿子的深度是总共O(n)O(n)O(n)的.所以加上外面的二分,总时间复杂度为O(nlog2n)O(nlog^2n)O(nlog2n).
Upd:Upd:Upd: 枚举轻儿子的深度的时间复杂度证明如下
时间复杂度O(n)O(n)O(n)。
分析如下:
每个点x只会暴力统计其所有轻儿子的信息,而每个轻儿子的信息大小为该轻儿子所在长链长度。
而当递归到x的父节点fa(x)时,若x不是fa(x)的重儿子,则fa(x)会暴力统计大小为x长链长度的信息。
故,每个长链只会对转移的复杂度做一次大小为其长度的贡献。
超强的长链剖分a!!! 证明摘自下面的博客.
学习长链剖分的看这里…博客传送门
CODE
#include<bits/stdc++.h>
using namespace std;
char cb[1<<15],*cs=cb,*ct=cb;
#define getc() (cs==ct&&(ct=(cs=cb)+fread(cb,1,1<<15,stdin),cs==ct)?0:*cs++)
template<class T>inline void read(T &res) {
char ch; int flg = 1; for(;!isdigit(ch=getc());)if(ch=='-')flg=-flg;
for(res=ch-'0';isdigit(ch=getc());res=res*10+ch-'0'); res*=flg;
}
template<class T>inline void chkmax(T &x, T y) { if(x < y) x = y; }
const int MAXN = 100005;
const double INF = 1e16;
int n, L, R, cnt, fir[MAXN], init_len[MAXN<<1];
struct edge { int to, nxt; double w; }e[MAXN<<1];
inline void add(int u, int v, int wt) {
e[cnt] = (edge) { v, fir[u], 0 }, init_len[cnt] = wt, fir[u] = cnt++;
e[cnt] = (edge) { u, fir[v], 0 }, init_len[cnt] = wt, fir[v] = cnt++;
}
int mxd[MAXN], dep[MAXN], son[MAXN], Eson[MAXN], dfn[MAXN], tmr;
void dfs1(int u, int ff) {
mxd[u] = dep[u] = dep[ff] + 1;
for(int i = fir[u], v; ~i; i = e[i].nxt)
if((v=e[i].to) != ff) {
dfs1(v, u); chkmax(mxd[u], mxd[v]);
if(mxd[v] > mxd[son[u]]) son[u] = v, Eson[u] = i;
}
}
void dfs2(int u, int ff) {
dfn[u] = ++tmr;
if(son[u]) dfs2(son[u], u);
for(int i = fir[u], v; ~i; i = e[i].nxt)
if((v=e[i].to) != ff && v != son[u])
dfs2(v, u);
}
int num[MAXN];
double dis[MAXN], mx[MAXN<<2];
void build(int i, int l, int r) {
mx[i] = -INF;
if(l == r) { num[l] = i; return; }
int mid = (l + r) >> 1;
build(i<<1, l, mid);
build(i<<1|1, mid+1, r);
}
void modify(int i, int l, int r, int x, double val) {
chkmax(mx[i], val);
if(l == r) return;
int mid = (l + r) >> 1;
if(x <= mid) modify(i<<1, l, mid, x, val);
else modify(i<<1|1, mid+1, r, x, val);
}
double query(int i, int l, int r, int x, int y) {
if(x > y) return -INF;
if(l == x && r == y) return mx[i];
int mid = (l + r) >> 1;
if(y <= mid) return query(i<<1, l, mid, x, y);
else if(x > mid) return query(i<<1|1, mid+1, r, x, y);
else return max(query(i<<1, l, mid, x, mid), query(i<<1|1, mid+1, r, mid+1, y));
}
double tmp[MAXN];
bool solve(int u, int ff) {
modify(1, 1, n, dfn[u], dis[u]);
if(son[u]) {
dis[son[u]] = dis[u] + e[Eson[u]].w;
if(solve(son[u], u)) return 1;
}
for(int i = fir[u], v; ~i; i = e[i].nxt)
if((v=e[i].to) != ff && v != son[u]) {
dis[v] = dis[u] + e[i].w;
if(solve(v, u)) return 1;
for(int j = 1; j <= mxd[v]-dep[u]; ++j) {
tmp[j] = mx[num[dfn[v]+j-1]];
if(j <= R) {
double temp = query(1, 1, n, dfn[u] + max(L-j, 0), dfn[u] + min(R-j, mxd[u]-dep[u]));
if(tmp[j] + temp - 2 * dis[u] >= 0) return 1;
}
}
for(int j = 1; j <= mxd[v]-dep[u]; ++j)
modify(1, 1, n, dfn[u]+j, tmp[j]);
}
return query(1, 1, n, dfn[u]+L, dfn[u]+min(R, mxd[u]-dep[u]))-dis[u] >= 0;
}
inline bool judge(double mid) {
for(int i = 0; i < cnt; ++i)
e[i].w = init_len[i] - mid;
build(1, 1, n);
return solve(1, 0);
}
int main() {
memset(fir, -1, sizeof fir);
read(n), read(L), read(R);
for(int i = 1, x, y, z; i < n; ++i)
read(x), read(y), read(z), add(x, y, z);
dfs1(1, 0), dfs2(1, 0);
double l = 0, r = 1e6, mid;
while(r - l > (1e-4)) {
mid = (l + r) / 2;
if(judge(mid)) l = mid;
else r = mid;
}
printf("%.3f\n", l);
}
CODE2
把二分转成迭代后直接洛谷rank 1了…
迭代就是随便取一个值作为mid,带进去算后得到更优的答案,然后把更优的答案作为mid继续迭代下去…直白点就是向答案逼近.
然后初值取一个适中的数就跑的超快了…感觉相当于没有了二分的那一个log
那分数规划都可以用迭代咯…
#include<bits/stdc++.h>
using namespace std;
char cb[1<<15],*cs=cb,*ct=cb;
#define getc() (cs==ct&&(ct=(cs=cb)+fread(cb,1,1<<15,stdin),cs==ct)?0:*cs++)
template<class T>inline void read(T &res) {
char ch; int flg = 1; for(;!isdigit(ch=getchar());)if(ch=='-')flg=-flg;
for(res=ch-'0';isdigit(ch=getchar());res=res*10+ch-'0'); res*=flg;
}
template<class T>inline void chkmax(T &x, T y) { if(x < y) x = y; }
const int MAXN = 100005;
const double INF = 1e16;
int n, L, R, cnt, fir[MAXN], init_len[MAXN<<1];
struct edge { int to, nxt; double w; }e[MAXN<<1];
inline void add(int u, int v, int wt) {
e[cnt] = (edge) { v, fir[u], 0 }, init_len[cnt] = wt, fir[u] = cnt++;
e[cnt] = (edge) { u, fir[v], 0 }, init_len[cnt] = wt, fir[v] = cnt++;
}
int mxd[MAXN], dep[MAXN], son[MAXN], Eson[MAXN], dfn[MAXN], tmr;
void dfs1(int u, int ff) {
mxd[u] = dep[u] = dep[ff] + 1;
for(int i = fir[u], v; ~i; i = e[i].nxt)
if((v=e[i].to) != ff) {
dfs1(v, u); chkmax(mxd[u], mxd[v]);
if(mxd[v] > mxd[son[u]]) son[u] = v, Eson[u] = i;
}
}
void dfs2(int u, int ff) {
dfn[u] = ++tmr;
if(son[u]) dfs2(son[u], u);
for(int i = fir[u], v; ~i; i = e[i].nxt)
if((v=e[i].to) != ff && v != son[u])
dfs2(v, u);
}
int num[MAXN];
double dis[MAXN];
struct node {
double x; int y; //x:最大路径和 y:边数
node(){}
node(double xx, int yy):x(xx), y(yy){}
inline bool operator <(const node &o)const { return x < o.x; }
}Ans, tmp[MAXN], mx[MAXN<<2];
void build(int i, int l, int r) {
mx[i] = node(-INF, 0);
if(l == r) { num[l] = i; return; }
int mid = (l + r) >> 1;
build(i<<1, l, mid);
build(i<<1|1, mid+1, r);
}
void modify(int i, int l, int r, int x, node val) {
chkmax(mx[i], val);
if(l == r) return;
int mid = (l + r) >> 1;
if(x <= mid) modify(i<<1, l, mid, x, val);
else modify(i<<1|1, mid+1, r, x, val);
}
node query(int i, int l, int r, int x, int y) {
if(x > y) return node(-INF, 0);
if(l == x && r == y) return mx[i];
int mid = (l + r) >> 1;
if(y <= mid) return query(i<<1, l, mid, x, y);
else if(x > mid) return query(i<<1|1, mid+1, r, x, y);
else return max(query(i<<1, l, mid, x, mid), query(i<<1|1, mid+1, r, mid+1, y));
}
void solve(int u, int ff) {
modify(1, 1, n, dfn[u], node(dis[u], dep[u]));
if(son[u]) dis[son[u]] = dis[u] + e[Eson[u]].w, solve(son[u], u);
for(int i = fir[u], v; ~i; i = e[i].nxt)
if((v=e[i].to) != ff && v != son[u]) {
dis[v] = dis[u] + e[i].w, solve(v, u);
for(int j = 1; j <= mxd[v]-dep[u]; ++j) {
tmp[j] = mx[num[dfn[v]+j-1]];
if(j <= R) {
node temp = query(1, 1, n, dfn[u] + max(L-j, 0), dfn[u] + min(R-j, mxd[u]-dep[u]));
chkmax(Ans, node(tmp[j].x + temp.x - 2 * dis[u], tmp[j].y + temp.y - 2*dep[u]));
}
}
for(int j = 1; j <= mxd[v]-dep[u]; ++j)
modify(1, 1, n, dfn[u]+j, tmp[j]);
}
node temp = query(1, 1, n, dfn[u]+L, dfn[u]+min(R, mxd[u]-dep[u]));
temp.x -= dis[u]; temp.y -= dep[u];
chkmax(Ans, temp);
}
inline node judge(double mid) {
for(int i = 0; i < cnt; ++i)
e[i].w = init_len[i] - mid;
build(1, 1, n);
Ans = node(-INF, 0);
solve(1, 0);
return Ans;
}
int main() {
memset(fir, -1, sizeof fir);
read(n), read(L), read(R);
for(int i = 1, x, y, z; i < n; ++i)
read(x), read(y), read(z), add(x, y, z);
dep[0] = -1; dfs1(1, 0), dfs2(1, 0);
double mid = 600000;
while(1) { //高级操作
node now = judge(mid);
double ans = (now.x + mid*now.y) / now.y; //把y条边减去的mid加回来,再算一次平均值
if(fabs(ans-mid) < 1e-3) break;
mid = ans;
}
printf("%.3f\n", mid);
}
BZOJ 1758 / Luogu P4292 [WC2010]重建计划 (分数规划(二分/迭代) + 长链剖分/点分治)的更多相关文章
- bzoj 1758 [Wc2010]重建计划 分数规划+树分治单调队列check
[Wc2010]重建计划 Time Limit: 40 Sec Memory Limit: 162 MBSubmit: 4345 Solved: 1054[Submit][Status][Disc ...
- BZOJ.1758.[WC2010]重建计划(分数规划 点分治 单调队列/长链剖分 线段树)
题目链接 BZOJ 洛谷 点分治 单调队列: 二分答案,然后判断是否存在一条长度在\([L,R]\)的路径满足权值和非负.可以点分治. 对于(距当前根节点)深度为\(d\)的一条路径,可以用其它子树深 ...
- BZOJ1758[Wc2010]重建计划——分数规划+长链剖分+线段树+二分答案+树形DP
题目描述 输入 第一行包含一个正整数N,表示X国的城市个数. 第二行包含两个正整数L和U,表示政策要求的第一期重建方案中修建道路数的上下限 接下来的N-1行描述重建小组的原有方案,每行三个正整数Ai, ...
- [WC2010]重建计划(分数规划+点分治+单调队列)
题目大意:给定一棵树,求一条长度在L到R的一条路径,使得边权的平均值最大. 题解 树上路径最优化问题,不难想到点分治. 如果没有长度限制,我们可以套上01分数规划的模型,让所有边权减去mid,求一条路 ...
- 洛谷 P4292 [WC2010]重建计划 解题报告
P4292 [WC2010]重建计划 题目描述 \(X\)国遭受了地震的重创, 导致全国的交通近乎瘫痪,重建家园的计划迫在眉睫.\(X\)国由\(N\)个城市组成, 重建小组提出,仅需建立\(N-1\ ...
- 「WC2010」重建计划(长链剖分/点分治)
「WC2010」重建计划(长链剖分/点分治) 题目描述 有一棵大小为 \(n\) 的树,给定 \(L, R\) ,要求找到一条长度在 \([L, R]\) 的路径,并且路径上边权的平均值最大 \(1 ...
- 洛谷 P4292 - [WC2010]重建计划(长链剖分+线段树)
题面传送门 我!竟!然!独!立!A!C!了!这!道!题!incredible! 首先看到这类最大化某个分式的题目,可以套路地想到分数规划,考虑二分答案 \(mid\) 并检验是否存在合法的 \(S\) ...
- 【BZOJ 1758】【WC 2010】重建计划 分数规划+点分治+单调队列
一开始看到$\frac{\sum_{}}{\sum_{}}$就想到了01分数规划但最终还是看了题解 二分完后的点分治,只需要维护一个由之前处理过的子树得出的$tb数组$,然后根据遍历每个当前的子树上的 ...
- P4292 [WC2010]重建计划 点分治+单调队列
题目描述 题目传送门 分析 看到比值的形式就想到 \(01分数规划\),二分答案 设当前的值为 \(mids\) 如果存在\(\frac{\sum _{e \in S} v(e)}{|S|} \geq ...
随机推荐
- 将本地镜像上传到Docker镜像仓库
⒈在Docker Hub(Docker官方镜像仓库)注册自己的用户名. ⒉在Docker中使用docker login命令登录自己的用户名密码. ⒊使用docker push 镜像名称 上传本地镜像到 ...
- java中的 |=、&=、^=
|= 关于 |= 运算符:|= 运算符和 += 这一类的运算符一样,拆解开就是 a = a | b: 代码如下: public static strictfp void main(String[] ...
- Jmeter之参数化(4种设置方法)
以多用户登录为例~~~ 参数化: 1.用户参数 2.CSV数据文件 3.函数助手CSVRead 4.用户自定义的变量 1.用户参数 脚本目录结构如下: 因为设置了2组账号密码,所以线程数设置为2(添加 ...
- 怎样在 Vue 里面使用自定义事件将子组件的数据传回给父组件?
首先, Vue 里面的组件之间的数据流动是 单向 的, 数据可以从父组件传递给子组件, 但不能从子组件传递给父组件, 因为组件和组件之间是 隔离 的. 就像两个嵌套的 黑盒子 . 能通过 props ...
- java 字节流与字符流的区别(转)
字节流与和字符流的使用非常相似,两者除了操作代码上的不同之外,是否还有其他的不同呢? 实际上字节流在操作时本身不会用到缓冲区(内存),是文件本身直接操作的,而字符流在操作时使用了缓冲区,通过缓冲区再操 ...
- 清除SQL日志文件
1.清除errorlog文件 MSSQL在 C:\Program Files\Microsoft SQL Server\MSSQL.1\MSSQL\LOG 目录下存放这一些日志文件,一共是7个,常常会 ...
- C#进阶之WebAPI(一)
最近出去面试,被问到关于WebAPI的知识,因为项目中没有单独写过WebAPI,使用的时候是和mvc结合在一起使用的,所以,在我的印象中WebAPI和mvc是差不多的,这种答案当然不能让人满意了,于是 ...
- js点击发送验证码 xx秒后重新发送
用于一些注册类的场景,点击发送验证码,xx秒后重新发送. 利用 setTimeout 方法,xx秒后执行指定的方法,修改button的属性值,disabled为true时为灰色,不可点击. <! ...
- day10 mysql常用操作
一. 目录 1.mysql的简介 2.增删改查操作 3.单表查询 4.多表查询常见的三种方式 5.pymysql模块操作数据库 二. 内容 一.mysql的简介 概述:mysql是一个关系型数据库, ...
- 8 java 笔记
1 import语句可以简化编程,可以导入指定包下面的某个类或者全部类 2 java.lang,Object类是所有类的父类,要么是其直接父类,要么是其间接父类 3 子类扩展了父类,子类是一种特殊的父 ...