bzoj5017: [Snoi2017]炸弹

最新推荐文章于 2021-08-29 14:00:28 发布
最新推荐文章于 2021-08-29 14:00:28 发布
阅读量1.2k 收藏 1
点赞数
分类专栏： 线段树（合并） 暴力/水题
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_36797743/article/details/77880242
版权

Description

在一条直线上有 N 个炸弹，每个炸弹的坐标是 Xi，爆炸半径是 Ri，当一个炸弹爆炸时，如果另一个炸弹所在位置 Xj 满足：
Xi−Ri≤Xj≤Xi+Ri,那么，该炸弹也会被引爆。
现在，请你帮忙计算一下，先把第 i 个炸弹引爆，将引爆多少个炸弹呢？
Input

第一行，一个数字 N，表示炸弹个数。
第 2∼N+1行，每行 2 个数字，表示 Xi，Ri，保证 Xi 严格递增。
N≤500000
−10^18≤Xi≤10^18
0≤Ri≤2×10^18
Output

一个数字，表示Sigma(i*炸弹i能引爆的炸弹个数),1<=i<=N mod10^9+7。
Sample Input

4

1 1

5 1

6 5

15 15
Sample Output

32

前言

这题是哈老师给的神题
然后他给了我一个有一点点复杂的做法
于是帅爷爷就又教了我一个方法
然后我在今天做操的时候仔细斟酌了一下，觉得很对，于是就A了

题解

首先，我们可以建图。。
就是每个点对他可以炸到的范围连一条边。。
然后对于每个点dfs一次，他可以从他连的点得到新的l和r
这里要注意的时，虽然说每个点必须要访问，更新答案，这样才能确保他的答案是对的
然后由于可能边数很大，于是要用线段树优化建图，就好了

CODE:

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const LL N=2000010;
int n;
LL X[N],R[N];
const int MOD=1e9+7;
struct qq
{
    int x,y,last;
}s[N*8];int Num,last[N];
int ll[N],rr[N];
void init (int x,int y)
{
    Num++;
    s[Num].x=x;s[Num].y=y;
    s[Num].last=last[x];
    last[x]=Num;
    return ;
}
void Ins ()
{
    Num=0;memset(last,-1,sizeof(last));
    scanf("%d",&n);
    for (int u=1;u<=n;u++)
        scanf("%lld%lld",&X[u],&R[u]);
}
struct qr
{
    int l,r;
    int s1,s2;
}tr[N];int num;
int pos[N];
void bt (int l,int r)
{
    int a=++num;
    tr[a].l=l;tr[a].r=r;
    ll[a]=l;rr[a]=r;
    if (l==r) 
    {
        pos[l]=a;
        return ;
    }
    int mid=(l+r)>>1;
    tr[a].s1=num+1;bt(l,mid);
    tr[a].s2=num+1;bt(mid+1,r);
    init(a,tr[a].s1);init(a,tr[a].s2);
}
void Link (int now,int l,int r,int x)
{
    if (tr[now].l==l&&tr[now].r==r)
    {
        init(x,now);
        return ;
    }
    int mid=(tr[now].l+tr[now].r)>>1;
    int s1=tr[now].s1,s2=tr[now].s2;
    if (r<=mid) Link(s1,l,r,x);
    else if (l>mid) Link(s2,l,r,x);
    else Link(s1,l,mid,x),Link(s2,mid+1,r,x);
}
void Bt ()
{
    num=0;bt(1,n);
    //for (int u=1;u<=n;u++) printf("%d ",pos[u]);
    for (int u=1;u<=n;u++)
    {
        int l=1,r=u;
        int p;
        while (l<=r)
        {
            int mid=(l+r)>>1;
            if (X[mid]+R[u]>=X[u]) {p=mid;r=mid-1;}
            else l=mid+1;
        }
        if (p!=u) Link(1,p,u-1,pos[u]);
    //  printf("%d %d %d\n",p,u-1,pos[u]);
        l=u;r=n;
        while (l<=r)
        {
            int mid=(l+r)>>1;
            if (X[mid]-R[u]<=X[u]) {p=mid;l=mid+1;}
            else r=mid-1;
        }
        if (p!=u) Link(1,u+1,p,pos[u]);
    //  printf("%d %d %d\n",u+1,u-1,pos[u]);
    }
}
bool vis[N];
int mymax (int x,int y){return x>y?x:y;}
int mymin (int x,int y){return x<y?x:y;}
void dfs (int x)
{
    vis[x]=true;
    for (int u=last[x];u!=-1;u=s[u].last)
    {
        int y=s[u].y;
        if (vis[y]==false) dfs(y);
        ll[x]=mymin(ll[x],ll[y]);
        rr[x]=mymax(rr[x],rr[y]);
    }
}
void solve ()
{
    memset(vis,false,sizeof(vis));
    for (int u=1;u<=num;u++) 
        dfs(u);
    long long ans=0;
    for (int u=1;u<=n;u++) 
    {
        //printf("YES:%d %d\n",ll[pos[u]],rr[pos[u]]);
        ans=(ans+(long long)u*(rr[pos[u]]-ll[pos[u]]+1)%MOD)%MOD;
    }
    printf("%lld\n",ans);
}
int main()
{
    Ins();
    Bt();
    solve();
    return 0;
}

Update 9.9 方法被人hack了。。
其实是bzoj的数据太水了。。
上面的方法其实有很大的漏洞，也跑不出来。。
其实有环的时候就会有很大的问题。。
这个自己对拍一下，就应该能找到了
所以我们要用强连通缩点，就可以跑过去了。。
我用的Claris的代码对拍

新的CODE:

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int N=2000010;
int n;
LL X[N],R[N];
const int MOD=1e9+7;
struct qq
{
    int x,y,last;
}s[N*8];int Num,last[N];
int ll[N],rr[N];
void init (int x,int y)
{
    Num++;
    s[Num].x=x;s[Num].y=y;
    s[Num].last=last[x];
    last[x]=Num;
    return ;
}
void Ins ()
{
    Num=0;memset(last,-1,sizeof(last));
    scanf("%d",&n);
    for (int u=1;u<=n;u++)
        scanf("%lld%lld",&X[u],&R[u]);
}
struct qr
{
    int l,r;
    int s1,s2;
}tr[N];int num;
int pos[N];
void bt (int l,int r)
{
    int a=++num;
    tr[a].l=l;tr[a].r=r;
    ll[a]=l;rr[a]=r;
    if (l==r) 
    {
        pos[l]=a;
        return ;
    }
    int mid=(l+r)>>1;
    tr[a].s1=num+1;bt(l,mid);
    tr[a].s2=num+1;bt(mid+1,r);
    init(a,tr[a].s1);init(a,tr[a].s2);
}
void Link (int now,int l,int r,int x)
{
    if (tr[now].l==l&&tr[now].r==r)
    {
        init(x,now);
        return ;
    }
    int mid=(tr[now].l+tr[now].r)>>1;
    int s1=tr[now].s1,s2=tr[now].s2;
    if (r<=mid) Link(s1,l,r,x);
    else if (l>mid) Link(s2,l,r,x);
    else Link(s1,l,mid,x),Link(s2,mid+1,r,x);
}
void Bt ()
{
    num=0;bt(1,n);
    //for (int u=1;u<=n;u++) printf("%d ",pos[u]);
    for (int u=1;u<=n;u++)
    {
        int l=1,r=u;
        int p;
        while (l<=r)
        {
            int mid=(l+r)>>1;
            if (X[mid]+R[u]>=X[u]) {p=mid;r=mid-1;}
            else l=mid+1;
        }
        if (p!=u) Link(1,p,u-1,pos[u]);
    //  printf("%d %d %d\n",p,u-1,pos[u]);
        l=u;r=n;
        while (l<=r)
        {
            int mid=(l+r)>>1;
            if (X[mid]-R[u]<=X[u]) {p=mid;l=mid+1;}
            else r=mid-1;
        }
        if (p!=u) Link(1,u+1,p,pos[u]);
    //  printf("%d %d %d\n",u+1,u-1,pos[u]);
    }
}
bool vis[N];
int mymax (int x,int y){return x>y?x:y;}
int mymin (int x,int y){return x<y?x:y;}
int L1[N],R1[N];
void dfs (int x)
{
    vis[x]=true;
    for (int u=last[x];u!=-1;u=s[u].last)
    {
        int y=s[u].y;
        if (vis[y]==false) dfs(y);
        L1[x]=mymin(L1[x],L1[y]);
        R1[x]=mymax(R1[x],R1[y]);
    }
}
int dfn[N],low[N],belong[N],cnt,sta[N],lalal,shen;
bool in[N];
void dfs1 (int x)
{
    dfn[x]=low[x]=++lalal;
    sta[++cnt]=x;
    in[x]=true;
    for (int u=last[x];u!=-1;u=s[u].last)
    {
        int y=s[u].y;
        if (dfn[y]==-1)
        {
            dfs1(y);
            low[x]=mymin(low[x],low[y]);
        }
        else if (in[y]) low[x]=mymin(dfn[y],low[x]);
    }
    if (low[x]==dfn[x])
    {
        shen++;
        L1[shen]=1<<30;
        R1[shen]=0;
        int now;
        do
        {
            now=sta[cnt--];
            belong[now]=shen;
            L1[shen]=mymin(ll[now],L1[shen]);
            R1[shen]=mymax(rr[now],R1[shen]);
            in[now]=false;
        }while (now!=x);
    }
}
void solve ()
{
    shen=lalal=cnt=0;
    memset(dfn,-1,sizeof(dfn));
    memset(in,false,sizeof(in));
    for (int u=1;u<=num;u++)
        if (dfn[u]==-1)
            dfs1(u);
   // for (int u=1;u<=num;u++) printf("%d %d\n",u,belong[u]);
    num=0;memset(last,-1,sizeof(last));
    for (int u=1;u<=Num;u++)
        if (belong[s[u].x]!=belong[s[u].y])
        {
            int X=belong[s[u].x],Y=belong[s[u].y];
            num++;
            s[num].x=X;s[num].y=Y;
            s[num].last=last[X];
            last[X]=num;
        }
    memset(vis,false,sizeof(vis));
    for (int u=1;u<=shen;u++) 
        dfs(u);
    long long ans=0;
    for (int u=1;u<=n;u++) 
    {
        ans=(ans+(long long)u*(R1[belong[pos[u]]]-L1[belong[pos[u]]]+1)%MOD)%MOD;
    }
    printf("%lld\n",ans);
}
int main()
{
    Ins();
    Bt();
    solve();
    return 0;
}
OI界第一麻瓜
关注
专栏目录
wangjunrui666#bzoj-problem##2017. [Usaco2009 Nov]硬币游戏1
07-25
初始时，一个有N(5 ,000)枚硬币的堆栈放在地上，从堆顶数起的第I枚硬币的币值为C_i (1 <= C_i ,000)。第1行
【SNOI2017】炸弹——并查集解法
袋鼠与蒟蒻
09-21 868
【SNOI2017】 炸弹 Description 在一条直线上有 N 个炸弹，每个炸弹的坐标是 Xi，爆炸半径是 Ri，当一个炸弹爆炸时，如果另一个炸弹所在位置 Xj 满足： Xi−Ri≤Xj≤Xi+Ri,那么，该炸弹也会被引爆。 现在，请你帮忙计算一下，先把第 i 个炸弹引爆，将引爆多少个炸弹呢？ Input 第一行，一个数字 N，表示炸弹个数。 第 2∼N+12∼N+12∼N+1行，每行...
5017: [Snoi2017]炸弹（线段树建图 强连通）
jkchen's Haven
08-13 322
original link - https://www.lydsy.com/JudgeOnline/problem.php?id=5017 题意： 给出位置升序的nnn个炸弹，每个有位置PiP_iPi​和爆炸范围RiR_iRi​，爆炸后波及的范围为：[P−R,P+R][P-R,P+R][P−R,P+R]。当然，会连锁爆炸。设FiF_iFi​为点燃炸弹iii后可以引爆的炸弹数，求∑Fi\sum F_...
[SNOI2017]炸弹（tarjan+DP+优化建图）
jwg2732的专栏
08-29 194
洛谷题目传送门 解题思路 嗯，我觉得看出是个有向图并不难 把每个炸弹向他能炸到的连一条边 然后缩点，那么每个炸弹一定能炸到他所属的强连通分量里的所有炸弹 缩完点后是一个DAG，那么对于新图上的每个点来说，他在DAG上能到达所有店都能炸到 并且我们能发现一个性质，一个炸弹能炸到的一定是连续的一段，所以我们建这个DAG的反图，进行DP，更新能炸到的左端点和右端点 但是数据有点大，所以一个一个连边就会T 一般情况，我们会选择线段树优化建图 但是众所周知，线段树优化建图细节很多，很难调 所以我们要另想一种方法优化
[bzoj5017][线段树优化建边]炸弹
Gzb1128的栖息地
03-27 257
5017: [Snoi2017]炸弹 Time Limit: 30 Sec Memory Limit: 512 MB Submit: 243 Solved: 95 [Submit][Status][Discuss] Description 在一条直线上有 N 个炸弹，每个炸弹的坐标是 Xi，爆炸半径是 Ri，当一个炸弹爆炸时，如果另一个炸弹所在位置 Xj 满足： Xi−Ri≤Xj...
SNOI 2019 题解
FSYo 的博客
11-05 469
[SNOI2019]纸牌 又是麻将 dpdpdp，令 fi,j,kf_{i,j,k}fi,j,k​ 表示 (i,i+1,i+2)(i,i+1,i+2)(i,i+1,i+2) 选了 jjj 个，(i−1,i,i+1)(i-1,i,i+1)(i−1,i,i+1) 选了 kkk 的方案数 fi,l,j<−fi−1,j,k∗(C−l−j−k3+1)f_{i,l,j}<-f_{i-1,j,k}*...
BZOJ_problem_rank:BZOJ问题等级按交流比排序
05-25
#BZOJ Problem Rankrank.cpp 程序文件data.dat bzoj题库数据done.dat AC过的题，初始可以把所有A过的题粘进去，正常退出的话自动维护。black.dat 黑名单。选题时会跳过。错题、神题、没题面、不想做等等。//Thank ...
BZOJ第一部分
11-03
【BZOJ第一部分】 BZOJ，全称为“北京邮电大学在线评测系统”（Beijing Zhiyuan Online Judge），是中国最早的一批在线编程竞赛平台之一，它为编程爱好者提供了一个练习和检验编程技能的环境。在这个系统中，用户...
BZOJ全代码
02-05
【标题】"BZOJ全代码"所指的是一份包含BZOJ（BestCoder Zhijiang Online Judge）平台所有编程题目解决方案的压缩文件。BZOJ是一个在线编程竞赛平台，它提供了大量的算法题目供参赛者练习和提交代码，以检验编程能力...
BZOJ3230相似子串测试数据
08-19
标题"BZOJ3230相似子串测试数据"指的是一个与在线判题平台BZOJ（BestCoder ZOJ，中国大学算法竞赛在线评测系统）相关的编程问题，具体问题编号为3230，涉及的主要知识点是"相似子串"的搜索和处理。在编程竞赛中，...
[bzoj5017][随机化]炸弹
Rose_max的博客
04-24 236
Description 在一条直线上有 N 个炸弹，每个炸弹的坐标是 Xi，爆炸半径是 Ri，当一个炸弹爆炸时，如果另一个炸弹所在位置 Xj 满足： Xi−Ri≤Xj≤Xi+Ri,那么，该炸弹也会被引爆。 现在，请你帮忙计算一下，先把第 i 个炸弹引爆，将引爆多少个炸弹呢？ Input 第一行，一个数字 N，表示炸弹个数。 第 2∼N+1行，每行 2 个数字，表示 ...
SNOI2018 暴力记
C_K_Y_的博客
04-23 3277
DAY-1：下午中考体考，早上就没去机房，十点才起的床，中午考试前和班里同学谈笑风生了一会，结果都奶我进队，之后体考也凉凉的，50m起跑慢了，1000m没安排好时间，实心球离满分差的贼远，最后47分滚粗，班里绝大多数人都49。DAY0：省选前一天，敲了一天模板，还打了三道SA题，已经为我省选翻车做好了铺垫。晚上机房人几乎都走完了，我一天也没打几个模板，感觉真的要凉了，就再也没管了颓圮去了，打了五局...
SNOI2018 滚粗记
Rising_shit的博客
04-23 958
DAY 0:      为了早上能多睡一会，在学校旁边定了酒店，结果好像有点不适应的原因，9：30睡觉1：30感觉才睡着，期间除了各种辗转反侧，还体验了学校旁边酒店隔音的差，就像在自己房间一样，听完了隔壁放的电影，好像很煽情的样子...隔音这么差，连毛巾都没有你是怎么收的我325的房费=_=DAY1:     早上被隔壁的闹铃叫醒...     T1：诶这题好像每次相当于中间画了一个分界线，那么只...
bzoj5017 [Snoi2017]炸弹（tarjan缩点+拓扑序dp+线段合并+线段树优化建图）
Icefox的博客
12-22 945
首先我们可以发现每个点能引爆的炸弹是一个连续的区间。相邻的点能引爆的区间肯定不会相离。每个点向他能引爆的点连边，建出一张有向图，Tarjan缩点（一个scc内的点都可以互达，合并他们能到达的区间即可），然后拓扑序倒序dp即可。转移相当于线段合并。每个点记录他最后能引爆的区间，答案就是区间大小。 可是我们最坏情况有O(n2)O(n^2)条边，无论是时间上还是空间上都是无法接受的。还好我们每个点能连向
[Snoi2017]炸弹
zhouyuheng2003的博客
10-25 1222
题目相关 题面： BZOJ的题目链接：https://www.lydsy.com/JudgeOnline/problem.php?id=5017 LOJ的数据下载链接：https://loj.ac/problem/2255/testdata 一句话题意？题面那么短，自己看吧 解法 首先看到题目，容易发现一个性质，一个炸弹能炸到的范围是一个区间，这个性质看起来很简单，但是事实上非常有用，方便解题 ...
bzoj5017 [Snoi2017]炸弹 O(n)递推
C_K_Y_的博客
04-17 952
5017: [Snoi2017]炸弹Time Limit: 30 Sec  Memory Limit: 512 MBSubmit: 320  Solved: 118[Submit][Status][Discuss]Description在一条直线上有 N 个炸弹，每个炸弹的坐标是 Xi，爆炸半径是 Ri，当一个炸弹爆炸时，如果另一个炸弹所在位置 Xj 满足： Xi−Ri≤Xj≤Xi+Ri,那么，该...
SNOI2017炸弹
bansi8227的博客
07-15 155
这个东西其实我是不太会的……但是勉强卡过去了。 首先肯定是建有向图，然后求每个节点能访问的节点个数，最裸的打法就是按照题意枚举建边然后tarjan缩点，用bitset记录一下访问节点，但是bitset开不了那么大，只能拿到50%的分，至于开数组枚举我没有试，目测高不了多少。 然后思考这样一个问题，我建的那么多边真的有用么。于是只建离他最近的两个点的边，然后直接topsort统计ans...
5017: [Snoi2017]炸弹
蒟蒻的Bzoj记录
02-10 374
题目链接题目大意：在一条直线上有 N 个炸弹，每个炸弹的坐标是 Xi，爆炸半径是 Ri，当一个炸弹爆炸时，如果另一个炸弹所在位置 Xj 满足： Xi−Ri≤Xj≤Xi+Ri,那么，该炸弹也会被引爆。 现在，请你帮忙计算一下，先把第 i 个炸弹引爆，将引爆多少个炸弹呢？ 题解：若x能引爆y，add(x,y)，易知i的答案即为i能到达的编号max-编号min+1考虑到i的连边是一个区间，可以用线
[bzoj5017/Snoi2017]炸弹
WorldWide_D的博客
10-23 736
题目大意在一条直线上有 N 个炸弹，每个炸弹的坐标是 Xi，爆炸半径是 Ri，当一个炸弹爆炸时，如果另一个炸弹所在位置 Xj 满足： Xi−Ri≤Xj≤Xi+Ri,那么，该炸弹也会被引爆。 对于i等于1到n，求把第i个炸弹引爆会导致多少个炸弹爆炸。n≤500000，|Xi|≤101810^{18}，0≤Ri≤2⋅10182·10^{18}，保证Xi严格递增分析首先引爆的炸弹肯定是一个区间。
BZOJ2839:集合计数(容斥,组合数学)
最新发布
06-08
这个问题是一个比较经典的容斥原理的应用题目。我们可以考虑使用容斥原理来求解这个问题。具体来说，我们可以考虑求出所有集合的个数，然后减去至少包含一个特定元素的集合的个数，再加上至少包含两个特定元素的集合的个数，以此类推，直到加上包含所有特定元素的集合的个数。 对于一个大小为n的集合，其子集的个数为2^n。因此，所有集合的个数为2^n。现在我们考虑如何计算至少包含一个特定元素的集合的个数。我们可以将这个特定元素看作是一个元素，然后考虑其余n-1个元素组成的集合的个数，即2^(n-1)。因此，至少包含一个特定元素的集合的个数为2^(n-1)。同理，至少包含两个特定元素的集合的个数为C(2,n)*2^(n-2)，其中C(2,n)表示从n个元素中选取2个元素的组合数。 将以上结果代入容斥原理的公式中，即可得到最终的计数公式： sum((-1)^k*C(k,n)*2^(n-k),k=0,1,...,n) 其中，sum表示对k从0到n进行求和，C(k,n)表示从n个元素中选取k个元素的组合数，2^(n-k)表示剩余元素组成集合的个数。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值