P5025 [SNOI2017]炸弹 题解

≯℡__Kan透↙ 提交于 2019-12-04 03:38:05

蒟蒻的第一篇黑题题解(学了这么长时间了才第一道也是没谁了。)

题目链接

Solution

朴素:

 根据题目描述,我们可以处理每一个x节点左右爆炸半径范围内的点,然后模拟一次爆炸 (for),遍历每一个点。每当我们遍历到一个点,我们就对这个点在进行一次处理半径+dfs直到没有能遍历的了,直接统计遍历的点数然后输出。

时间复杂度:O(n^2)级别。

于是乎:

  

 

 是不是可以暴力踩标算了呢?

显然完全接受不了。

因为这个算法要考虑到连边操作。对于每一个点,我们要将它自己连向能够遍历到的边,那么意味着,一个边要被连很多次。

就如下图:

 

 其中编号为1的节点能炸到1 2 3 4 5,2能炸到1 2 3,3能炸到1 2 3 4 5,其中2 3 4挺惨的。

 于是有了这么一个线段树建图的思路:

考虑线段树的性质,上层节点可以代表下层节点的一段区间,那么我们是不是也可以建这类的上层边,表示通过这条边能遍历到它对应的下层边的所有点(就是这句精华)

这样我们先把所有点进行从1到n的编号,然后跑线段树的build建边,然后依次对每一个点进行从底层(只能这么理解)到上层节点的连边表示能遍历到上层节点能遍历到的节点。

这样的话,会形成很多环。在对每一个点都进行范围连边以后,对于一个环(或者强连通分量)内的点,自然能互相遍历到,为了方便统计,我们需要跑一个tarjan记录强连通分量内的点能遍历到的左右区间范围。

跑完tarjan就可以O(n)遍历一遍按照题目所说进行统计了呢。

部分代码

build:

void build(int now,int l,int r)
{
    canl[now]=l;canr[now]=r;//can表示当前节点能遍历到的左右区间范围 
    if(l==r)//线段树日常操作 
    {
        id[r]=now;
        return;
    }
    int mid=(l+r)>>1;
    build(now<<1,l,mid);build(now<<1|1,mid+1,r);
    add(now,now<<1);add(now,now<<1|1);//不同于线段树的地方,向左右两个节点连边,往下层伸展 
}

下层往上连边部分:

void lb(int mb,int l,int r,int L,int R,int now)
{
    if(l>=L&&r<=R)
    {
        if(mb==now)//mb表示目标节点,也就是当前初始节点 
            return;
        add(mb,now);//now表示线段树上当前节点,从mb到now连边表示mb这个节点可以爆炸到now节点代表的区间 
        return;
    }
    int mid=(l+r)>>1;
    if(L<=mid)lb(mb,l,mid,L,R,now<<1);//左右区间连边 
    if(R>mid)lb(mb,mid+1,r,L,R,now<<1|1);
}

tarjan部分:

void tarjan(int now)
{
    dfn[now]=++tarjannum;//dfn都知道是什么吧 
    low[now]=tarjannum;
    vis[now]=1;//设置入栈 
    st[++top]=now;
    for(int i=hea[now];i;i=edge[i].nex)
    {
        int v=edge[i].to;
        if(!dfn[v])
        {
            tarjan(v);low[now]=min(low[now],low[v]);
        }
        else
            if(vis[v]) //在栈中 
                low[now]=min(low[now],dfn[v]);//以上正常tarjan操作 
    }
    if(low[now]==dfn[now])
    {
        numdag++;//标记为同一个dag 
        do
        {
            int topp=st[top--];
            dagbh[topp]=numdag;
               leff[numdag]=min(leff[numdag],canl[topp]);
//在统计tarjan的同时,把这个强联通分量里面所有的点能到的左右区间范围统计一下 
            rigg[numdag]=max(rigg[numdag],canr[topp]);
            vis[topp]=0;//出栈标记 
        }while(st[top+1]!=now);//dowhile保证能够便历到,因为上面top已经--了所以要判断top+1 
    }
}

高清无码完整代码:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#define N 500003
#define modd 1000000007 
using namespace std;
long long read()
{
    long long ans=0;
    char ch=getchar(),last=' ';
    while(ch<'0'||ch>'9')last=ch,ch=getchar();
    while(ch>='0'&&ch<='9')ans=(ans<<3)+(ans<<1)+ch-'0',ch=getchar();
    return last=='-'?-ans:ans;
}
int n,id[N],hea[N*4],num,vis[N*4],dagbh[N*4],lef,rig,top,canl[N*4],canr[N*4];
int st[N*4],dfn[N*4],low[N*4],leff[N*4],rigg[N*4],tarjannum,numdag;
long long r[N],x[N],ans;
struct edg{
    int nex,to;
}edge[N*12];
inline void add(int from,int to)
{
    num++;
    edge[num]={hea[from],to};
    hea[from]=num;
}
void build(int now,int l,int r)
{
    canl[now]=l;canr[now]=r;//can表示当前节点能遍历到的左右区间范围 
    if(l==r)//线段树日常操作 
    {
        id[r]=now;
        return;
    }
    int mid=(l+r)>>1;
    build(now<<1,l,mid);build(now<<1|1,mid+1,r);
    add(now,now<<1);add(now,now<<1|1);//不同于线段树的地方,向左右两个节点连边,往下层伸展 

void lb(int mb,int l,int r,int L,int R,int now)
{
    if(l>=L&&r<=R)
    {
        if(mb==now)//mb表示目标节点,也就是当前初始节点 
            return;
        add(mb,now);//now表示线段树上当前节点 
        return;
    }
    int mid=(l+r)>>1;
    if(L<=mid)lb(mb,l,mid,L,R,now<<1);//左右区间连边 
    if(R>mid)lb(mb,mid+1,r,L,R,now<<1|1);
}
void tarjan(int now)
{
    dfn[now]=++tarjannum;//dfn都知道是什么吧 
    low[now]=tarjannum;
    vis[now]=1;//设置入栈 
    st[++top]=now;
    for(int i=hea[now];i;i=edge[i].nex)
    {
        int v=edge[i].to;
        if(!dfn[v])
        {
            tarjan(v);low[now]=min(low[now],low[v]);
        }
        else
            if(vis[v]) //在栈中 
                low[now]=min(low[now],dfn[v]);//以上正常tarjan操作 
    }
    if(low[now]==dfn[now])
    {
        numdag++;//标记为同一个dag 
        do
        {
            int topp=st[top--];
            dagbh[topp]=numdag;
               leff[numdag]=min(leff[numdag],canl[topp]);
//在统计tarjan的同时,把这个强联通分量里面所有的点能到的左右区间范围统计一下 
            rigg[numdag]=max(rigg[numdag],canr[topp]);
            vis[topp]=0;//出栈标记 
        }while(st[top+1]!=now);//dowhile保证能够便历到,因为上面top已经--了所以要判断top+1 
    }
}
int main()
{
    memset(leff,0x3f,sizeof(leff));
    n=read();
    for(int i=1;i<=n;i++)
    {
        x[i]=read(),r[i]=read();//x为横轴坐标,r为爆炸半径 
    }
    x[n+1]=0x3f3f3f3f3f3f3f3fll;
    build(1,1,n);
    for(int i=1;i<=n;i++)
    {
        if(!r[i]) continue;
        rig=upper_bound(x+1,x+1+n,x[i]+r[i])-x-1;
        lef=lower_bound(x+1,x+1+n,x[i]-r[i])-x;
        lb(id[i],1,n,lef,rig,1);
    }    
    tarjan(1);//从1号节点往下跑tarjan缩点跑成dag
    for(int i=1;i<=n;i++)
    {
        int fromm=dagbh[id[i]];
        ans=(ans+((long long)i*(rigg[fromm]-leff[fromm]+1))%modd)%modd;//ll
    }
    printf("%lld",ans%modd);
}

完结。。

可能是最水的黑题了其他的都不会做。

 希望讲解过后能对泥萌有所帮助qwq

易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!