CCF认证 2017-12 行车路线

最短路的变形,有两种道路,需要分开计算最短距离,要不然得不到全局最优解

考试得分80,后来看了好久发现,竟然是int爆掉了,导致没有满分(囧)

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
const int N=500+10;
const ll inf=1e18; //需要开long long,虽然结果不大于1e6,但中间值依然可能爆int,就因为这个80分,天真了 
bool vis[N];
int que[N];
ll dist[N],dist0[N]; //前驱为大路和前驱为小路分别计算距离 
ll g[N][N],g0[N][N]; //大路邻接矩阵和小路邻接矩阵 
void SPFA(int start,int n)
{
    int front=0,rear=0;
    for(int v=1;v<=n;v++)
    {
        if(v==start)
        {
            que[rear++]=v;
            vis[v]=true;
            dist[v]=dist0[v]=0;
        }
        else
        {
            vis[v]=false;
            dist[v]=dist0[v]=inf;
        }
    }
    while(front!=rear)
    {
        int u=que[front++];
        vis[u]=false;
        if(front>=N)front=0; 
        for(int i=1;i<=n;i++)
        {
            ll v=g[u][i];
            if(dist[i]>dist[u]+v) //走大路 
            {
                dist[i]=dist[u]+v;
                if(!vis[i])
                {
                    vis[i]=true;
                    que[rear++]=i; 
                    if(rear>=N)rear=0;
                }
            }
            if(dist[i]>dist0[u]+v) //走大路 
            {
                dist[i]=dist0[u]+v;
                if(!vis[i])
                {
                    vis[i]=true;
                    que[rear++]=i; 
                    if(rear>=N)rear=0;
                }
            }
			if(g0[u][i]!=inf) //如果小路能到达 
			{
				v=g0[u][i]*g0[u][i];
				if(dist0[i]>dist[u]+v) //走小路 ,前驱为小路时不能走小路 
            	{
               		dist0[i]=dist[u]+v;
               	 	if(!vis[i])
                	{
                  	 	vis[i]=true;
                  		que[rear++]=i; 
                   	 	if(rear>=N)rear=0;
               	 	}
           	 	}
			}
        }
    }
}
int main()
{
	int n,m;
	scanf("%d%d",&n,&m);
	ll t,a,b,c;
	for(int i=0;i<=n;i++)
	for(int j=0;j<=n;j++)
	g[i][j]=g0[i][j]=inf;
	for(int i=0;i<m;i++)
	{
		scanf("%lld%lld%lld%lld",&t,&a,&b,&c);
		if(t&&g0[a][b]>c) g0[a][b]=g0[b][a]=c;
		if(!t&&g[a][b]>c) g[a][b]=g[b][a]=c;
	}
	for(int i=1;i<=n;i++) //这里用floyd事先计算好只走小路时两两点之间的最短距离
	for(int j=i+1;j<=n;j++)
	{
		for(int k=1;k<=n;k++)
		{
			if(k==i||k==j) continue;
			if(g0[i][j]>g0[i][k]+g0[k][j]) g0[i][j]=g0[j][i]=g0[i][k]+g0[k][j];
		}
	}
	SPFA(1,n); //spfa 
	printf("%lld\n",min(dist[n],dist0[n]));
	return 0;
} 


  • 8
    点赞
  • 18
    收藏
    觉得还不错? 一键收藏
  • 15
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值