挖隧道 最短路

NKOJ3543，题面

首先容易反应出最短路算法。

如果直接乱连边，那么线段上每一个点都要讨论是否能够挖隧道。由于坐标都是实数，那么几乎有无数个点需要讨论，这显然是不行的。因此，我们需要得到一些结论，很容易猜想只有折点才可以作为隧道的入口或出口。

接下来想办法证明一下。首先考虑隧道只会穿过一个“山峰”的情形，如下图所示。

考虑如果不在折点A打隧道，比如图中的黄色部分，由于无论如何都会经过B点，那么可以通过数学的推导（如相似三角形）得出从A到B的时间与隧道长度的关系，发现是一个一次函数的形式。根据一次函数的性质，最值只会在端点取到，即要么在折点打隧道，要么就在这一段路都不打隧道。

接下来就考虑比较复杂的情形，如果隧道要穿过多个“山峰”又怎么办呢？

其实结论是一样的。

跨过多个山峰，无非有两种方法：

1.如果能直接打一条隧道穿过所有，打一条直接穿过所有山峰的隧道；
2.对于每个山峰，讨论上面“只会超过一个山峰”的两种情况。

那么只需要讨论方法1的实现。

因为是在绿色虚线（表示中间的山峰中最低山谷所在的水平线）以下打隧道，而实际上绿色虚线以上部分是每次都一样的，考虑的只是两侧的部分。对黑色虚线补成的三角形运用相似三角形推导，又可以得到一个一次函数关系。即在情况2下，要么打绿色虚线的隧道，要么打蓝色隧道。而绿色虚线的隧道会在之后的加边中体现为两个或多个（“多个”是因为可能最低山谷有多个）隧道连接的形式，所以实现起来还是比较简单的。

接下来算法就明确了：

对于每一个折点，讨论它是否可以作为一个隧道的起点或终点，如果可以就打一条隧道，最后用最短路算法得出答案。

打隧道可能实现起来觉得比较恶心，由于打隧道与不打隧道存在重复的道路，每次打隧道只需要连两条边，具体实现参照代码。

#include<stdio.h>
#include<queue>
#include<cmath>
#define MAXN 4005
#define MAXM 16005
#define dd double
using namespace std;

int N,tot;
dd W,C,X[MAXN],Y[MAXN];

dd Dis(int x,int y){return sqrt((X[x]-X[y])*(X[x]-X[y])+(Y[x]-Y[y])*(Y[x]-Y[y]));}

int Tot,nex[MAXM],las[MAXN],en[MAXM];
dd len[MAXM];
void ADD(int x,int y,dd z)
{
    en[++Tot]=y;
    nex[Tot]=las[x];
    las[x]=Tot;
    len[Tot]=z;
}

bool mark[MAXN];
dd dis[MAXN];
void SPFA()
{
    int i,x,y;
    queue<int>Q;

    for(i=1;i<=tot;i++)dis[i]=1e9;
    dis[1]=0;
    Q.push(1);
    while(Q.size())
    {
        x=Q.front();Q.pop();mark[x]=false;
        for(i=las[x];i;i=nex[i])
        {
            y=en[i];
            if(dis[y]>dis[x]+len[i])
            {
                dis[y]=dis[x]+len[i];
                if(!mark[y])mark[y]=true,Q.push(y);
            }
        }
    }
}

int main()
{
    int i,j;
    dd tmp;

    scanf("%d%lf%lf",&N,&W,&C);
    for(i=1;i<=N;i++)scanf("%lf%lf",&X[i],&Y[i]);
    tot=N;

    for(i=1;i<N;i++)
    {
        for(j=i+1;j<=N&&Y[j]>Y[i];j++);
        if(j==N+1||j==i+1)continue; //越出边界或不存在山峰，特判
        tot++;
        tmp=X[j-1]-X[i]+(X[j]-X[j-1])*(Y[j-1]-Y[i])/(Y[j-1]-Y[j]);//相似三角形
        ADD(i,tot,tmp/C);//加一条隧道
        tmp=Dis(j-1,j)*((Y[j-1]-Y[j])-(Y[j-1]-Y[i]))/(Y[j-1]-Y[j]);//相似三角形
        ADD(tot,j,tmp/W);//隧道出口与下一个折点连边
    }
    //讨论是否可以作为起点，下面的情况同理
    for(i=N;i>1;i--)
    {
        for(j=i-1;j&&Y[j]>Y[i];j--);
        if(j==1||j==i-1)continue;
        tot++;
        tmp=X[i]-X[j+1]+(X[j+1]-X[j])*(Y[j+1]-Y[i])/(Y[j+1]-Y[j]);
        ADD(tot,i,tmp/C);
        tmp=Dis(j+1,j)*((Y[j+1]-Y[j])-(Y[j+1]-Y[i]))/(Y[j+1]-Y[j]);
        ADD(j,tot,tmp/W);
    }
    //讨论是否可以作为终点

    for(i=1;i<N;i++)ADD(i,i+1,Dis(i,i+1)/W),ADD(i+1,i,Dis(i,i+1)/W);
    SPFA();
    printf("%.6lf",dis[N]);
}