Codeforces Round #531 (Div. 3) F.Elongated Matrix(状压dp)

题意

n<=16,m<=1e4

给你一个n*m的矩阵，要求你交换其中一些行，行只能完整地被交换

使得S型按列遍历，

即第一行第一个到第二行第一个到最后一行第一个

再到第一行第二个……

这样形成一个遍历序列后

使得相邻元素之差的绝对值最大

思路来源

https://blog.csdn.net/weixin_38686780/article/details/86192979

题解

定义行与行之间的距离为这两行对应元素之间距离的最小值

因为相邻行一旦固定，肯定要在这个序列中出现每一个对应位置的距离访问

那么就把行缩成了一个点

而从第一行到最后一行的安排，

若令第一个访问的是第一行，最后一个访问的是最后一行

就是我们怎么安排这个访问顺序，使得这个值最大

那就变成了安排顺序的哈密顿开路径1->2->…->n

枚举第一行，状压dp到所有都访问完，再枚举最后一行

因为还要考虑最后一行到第一行的下一个元素的转移，使答案变小

复杂度，

但考虑到不是严格的0到(1<<n-1)，而且每次只有合法位置才会转移，

常数比较小吧，再加上4s的时限，就没什么阻力了QAQ

代码

#include <iostream>
#include <algorithm> 
#include <cstring>
#include <cstdio>
#include <cmath>
#include <set>
#include <map>
#include <vector>
#include <stack>
#include <queue>
#include <functional>
const int INF=0x3f3f3f3f;
const int maxn=17; 
const int maxm=1e4+10;
const int mod=1e9+7;
const int MOD=998244353;
const double eps=1e-7;
typedef long long ll;
#define vi vector<int> 
#define si set<int>
#define pii pair<int,int> 
#define pi acos(-1.0)
#define pb push_back
#define mp make_pair
#define lowbit(x) (x&(-x))
#define sci(x) scanf("%d",&(x))
#define scll(x) scanf("%I64d",&(x))
#define sclf(x) scanf("%lf",&(x))
#define pri(x) printf("%d",(x))
#define rep(i,j,k) for(int i=j;i<=k;++i)
#define per(i,j,k) for(int i=j;i>=k;--i)
#define mem(a,b) memset(a,b,sizeof(a)) 
using namespace std;
int n,m,a[maxn][maxm],dis[maxn][maxn],dp[1<<maxn][maxn],res;
//dp[i][j]代表当前状态为i接下来从j出发 
int main()
{ 
   sci(n),sci(m);
   rep(i,0,n-1)
   {
   	rep(j,0,m-1)
   	{
   		sci(a[i][j]);
   	}
   }
   rep(i,0,n-1)
   {
   	rep(j,i+1,n-1)
   	{
   		dis[i][j]=INF;
   		rep(k,0,m-1)
   		{
   			dis[i][j]=min(dis[i][j],abs(a[i][k]-a[j][k]));
   		}
   		dis[j][i]=dis[i][j];
   	}
   }
   rep(s,0,n-1)//枚举第一行 
   {
   	mem(dp,0);
   	dp[(1<<s)][s]=INF; 
   	rep(i,(1<<s),(1<<n)-1)//显然只能转移到比1<<s大的地方去 
   	{
   	 rep(from,0,n-1)//当前行是from 
   	 { 
   	  if(dp[i][from])//如果from的1在i中即状态合法,合法状态->下一个合法状态 
   	  {
   	    rep(to,0,n-1)//下一行是to 
	    { 
	   	 if(!(i&(1<<to)))
	   	 {
	   	 	dp[i|(1<<to)][to]=max(dp[i|(1<<to)][to],min(dp[i][from],dis[from][to]));
	   	 }
	    }
          }
	 }
        }
	 rep(t,0,n-1)//枚举最后一行 显然包含了只有一行的情况 
	 {
		int ans=dp[(1<<n)-1][t];
		rep(i,0,m-2)//从最后一行转到下一列的第一行 
		{
			ans=min(ans,abs(a[t][i]-a[s][i+1])); 
		}
		res=max(res,ans);
	 } 
    }
    printf("%d\n",res);
   return 0;
}