【codevs4632】【BZOJ4326】运输计划，链剖+二分+差分

最新推荐文章于 2019-08-11 16:28:42 发布
最新推荐文章于 2019-08-11 16:28:42 发布
阅读量1k 收藏 1
点赞数 2
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xym_CSDN/article/details/51174138
版权

传送门1
传送门2
写在前面：
当初洗澡的时候脑补了下这个题，觉得很简单直到看到题面，才发现自己记错了，自此开始不归路
思路：这里的链剖实际上就是求每个操作的路径长度（也可以LCA求），然后我想到了二分答案（最大求最小谁不会啊！）进行判断，也想到求交集了，但是卡在交集这里了，发现好几个问题，比如究竟删掉的边是交最多的还是权值最大的还是另有他法，而且这也用线段树来搞的话会被卡常……想的方法也是实现了好久……
二分总路径长度x，把大于x的操作在线段树上打一下标记并sum++，最后看一下每条边的标记个数，如果标记个数=sum那么这条边是可以删去的（如果删去那些标记个数＜sum的边，肯定有长度>x的路径没受到影响），最后判断操作中最长的路径长度-标记=sum的最长边是否小于等于x就可以了，复杂度 O(mlog2n)
95分线段树版本

#include<bits/stdc++.h>
#define M 300002
using namespace std;
int in()
{
    int t=0;char ch=getchar();
    while (ch>'9'||ch<'0') ch=getchar();
    while (ch>='0'&&ch<='9') t=t*10+ch-'0',ch=getchar();
    return t;
} 
int ans,n,m,tot,cnt,sum;
int num[M],fa[M],son[M],t[M],siz[M],L[M],pre[M],dep[M],first[M],top[M],tree[M<<2],a[M<<2];
bool flag[M];
struct edge
{
    int u,v,w,next;
}e[M<<1];
struct os
{
    int x,y,dis;
}q[M];
void add(int x,int y,int z)
{
    e[++tot].v=y;
    e[tot].u=x;
    e[tot].w=z;
    e[tot].next=first[x];
    first[x]=tot;
}
void dfs1(int x)
{
    siz[x]=1;
    for (int i=first[x];i;i=e[i].next)
    if (e[i].v!=fa[x])
    {
        dep[e[i].v]=dep[x]+1;
        fa[e[i].v]=x;
        dfs1(e[i].v);
        if (siz[e[i].v]>siz[son[x]]) son[x]=e[i].v;
        siz[x]+=siz[e[i].v]; 
    }
}
void dfs2(int x,int tp)
{
    top[x]=tp;
    L[x]=++cnt;
    pre[cnt]=x;
    if (son[x]) dfs2(son[x],tp);
    for (int i=first[x];i;i=e[i].next)
    if (e[i].v!=fa[x]&&e[i].v!=son[x]) dfs2(e[i].v,e[i].v);
}
void build(int now,int begin,int end)
{
    if (begin==end){tree[now]=t[pre[end]];return;}
    int mid=(begin+end)>>1;
    build(now<<1,begin,mid);
    build(now<<1|1,mid+1,end);
    tree[now]=tree[now<<1]+tree[now<<1|1];
}
int get(int now,int begin,int end,int l,int r)
{
    if (l<=begin&&end<=r) return tree[now];
    int mid=(begin+end)>>1,ans=0;
    if (mid>=l) ans+=get(now<<1,begin,mid,l,r);
    if (mid<r) ans+=get(now<<1|1,mid+1,end,l,r);
    return ans;
}
void pushdown(int now,int begin,int end)
{
    if (!a[now]) return;
    int mid=(begin+end)>>1;
    a[now<<1]+=a[now];
    a[now<<1|1]+=a[now];
    a[now]=0;
} 
void update(int now,int begin,int end,int l,int r,int mk)
{
    if (l<=begin&&end<=r) {a[now]+=mk;return;}
    pushdown(now,begin,end);
    int mid=(begin+end)>>1;
    if (mid>=l) update(now<<1,begin,mid,l,r,mk);
    if (mid<r) update(now<<1|1,mid+1,end,l,r,mk);
}
void fills(int x,int y,int mk)
{
    int f1=top[x],f2=top[y];
    while (f1!=f2)
    {
        if (dep[f1]<dep[f2]) swap(f1,f2),swap(x,y);
        update(1,1,cnt,L[f1],L[x],mk);
        x=fa[f1];f1=top[x];
     }
    if (dep[x]>dep[y]) swap(x,y);
    if (x!=y) update(1,1,cnt,L[son[x]],L[y],mk);
}
void down(int now,int begin,int end)
{
    if (begin==end) {num[pre[end]]=a[now];return;}
    pushdown(now,begin,end);
    int mid=(begin+end)>>1;
    down(now<<1,begin,mid);
    down(now<<1|1,mid+1,end);
}
bool check(int x)
{
    int maxn_q=0,maxn_e=0;
    for (int i=1;i<=m;i++)
    if (q[i].dis>x)
    {
        maxn_q=max(maxn_q,q[i].dis);
        if (!flag[i]) flag[i]=1,sum++,fills(q[i].x,q[i].y,1);
    }
    else
        if (flag[i]) flag[i]=0,sum--,fills(q[i].x,q[i].y,-1);
    down(1,1,cnt);//把标记放下来，查询每个线段被覆盖的次数
    for (int i=1;i<=cnt;i++)
    if (num[i]==sum) maxn_e=max(maxn_e,t[i]);
    return maxn_q-maxn_e<=x;
}
main()
{
    n=in();m=in();
    int x,y,z,f1,f2;
    int l=0,r=0;
    for (int i=1;i<n;i++)
    x=in(),y=in(),z=in(),
    add(x,y,z),add(y,x,z);
    dfs1(1);
    dfs2(1,1);
    for (int i=1;i<=tot;i++)
    {
        if (dep[e[i].v]<dep[e[i].u]) swap(e[i].v,e[i].u);
        t[e[i].v]=e[i].w;
    }
    build(1,1,cnt);
    for (int i=1;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        q[i].x=x;q[i].y=y;
        f1=top[x];f2=top[y];
        while (f1!=f2)
        {
            if (dep[f1]<dep[f2]) swap(f1,f2),swap(x,y);
            q[i].dis+=get(1,1,cnt,L[f1],L[x]);
            x=fa[f1];f1=top[x];
        }
        if (dep[x]>dep[y]) swap(x,y);
        if (x!=y) q[i].dis+=get(1,1,cnt,L[son[x]],L[y]);
        r=max(r,q[i].dis);
    }
    while (l<=r)
    {
        int mid=(l+r)>>1;
        if (check(mid)) ans=mid,r=mid-1;
        else l=mid+1;
    }
    printf("%d",ans);
}

后来又调了好久的差分……还是不太能理解
100分差分版本

#include<bits/stdc++.h>
#define M 300002
using namespace std;
int in()
{
    int t=0;char ch=getchar();
    while (ch>'9'||ch<'0') ch=getchar();
    while (ch>='0'&&ch<='9') t=t*10+ch-'0',ch=getchar();
    return t;
} 
int ans,n,m,tot,cnt;
int t[M],num[M],fa[M],son[M],siz[M],L[M],pre[M],dep[M],first[M],top[M],tree[M<<2];
struct edge
{
    int u,v,w,next;
}e[M<<1];
struct os
{
    int x,y,dis;
}q[M];
void add(int x,int y,int z)
{
    e[++tot].u=x;
    e[tot].v=y;
    e[tot].w=z;
    e[tot].next=first[x];
    first[x]=tot;
}
void dfs1(int x)
{
    siz[x]=1;
    for (int i=first[x];i;i=e[i].next)
    if (e[i].v!=fa[x])
    {
        dep[e[i].v]=dep[x]+1;
        fa[e[i].v]=x;
        dfs1(e[i].v);
        if (siz[e[i].v]>siz[son[x]]) son[x]=e[i].v;
        siz[x]+=siz[e[i].v]; 
    }
}
void dfs2(int x,int tp)
{
    top[x]=tp;
    L[x]=++cnt;
    pre[cnt]=x;
    if (son[x]) dfs2(son[x],tp);
    for (int i=first[x];i;i=e[i].next)
    if (e[i].v!=fa[x]&&e[i].v!=son[x]) dfs2(e[i].v,e[i].v);
}
void build(int now,int begin,int end)
{
    if (begin==end){tree[now]=t[pre[end]];return;}
    int mid=(begin+end)>>1;
    build(now<<1,begin,mid);
    build(now<<1|1,mid+1,end);
    tree[now]=tree[now<<1]+tree[now<<1|1];
}
int get(int now,int begin,int end,int l,int r)
{
    if (l<=begin&&end<=r) return tree[now];
    int mid=(begin+end)>>1,ans=0;
    if (mid>=l) ans+=get(now<<1,begin,mid,l,r);
    if (mid<r) ans+=get(now<<1|1,mid+1,end,l,r);
    return ans;
}
bool check(int x)
{
    memset(num,0,sizeof(num));
    int sum=0,mp=0,l,r,f1,f2,maxn_q=0,maxn_e=0;
    for (int i=1;i<=m;i++)
    if (q[i].dis>x)
    {
        sum++;
        maxn_q=max(maxn_q,q[i].dis);
        l=q[i].x;r=q[i].y;
        f1=top[l];f2=top[r];
        while (f2!=f1)
        {
            if (dep[f1]<dep[f2]) swap(f1,f2),swap(l,r);
            num[L[f1]]++;num[L[l]+1]--;
            l=fa[f1];f1=top[l];
        }
        if (l==r)continue;
        if (dep[l]>dep[r]) swap(l,r);
        num[L[son[l]]]++;num[L[r]+1]--; 
    }
    for (int i=1;i<=cnt;i++)
    {
        mp+=num[i];
        if (mp==sum) maxn_e=max(maxn_e,t[pre[i]]);
    }
    return maxn_q-maxn_e<=x;
}
main()
{
    n=in();m=in();
    int x,y,z,f1,f2;
    int l=0,r=0;
    for (int i=1;i<n;i++)
    x=in(),y=in(),z=in(),
    add(x,y,z),add(y,x,z);
    dfs1(1);
    dfs2(1,1);
    for (int i=1;i<=tot;i++)
    {
        if (dep[e[i].v]<dep[e[i].u]) swap(e[i].v,e[i].u);
        t[e[i].v]=e[i].w;
    }
    build(1,1,cnt);
    for (int i=1;i<=m;i++)
    {
        x=in();y=in();
        q[i].x=x;q[i].y=y;
        f1=top[x];f2=top[y];
        while (f1!=f2)
        {
            if (dep[f1]<dep[f2]) swap(f1,f2),swap(x,y);
            q[i].dis+=get(1,1,cnt,L[f1],L[x]);
            x=fa[f1];f1=top[x];
        }
        if (dep[x]>dep[y]) swap(x,y);
        if (x!=y) q[i].dis+=get(1,1,cnt,L[son[x]],L[y]);
        r=max(r,q[i].dis);
    }
    while (l<=r)
    {
        int mid=(l+r)>>1;
        if (check(mid)) ans=mid,r=mid-1;
        else l=mid+1;
    }
    printf("%d",ans);
}

这道题给我的启发还是挺深的。，别人的code是别人思想方法的结晶，即使是偶尔想偷懒地模仿也是不可能成功的，更是白白浪费时间，深入的思考才能真正体会到题目的内涵与解决方法；小细节也是要注意的，而且不要去尝试从code上体会别人的思想，而是自己思考后结合别人的文字思路来加深印象与体会

iamxym
关注
【BZOJ 4326】运输计划【树链剖分+差分+二分答案】
skysys的研究小屋
11-08 494
首先，此题要求最大值最小，容易可以想到二分。现在要想如何check显然，要变为虫洞的边应该在所有比mid值大的路径上，于是我们需要快速求出所有这类路径的交集上的最大边。我们可以用差分来进行求交集，之后乱搞一下就好了。上面是Cycleke的原话。#include<bits/stdc++.h> using namespace std; #define g() getchar() #define d(x)
BZOJ 2135 刷题计划（贪心，求导，二分）【BZOJ 修复工程】
繁凡さん的博客
09-07 451
【BZOJ修复计划 #5】BZOJ 2135 刷题计划 【国家集训队作业】
BZOJ_P4326[NOIP]2015 运输计划(LCA+Tarjan+二分)
BeiYu
01-12 922
NOIP2015 运输计划Time Limit: 30 Sec Memory Limit: 128 MB Submit: 248 Solved: 166 [Submit][Status][Discuss] Description公元 2044 年，人类进入了宇宙纪元。L 国有 n 个星球，还有 n−1 条双向航道，每条航道建立在两个星球之间，这 n−1 条航道连通了 L 国的所有星球。小
[BZOJ4326][NOIP2015]运输计划（二分答案+树上差分）
女装的你，如此好看！
08-25 710
欢迎使用Markdown编辑器写博客本Markdown编辑器使用StackEdit修改而来，用它写博客，将会带来全新的体验哦： Markdown和扩展Markdown简洁的语法 代码块高亮 图片链接和图片上传 LaTex数学公式 UML序列图和流程图 离线写博客 导入导出Markdown文件 丰富的快捷键 快捷键 加粗 Ctrl + B 斜体 Ctrl + I 引用 Ctrl
BZOJ4326(NOIP2015)运输计划--二分+LCA+差分
Greninja的博客
10-16 696
bzoj4326
[BZOJ4326][NOIP2015]运输计划（二分+dfs序+树上差分）
zyf2000
10-18 2002
“过了这么多年，我终于发现了那些让我喜欢的人的特质，他们总是生机勃勃的，有自己喜欢的事情……容得下生命的不完美，也经得起世事的颠簸。”——扶南
bzoj4326
weixin_30378623的博客
09-15 72
二分+树剖+差分 之前的做法naive，莫名其妙的wa，明明uoj95分 看到最小最大上二分，树上路径问题直接剖，然后问题就转化成了一个判定问题，每次二分出最长路径长度，问能不能达到。那么我们就把所有长度大于二分出的d的路径拉出来，求出他们公共路径的最大长度，看减去能不能满足。那么现在的问题就转化成了求路径的交。我们可以利用bzoj4390的树上差分的方法解决，我们先把边转化为点，即...
BZOJ 2144 跳跳棋（LCA+欧几里德+二分答案）
Mogician_Evian的博客
10-28 573
跳跳棋问题描述 跳跳棋是在一条数轴上进行的。棋子只能摆在整点上。每个点不能摆超过一个棋子。我们用跳跳棋来做一个简单的游戏：棋盘上有3颗棋子，分别在a，b，c这三个位置。我们要通过最少的跳动把他们的位置移动成x，y，z。（棋子是没有区别的）跳动的规则很简单，任意选一颗棋子，对一颗中轴棋子跳动。跳动后两颗棋子距离不变。一次只允许跳过1颗棋子。 写一个程序，首先判断是否可以完成任务。如
bzoj2144 跳跳棋 二分+lca
老臣
08-15 739
这道题刚开始看以为是DP，，，我真是智障了。该DP的看成搜索，该二分的看成DP。。这状态有点危险。。 题意:现在有三个人a,b,c，每一个人可以跳过他左/右边的人，比如说一个人在x，一个在y，那么x可以跳到y+(y-x)（y>x），左边同理。。但是两个人之间不能有另一个人。 现在有三个目标位置和起始位置，让你计算最少的步数让三个人到达目标位置，并不要求人的顺序（即无论哪一个人都可以去任意的终点）
BZOJ2654——WQS二分+最小生成树
stevensonson的博客
05-08 854
Description 给你一个无向带权连通图，每条边是黑色或白色。让你求一棵最小权的恰好有need条白色边的生成树。 题目保证有解。 Input 第一行V,E,need分别表示点数，边数和需要的白色边数。 接下来E行,每行s,t,c,col表示这边的端点(点从0开始标号)，边权，颜色(0白色1黑色)。 Output 一行表示所求生成树的边权和。 V&lt;=50000,E&lt;...
[二分答案 树链剖分] BZOJ 4326 NOIP2015 运输计划
雯舞
04-22 641
考虑二分 把所有超过mid的链取交集，然后取交集上最长的边权赋0，判断 #include #include #include #define V G[p].v using namespace std; inline char nc() { static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+frea
bzoj 4326 运输计划 （树链剖分 + 树上差分 + 二分）
阿柟喜欢吃小星星？
10-26 289
不想贴题面系列…… 这道题就是先二分一个答案，然后check，如果路线的总长度比二分的答案大，就把它们拎出来进行树上差分（求它们所有路线都经过的边），然后把这条边的权值减掉（贪心），反复check即可。#include <bits/stdc++.h> using namespace std;const int N = 300010; int n, m;struct E { int v, w
【二分+LCA+差分】BZOJ4326(NOIP2015)[运输计划]题解
ZigZagK的博客
10-16 470
题目概述给 nn 个点的树和 mm 个运输计划( x→yx\to y )，完成所有运输计划的时间是运输计划中的最大时间。现在可以把一条边的边权改为 00 ，求完成所有运输计划的最快时间。解题报告由于求最大值的最小值，所以想到二分枚举答案 midmid 。然后预先求出LCA算出每个运输计划的时间，我们就知道有多少个运输计划超出了 midmid 。对于每个超出 midmid 的运输计划，我们肯定需要把这
BZOJ 4326: NOIP2015 运输计划
Doleful Aria
09-01 597
去年这道水题没做出来 然而今儿看到正解也是崩溃的 正解是个垃圾树剖。。 noip什么时候开始考树剖了？？？ 好吧去年联赛的时候还没学树剖。。。第一遍dfsdfs除了处理那些个节点的重量外顺便处理下每个节点到根的距离 其他和普通树剖没什么区别
学习记录——树上差分（BZOJ4326，BZOJ3631）
98kaiii的博客
08-11 140
前言 最近做的很多题都与树上差分有点儿关系，太久没用有点儿忘了，用这篇博客总结一下。 什么是差分 可以先考虑差分数组，即可以在线性时间内对一个数组的进行次区间修改（区间加减）。即额外维护b数组，初始化全为0，每一次区间修改，如把数组中l~r的数都加上a，那么就把b[l]+=a,b[r+1]-=a，然后多次修改之后，b[i]数组的前缀和就是原数组第i个节点的修改的值。 对应到树上，就是把从节点x到节...
bzoj4326: NOIP2015 运输计划
xumingyang0的博客
05-31 324
题目 题解： 大佬的博客 思路其实很简单，就是tarjan求lca+二分答案+树上差分（貌似不需要解释，看到这个思路就应该基本上会写了吧），实现起来也听简单的 #include&lt;bits/stdc++.h&gt; using namespace std; const int N=300003; struct node{ int to,w,ne; }e[N&lt;&lt;1...
bzoj 4453: cys就是要拿英魂！（后缀数组+二分）
clover_hxy的博客
08-21 728
4453: cys就是要拿英魂！ Time Limit: 3 Sec  Memory Limit: 128 MB Submit: 77  Solved: 39 [Submit][Status][Discuss] Description pps又开始dota视频直播了！一群每天被pps虐的蒟蒻决定学习pps的操作技术，他们把pps在这局放的技能记录了下 来，每个技能用一个字符表示
【BZOJ4370】【IOI2015】horses 数据结构 平衡树+线段树
function2的博客
02-04 1260
4370: [IOI2015]horses马Time Limit: 30 Sec Memory Limit: 1500 MB Description像他的祖先一样，Mansur喜欢繁殖马匹。目前，他拥有哈萨克斯坦最大的马场。以前情况可不是这样，N年前Mansur年轻时，他只拥有一匹马，但他一直梦想着成为富豪，最终,他美梦成真。 按照时间的先后顺序将年份编号为0到N-1（即N-1年是最近的一年）
【NOIP2015】【BZOJ4326】运输计划
function2的博客
11-11 6009
4326: NOIP2015 运输计划Time Limit: 20 Sec Memory Limit: 128 MB Submit: 30 Solved: 18 [Submit][Status][Discuss] Description公元 2044 年，人类进入了宇宙纪元。L 国有 n 个星球，还有 n−1 条双向航道，每条航道建立在两个星球之间，这 n−1 条航道连通了 L 国的所有星
bzoj做题计划 excel
最新发布
07-16
### 回答1： bzoj作为一个计算机竞赛的在线评测系统，不仅可以提供大量的题目供程序员练习和学习，还可以帮助程序员提升算法和编程能力。为了更好地利用bzoj进行题目的学习和刷题，制定一个bzoj做题计划是非常有必要的。 首先，我们需要合理安排时间，每天留出一定的时间来做bzoj的题目。可以根据自己的时间安排，每天挑选适量的题目进行解答。可以先从难度较低的题目开始，逐渐提高难度，这样既能巩固基础知识，又能挑战自己的思维能力。 其次，要有一个计划和目标。可以规划一个每周或每月的题目数量目标，以及每个阶段要学习和掌握的算法知识点。可以根据bzoj的题目分类，如动态规划、图论、贪心算法等，结合自己的实际情况，有针对性地选择题目进行学习。 此外，要充分利用bzoj提供的资源。bzoj网站上有很多高质量的题解和优秀的解题代码，可以参考和学习。还有相关的讨论区，可以与其他程序员交流和讨论，共同进步。 最后，要坚持并保持思考。做题不是单纯为了刷数量，更重要的是学会思考和总结。遇到难题时，要有耐心，多思考，多尝试不同的解法。即使不能一次性解出来，也要学会思考和分析解题过程，以及可能出现的错误和优化。 总之，bzoj做题计划的关键在于合理安排时间、制定目标、利用资源、坚持思考。通过有计划的刷题，可以提高算法和编程能力，并培养解决问题的思维习惯，在计算机竞赛中取得更好的成绩。 ### 回答2： bzoj做题计划是指在bzoj这个在线测评系统上制定一套学习和刷题的计划，并且将计划记录在excel表格中。该计划主要包括以下几个方面的内容。 首先是学习目标的设定。通过分析自己的水平和知识缺口，可以设定一个合理的目标，比如每天解决一定数量的题目或者提高特定的算法掌握程度。 其次是题目选择的策略。在excel表格中可以记录下自己选择的题目编号、题目类型和难度等信息。可以根据题目的类型和难度来安排每天的刷题计划，确保自己可以逐步提高技巧和解题能力。 然后是学习进度的记录和管理。将每天的完成情况记录在excel表格中，可以清晰地看到自己的学习进度和任务完成情况。可以使用图表等功能来对学习进度进行可视化展示，更好地管理自己的学习计划。 同时，可以在excel表格的备注栏中记录下每道题目的解题思路、关键点和需要复习的知识点等信息。这样可以方便自己回顾和总结，巩固所学的知识。 最后，可以将excel表格与其他相关资料进行整合，比如算法教材、题目解析和学习笔记等。这样可以形成一个完整的学习档案，方便自己进行系统的学习和复习。 总之，bzoj做题计划excel的制定和记录可以帮助我们更加有条理和高效地进行学习和刷题。通过合理安排学习目标和题目选择策略，记录学习进度和思路，并整合其他学习资料，我们可以提高自己的解题能力，并在bzoj上取得更好的成绩。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值