bzoj 3672 利用点分治将CDQ分治推广到树型结构上
最大的收获就是题目所说。
deal(s) : 处理节点s所在块的问题,并保证:
1、s是该块中最靠近根节点的点,没有之一。
2、s所在块到根节点的路径上的点全都用来更新过了s所在块的所有节点。
然后步骤是:
1、找s所在块的重心c。
2、如果s就是c,那么用c更新当前块的所有节点,然后“删除c”,递归处理新产生的子块。
3、否则,删除c,deal(s),用c到s的路径(不包括c,包括s)更新c除了s子块的其他子块以及c,然后再用c去更新一次。
4、递归处理其它子块。
#include <cstdio>
#include <cassert>
#include <cstring>
#include <algorithm>
#define N 200010
#define M N<<1
#define fill(arr,lf,rg,v) memset(arr+lf,v,sizeof(arr[0])*(rg-lf+1))
using namespace std; typedef long long dnt;
struct Vector {
dnt x, y;
Vector(){}
Vector( dnt x, dnt y ):x(x),y(y){}
Vector operator+( const Vector &o ) const { return Vector(x+o.x,y+o.y); }
Vector operator-( const Vector &o ) const { return Vector(x-o.x,y-o.y); }
double operator^( const Vector &o ) const { return (double)x*o.y-(double)y*o.x; }
};
typedef Vector Point;
bool onleft( const Point &a, const Point &b, const Point &c ) {
return ((b-a)^(c-a)) >= 0.0;
}
struct Convex {
Point stk[N];
int top;
void init() {
top=-;
}
inline void append( const Point &p ) {
while( top> && onleft(stk[top-],stk[top],p) ) top--;
stk[++top] = p;
}
const Point& query( dnt k ) {
int lf=;
int rg=top;
if( lf==rg ) return stk[top];
assert(k*(stk[lf].x-stk[lf+].x)>=);
if( k*(stk[lf].x-stk[lf+].x) >= (stk[lf].y-stk[lf+].y) )
return stk[lf];
assert(k*(stk[rg-].x-stk[rg].x)>=);
if( k*(stk[rg-].x-stk[rg].x) <= (stk[rg-].y-stk[rg].y) )
return stk[rg];
lf++;
rg--;
while( lf<rg ) {
int mid=(lf+rg)>>;
assert(k*(stk[mid].x-stk[mid+].x)>=);
if( k*(stk[mid].x-stk[mid+].x) > (stk[mid].y-stk[mid+].y) )
rg=mid;
else
lf=mid+;
}
return stk[lf];
}
}; int n, case_type;
int head[N], next[M], dest[M], etot;
dnt wp[N], wq[N], lim[N], ws[M];
int anc[N], fat[N], vis[N], siz[N], bac[N];
dnt dep[N], dp[N];
int qu[N], bg, ed;
Convex convex; void adde( int u, int v, dnt s ) {
etot++;
dest[etot] = v;
next[etot] = head[u];
ws[etot] = s;
head[u] = etot;
}
void bfs( int s ) {
qu[bg=ed=] = s;
anc[s] = ;
dep[s] = ;
while( bg<=ed ) {
int u=qu[bg++];
for( register int t=head[u]; t; t=next[t] ) {
int v=dest[t];
if( v==anc[u] ) continue;
qu[++ed] = v;
anc[v] = u;
dep[v] = dep[u]+ws[t];
}
}
}
int getc( int s ) {
qu[bg=ed=] = s;
fat[s] = ;
bac[s] = ;
siz[s] = ;
while( bg<=ed ) {
int u=qu[bg++];
for( register int t=head[u]; t; t=next[t] ) {
int v=dest[t];
if( vis[v] || v==fat[u] ) continue;
qu[++ed] = v;
fat[v] = u;
bac[v] = ;
siz[v] = ;
}
}
int c = ;
for( register int i=ed; i>=; i-- ) {
int u=qu[i];
siz[u]++;
if( fat[u] ) {
int f=fat[u];
siz[f]+=siz[u];
if( siz[u]>bac[f] ) bac[f]=siz[u];
}
}
for( register int i=; i<=ed; i++ ) {
int u=qu[i];
if( bac[u]<siz[s]-siz[u] ) bac[u]=siz[s]-siz[u];
if( !c || bac[u]<bac[c] ) c=u;
}
//fprintf( stderr, "%d\n", c );
return c;
}
void flood( int s ) {
qu[bg=ed=] = s;
fat[s] = ;
while( bg<=ed ) {
int u=qu[bg++];
for( register int t=head[u]; t; t=next[t] ) {
int v=dest[t];
if( vis[v] || v==fat[u] ) continue;
qu[++ed] = v;
fat[v] = u;
}
}
}
inline void update( const Point &p, int u ) {
dnt v = p.y + wp[u]*(dep[u]-p.x) + wq[u];
if( dp[u]>v ) dp[u]=v;
}
bool cmp( int a, int b ) {
return dep[a]-lim[a]>dep[b]-lim[b];
}
void vdcp( int s ) {
int c=getc(s);
vis[c] = true;
if( c==s ) {
Point pc = Point(dep[c],dp[c]);
flood(s);
for( int i=; i<=ed; i++ ) {
int u=qu[i];
if( u!=s && dep[u]-lim[u]<=pc.x ) update(pc,u);
}
} else {
vdcp(s);
flood(c);
sort( qu+, qu++ed, cmp );
convex.init();
int cur = c;
for( register int i=; i<=ed; i++ ) {
int u=qu[i];
while( cur!=s && dep[anc[cur]]>=dep[u]-lim[u] ) {
cur=anc[cur];
convex.append( Point(dep[cur],dp[cur]) );
}
if( convex.top>= )
update( convex.query(wp[u]), u );
}
Point cp = Point(dep[c],dp[c]);
for( register int i=ed; i>=; i-- ) {
int u=qu[i];
if( u==c ) continue;
if( dep[u]-lim[u]>dep[c] ) break;
update( cp, u );
}
}
for( int t=head[c]; t; t=next[t] ) {
int v=dest[t];
if( vis[v] ) continue;
vdcp(v);
}
}
int main() {
scanf( "%d%d", &n, &case_type );
for( int i=; i<=n; i++ ) {
int f;
dnt s;
scanf( "%d%lld%lld%lld%lld", &f, &s, wp+i, wq+i, lim+i );
adde( f, i, s );
adde( i, f, s );
}
fill( dp, , n, 0x3f );
dp[] = ;
bfs();
vdcp();
for( int i=; i<=n; i++ )
printf( "%lld\n", dp[i] );
}
bzoj 3672 利用点分治将CDQ分治推广到树型结构上的更多相关文章
- COGS 2479. [HZOI 2016]偏序 [CDQ分治套CDQ分治 四维偏序]
传送门 给定一个有n个元素的序列,元素编号为1~n,每个元素有三个属性a,b,c,求序列中满足i<j且ai<aj且bi<bj且ci<cj的数对(i,j)的个数. 对于100%的 ...
- NOI 2007 货币兑换Cash (bzoj 1492) - 斜率优化 - 动态规划 - CDQ分治
Description 小Y最近在一家金券交易所工作.该金券交易所只发行交易两种金券:A纪念券(以下简称A券)和 B纪念券(以下 简称B券).每个持有金券的顾客都有一个自己的帐户.金券的数目可以是一个 ...
- 【BZOJ】1176: [Balkan2007]Mokia(cdq分治)
http://www.lydsy.com/JudgeOnline/problem.php?id=1176 在写这题的时候思维非常逗啊........2333................... 最后 ...
- BZOJ 2726: [SDOI2012]任务安排( dp + cdq分治 )
考虑每批任务对后面任务都有贡献, dp(i) = min( dp(j) + F(i) * (T(i) - T(j) + S) ) (i < j <= N) F, T均为后缀和. 与j有关 ...
- BZOJ 2716 Violet 3 天使玩偶 CDQ分治
题目大意:初始给定平面上的一个点集.提供两种操作: 1.将一个点增加点集 2.查询距离一个点最小的曼哈顿距离 K-D树是啥...不会写... 我仅仅会CDQ分治 对于一个询问,查询的点与这个点的位置关 ...
- BZOJ 2716: [Violet 3]天使玩偶( CDQ分治 + 树状数组 )
先cdq分治, 然后要处理点对答案的贡献, 可以以询问点为中心分成4个区域, 然后去掉绝对值(4种情况讨论), 用BIT维护就行了. --------------------------------- ...
- 【BZOJ 3262】 3262: 陌上花开 (CDQ分治)
3262: 陌上花开 Description 有n朵花,每朵花有三个属性:花形(s).颜色(c).气味(m),又三个整数表示.现要对每朵花评级,一朵花的级别是它拥有的美丽能超过的花的数量.定义一朵花A ...
- BZOJ 2716 [Violet 3]天使玩偶 (CDQ分治、树状数组)
题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=2716 怎么KD树跑得都那么快啊..我写的CDQ分治被暴虐 做四遍CDQ分治,每次求一个 ...
- BZOJ 1935 Tree 园丁的烦恼 CDQ分治/主席树
CDQ分治版本 我们把询问拆成四个前缀和,也就是二维前缀和的表达式, 我们把所有操作放入一个序列中 操作1代表在x,y出现一个树 操作2代表加上在x,y内部树的个数 操作3代表减去在x,y内部树的个数 ...
随机推荐
- springCloud全实战超详细代码demo+笔记
码云: https://gitee.com/houzheng1216/springcloud
- v8-su-root
1.下载userdebug版本 2.设置模块打开develop options 3.勾选usb debugging 4.adb remount 5.解压SuperSU_N.7z(联系我索取)并push ...
- DirectFB简介以及移植[一]【转】
转自:https://blog.csdn.net/wavemcu/article/details/39251805 ****************************************** ...
- MySQL异步复制、半同步复制详解
MySQL数据复制的原理图大致如下: 从上图我们可以看出MySQL数据库的复制需要启动三个线程来实现: 其中1个在主服务器上,另两个在从服务器上.当发出START SLAVE时,从服务器创建一个I/O ...
- 个性化你的Git Log的输出格式
git已经变成了很多程序员日常工具之一. git log是查看git历史的好工具,不过默认的格式并不是特别的直观. 很多时候想要更简便的输出更多或者更少的信息,这里列出几个git log的format ...
- 『实践』百度地图给多个marker添加右键菜单(删除、更新)
js: $.getJSON("./GetStationPlaceServlet",function(json){ for(var i=0;i<json.length;i++) ...
- Flask:操作SQLite3(0.1)
Windows 10家庭中文版,Python 3.6.4,Flask 1.0.2 本文介绍了第一次在Flask框架中操作SQLite3数据库的测试,参考了官网的文档Using SQLite 3 wit ...
- sizeof求结构体大小
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 3 ...
- es6遍历数组forof
- java算法小例子
作为一个程序员,有时候我觉得自己都不适合,因为自己数学不好,算法不好,脑子不够灵活.而对于算法,感觉就像是数学题,总觉得很难.以前上学,在班里总有几个什么都不好,但唯独数学很厉害,真气人!面对难题时, ...