description

题面

solution

线段树分治+斜率优化毒瘤题

题目可以简化为:

你要维护一个包含元素\((x,c)\)的集合

修改操作为从以前的一个版本更新,修改内容为添加或删除一个元素

查询操作给出\(x_0\),查询某个版本中的\(min\{(x-x_0)^2+c\}\)

可以知道版本之间的时间关系形成一颗树

如果在一个版本删除了某个元素,那么在这个版本的子树中都不会再有这个版本

由于子树的\(dfn\)是连续的,因此操作可以简化为在序列上进行,总共有\(O(m)\)个区间

最简单的方法是对于线段树上的每个节点开\(vector\)<\(query\)>和\(vector\)<\(modify\)>,

注意区间查询和单点修改的插入方式有所不同

斜率优化的式子就是\(max_i\{x_i^2+c_i-2x_0x_i\}+x_0^2\)

然而这题比较丧,它会给出两个横坐标相等的点

因此一定要判断有向斜率

我们先将修改和询问按照\(x\)和\(k\)外部排序,这样可以保证线段树分治时\(x\)和\(k\)都单调

于是在每个线段树节点里可以单调队列直接做

最后这题卡空间

于是我们只插入修改,对于询问直接递归求解即可

总复杂度为\(O((n+m)logn)\)

code

#include<bits/stdc++.h>
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<iomanip>
#include<cstring>
#include<complex>
#include<vector>
#include<cstdio>
#include<string>
#include<bitset>
#include<ctime>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<set>
#define FILE "a"
#define mp make_pair
#define pb push_back
#define RG register
#define il inline
using namespace std;
typedef unsigned long long ull;
typedef vector<int>VI;
typedef long long ll;
typedef double dd;
const dd eps=1e-10;
const int mod=1e9+7;
const int N=5e5+10;
const int T=5e6+10;
const dd pi=acos(-1);
const int inf=2147483647;
const ll INF=1e18+1;
il ll read(){
RG ll data=0,w=1;RG char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
if(ch=='-')w=-1,ch=getchar();
while(ch<='9'&&ch>='0')data=data*10+ch-48,ch=getchar();
return data*w;
}
il void file(){
srand(time(NULL)+rand());
freopen(FILE".in","r",stdin);
freopen(FILE".out","w",stdout);
} int n,m,mt,b[N],s[N];ll x[N],y[N],ans[N];
struct modify{int l,r;ll x,y;}M[N];
bool cmpm(modify a,modify b){return a.x<b.x;}
struct query{int id,s;ll x;}Q[N];
bool cmpq(query a,query b){return a.x<b.x;}
int head[N],nxt[N<<1],to[N<<1],cnt;
il void add(int u,int v){
to[++cnt]=v;
nxt[cnt]=head[u];
head[u]=cnt;
} VI ML[N],MR[N];int w[N],cntw,sz[N];
void dfs(int u,int fa){
w[u]=++cntw;sz[u]=1;
if(!b[u])ML[s[u]].pb(w[u]);
else MR[s[u]].pb(w[u]-1);
for(RG int i=head[u];i;i=nxt[i]){
RG int v=to[i];if(v==fa)continue;
dfs(v,u);sz[u]+=sz[v];
}
if(!b[u])MR[s[u]].pb(w[u]+sz[u]-1);
else ML[s[u]].pb(w[u]+sz[u]);
} VI f[N<<2];int top[N<<2];
#define mid ((l+r)>>1)
#define ls (i<<1)
#define rs (i<<1|1)
il double getk(ll ax,ll ay,ll bx,ll by){
if(ax==bx)return ay>by?-inf:inf;
return (ay-by)*1.0/(ax-bx);
} void insertmodify(int i,int l,int r,int xl,int yr,int id){
if(xl<=l&&r<=yr){
RG int R=f[i].size()-1;RG ll qx=x[id],qy=y[id];
while(R>0&&getk(x[f[i][R-1]],y[f[i][R-1]],x[f[i][R]],y[f[i][R]])>getk(x[f[i][R]],y[f[i][R]],qx,qy))R--,f[i].pop_back();
R++;f[i].pb(id);return;
}
if(xl<=mid)insertmodify(ls,l,mid,xl,yr,id);
if(mid<yr)insertmodify(rs,mid+1,r,xl,yr,id);
} ll ask(int i,int l,int r,int p,query q){
RG int &L=top[i],R=f[i].size()-1;RG ll ans=INF;
while(L<R&&getk(x[f[i][L]],y[f[i][L]],x[f[i][L+1]],y[f[i][L+1]])<q.x)
L++;
if(L<=R)ans=y[f[i][L]]-q.x*x[f[i][L]]+q.x*q.x;
if(l==r)return ans;
if(p<=mid)ans=min(ans,ask(ls,l,mid,p,q));
else ans=min(ans,ask(rs,mid+1,r,p,q));
return ans;
} int main()
{
n=read();m=read();s[1]=1;y[1]=read();
for(RG int i=2,o,f,id;i<=n;i++){
o=read();f=read()+1;id=read()+1;add(f,i);add(i,f);b[i]=o;s[i]=id;
if(!o){x[id]=read();read();read();y[id]=read();}
}
dfs(1,0); for(RG int i=1;i<=n;i++)
for(RG int j=0,sz=ML[i].size();j<sz;j++)
M[++mt]=(modify){ML[i][j],MR[i][j],2ll*x[i],1ll*x[i]*x[i]+y[i]};
sort(M+1,M+mt+1,cmpm);
for(RG int i=1;i<=mt;i++)x[i]=M[i].x,y[i]=M[i].y;
for(RG int i=1;i<=m;i++)
Q[i].id=i,Q[i].s=read()+1,Q[i].x=read(),ans[i]=INF;
sort(Q+1,Q+m+1,cmpq); for(RG int i=1;i<=mt;i++)insertmodify(1,1,n,M[i].l,M[i].r,i);
for(RG int i=1;i<=m;i++)ans[Q[i].id]=ask(1,1,n,w[Q[i].s],Q[i]);
for(RG int i=1;i<=m;i++)printf("%lld\n",ans[i]);
return 0;
}

[CTSC2016]时空旅行的更多相关文章

  1. Luogu P5416 [CTSC2016]时空旅行

    第一次写线段树分治的题目,没想到是道这么毒的题233 首先发现题目里的\((x,y,z,c)\)就是在放屁,只有\((x,c)\)是有用的 因此我们可以把题意转化为,在某一个时间节点上,求出所有元素的 ...

  2. [UOJ198][CTSC2016]时空旅行

    uoj description 你要维护若干个集合,每个集合都是有一个编号比他小的集合扩展而来,扩展内容为加入一个新的元素\((x,c)\)或者删除一个已有元素.集合的扩展关系之间构成一个树形结构. ...

  3. [CTSC2016]时空旅行(线段树+凸包)

    应该是比较套路的,但是要A掉仍然不容易. 下面理一下思路,思路清楚了也就不难写出来了. 0.显然y,z坐标是搞笑的,忽略即可. 1.如果x不变,那么直接set即可解决. 2.考虑一个空间和询问x0,通 ...

  4. CTSC2016时空旅行

    当时看这道题AC的人数比较多,就开了这道题. 很容易发现是这是一个有关凸包的题. 然后不知道怎么维护凸包,一直在想cdq,感觉复杂度不行,于是被这玩意难住了…… 幸好有亲学长yyh造福人类的题解:ht ...

  5. @loj - 2987@ 「CTSC2016」时空旅行

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 2045 年,人类的技术突飞猛进,已经找到了进行时空旅行的方法. ...

  6. 【CTSC2016】时空旅行

    链接 http://uoj.ac/problem/198 题解 首先要发现答案要我们求这个式子: \[ ans=min\bigl((x_i-x)^2+c_i\bigr) \] 显而易见的是这种时空嫁接 ...

  7. uoj198【CTSC2016】时空旅行

    传送门:http://uoj.ac/problem/198 [题解] 首先y.z是没有用的.. 然后式子就是w = (x0-xi)^2+ci的最小值,化出来可以变成一个直线的形式. 然后我们可以用线段 ...

  8. 【UOJ #198】【CTSC 2016】时空旅行

    http://uoj.ac/problem/198 (先补一下以前的题解) 这道题5分暴力好写好调,链上部分分可以用可持久化线段树,每次旅行\(x\)值相同的可以用标记永久化线段树.我还听到某些神犇说 ...

  9. [UOJ198]时空旅行

    看懂题目就知道$y,z$是没用的,这题相当于是给一堆$(x_i,c_i)$和询问$x_q$,要求$(x_q-x_i)^2+c_i$的最大值 先把这个式子拆开:$-2x_ix_q+x_i^2+c_i+x ...

随机推荐

  1. 获取Java线程返回值的几种方式

    在实际开发过程中,我们有时候会遇到主线程调用子线程,要等待子线程返回的结果来进行下一步动作的业务. 那么怎么获取子线程返回的值呢,我这里总结了三种方式: 主线程等待. Join方法等待. 实现Call ...

  2. 《More Effective C++ 》读书笔记(二)Exception 异常

    这事篇读书笔记,只记录自己的理解和总结,一般情况不对其举例子具体说明,因为那正是书本身做的事情,我的笔记作为梳理和复习之用,划重点.我推荐学C++的人都好好读一遍Effective C++ 系列,真是 ...

  3. 高可用Kubernetes集群-5. 部署flannel网络

    七.部署flannel网络 kubernetes支持基于vxlan方式的flannel与weave网络,基于BGP路由的Calico网络,本节采用flannel网络. Flannel网络采用etcd等 ...

  4. 华为云分布式缓存服务DCS与开源服务差异对比

    华为云分布式缓存DCS提供单机.主备.集群等丰富的实例类型,满足用户高读写性能及快速数据访问的业务诉求.支持丰富的实例管理操作,帮助用户省去运维烦恼.用户可以聚焦于业务逻辑本身,而无需过多考虑部署.监 ...

  5. JAVA学习笔记--数组初始化

    JAVA中,数组只是相同类型的.用一个标识符名称封装到一起的一个对象序列或基本类型数据序列.数组通过方括号下标操作符[]来定义和使用,要定义一个数组只需在类型名后面加上一个方括号即可,如: int[] ...

  6. Just a Hook:线段树+区间修改

    E - Just a Hook In the game of DotA, Pudge’s meat hook is actually the most horrible thing for most ...

  7. day-19 多种优化模型下的简单神经网络tensorflow示例

    如下样例基于tensorflow实现了一个简单的3层深度学习入门框架程序,程序主要有如下特性: 1.  基于著名的MNIST手写数字集样例数据:http://yann.lecun.com/exdb/m ...

  8. UI优秀框架(库)

    1.vux 官网:https://doc.vux.li/zh-CN/ Github:https://github.com/airyland/vux 13818  Stars  3064 Forks   ...

  9. CSS3在线实战

    作者声明:本博客中所写的文章,都是博主自学过程的笔记,参考了很多的学习资料,学习资料和笔记会注明出处,所有的内容都以交流学习为主.有不正确的地方,欢迎批评指正. 本节课视频网站:https://www ...

  10. js操作对象属性值为字符串

    今天在项目开发中遇到一个没遇到过的问题,这个问题是需要对比两个对象a和b,a是一个只有一个属性的对象,b是一个含有多个属性对象,如果b中包含和a一模一样的属性名和值,则把这个一样的属性和值从b中删除了 ...