「题解」:e

问题 B: e

时间限制: 2 Sec  内存限制: 512 MB

题面



题面谢绝公开。

题解



话说一天考两个主席树这回事……

正解可以叫树上主席树??(脸哥说也叫主席树上树???)

对于树上的每一条链建主席树,支持链上查询前驱和后继。

对于所有的$p[i]$,他说怎么得到就按他说的做就好,然后求所有$p[i]$的$LCA$。

对于每个$p[i]$到$LCA$的链上查一次$r$的前驱和后继更新答案即可。

注意:参数不要传反、别一个特判把自己判掉、pre和nxt的代码不要粘贴,粘贴了不要忘记改掉内部递归函数……

(昨天下午大概帮3、4个人调这题代码???RP++)

#include<bits/stdc++.h>
#define rint register int
using namespace std;
const int N=100005,M=1000000000;
inline void read(int &A)
{
    A=0;int B=1;char ch=getchar();
    while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)B=-1;ch=getchar();}
    while(ch>=‘0‘&&ch<=‘9‘){A=(A<<3)+(A<<1)+ch-‘0‘;ch=getchar();}
    A=A*B;
}
int n,q,tp,ans,a[N],p[N];
int v[N<<1],nxt[N<<1],first[N<<1],tot;
int root[N],ls[N<<6],rs[N<<6],val[N<<6],cnt;
int dep[N],f[N][23];
inline void build_line(int uu,int vv)
{
    v[++tot]=vv,nxt[tot]=first[uu];
    first[uu]=tot;return ;
}
inline void insert(int &k,int l,int r,int his,int dat)
{
    k=++cnt;val[k]=val[his]+1;
    if(l==r)return ;
    ls[k]=ls[his],rs[k]=rs[his];
    int mid=(l+r)>>1;
    if(dat<=mid)insert(ls[k],l,mid,ls[his],dat);
    else insert(rs[k],mid+1,r,rs[his],dat);
}
inline int get_pre(int k1,int k2,int l,int r,int dat)
{
    if(!(val[k2]-val[k1]))return 0;
    if(l==r)return l;int mid=(l+r)>>1;
    if(dat<=mid)return get_pre(ls[k1],ls[k2],l,mid,dat);
    int lin=get_pre(rs[k1],rs[k2],mid+1,r,dat);
    if(!lin)return get_pre(ls[k1],ls[k2],l,mid,dat);
    return lin;
}
inline int get_nxt(int k1,int k2,int l,int r,int dat)
{
    if(!(val[k2]-val[k1]))return 0;
    if(l==r)return l;int mid=(l+r)>>1;
    if(dat>mid)return get_nxt(rs[k1],rs[k2],mid+1,r,dat);
    int lin=get_nxt(ls[k1],ls[k2],l,mid,dat);
    if(!lin)return get_nxt(rs[k1],rs[k2],mid+1,r,dat);
    return lin;
}
inline void bfs()
{
    queue <int> q;
    q.push(1);dep[1]=1;
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(rint i=first[x];i;i=nxt[i])
        {
            int y=v[i];if(dep[y])continue;
            dep[y]=dep[x]+1;f[y][0]=x;
            for(rint i=1;i<=20;++i)
                f[y][i]=f[f[y][i-1]][i-1];
            q.push(y);
        }
    }
    return ;
}
inline int get_lca(int x,int y)
{
    if(dep[x]<dep[y])swap(x,y);
    for(rint i=20;i>=0;--i)
        if(dep[f[x][i]]>=dep[y])x=f[x][i];
    if(x==y)return x;
    for(rint i=20;i>=0;--i)
        if(f[x][i]!=f[y][i])
            x=f[x][i],y=f[y][i];
    return f[x][0];
}
inline void dfs_build(int x,int fa)
{
    insert(root[x],1,M,root[fa],a[x]);
    for(rint i=first[x];i;i=nxt[i])
        if(v[i]!=fa)dfs_build(v[i],x);
    return ;
}
int main()
{
    read(n),read(q),read(tp);
    for(rint i=1;i<=n;++i)read(a[i]);
    for(rint i=1,ST,EN;i<n;++i)
        read(ST),read(EN),build_line(ST,EN),build_line(EN,ST);
    dfs_build(1,0);bfs();
    for(rint i=1,r,k;i<=q;++i)
    {
        read(r),read(k);
        for(rint j=1,x;j<=k;++j)
            read(x),p[j]=(x-1+ans*tp)%n+1;
        int lca=p[1];
        for(rint j=2;j<=k;++j)
            lca=get_lca(lca,p[j]);
        ans=0x7fffffff;
        for(rint j=1;j<=k;++j)
        {
            int pre=get_pre(root[f[lca][0]],root[p[j]],1,M,r);
            int nxt=get_nxt(root[f[lca][0]],root[p[j]],1,M,r);
            if(pre)ans=min(r-pre,ans);
            if(nxt)ans=min(nxt-r,ans);
        }
        printf("%d\n",ans);
    }
}

原文地址:https://www.cnblogs.com/xingmi-weiyouni/p/11655800.html

时间: 2024-11-02 23:33:45

「题解」:e的相关文章

「题解」kuangbin 最小生成树

POJ-1251 Jungle Roads (水题,%c) POJ-1287 Networking (水) POJ-2031 Building a Space Station (%f浮点数尴尬精度,两球间距离) POJ-2421 Constructing Roads (一些边已建好,简单处理一下) ZOJ-1586 QS Network (处理一下边权) HDU-1233 还是畅通工程 (水) HDU-1875 畅通工程再续 (浮点数,条件连边) HDU-1301 Jungle Roads (重

「题解」「美团 CodeM 资格赛」跳格子

目录 「题解」「美团 CodeM 资格赛」跳格子 题目描述 考场思路 思路分析及正解代码 「题解」「美团 CodeM 资格赛」跳格子 今天真的考自闭了... \(T1\) 花了 \(2h\) 都没有搞定,最后无奈 \(90pts\) . 然而 \(T2\) 想到很多很奇怪的做法,结果正解在 \(28min\) 之内做出... 结果 \(T3\) 是本人最不擅长的伪期望,直接跳过,啥都没得. 来水一发 \(T1\) 的题解... 题目描述 点这里 考场思路 其实并没有什么十分特别的思路,就是一通乱

「题解」:[loj2763][JOI2013]现代豪宅

问题 A: 现代豪宅 时间限制: 1 Sec  内存限制: 256 MB 题面 题目描述 (题目译自 $JOI 2013 Final T3$「現代的な屋敷」) 你在某个很大的豪宅里迷路了.这个豪宅由东西方向$M$列,南北方向$N$行的正方形房间组成. 从西面开始第$x$列,从南面开始第y行的房间用$(x,y)$表示. 相邻的两个房间之间都有一扇门.对于每扇门,门关上表示不可通行,门打开表示可以通行. 当门打开时,从门一边的房间走到另一边的房间需要$1$分钟. 另外,一些房间中有一个开关,如果连续

「题解」:世界线

问题 A: 世界线 时间限制: 1 Sec  内存限制: 256 MB 题面 题面谢绝公开. 题解 我刚学完bitset就考了???然而我赛时并没有想到是bitset…… 话说我bitset还是没有颓完啊赶紧去补一波时空复杂度…… bitset优化暴力(据某诺神说这是bitset果题??) 注意dfs的时侯如果目标节点已经搜过了直接用当前节点的bitset或上目标节点的bitset值即可. 还有$6e4×6e4$会炸空间.开一个$6e4×3e2$的然后一半一半跑即可. 至于怎么跑一半,直接判断和

「题解」:联

问题 A: 联 时间限制: 2 Sec  内存限制: 256 MB 题面 题面谢绝公开. 题解 解法1:离散化+线段树. 1e18的数据范围直接离散化掉所有的l和r,加一个映射数组表示间距即可. 线段树维护区间和,扫0的时候判定子树和等不等于子树大小. 维护两个标记:lazy(将一个区间赋值成什么).xr(是否异或整个区间)跑就完了. (我会说我沉迷解法二而解法一是水的么 感谢Larry提供的优质讲解及代码) #include<bits/stdc++.h> #define int long l

「题解」:X国的军队

问题 A: X国的军队 时间限制: 1 Sec  内存限制: 256 MB 题面 题面谢绝公开. 题解 简单贪心. 按照存活的士兵数量(即参加战斗的士兵数量减去阵亡的士兵数量)排序. 若存活士兵数量相同则按照参与战斗的士兵数量排序. 顺序扫一遍统计答案. #include<bits/stdc++.h> #define int long long #define rint register int #define read(A) A=init() using namespace std; inl

「题解」:Simple

问题 A: Simple 时间限制: 1 Sec  内存限制: 256 MB 题面 题面谢绝公开. 题解 不算数学的数学题?? 直接枚举会重.$60%$两种算法:1.无脑$vis$数组记录.2.$exgcd$解方程判定是否有解. $100%$:首先考虑特殊情况:$n$.$m$互质. 我们设$n*x+m*y=z$,考虑枚举$y$和$x$,不难发现,当$y>=x$的时候均能找到一个$y'$使得$n|(y-y')$. 于是会出现重复.因此只需枚举$y([0,n-1])$,计算贡献即可. 对于一般情况,

「题解」:毛一琛/$cow$ $subsets$

问题 A: 毛一琛/$cow$ $subsets$ 时间限制: 1 Sec  内存限制: 512 MB 题面 题面谢绝公开. 题解 题名貌似是个大神??看起来像是签到题然后就死了. 首先$O(3^n)$算法显然.也显然过不去$20$的测试点. 正解是赫赫有名的$meet$ $in$ $the$ $middle$算法.数据在40以内的都能用$meet$ $in$ $the$ $middle$?? 对于两半路径,可以拼起来并且构成合法答案的条件是两人获得的分数相同. 所以一个比较聪明的办法是,不去记

「题解」:毛三琛

问题 C: 毛三琛subset 时间限制: 1 Sec  内存限制: 512 MB 题面 题面谢绝公开. 题解 一眼题解随机化,吓够呛.一句话题解:二分答案加剪枝. 外层枚举$x$,然后二分答案暴力$check$.如果当前答案对于x的check失败就continue, 因为在当前的x中不可能找到比当前答案更优秀的解.加clock卡常可以A. 貌似不需要看脸.毕竟我这个非洲人都一遍A了.复杂度$O(np+nlognlogp)$. 代码:(ps.$¥$神指出了我代码的缺陷:其实在外面直接赋值可以少一