码迷,mamicode.com
首页 > 其他好文 > 详细

[tjoi2012]桥

时间:2018-10-05 18:54:42      阅读:141      评论:0      收藏:0      [点我收藏+]

标签:get   code   math   n+1   esc   long   而且   rsa   algorithm   

Description
有n个岛屿,m座桥,每座桥连通两座岛屿,桥上会有一些敌人,玩家只有消灭了桥上的敌人才能通过,与此同时桥上的敌人会对玩家造成一定伤害。而且会有一个大Boss镇守一座桥,以玩家目前的能力,是不可能通过的。而Boss是邪恶的,Boss会镇守某一座使得玩家受到最多的伤害才能从岛屿1到达岛屿n(当然玩家会选择伤害最小的路径)。
问,Boss可能镇守桥有哪些。

Input
第一行两个整数n,m
接下来m行,每行三个整数s,t,c,表示一座连接岛屿s和t的桥上的敌人会对玩家造成c的伤害。
1<=n<=100000,1<=m<=200000,1<=c<=10000,数据保证玩家可以从岛屿1到达岛屿n

Output
一行,两个整数d,cnt,d表示有Boss的情况下,玩家要受到的伤害,cnt表示Boss有可能镇守的桥的数目。

Sample Input
3 4
1 2 1
1 2 2
2 3 1
2 3 1

Sample Output
3 1

首先简化题面,问去掉原图的一条边之后,令新图S到T的最短路长度最大的方案数

我们要明白一件事,割的边肯定是最短路径上的边,那么我们就找一条最短路径,考虑经过其他的边的贡献。如果经过其他路径的话,我们一定是走形如\(S\rightarrow A\rightarrow B\rightarrow T\)的形式,并且\(A\rightarrow B\)一定是一条边,是路径我就在其中找一条边出来充当\(A\rightarrow B\)

那么\(A\rightarrow B\)会对哪些点造成贡献呢?首先明白\(S\rightarrow A\)的路线,一定是\(S\rightarrow lca(A,T)\rightarrow A\),那么有贡献的肯定是\(lca(A,T)\)下面的点。上界找完了,下界呢?

下界其实和上界类似,就是\(lca(T,v)\),至于为什么?我不会。。。你们就感性理解一下吧。。。(挖坑,如果我之后会证明了我会来填坑的)

这样弄好之后,我们相当于枚举一条边,在一段区间进行更新,这样我们就可以用线段树来进行维护了

/*program from Wolfycz*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline int read(){
    int x=0,f=1;char ch=getchar();
    for (;ch<'0'||ch>'9';ch=getchar())  if (ch=='-')    f=-1;
    for (;ch>='0'&&ch<='9';ch=getchar())    x=(x<<1)+(x<<3)+ch-'0';
    return x*f;
}
inline void print(int x){
    if (x>=10)  print(x/10);
    putchar(x%10+'0');
}
const int N=1e5,M=2e5;
int pre[(M<<1)+10],now[N+10],child[(M<<1)+10],val[(M<<1)+10];
int Fa[N+10],Lca[N+10],stack[N+10],dfn[N+10];
bool vis[N+10],flag[(M<<1)+10];
int n,m,tot,num,Ans,Max,Dis;
struct S1{
    #define ls (p<<1)
    #define rs (p<<1|1)
    #define fa (p>>1)
    struct node{
        int x,v;
        bool operator <(const node &a)const{return v<a.v;}
    }Q[N+10];
    int tot;
    void insert(int x,int v){
        Q[++tot]=(node){x,v};
        int p=tot;
        while (p!=1&&Q[p]<Q[fa])    swap(Q[p],Q[fa]),p=fa;
    }
    void Delete(){
        Q[1]=Q[tot--];
        int p=1,son;
        while (ls<=tot){
            if (rs>tot||Q[ls]<Q[rs])    son=ls;
            else    son=rs;
            if (Q[son]<Q[p])    swap(Q[p],Q[son]),p=son;
            else    break;
        }
    }
}Heap;
struct Segment{
    int tree[(N<<2)+10];
    Segment(){memset(tree,127,sizeof(tree));}
    void Add_tag(int p,int v){tree[p]=min(tree[p],v);}
    void pushdown(int p){
        if (tree[p]==inf)   return;
        Add_tag(ls,tree[p]);
        Add_tag(rs,tree[p]);
//      Lazy[p]=inf;
    }
    void Modify(int p,int l,int r,int x,int y,int v){
        if (x<=l&&r<=y){
            Add_tag(p,v);
            return;
        }
//      pushdown(p);
        int mid=(l+r)>>1;
        if (x<=mid) Modify(ls,l,mid,x,y,v);
        if (y>mid)  Modify(rs,mid+1,r,x,y,v);
//      tree[p]=max(tree[ls],tree[rs]);
    }
    void traversals(int p,int l,int r){
        if (l==r){
            if (tree[p]>Max)    Max=tree[p],Ans=0;
            Ans+=tree[p]==Max;
            return;
        }
        pushdown(p);
        int mid=(l+r)>>1;
        traversals(ls,l,mid);
        traversals(rs,mid+1,r);
    }
}Tree;
struct S2{
    int dis[N+10];
    S2(){memset(dis,63,sizeof(dis));}
}Frw,Bck;
void join(int x,int y,int z){pre[++tot]=now[x],now[x]=tot,child[tot]=y,val[tot]=z;}
void insert(int x,int y,int z){join(x,y,z),join(y,x,z);}
void dijkstra(int x,int *dis){
    memset(vis,0,sizeof(vis));
    Heap.insert(x,dis[x]=0);
    while (Heap.tot){
        int Now=Heap.Q[1].x;
        Heap.Delete();
        if (vis[Now])   continue;
        vis[Now]=1;
        for (int p=now[Now],son=child[p];p;p=pre[p],son=child[p]){
            if (dis[son]>dis[Now]+val[p]){
                dis[son]=dis[Now]+val[p];
                Heap.insert(son,dis[son]);
            }
        }
    }
}
void put_flag(int x){
    for (int p=now[x],son=child[p];p;p=pre[p],son=child[p]){
        if (Frw.dis[x]+val[p]+Bck.dis[son]==Dis){
            flag[p]=flag[p+(p&1?1:-1)]=1;
            put_flag(son);
            return;
        }
    }
}
void Find_fa(int x){
    for (int p=now[x],son=child[p];p;p=pre[p],son=child[p])
        if (Frw.dis[x]+val[p]==Frw.dis[son]&&!Fa[son])  Fa[son]=x,Find_fa(son);
}
void Fill(int x,int v){
    if (Lca[x]) return;
    Lca[x]=v;
    for (int p=now[x],son=child[p];p;p=pre[p],son=child[p]) if (Fa[son]==x) Fill(son,v);
}
void Get_Lca(int x){
    int top=0;
    while (x){
        Fill(x,x);
        stack[++top]=x;
        x=Fa[x];
    }
    for (int i=top;i;i--)   dfn[stack[i]]=num++;
    num--;
}
int main(){
    int n=read(),m=read();
    for (int i=1;i<=m;i++){
        int x=read(),y=read(),z=read();
        insert(x,y,z);
    }
    dijkstra(1,Frw.dis);
    dijkstra(n,Bck.dis);
    Dis=Frw.dis[n];
    put_flag(1);
    Find_fa(1);
    Get_Lca(n);
    for (int i=1;i<=tot;i+=2){
        if (flag[i])    continue;
        int x=child[i],y=child[i+1],v=val[i];
        if (Lca[x]==Lca[y]) continue;
        if (dfn[Lca[x]]>dfn[Lca[y]])    swap(x,y);
        int tmp=Frw.dis[x]+v+Bck.dis[y];
        Tree.Modify(1,1,num,dfn[Lca[x]]+1,dfn[Lca[y]],tmp);
    }
    Tree.traversals(1,1,num);
    printf("%d %d\n",Max,Max==Dis?m:Ans);
    return 0;
}

[tjoi2012]桥

标签:get   code   math   n+1   esc   long   而且   rsa   algorithm   

原文地址:https://www.cnblogs.com/Wolfycz/p/9745288.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!