洛谷——P1907 设计道路-CSDN博客

P1907 设计道路

题目描述

Caesar远征高卢回来后，对你大加赞赏，他亲自来到Genoa视察。

Genoa在你的建设下变得无比繁荣，由于财政收入的增加，你为城市修建了交通系统。古罗马的交通系统由两部分组成——Dirt Road和Rome Road。两个路口间只可能是其中一种道路。在Rome Road上可以驾驶马车，而在Dirt Road上则不行。由于修建道路是一项浩大的工程，使得你无法将整个城市用Rome Road连接起来。

现在Caesar已经到达码头，他要求去你家参观。Caesar由一个癖好，喜欢坐车而不喜欢走路。所以Caesar走Dirt Road时的不满值要比走Rome Road时大。

为了不让Caesar过于不满而罢免你的职位，请设计路线使得Caesar的不满值最小。

输入输出格式

输入格式：

输入数据第一行有两个实数，分别表示走Dirt Road和Rome Road一个单位长度时Caesar的不满值。接下来是一个整数N（N<=1000），代表路口总数。接下来有N行，每行一组实数(x,y)分别描述这N个路口的坐标。接下来有若干行，每行一组整数（i,j），表示第i个路口与第j个路口间为Rome Road，以0 0结束。最后两行，每行一对实数，分别描述码头和你家的坐标。

输出格式：

输出Caesar从码头到你家的最小不满值，保留4位小数

输入输出样例

输入样例#1：复制

输出样例#1：复制

202.8284


最短路
什么鬼畜的操作、、、

再跑这一步的时候必须要分开写，否则出错？！

++n;scanf("%lf%lf",&x[n],&y[n]);

注意数组大小！

#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define N 5010
#define maxn 9999999
using namespace std;
queue<int>q;
bool vis[N];
int n,m,t,xx,yy,tot,head[N*10];
double a,b,z,dis[N],x[N],y[N],f[N][N];
int read()
{
    int x=0,f=1; char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
struct Edge
{
    int to,next;
}edge[N*10];
int add(int x,int y)
{
    tot++;
    edge[tot].to=y;
    edge[tot].next=head[x];
    head[x]=tot;
}
int spfa(int s)
{
    for(int i=1;i<=n;i++) dis[i]=maxn,vis[i]=0;
    vis[s]=true,dis[s]=0,q.push(s);
    while(!q.empty())
    {
        xx=q.front();q.pop();vis[xx]=false;
        for(int i=head[xx];i;i=edge[i].next)
        {
            int to=edge[i].to;
            if(dis[to]>dis[xx]+f[xx][to])
            {
                dis[to]=dis[xx]+f[xx][to];
                if(!vis[to]) vis[to]=true,q.push(to);
            }
        }
    }
}
int main()
{
    scanf("%lf%lf",&a,&b);
    n=read();int m=2;
    for(int i=1;i<=n;i++) 
     scanf("%lf%lf",&x[i],&y[i]);
    while(1)
    {
        xx=read(),yy=read();
        add(xx,yy),add(yy,xx);
        if(xx==0&&yy==0) break;
        z=sqrt(pow(x[xx]-x[yy],2)+pow(y[xx]-y[yy],2));
        f[xx][yy]=f[yy][xx]=z*b;
    }
    for(int i=1;i<=n;i++)
     for(int j=i+1;j<=n;j++)
      if(!f[i][j])
      {
          add(i,j),add(j,i);
          z=sqrt(pow(x[i]-x[j],2)+pow(y[i]-y[j],2));
          f[i][j]=f[j][i]=z*a;
      }
    while(m--)
    {
        ++n;
        scanf("%lf%lf",&x[n],&y[n]);
        for(int i=1;i<n;i++)
        {
            add(n,i),add(i,n);
            z=sqrt(pow(x[i]-x[n],2)+pow(y[i]-y[n],2));
              f[n][i]=f[i][n]=z*a;
        } 
    }
    spfa(n-1);
    printf("%.4lf",dis[n]);
    return 0;
}