CCF认证 2017-12 行车路线
2018-02-05 22:19
288 查看
最短路的变形,有两种道路,需要分开计算最短距离,要不然得不到全局最优解
考试得分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
;
int que
;
ll dist
,dist0
; //前驱为大路和前驱为小路分别计算距离
ll g
,g0
; //大路邻接矩阵和小路邻接矩阵
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
,dist0
));
return 0;
}
考试得分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
;
int que
;
ll dist
,dist0
; //前驱为大路和前驱为小路分别计算距离
ll g
,g0
; //大路邻接矩阵和小路邻接矩阵
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
,dist0
));
return 0;
}
相关文章推荐
- CCF认证 2017-12 行车路线
- CCF认证 2017-12 行车路线
- CSP 2017_12_4行车路线
- CCF认证 2017-12 Crontab
- CCF认证 行车路线
- CCF认证 2017-12 最小差值
- CCF认证 2017-12 游戏
- CCF认证 2017-12 商路
- CCF认证 2015-12 画图
- C++基础 多态与虚函数 2017_12_03
- 12/04/2017 - 12/10/2017
- (2017-12)最新12306爬虫
- 白盒测试实践--Day3 12/19/2017
- C语言趣味一百道 第19题 2017_12_21
- 安徽科技学院2017-2018-1学期15电信12《Java编程技术》期末考试
- JZOJ5453. 【NOIP2017提高A组冲刺11.5】好路线
- CCF认证 2013-12 有趣的数
- 2017_12_05学习Java心得笔记
- fundamental_2017_12_9
- 2017-12-15 1周5次课 Linux目录结构、ls命令、文件类型、别名