【题】【二分答案&倍增(lca)】NKOJ 3560 运输计划

9 篇文章 0 订阅
4 篇文章 0 订阅

NKOJ 3560 运输计划
时间限制 : 20000 MS 空间限制 : 262144 KB

问题描述
公元 2044 年,人类进入了宇宙纪元。
L 国有 n 个星球,还有 n−1 条双向航道,每条航道建立在两个星球之间,这 n−1 条航道连通了 L 国的所有星球。
小 P 掌管一家物流公司, 该公司有很多个运输计划,每个运输计划形如:有一艘物流飞船需要从 ui 号星球沿最快的宇航路径飞行到 vi 号星球去。显然,飞船驶过一条航道是需要时间的,对于航道 j,任意飞船驶过它所花费的时间为 tj,并且任意两艘飞船之间不会产生任何干扰。
为了鼓励科技创新, L 国国王同意小 P 的物流公司参与 L 国的航道建设,即允许小P 把某一条航道改造成虫洞,飞船驶过虫洞不消耗时间。
在虫洞的建设完成前小 P 的物流公司就预接了 m 个运输计划。在虫洞建设完成后,这 m 个运输计划会同时开始,所有飞船一起出发。当这 m 个运输计划都完成时,小 P 的物流公司的阶段性工作就完成了。
如果小 P 可以自由选择将哪一条航道改造成虫洞, 试求出小 P 的物流公司完成阶段性工作所需要的最短时间是多少?

输入格式
/输入文件名为transport.in。/
第一行包括两个正整数n,m,表示L国中星球的数量及小P公司预接的运输计划的数量,星球从1到n编号。
接下来n-1行描述航道的建设情况,其中第i行包含三个整数ai,bi和ti,表示第i条双向航道修建在ai与bi两个星球之间,任意飞船驶过它所花费的时间为ti。
接下来m行描述运输计划的情况,其中第j行包含两个正整数ui和vi,表示第j个运输计划是从uj号星球飞往vj号星球。

输出格式
/输出文件名为transport.out。/
共1行,包含1个整数,表示小P的物流公司完成阶段性工作所需要的最短时间。

样例输入
样例输入1:
6 3
1 2 3
1 6 4
3 1 7
4 3 6
3 5 5
3 6
2 5
4 5

样例输出
样例输出1:
11

提示
输入输出样例1说明:
将第1条航道改造成虫洞,则三个计划耗时分别为11、12、11,故需要花费的时间为12;
将第2条航道改造成虫洞,则三个计划耗时分别为7、15、11,故需要花费的时间为15;
将第3条航道改造成虫洞,则三个计划耗时分别为4、8、11,故需要花费的时间为11;
将第4条航道改造成虫洞,则三个计划耗时分别为11、15、5,故需要花费的时间为15;
将第5条航道改造成虫洞,则三个计划耗时分别为11、10、6,故需要花费的时间为11。
故将第3条或第5条航道改造成虫洞均可使得阶段性工作的耗时最短,需要花费的时间为11。

数据规模于约定:
这里写图片描述

来源
感谢nodgd手打题目
NOIP 2015

思路:
该题即为删除一条边使得所有航线中最长的一条最短
考虑二分答案
judge:
找出所有路径长度大于当前枚举答案的路径,删除被所有路径覆盖的边中最长的一条,如果没有被所有路径覆盖的边,false。如果删除后还有路径大于枚举的答案,false;

结论:求被所有路径覆盖的边:http://blog.csdn.net/Y__XV/article/details/52473888

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
using namespace std;
const double lg2=log(2);
const int need=300003;

struct fy{int a,b,z,len;} e[need];//a、b表示端点,z表示lca,len表示a到b航线距离 

int dis[need],dis_[need],fa[need][19];
//dis表示该点到根节点的距离,dis_[i]表示i点到父亲的距离 
int n,m;
//.......................................................
int la[need*2],len[need*2],en[need*2],fi[need];
int tot;

void add(int a,int b,int c)
{
    tot++;
    la[tot]=fi[a];
    fi[a]=tot;
    en[tot]=b;
    len[tot]=c;
}
//.......................................................
int la1[need*2],len1[need*2],en1[need*2],fi1[need];
int tot1;

void add1(int a,int b)
{
    tot1++;
    la1[tot1]=fi1[a];
    fi1[a]=tot1;
    en1[tot1]=b;
}
//.......................................................
bool visit[need];
int dep[need],dm;

void dfs(int s)
{
    int k=ceil(log(dep[s])/lg2);
    dm=max(k,dm);
    for(int i=1;i<=k;i++) 
    {
        fa[s][i]=fa[fa[s][i-1]][i-1];
    }
    for(int t=fi[s],y;t;t=la[t])
    {
        y=en[t];
        if(visit[y]) continue;
        visit[y]=true;
        dep[y]=dep[s]+1;
        dis[y]=dis[s]+len[t];
        dis_[y]=len[t];
        fa[y][0]=s;
        add1(s,y);
        dfs(y);
    }
}
//.......................................................
int lca(int v,int u)
{
    if(dep[u]>dep[v]) swap(v,u);
    int p=dep[v]-dep[u];
    for(int i=0;i<=dm;i++) if((p>>i)&1) v=fa[v][i];
    if(v==u) return v;
    for(int i=dm;i>=0;i--) if(fa[v][i]!=fa[u][i]) v=fa[v][i],u=fa[u][i];
    return fa[v][0];
}
//.......................................................
int w[need],j[need],tt1;

void dfs1(int s)
{
    for(int t=fi1[s],y;t;t=la1[t])
    {
        y=en1[t];
        dfs1(y);
        j[s]+=j[y]; 
    }
}

bool judge(int s)
{
    tt1=0;
    memset(j,0,4*(n+1));
    for(int i=1;i<=m;i++) 
    {
        if(e[i].len>s) 
        {
            j[e[i].a]++,j[e[i].b]++;
            j[e[i].z]-=2;
            tt1++;
        }
    }
    dfs1(1);
    int flag=0;
    for(int i=1;i<=n;i++) if(j[i]==tt1) flag=max(flag,dis_[i]);
    if(flag==0) return false;
    for(int i=1;i<=m;i++) if(e[i].len-flag>s) return false;
    return true;
}
//.......................................................

int main_()
{
    scanf("%d%d",&n,&m);
    for(int i=1,a,b,c;i<n;i++)
    {
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c),add(b,a,c);
    }
    visit[1]=true;
    dfs(1);
    int l=0,r=0,mid;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&e[i].a,&e[i].b);
        e[i].z=lca(e[i].a,e[i].b);
        e[i].len=dis[e[i].a]+dis[e[i].b]-(dis[e[i].z]<<1);
        r=max(r,e[i].len);
        //cout<<e[i].z<<endl;
    }
    while(l<=r)
    {
        mid=(l+r)>>1;
        if(judge(mid)) r=mid-1;
        else l=mid+1;
    }
    printf("%d",l);
    return 0;
}

const int main_stack=16;
char my_stack[128<<20];

int main(){
     __asm__("movl %%esp, (%%eax);\n"::"a"(my_stack):"memory");
     __asm__("movl %%eax, %%esp;\n"::"a"(my_stack+sizeof(my_stack)-main_stack):"%esp");
     main_();
     __asm__("movl (%%eax), %%esp;\n"::"a"(my_stack):"%esp");
     return 0;
} 
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值