Codeforces Round #355 (Div. 2) D. Vanya and Treasure 分治暴力

D. Vanya and Treasure

Vanya is in the palace that can be represented as a grid n × m. Each room contains a single chest, an the room located in the i-th row and j-th columns contains the chest of type aij. Each chest of type x ≤ p - 1 contains a key that can open any chest of type x + 1, and all chests of type 1 are not locked. There is exactly one chest of type p and it contains a treasure.

Vanya starts in cell (1, 1) (top left corner). What is the minimum total distance Vanya has to walk in order to get the treasure? Consider the distance between cell (r1, c1) (the cell in the row r1 and column c1) and (r2, c2) is equal to |r1 - r2| + |c1 - c2|.

Input

The first line of the input contains three integers n, m and p (1 ≤ n, m ≤ 300, 1 ≤ p ≤ n·m) — the number of rows and columns in the table representing the palace and the number of different types of the chests, respectively.

Each of the following n lines contains m integers aij (1 ≤ aij ≤ p) — the types of the chests in corresponding rooms. It's guaranteed that for each x from 1 to p there is at least one chest of this type (that is, there exists a pair of r and c, such that arc = x). Also, it's guaranteed that there is exactly one chest of type p.

Output

Print one integer — the minimum possible total distance Vanya has to walk in order to get the treasure from the chest of type p.

Examples

input

3 4 3
2 1 1 1
1 1 1 1
2 1 1 3

output

5

input

3 3 9
1 3 5
8 9 7
4 6 2

output

22

input

3 4 12
1 2 3 4
8 7 6 5
9 10 11 12

output

11

题意:给你一个n*m的方格，你需要把权值为p的那个门打开，前提是你必须打开过权值为p-1的门。然后问你打开权值为p的门，你最少走的距离是多少，一开始你在(1,1)注意：一个点的门即使没有被打开，也是可以走的，只是门没开而已。
我们对于这个题进行优雅的暴力.
假设我们现在走的点数为i,那么点数为i的格子的个数如果小于sqrt(n*m)我们直接暴力更新
如果大于sqrt(n*m)我们就跑一遍bfs,这样做均摊了复杂度n^3
代码如下:

#include <bits/stdc++.h>

using namespace std;
const int maxn = 330;
int a[maxn][maxn];
int dp[maxn][maxn];
int inq[maxn][maxn];
int dis[maxn][maxn];
int dx[4]={1,-1,0,0};
int dy[4]={0,0,1,-1};
vector < pair<int,int> > vec[maxn*maxn];
int n,m,p;
int main()
{
    //freopen("de.txt","r",stdin);
    while (~scanf("%d%d%d",&n,&m,&p)){
        for (int i=0;i<maxn*maxn;++i)
        vec[i].clear();
        memset(inq,0,sizeof inq);
        memset(dp,127,sizeof dp);
        for (int i=1;i<=n;++i){
            for (int j=1;j<=m;++j){
                scanf("%d",&a[i][j]);
                if (a[i][j]==1) dp[i][j] = i+j-2;
                vec[a[i][j]].push_back(make_pair(i,j));
            }
        }
        int level = (int) sqrt(n*m);
        int idx = 1;//idx这样处理可以避免memset浪费时间
        queue <pair<int,int> > q;
        for (int i=1;i<p;++i){
            int sz = vec[i].size();
            if (sz<=level){
                for (auto u:vec[i]){
                    for (auto v:vec[i+1]){
                        dp[v.first][v.second] = min(dp[v.first][v.second],dp[u.first][u.second]+abs(u.first-v.first)+abs(u.second-v.second));
                    }
                }
            }
            else{
                idx++;
                for (auto v:vec[i])
                    q.push(v);
                memset(dis,127,sizeof dis);
                for (auto v:vec[i])
                    dis[v.first][v.second] = dp[v.first][v.second];
                while (!q.empty()){
                    pair<int,int > now = q.front();
                    q.pop();
                    inq[now.first][now.second]=0;
                    for (int k=0;k<4;++k){
                        pair<int,int> nxt = now;
                        nxt.first+=dx[k];
                        nxt.second+=dy[k];
                        if (nxt.first<1||nxt.first>n) continue;
                        if (nxt.second<1||nxt.second>m) continue;
                        if (dis[nxt.first][nxt.second]>dis[now.first][now.second]+1){
                            dis[nxt.first][nxt.second]=dis[now.first][now.second]+1;
                            if (inq[nxt.first][nxt.second]<idx){
                                inq[nxt.first][nxt.second]=idx;
                                q.push(nxt);
                            }
                        }
                    }
                }
                for (auto v:vec[i+1])
                    dp[v.first][v.second]=dis[v.first][v.second];
            }
        }
        printf("%d
",dp[vec[p][0].first][vec[p][0].second]);
    }
    return 0;
}