[tjoi2012]桥

标签：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