10.15T2 生成树+非树边暴力

Description

在这一切发生之前，学园都市表面看起来其实是非常平静的。但仅仅只是表面而已。例如第五市场下面的地下街又传来了哔哩哔哩的声音，很显然，这是电击使在使用能力，而她的面前站着一个抬起右臂的少年。为什么放电妹要去攻击刺猬头？通过御坂美琴，我们可以知道，刺猬头——上条当麻在大霸星祭的时候输掉了和御坂美琴的赌约。
作为惩罚游戏，上条当麻这一天要无条件接受美琴的所有条件。而上条当麻趁御坂美琴在办理他们两人的双人通话套餐的一小会儿时间，又“勾搭”上了御坂妹。
在御坂妹和最后之作离开过后，御坂美琴拉着上条当麻去了地下街的游戏厅。现在的美琴小姐十分气氛，因为总有人来打扰她和上条当麻的私会。上条当麻现在也发现了这个问题，正在（被迫）想如何让美琴开心。这时，一台崭新的街机出现在他们面前，而最终大奖是一个3米长的呱太玩偶。此时的美琴，早已准备跃跃欲试，可当麻却突然挡在了他的面前，“这次就让我来吧！（被迫）”
这台游戏机里面的游戏是一个十分无聊的游戏：首先，屏幕上面会出现一张多个点多条无向边的图，两点之间最多只会有一条边，且这张地图是不会有自环的。每条边的权值也标在了图上。每一关，屏幕上会给出两点，让你在三秒之内答出这两点之间的最短距离。游戏保证两点之间一定有边。
作为无能力者（bug 者），当麻的计算能力肯定没有美琴强，但他也不像在美琴面前出丑。所以他悄悄地找到了来自于BSKingSchool 的你，让你帮他计算一下游戏的答案。

Input

第一行为两个整数 n，m。表示地图的点数和边数。
第二行到第m+1 行，每行三个整数 x,y ,z ，表示 x，y 两点之间存在一条直接路径，长度为z 。
第m+2 行为一个整数 q，表示机器给出了q 个提问。
第m+3 行到第m+3+q 行，每行二个整数 a，b 。表示机器所给出的两点。

Output

共q行，每行一个数表示答案。

Sample Input

3 3 1 2 3 2 3 1 3 1 5 3 1 2 1 3 2 3

Sample Output

3 4 1 样例解释 从1号节点直接走到至2号节点，距离为3。 从1号节点经2号节点走到至3号节点，距离为4。 从2号节点直接走到至3号节点，距离为1。

Hint

数据范围及约定

Source

from LH&LYT

Hide Information »

 

Case Time Limit:	2000 ms
Memory Limit:	256536 KB
Case score:	10
Comparison:	Traditional
Level:	4

Congratulation !
Users Submitted:	9
Users Solved:	8
Total Submits:	17
Accepted:	8
Time Out:	3
Wrong Answer:	1
Runtime Error:	3
Compile Error:	2

Submit         Status          Discuss

 

 

Add notes and tags...

 

 

 Mark

 

 

题解：由于m只比n大一点，我们可以先找个树出来，然后两点之间最短路要么经过树上，要么经过非树边

所以我们对非树边每个点跑个spfa，然后询问时每个依次暴力比较，因为只有21，所以不会T

code：

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
#define N 300006
using namespace std;
struct node{
    long long u,v,w;
}e[N],e1[N],a[N],temp[N];
long long cnt1,first1[N],nxt1[N];
void add1(long long u,long long v,long long w){
    e1[++cnt1].u=u;
    e1[cnt1].v=v;
    e1[cnt1].w=w;
    nxt1[cnt1]=first1[u];
    first1[u]=cnt1;
}
long long first[N],nxt[N],cnt;
void add(long long u,long long v,long long w){
    e[++cnt].u=u;
    e[cnt].v=v;
    e[cnt].w=w;
    nxt[cnt]=first[u];
    first[u]=cnt;
}
long long fa[N];
long long find(long long x){
    if(x!=fa[x])return fa[x]=find(fa[x]);
    return fa[x];
}
void merge(long long x,long long y){
    long long f1=find(x),f2=find(y);
    if(f1!=f2){
        fa[f1]=f2;
    }
}
long long dep[N],f[N][32];
long long lca(long long x,long long y){
    if(dep[y]>dep[x])swap(x,y);
    for(long long i=30;i>=0;i--){
        if(dep[f[x][i]]>=dep[y])
        x=f[x][i];
        if(x==y)return x;
    }
    for(long long i=30;i>=0;i--){
        if(f[x][i]==f[y][i])continue;
        x=f[x][i],y=f[y][i];
    }
    return f[x][0];
}
long long check[N],n,m,cnt_;
void kruskal(){
    long long cnt_=0;
    for(long long i=1;i<=m;i++){
        long long u=a[i].u,v=a[i].v,w=a[i].w;
        if(find(u)!=find(v)){
            merge(u,v);
            check[i]=1;
            add1(u,v,w);
            add1(v,u,w);
            cnt_++;
            if(cnt_==n-1)return;
        }
    }
}
long long vis[N],dis[50][N];
void spfa(long long s,long long now){
    queue<long long>q;
    memset(vis,0,sizeof vis);
    dis[s][now]=0;
    vis[now]=1;
    q.push(now);
    while(!q.empty()){
        long long u=q.front();
        q.pop();
        vis[u]=0;
        for(long long i=first[u];i;i=nxt[i]){
            long long v=e[i].v;
            if(dis[s][v]>dis[s][u]+e[i].w){
                dis[s][v]=dis[s][u]+e[i].w;
                if(!vis[v]){
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
}
bool cmp(const node &a,const node &b){
    return a.w<b.w; 
}
long long distree[N];
void dfs(long long x){
    for(long long i=first1[x];i;i=nxt1[i]){
        long long v=e1[i].v;
        if(v==f[x][0])continue;
        f[v][0]=x;
        dep[v]=dep[x]+1;
        distree[v]=distree[x]+e1[i].w;
        dfs(v);
    }
}
long long tot,tot_point,dian[N];
long long read(){
    long long x=0,f=1;
    char c=getchar();
    while(!isdigit(c)){
        if(c=='-')f=-1;
        c=getchar();
    }
    while(isdigit(c)){
        x=(x<<3)+(x<<1)+c-'0';
        c=getchar();
    }
    return x*f;
}
int main(){
//    freopen("game.out","w",stdout);
    memset(dis,0x3f3f3f3f,sizeof dis);
    n=read(),m=read();
    for(long long i=1;i<=n;i++)fa[i]=i;
    for(long long i=1;i<=m;i++){
        a[i].u=read(),a[i].v=read(),a[i].w=read();
        add(a[i].u,a[i].v,a[i].w);
        add(a[i].v,a[i].u,a[i].w);
    }
    sort(a+1,a+m+1,cmp);
    kruskal();
    for(long long i=1;i<=m;i++){
        if(check[i])continue;
//        cout<<a[i].u<<" ||| "<<a[i].v<<" ||| "<<a[i].w<<'
';
        temp[++tot].u=a[i].u;
        temp[tot].v=a[i].v;
        temp[tot].w=a[i].w;
        dian[++tot_point]=a[i].u;
        dian[++tot_point]=a[i].v;
        spfa(tot_point-1,a[i].u);
        spfa(tot_point,a[i].v);
    }
    //cout<<dian[3]<<" &&&& "<<dis[3][2]<<'
';
    f[1][0]=1;
    dfs(1);
//    for(int i=1;i<=tot_point;i++){
//        for(int j=1;j<=n;j++){
//            cout<<"st->"<<dian[i]<<" ed->"<<j<<" dist->"<<dis[i][j]<<'
';
//        }
//    }
//    for(int i=1;i<=n;i++){
//        cout<<"father->"<<f[i][0]<<" dis->"<<distree[i]<<" dep->"<<dep[i]<<endl;
//    }
    for(long long i=1;i<=30;i++){
        for(long long j=1;j<=n;j++){
            f[j][i]=f[f[j][i-1]][i-1];
        }
    }
//    cout<<"lca->"<<lca(5,2)<<endl;
//    for(int i=1;i<=n;i++){
//        cout<<"father->"<<f[i][1]<<" dis->"<<distree[i]<<" dep->"<<dep[i]<<endl;
//    }
    long long Q;
    Q=read();//dp，数论，图论，数据结构高级算法 
    while(Q--){
        long long u,v;
        u=read();v=read();
        long long ans=distree[u]+distree[v]-2*distree[lca(u,v)];
        for(long long i=1;i<=tot;i++){
            ans=min(ans,dis[i*2-1][u]+dis[i*2][v]+temp[i].w);
            ans=min(ans,dis[i*2-1][v]+dis[i*2][u]+temp[i].w);
        }
        cout<<ans<<'
';
    }
    return 0;
}
/*
5 6
1 2 1
2 3 2
3 4 3
4 5 4
5 1 5
1 3 1
5
1 5
2 4
3 4
2 1
3 5
*/

over