[bzoj 2709]迷宫花园

普通得甚至有些二逼的矮穷挫少年——Dios，不可避免地遇到了他生命中的劫数，白富美少女 Nyution。但是按照正常的校园故事的发展，Nyution 是无论如何不会喜欢上各方面条件都差到不行的 Dios 的。不过，Dios 还是面对 Nyution 颤抖着说出了那三个字。Nyution既不想过分地让 Dios 伤心，又不想接受她根本看不上的 Dios，于是决定让 Dios 走一个建在她家后院里的迷宫花园——如果 Dios 能很快地从起点走到终点，证明他的聪明才智，Nyution就答应他的表白。
当然 Nyution 敢这么说肯定是有准备的。Nyution 的花园可以看做一个迷宫，在迷宫内部有起点和终点。Dios 要从起点走到终点，并且他只能选择前后左右四个方向行走，而且显然不能走到篱笆上，也不能走出迷宫的边界。Nyution 经过仔细的调查发现，Dios 移动到相邻格子的耗时肯定是 1。同时，Nyution 将在 Dios 的挑战开始前，通过进行适当的路面调 整，使 Dios 在南北方向（数据中的上下方向）的移动时间由 1 变成实数v 。首先，Nyution 不能让 Dios 过快地到达终点，这样她就得接受表白；其次，Nyution 也不想让 Dios 开了小宇宙之后还是过慢地到达终点，这样显得她在刁难 Dios。最后她确定了一个实数 L ——就 是最坏情况（也就是 Dios 最神勇威武耗时最短的情况）下，Dios 将花费 L 的时间由起点到达终点。但是 Nyution 显然不会求此时的v 值，于是她找到了一向以算法达人著称的你。你当然不会拒绝白富美 Nyution 的请求，决定帮她算出此时的v 。由于 Nyution 不仅是白富美同时也是三好学生，所以她肯定不会给你一个无解的任务。并且，Nyution 的迷宫中一定没有水平的从起点到终点的通路。

这道题看到它给出了答案范围，又发现它满足二分性，那我们就往二分那里想，那这道题转换为一道spfa大水题，注意精度问题，用1e-7。（然而身为蒟蒻的我认为这题难点实际上是输入）

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;
double L;
int R,C,stx,sty,edx,edy;
struct node
{
    int x,y;
}list[11000];
int dx[5]={1,0,-1,0};
int dy[5]={0,1,0,-1};
double ss[110][110];
char s[110][110];
bool v[110][110];
bool check(double k)
{
    memset(v,true,sizeof(v));
    for(int i=1;i<=100;i++)
    {
        for(int j=1;j<=100;j++)ss[i][j]=999999999.0;
    }   
    int head=1,tail=2;double ans=999999999.0;
    list[1].x=stx,list[1].y=sty;
    v[stx][sty]=false;
    ss[stx][sty]=0.0;
    while(head!=tail)
    {
        node tno=list[head];
        for(int i=0;i<4;i++)
        {
            tno.x=list[head].x+dx[i],tno.y=list[head].y+dy[i];
            if(tno.x>=1 && tno.x<=R && tno.y>=1 && tno.y<=C && s[tno.x][tno.y]!='#')
            {
                double he=1.0;
                if(i==0 || i==2)he=k;
                if(ss[tno.x][tno.y]>ss[list[head].x][list[head].y]+he)
                {
                    ss[tno.x][tno.y]=ss[list[head].x][list[head].y]+he;
                    if(v[tno.x][tno.y]==true)
                    {
                        v[tno.x][tno.y]=false;
                        list[tail++]=tno;
                        if(tail==10001)tail=1;
                        if(tno.x==edx && tno.y==edy)
                        {
                            if(ans>ss[tno.x][tno.y])ans=ss[tno.x][tno.y];
                        }
                    }
                }
            }
        }
        v[list[head].x][list[head].y]=true;
        head++;
        if(head==10001)head=1;
    }
    if(ans>=L)return true;
    else return false;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%lf%d%d",&L,&R,&C);
        for(int i=1;i<=R;i++)
        {
            char c;
            scanf("%c",&c);
            while (c!='\n') scanf("%c",&c);  
            for(int j=1;j<=C;j++)
            {
                scanf("%c",&c);  
                while (c=='\n') scanf("%c",&c);                 
                s[i][j]=c;
                if(s[i][j]=='S')stx=i,sty=j;
                if(s[i][j]=='E')edx=i,edy=j;
            }
        }
        double l=0.0,r=10.0,mid,ans;
        while(r-l>=1e-7)
        {
            mid=(l+r)/2.0;
            if(check(mid)==true)
            {
                ans=mid;r=mid;
            }
            else l=mid;
        }
        printf("%.5lf\n",ans);
    }
    return 0;
}