货车运输题解

https://i.cnblogs.com/EditPosts.aspx?opt=1

题目要求的是求出点x到点y能运输的最大货物量,换句话说就是找到x到y的一条路径,使最小边权值尽可能大。

第一眼想到用类似最短路的方法,但是数据范围太大,如果个x,y都跑一遍最短路显然会超时。

其实不难发现,我们只需要选择尽可能大的边,让他们组成一颗颗最大生成树,保证了原来图的联通性,又把一张图退化成树。

边权最小的值在倍增的过程中顺便处理,时间复杂度O(q*logn)

代码有点长,但都是模板

#include<cstdio>
#include<iostream>
#include<algorithm>
#define N 1000000
#define maxint 0x7f7f7f7f
using namespace std;
int n,m,q,tot;
int fir[N],to[N],nex[N],dis[N];
int fa[N],f[N][25],deep[N],d[N][25];
bool vis[N];
struct node{
    int x,y,k;
}e[N];
inline void r(int &x){
    x=0;
    char ch=getchar();
    while(ch<'0'||ch>'9') ch=getchar();
    while(ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
}
inline void add(int x,int y,int k){//树上的无向边建立
    to[++tot]=y,nex[tot]=fir[x],dis[tot]=k,fir[x]=tot;
    to[++tot]=x,nex[tot]=fir[y],dis[tot]=k,fir[y]=tot;
}
inline bool cmp(node x,node y){//边权按大到小排序
    return x.k>y.k;
}
inline int anc(int x){//路径压缩
    if(fa[x]!=x) fa[x]=anc(fa[x]);
    return fa[x];
}
inline int dfs(int x){//搜索树,记录f[][0]
    vis[x]=1;
    int i,v;
    for(i=fir[x];v=to[i],i;i=nex[i]){
        if(!vis[v]){
            deep[v]=deep[x]+1;
            f[v][0]=x;
            d[v][0]=dis[i];
            dfs(v);
        }
    }
}
inline void cal(){//倍增,求出每个节点的祖先,以及每个节点到祖先的路径上最小边权
    for(int i=1;i<=19;i++)
        for(int j=1;j<=n;j++){
            f[j][i]=f[f[j][i-1]][i-1];
            d[j][i]=min(d[j][i-1],d[f[j][i-1]][i-1]);//最小边权用d记录
        }
}
inline int lca(int x,int y){//求答案
    if(f[x][19]!=f[y][19])//如果不等,那么他们不在一棵树上
        return -1;
    int i,ans=maxint;
    if(deep[x]<deep[y]) swap(x,y);//以下lca模板,在求lca过程顺便更新an最小值
    for(i=19;i>=0;i--){
        if(deep[f[x][i]]>=deep[y]){
            ans=min(ans,d[x][i]);
            x=f[x][i];
        }
        if(x==y) return ans;
    }
    for(i=19;i>=0;i--)
        if(f[x][i]!=f[y][i]){
            ans=min(ans,min(d[x][i],d[y][i]));
            x=f[x][i];
            y=f[y][i];
        }
    ans=min(ans,min(d[x][0],d[y][0]));
    return ans;
}
int main()
{
    int i,j,x,y,k,a,b;
    r(n),r(m);
    for(i=1;i<=m;i++)
        r(e[i].x),r(e[i].y),r(e[i].k);
    sort(e+1,e+1+m,cmp);
    for(i=1;i<=n;i++) fa[i]=i;
    for(i=1;i<=m;i++){//最大生成树模板
        a=anc(e[i].x),b=anc(e[i].y);
        if(a==b) continue;
        add(e[i].x,e[i].y,e[i].k);
        fa[a]=b;
    }
    for(i=1;i<=n;i++){
        if(!vis[i]){//因为可能有多棵树,因此要遍历每个点,看看它是否在已遍历的树上
            f[i][0]=i;
            d[i][0]=maxint;
            deep[i]=1;
            dfs(i);
        }
    }
    cal();
    r(q);
    while(q--){
        r(x),r(y);
        printf("%d\n",lca(x,y));
    }
}

 

转载于:https://www.cnblogs.com/quitter/p/11497017.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
06-01
这道是一道典型的费用限制最短路目,可以使用 Dijkstra 算法或者 SPFA 算法来解决。 具体思路如下: 1. 首先,我们需要读入输入数据。输入数据中包含了道路的数量、起点和终点,以及每条道路的起点、终点、长度和限制费用。 2. 接着,我们需要使用邻接表或邻接矩阵来存储图的信息。对于每条道路,我们可以将其起点和终点作为一个有向边的起点和终点,长度作为边权,限制费用作为边权的上界。 3. 然后,我们可以使用 Dijkstra 算法或 SPFA 算法求解从起点到终点的最短路径。在这个过程中,我们需要记录到每个点的最小费用和最小长度,以及更新每条边的最小费用和最小长度。 4. 最后,我们输出从起点到终点的最短路径长度即可。 需要注意的是,在使用 Dijkstra 算法或 SPFA 算法时,需要对每个点的最小费用和最小长度进行松弛操作。具体来说,当我们从一个点 u 经过一条边 (u,v) 到达另一个点 v 时,如果新的费用和长度比原来的小,则需要更新到达 v 的最小费用和最小长度,并将 v 加入到优先队列(Dijkstra 算法)或队列(SPFA 算法)中。 此外,还需要注意处理边权为 0 或负数的情况,以及处理无法到达终点的情况。 代码实现可以参考以下样例代码: ```c++ #include <cstdio> #include <cstring> #include <queue> #include <vector> using namespace std; const int MAXN = 1005, MAXM = 20005, INF = 0x3f3f3f3f; int n, m, s, t, cnt; int head[MAXN], dis[MAXN], vis[MAXN]; struct Edge { int v, w, c, nxt; } e[MAXM]; void addEdge(int u, int v, int w, int c) { e[++cnt].v = v, e[cnt].w = w, e[cnt].c = c, e[cnt].nxt = head[u], head[u] = cnt; } void dijkstra() { priority_queue<pair<int, int>, vector<pair<int, int>>, greater<pair<int, int>>> q; memset(dis, 0x3f, sizeof(dis)); memset(vis, 0, sizeof(vis)); dis[s] = 0; q.push(make_pair(0, s)); while (!q.empty()) { int u = q.top().second; q.pop(); if (vis[u]) continue; vis[u] = 1; for (int i = head[u]; i != -1; i = e[i].nxt) { int v = e[i].v, w = e[i].w, c = e[i].c; if (dis[u] + w < dis[v] && c >= dis[u] + w) { dis[v] = dis[u] + w; q.push(make_pair(dis[v], v)); } } } } int main() { memset(head, -1, sizeof(head)); scanf("%d %d %d %d", &n, &m, &s, &t); for (int i = 1; i <= m; i++) { int u, v, w, c; scanf("%d %d %d %d", &u, &v, &w, &c); addEdge(u, v, w, c); addEdge(v, u, w, c); } dijkstra(); if (dis[t] == INF) printf("-1\n"); else printf("%d\n", dis[t]); return 0; } ```

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值