单源最短路的综合应用

最新推荐文章于 2024-10-09 22:28:58 发布

没头发的年轻人

最新推荐文章于 2024-10-09 22:28:58 发布

阅读量218

点赞数

文章标签：算法

本文链接：https://blog.csdn.net/lee_14141/article/details/126911099

版权

1.信息学奥赛一本通（C++版）在线评测系统

（1）先预处理出从1 a b c d e出发到其他所有点的最短路

（2）DFS所有拜访顺序 5! 对于每种拜访顺序可以通过查表的方式算出最短距离

O（ k*6*m + 5！）

#include <bits/stdc++.h>
using namespace std;
const int N=50010,M=200010,INF=0x3f3f3f3f;
int n,m;
int source[6];
int h[N];
int e[M],w[M],ne[M],idx;
int q[N],dist[6][N];
bool st[N];
void add(int a,int b,int c)
{
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

void spfa(int start,int dist[])
{
    memset(dist,0x3f,4*N);
    dist[start]=0;
    int hh=0,tt=1;
    q[0]=start;
    while(hh!=tt)
    {
        int t=q[hh++];
        if(hh==N) hh=0;
        st[t]=false;
        for(int i=h[t];~i;i=ne[i])
        {
            int j=e[i];
            if(dist[j]>dist[t]+w[i])
            {
                dist[j]=dist[t]+w[i];
                if(!st[j])
                {
                    q[tt++]=j;
                    if(tt==N) tt=0;
                    st[j]=true;
                }
            }
        }
    }
}
int dfs(int u,int start,int distance)
{
    if(u==6) return distance;
    int res=INF;
    for(int i=1;i<=5;i++)
    {
        if(!st[i])
        {
            int next=source[i];
            st[i]=true;
            res=min(res,dfs(u+1,i,distance+dist[start][next]));
            st[i]=false;
        }
    }
    return res;
}
int main()
{
    scanf("%d%d",&n,&m);
    source[0]=1;
    for(int i=1;i<=5;i++)
        scanf("%d",&source[i]);
    memset(h,-1,sizeof h);
    while(m--)
    {
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
        add(b,a,c);
    }
    for(int i=0;i<=5;i++)
    {
        spfa(source[i],dist[i]);
    }
    printf("%d\n",dfs(1,0,0));
    return 0;
}

2.340. 通信线路 - AcWing题库

二分+最短路

#include <bits/stdc++.h>
using namespace std;
const int N=1010,M=20010;
int n,m,k;
int h[N],e[M],ne[M],w[M],idx;
deque<int> q;
int dis[N];
bool st[N];
void add(int a,int b,int c)
{
    e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}
bool check(int bound)
{
    memset(st,false,sizeof st);
    memset(dis,0x3f,sizeof dis);
    dis[1]=0;
    q.push_back(1);
     while(q.size())
     {
         int t=q.front();
         q.pop_front();
         if(st[t]) continue;
         st[t]=true;
         for(int i=h[t];~i;i=ne[i])
         {
             int j=e[i],v=w[i]>bound;
             if(dis[j]>dis[t]+v)
             {
                 dis[j]=dis[t]+v;
                 if(!v)
                 {
                     q.push_front(j);
                 }
                 else q.push_back(j);
             }
         }
     }
     return dis[n]<=k;
}
int main()
{
    cin>>n>>m>>k;
    memset(h,-1,sizeof h);
    while(m--)
    {
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
        add(b,a,c);

    }
    int l=0,r=1e6+1;
    while(l<r)
    {
        int mid=l+r>>1;
        if(check(mid)) r=mid;
        else l=mid+1;
    }
    if(r==1e6+1) r=-1;
    cout<<r<<endl;
    return 0;
}

3.342. 道路与航线 - AcWing题库

#include <bits/stdc++.h>
#define x first
#define y second
using namespace std;
typedef pair<int,int> PII;
const int N=25010,M=150010,INF=0x3f3f3f3f;
int n,mr,mp,s;
int h[N],e[M],ne[M],w[M],idx;
int id[N];
bool st[N];
vector<int> block[N];
int dis[N],din[N];
int bcnt;
queue<int> q;
void add(int a,int b,int c)
{
    e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}
void dijkstra(int bid)
{
    priority_queue<PII,vector<PII>,greater<PII>> heap;
    for(auto ver:block[bid]) heap.push({dis[ver],ver});
    while(heap.size())
    {
        auto t=heap.top();
        heap.pop();
        int ver=t.y,distance=t.x;
        if(st[ver]) continue;
        st[ver]=true;
        for(int i=h[ver];~i;i=ne[i])
        {
            int j=e[i];
            if(dis[j]>dis[ver]+w[i])
            {
                dis[j]=dis[ver]+w[i];
                if(id[j]==id[ver])
                    heap.push({dis[j],j});
            }
            if(id[j]!=id[ver]&&--din[id[j]]==0)
            {
                q.push(id[j]);
            }
        }
    }
}
void dfs(int u,int bid)
{
    id[u]=bid;
    block[bid].push_back(u);
    for(int i=h[u];~i;i=ne[i])
    {
        int j=e[i];
        if(!id[j])
        {
            dfs(j,bid);
        }
    }
}
void topsort()
{
    memset(dis,0x3f,sizeof dis);
    dis[s]=0;
    for(int i=1;i<=bcnt;i++)
    {
        if(!din[i]) q.push(i);
    }
    while(q.size())
    {
        int t=q.front();
        q.pop();
        dijkstra(t);
    }
}
int main()
{
    scanf("%d%d%d%d",&n,&mr,&mp,&s);
    memset(h,-1,sizeof h);
    while(mr--)
    {
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c),add(b,a,c);
    }
    for(int i=1;i<=n;i++)
    {
        if(!id[i])
        {
            dfs(i,++bcnt);
        }
    }
    while(mp--)
    {
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
        din[id[b]]++;
    }
    topsort();
    for(int i=1;i<=n;i++)
    {
        if(dis[i]>INF/2) puts("NO PATH");
        else printf("%d\n",dis[i]);
    }
    return 0;
}

4.341. 最优贸易 - AcWing题库

dp+最短路

为什么要建反向图

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10,M=1e6+10;
int hs[N],ne[M],e[M],idx,w[M],ht[N];
int n,m;
int dmin[N],dmax[N];
int q[N];
bool st[N];
void add(int h[],int a,int b)
{
    e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void spfa(int h[],int dist[],int type)
{
    int hh=0,tt=1;
    if(type==0)
    {
        memset(dist,0x3f,sizeof dmin);
        dist[1]=w[1];
        q[0]=1;
    }
    else
    {
        memset(dist,-0x3f,sizeof dmax);
        dist[n]=w[n];
        q[0]=n;
    }
    while(hh!=tt)
    {
        int t=q[hh++];
        if(hh==N) hh=0;
        st[t]=false;
        for(int i=h[t];~i;i=ne[i])
        {
            int j=e[i];
            if((type==0&&dist[j]>min(dist[t],w[j]))||(type==1&&dist[j]<max(dist[t],w[j])))
            {
                if(type==0) dist[j]=min(dist[t],w[j]);
                else dist[j]=max(dist[t],w[j]);
                if(!st[j])
                {
                    q[tt++]=j;
                    if(tt==N) tt=0;
                    st[j]=true;
                }
            }
        }
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        scanf("%d",&w[i]);
    memset(hs,-1,sizeof hs);
    memset(ht,-1,sizeof ht);
    while(m--)
    {
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        add(hs,a,b);
        add(ht,b,a);
        if(c==2)
        {
            add(hs,b,a),add(ht,a,b);
        }
    }
    spfa(hs,dmin,0);
    spfa(ht,dmax,1);
    int res=0;
    for(int i=1;i<=n;i++)
        res=max(res,dmax[i]-dmin[i]);
    printf("%d\n",res);
    return 0;
}