HDU 1428 漫步校园（spfa最短路+记忆化搜索）

题目

LL最近沉迷于AC不能自拔，每天寝室、机房两点一线。由于长时间坐在电脑边，缺乏运动。他决定充分利用每次从寝室到机房的时间，在校园里散散步。整个HDU校园呈方形布局，可划分为n*n个小方格，代表各个区域。例如LL居住的18号宿舍位于校园的西北角，即方格(1,1)代表的地方，而机房所在的第三实验楼处于东南端的(n,n)。因有多条路线可以选择，LL希望每次的散步路线都不一样。另外，他考虑从A区域到B区域仅当存在一条从B到机房的路线比任何一条从A到机房的路线更近(否则可能永远都到不了机房了…)。现在他想知道的是，所有满足要求的路线一共有多少条。你能告诉他吗?

Input
每组测试数据的第一行为n(2<=n<=50)，接下来的n行每行有n个数，代表经过每个区域所花的时间t(1<=t<=50)(由于寝室与机房均在三楼，故起点与终点也得费时)。

Output
针对每组测试数据，输出总的路线数(小于2^63)。

Sample Input
3
1 2 3
1 2 3
1 2 3
3
1 1 1
1 1 1
1 1 1

Sample Output
1
6

解题

问你从（1，1）到（n，n）有多少条不同的路径。
题目有个比较含蓄的说法“他考虑从A区域到B区域仅当存在一条从B到机房的路线比任何一条从A到机房的路线更近(否则可能永远都到不了机房了…)”。
这个意思就是你只能越走离（n，n）越近，中途不能走远方去。
先求出以（n，n）为源点的单源最短路d[i][j]。
然后问题就转化成了在一个二维的格子里，从（1，1）到（n，n）的走法数，要求只能往上下左右四个方向走且下一步到达的方格里的值d[fi][fj]要比当前的d值小。
这样，记忆化搜一遍就好了。

AC代码

#include <cstdio>
#include <algorithm>
#include <iostream>
#include <vector>
#include <cstring>
#include <queue>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int maxn=55;
int n;
int g[maxn][maxn];
int d[maxn][maxn];//d[i][j]表示机房到（i，j）的最短路
bool inque[maxn][maxn];
int dir[4][2]={1,0,-1,0,0,1,0,-1};
struct node
{
    int x,y,w;
    node(int x,int y,int w):x(x),y(y),w(w){}
    node(){}
};

bool valid(int x,int y)
{
    return x>=1&&x<=n&&y>=1&&y<=n;
}
void spfa()
{
    memset(d,0x3f,sizeof(d));
    memset(inque,false,sizeof(inque));
    queue<node> Q;
    d[n][n]=g[n][n];
    Q.push(node(n,n,d[n][n]));
    inque[n][n]=true;
    while(!Q.empty())
    {
        node now=Q.front();Q.pop();
        int x=now.x,y=now.y,w=now.w;
        inque[x][y]=false;
        for(int i=0;i<4;i++)
        {
            int fx=x+dir[i][0],fy=y+dir[i][1];
            if(!valid(fx,fy)) continue;
            if(d[fx][fy]>d[x][y]+g[fx][fy])
            {
                d[fx][fy]=d[x][y]+g[fx][fy];
                if(inque[fx][fy]) continue;
                Q.push(node(fx,fy,d[fx][fy]));
            }
        }
    }
}
ll sum[maxn][maxn];
ll dp(int x,int y)
{
    if(sum[x][y]) return sum[x][y];
    ll &ans=sum[x][y];
    for(int i=0;i<4;i++)
    {
        int fx=x+dir[i][0],fy=y+dir[i][1];
        if(!valid(fx,fy)) continue;
        if(d[fx][fy]>=d[x][y]) continue;
        ans+=dp(fx,fy);
    }
    return ans;
}
int main()
{
    while(~scanf("%d",&n))
    {
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
            scanf("%d",&g[i][j]);
        spfa();
        memset(sum,0,sizeof(sum));
        sum[n][n]=1;
        ll ans=dp(1,1);//记忆化搜索
        printf("%lld\n",ans);
    }
    return 0;
}