[2016北京集训测试赛15]statement-[线段树+拆环]

Description

Solution

由于题目要求,将a[i]->b[i](边权为i)后所得的图应该是由森林和环套树组合而成。

假如是树形结构,所有的t[i]就直接在线段树t[i]点的dfs序(即in[t[i]],out[t[i]]区间)处记录t[i]点的深度。

这样,针对所有的f[i],在线段树上查找所有包含in[f[i]]点的区间所记录的最大深度d。(这个深度就是在离f[i]最近并且已经验证了是真命题的祖先的深度)

然后用倍增算出f[i]向上到深度d,所经过的编号最大值c。ans=min(ans,c)。

原因:ans是指,图中存在a[i]->b[i](1<=i<=ans)时该询问刚好出现矛盾。如果ans-=1,则所有的f[i]都无法到达离它最近的并且已经验证了是真命题的祖先点。

环套树结构:我们把环拆成树

其中Ca->Cb的边权还是a->b边的边权。以此类推其他都是这样。(原本环套树上以a,b,c,d,e为根的树还是以a,b,c,d,e为根)如此,假如说t[i]在环上,我们除了在in[t[i]],out[t[i]]区间记录t[i]点的深度,还要在in[C(t[i])],out[C(t[i])]区间记录C(t[i])点的深度。这样,不论是原本以t[i]为根的子树,还是环上的以其他点为根的树的信息都可以更新完毕。

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int N=2e5+10;
int n,_n,m,Q,a,b;
int fa1[N<<1],val[N<<1],cir[N];
int tag[N];//0-not_visit 1-tag the stack 2-out the stack
int rt[N],cnt;

int y[N<<1],nxt[N<<1],h[N<<1],tot;
void link(int _x,int _y){y[++tot]=_y;nxt[tot]=h[_x];h[_x]=tot;}
void build(int u)
{
    if (!fa1[u]){rt[++cnt]=u;tag[u]=2;return;}
    if (tag[fa1[u]]==1)
    {
        int x=u;
        for(int t=fa1[x];t!=u;t=fa1[t])
        {
            cir[t]=++n;val[n]=val[t];link(n,x);
            fa1[x]=n;x=fa1[x];
        }
        cir[u]=++n;link(n,x);rt[++cnt]=n;
        tag[u]=2;
        return;
    }
    tag[u]=1;
    if (!tag[fa1[u]]) build(fa1[u]);link(fa1[u],u);
    tag[u]=2;
}
int dfn,in[N<<1],out[N<<1],dep[N<<1],fa[N<<1][20],mxe[N<<1][20];
void dfs(int x)
{
    in[x]=++dfn;mxe[x][0]=val[x];dep[x]=dep[fa[x][0]]+1;
    for (int i=1;i<=19;i++)
    {
        fa[x][i]=fa[fa[x][i-1]][i-1];
        mxe[x][i]=max(mxe[x][i-1],mxe[fa[x][i-1]][i-1]);
    }
    for (int i=h[x];i;i=nxt[i]) fa[y[i]][0]=x,dfs(y[i]);
    out[x]=dfn;
}
int cur[N<<3],mxd[N<<3];
void modify(int k,int l,int r,int askx,int asky,int d)
{
    if (askx<=l&&r<=asky){if (Q<cur[k]||mxd[k]<d) cur[k]=Q,mxd[k]=d;return;}
    int mid=(l+r)/2;
    if (askx<=mid) modify(k<<1,l,mid,askx,asky,d);
    if (asky>mid) modify(k<<1|1,mid+1,r,askx,asky,d);
}
int query(int k,int l,int r,int x)
{
    int mid=(l+r)/2,re=cur[k]==Q?mxd[k]:0;
    if (l==r) return re;
    if (x<=mid) re=max(re,query(k<<1,l,mid,x));else re=max(re,query(k<<1|1,mid+1,r,x));
    return re;
}
int t,f,ret,ref;
int main()
{
    scanf("%d%d",&n,&m);_n=n;
    for (int i=1;i<=m;i++) {scanf("%d%d",&a,&b);fa1[b]=a;val[b]=i;}
    for (int i=1;i<=_n;i++) if (!tag[i]) build(i);
    for (int i=1;i<=cnt;i++) dfs(rt[i]);
    scanf("%d",&Q);
    while (Q--)
    {
        scanf("%d",&ret);
        for (int i=1;i<=ret;i++)
        {
            scanf("%d",&t),modify(1,1,n,in[t],out[t],dep[t]);
            if (cir[t]) modify(1,1,n,in[cir[t]],out[cir[t]],dep[cir[t]]);
        }
        int ans=m+1,d,c;
        scanf("%d",&ref);
        for (int i=1;i<=ref;i++)
        {
            scanf("%d",&f);
            d=query(1,1,n,in[f]);c=0;
            if (d==0) continue;else d=dep[f]-d;
            for(int j=0;d;d>>=1,j++)
            if (d&1) c=max(c,mxe[f][j]),f=fa[f][j];
            ans=min(ans,c);
        }
        if (ans==m+1) printf("OK\n");else printf("%d\n",ans);
    }
}

原文地址:https://www.cnblogs.com/coco-night/p/9715571.html

时间: 2024-11-07 17:35:31

[2016北京集训测试赛15]statement-[线段树+拆环]的相关文章

【2016北京集训测试赛(八)】 crash的数列

Description 题解 题目说这是一个具有神奇特性的数列!这句话是非常有用的因为我们发现,如果套着这个数列的定义再从原数列引出一个新数列,它居然还是一样的...... 于是我们就想到了能不能用多点数列套着来加速转移呢? 但是发现好像太多数列套起来是可以烦死人的...... 我们就采用嵌套两次吧,记原数列为A,第一层嵌套为B,第二层嵌套为C. 我们其实可以发现一些规律,对于Ci,它对应了B中i的个数:对于Bi,它对应了A中i的个数. 稍加处理即可,我们一边计算一边模拟数列的运算,同时可以计算

【2016北京集训测试赛(八)】直径

注意:时限更改为4s 题解 考虑最原始的直径求法:找到离根节点(或离其他任意一点)最远的节点far1,再从far1出发找到离far1最远的节点far2,far1至far2的距离即为直径. 题目中提到了将原树的子树复制成新子树这一操作,显然如果我们将子树真正复制出来是会爆炸的. 实际上我们可以将每棵新子树中,真正有用的节点提取出来,以简化每个子树的结构,再按照题目的要求连接各个新子树. 我们用虚树来重构每一棵子树.每棵子树的虚树的关键点应至少包含: 子树的根节点. 这棵子树内部的直径的两端节点.

【2016北京集训测试赛(七)】自动机 (思考题)

Time Limit: 1000 ms Memory Limit: 256 MB Description Solution 这是一道看起来令人毫无头绪的题,然而确实十分简单巧妙TAT. 题目要求所有点执行相同指令后都回到初始状态. 我们先来考虑只有两种状态的情况:初始状态$T_0$与另一个状态$T_x$. 这样,我们可以通过一个二元记忆化深搜,来得到一种方案A,使得$T_0$回到$T_0$,且$T_x$回到$T_0$.如果这个方案都不存在,那么此时无解. 现在我们知道,执行方案A后,$T_x$与

【2016北京集训测试赛】river

HINT 注意是全程不能经过两个相同的景点,并且一天的开始和结束不能用同样的交通方式. [吐槽] 嗯..看到这题的想法的话..先想到了每个点的度为2,然后就有点不知所措了 隐隐约约想到了网络流,但并没有继续往下想了... 听完学长的讲评之后(%xj)个人觉得建图还是很有意思的ovo [题解] 因为每个点到对面都有k种方式,那就想到每个点原来的点$x_0$拆成k个点$x_1$, $x_2$, $x_3$... $x_k$ 然后很自然地$x_0$和拆成的点之间要连边 容量的话,因为hint里面的限制

【2016北京集训测试赛】azelso(unfinished)

[吐槽] 首先当然是要orzyww啦 以及orzyxq奇妙顺推很强qwq 嗯..怎么说呢虽然说之前零零散散做了一些概d的题目但是总感觉好像并没有弄得比较明白啊..(我的妈果然蒟蒻) 这题的话可以说是难得的一道搞得比较清楚的概d题目吧记录一下还是挺有意思的ovo 当然咯..显然考场上并没有推出来..嗯qwq [题解] 看到说要求期望的距离,然后总的长度又被分成了一段一段的(各个事件) 所以就有一个比较直接的想法:将每一段期望走的次数算出来然后再乘上每一段的距离,加起来就是答案啦 那么现在问题来了怎

2016北京集训测试赛(十七)- 小结

先说一下总体的情况. 场上期望得分 50 + 40 + 30 = 120 , 最后得分 50 + 0 + 30 = 80 , 实际上自己能力所及能做到的 50 + 65 + 30 = 145 分. 第二题爆零是因为我开始写了一个做法, 后来发现这个做法是错的, 然后开始随机化, 调着调着突然发现只有一分钟了, 然后自己把之前调的改回来, 然后发现怎么全都输出 0 ??? Excuse me ?? 原本不用随机化的做法可以拿 40 分, 如果结合上暴力就有 65 了. 这几天打起比赛来还是暴露了许

[2016北京集训测试赛(一)]奇怪的树-[树链剖分]

Description Solution 对于操作1,不论选了哪个点为a,最后反转颜色的点集都只有两种类型(显然啦). 暴力解法:对每个操作3,从a向上直到根节点,每到一个节点记录(它父亲的黑点数减去自己的黑点数)*父亲节点的编号.另外,还要记录a子树内的黑点.. 这种O(n2)的做法肯定会爆,考虑优化.由于这是一棵静态树,考虑树链剖分. 需要记录一个数组re[x][0/1][0/1].第2维表示深度的奇偶,第3维表示点的颜色.例如:re[x][0][0]记录的是初始情况下以x为根的子树中深度为

[2016北京集训测试赛(五)]打地鼠-[思考题]

Description Solution 我们先只考虑一只地鼠的情况,依题意得,在某一个时刻该地鼠的可能停留位置是一个公差为2的等差数列.我们设这个等差数列的两端为[L,R].则如果区间[L+1,R-1]的格子被打实际上是不会影响L和R的(列一个等差数列实际模拟一下就发现啦).而如果格子L被打,则L+2:如果格子R被打,则R-2.打了格子后,别忘了L--,R++. 嗯根据以上性质,我们可以知道,地鼠1,3,5,7,9...的L是非递减的,地鼠2,4,6,8,10...的L也是非递减的. 然后看一

[2016北京集训测试赛17]crash的游戏-[组合数+斯特林数+拉格朗日插值]

Description Solution 核心思想是把组合数当成一个奇怪的多项式,然后拉格朗日插值..:哦对了,还要用到第二类斯特林数(就是把若干个球放到若干个盒子)的一个公式: $x^{n}=\sum _{i=0}^{n}C(n,i)*i!*S(i,x)$ 围观大佬博客(qaq公式太难打了) Code #include<iostream> #include<cstdio> #include<cstring> #include<cmath> using na