本题思路参考自大佬https://blog.csdn.net/stl112514/article/details/51563993
我们可以发现,对于每个当前点到下一个点的最优路径一定是一段横(或没有)+一段竖(或没有),因为每个转折点一定位于上一步的同一行和下一步的同一列,所以我们可以枚举转折点,利用每个已经更新过的点来更新其所属横的每个点作为前往下一个点的转折点,利用vis将其标记用对应钥匙,并更新其上一个出发点到这个转折点的最小贡献,对于每个下一步,我们可以枚举同一列的点作为转折点,通过标记判断是否合法(即可否为转折点),若合法则更新转折点即可,总复杂度O(n^3)
细节见代码注释
#include<iostream>
#include<stack>
#include<list>
#include<set>
#include<vector>
#include<algorithm>
#include<math.h>
#include<numeric>
#include<map>
#include<cstring>
#include<queue>
#include<iomanip>
#include<cmath>
#include<queue>
#include <bitset>
#include<unordered_map>
#ifndef local
#define endl '\n'
#endif */
#define mkp make_pair
using namespace std;
using std::bitset;
typedef long long ll;
typedef long double ld;
const int inf=0x3f3f3f3f;
const ll MAXN=2e6+10;
const ll N=1e5+100;
const ll mod=1e9+7;
const ll hash_p1=1610612741;
const ll hash_p2=805306457;
const ll hash_p3=402653189;
//-----------------------------------------------------------------------------------------------------------------*/
// ll head[MAXN],net[MAXN],to[MAXN],edge[MAXN]/*流量*/,cost[MAXN]//费用;
/*
void add(ll u,ll v,ll w,ll s){
to[++cnt]=v;net[cnt]=head[u];edge[cnt]=w;cost[cnt]=s;head[u]=cnt;
to[++cnt]=u;net[cnt]=head[v];edge[cnt]=0;cost[cnt]=-s;head[v]=cnt;
}
struct elemt{
int p,v;
};
struct comp{
public:
bool operator()(elemt v1,elemt v2){
return v1.v<v2.v;
}
};
-----------------------------------
求[1,MAXN]组合式和逆元
ll mi(ll a,ll b){
ll res=1;
while(b){
if(b%2){
res=res*a%mod;
}
a=a*a%mod;
}
return res;
}
ll fac[MAXN],inv[MAXN]
fac[0]=1;inv[0]=1;
for(int i=1;i<=MAXN;i){
fac[i]=(fac[i-1]*i)%mod;
inv[i]=mi(fac[i],mod-2);
}
ll C(int m,int n){//组合式C(m,n);
if(!n){
return 1;
}
return fac[m]*(inv[n]*inv[m*-n]%mod)%mod;
}
---------------------------------
unordered_map<int,int>mp;
//优先队列默认小顶堆 , greater<int> --小顶堆 less<int> --大顶堆
priority_queue<elemt,vector<elemt>,comp>q;
set<int>::iterator it=st.begin();
*/
// vector<vector<int>>edge; 二维虚拟储存坐标
//-----------------------------------------------------------------------------------------------------------------*/
//map<int,bool>mp[N];
vector<pair<int,int>>a[N];
int ans[310][310];//打开当前点宝箱所需总步数
int dis[310][310];//从上一点到该转折点所需最小步数
int vis[310][310];//转折点对应的可达钥匙
int main(){
/*cout<<setiosflags(ios::fixed)<<setprecision(8)<<ans<<endl;//输出ans(float)格式控制为8位小数(不含整数部分)*/
/*cout<<setprecision(8)<<ans<<endl;//输出ans(float)格式控制为8位小数(含整数部分)*/
ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);//同步流
int n,m,p;
memset(ans,inf,sizeof(ans));
cin>>n>>m>>p;
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
int d;
cin>>d;
a[d].push_back({i,j});
}
}
a[0].push_back({1,1});//从[1][1]出发,但可能会导致后面使得ans[1][1]被更新为0,和实际不符,需要特判
cout<<ans[1][1]<<endl;
for(int i=0;i<=p;i++){//枚举要打开的宝箱
for(int j=0;j<a[i].size();j++){//更新当前ans
int x=a[i][j].first,y=a[i][j].second;
if(x==1&&y==1&&ans[x][y]==0){//特判
ans[x][y]=inf;
}
for(int k=1;k<=n;k++){
if(vis[k][y]==i){//枚举转折点来更新答案
ans[x][y]=min(ans[x][y],dis[k][y]+abs(k-x));//
}
}
}
for(int j=0;j<a[i].size();j++){//更新下一次的dis
int x=a[i][j].first,y=a[i][j].second;
for(int k=1;k<=m;k++){//更新到达转折点所需最小步数
if(vis[x][k]!=i+1){
vis[x][k]=i+1;
dis[x][k]=ans[x][y]+abs(k-y);
}
else{
dis[x][k]=min(dis[x][k],ans[x][y]+abs(k-y));
}
}
}
}
cout<<ans[a[p][0].first][a[p][0].second]<<endl;;//输出答案
return 0;
}