传送门
题意:给定了一个n*m的矩阵,(1<=n<=16,1<=m<=10000) 。我们可以任意交换行的顺序,交换完成后。我们定义k为:从第一列开始向下遍历元素,完成后从第二列开始向下遍历… 得到a11 a21…an1 a12 a22…an2…anm . k为相邻元素差值绝对值的最小值。
求k最大可以是多少?
分析:我们看到n的范围,大概可以猜到是状压dp 。
-
我们具体分析遍历方式可以发现,相邻元素的差值最小值其实就是i j行元素对应位置相差的最小值,还有首末行是错开一个位置相差最小值。 我们可以预处理得出。
-
dp[i][j][k]
:状态为i(1表示这行目前选了 0表示还没),起始行为j,当前选了第k行时,任意两行的相差最小值(我们先不考虑首末行,最后再考虑) -
dp[1<<(i-1)][i][i]
作为边界可以求出 -
dp[sta][st][to]=max(dp[sta][st][to],min(dp[st][ed][last],dis[ed][to])
-
具体的含义就看代码的注释吧 ,有个注意的地方是我们dp状态定义的是不考虑首末行的最优解,如果上一个状态只有1行,那么直接拿
dis[ed][to]
更新当前状态就好了。
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pii pair<int,int>
const int maxn = 1e4+10;
const int mx = 40;
const int mod = 1e9+7;
const ll inf = 34359738370;
const int INF = 1e9+7;
int dis[17][17],tr[17][17];//i j行相差的最小值 / 首行是i 末行是j 相差最小值
int a[17][maxn];
int n,m;
int dp[1<<16][17][17];//状态是k 首行是i 最后是j行 不考虑首末行时的最优解
int main()
{
scanf("%d %d",&n,&m);
memset(dp,-1,sizeof dp);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++) scanf("%d",&a[i][j]);
}
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
int mmin=INF;
for(int k=1;k<=m;k++) mmin=min(mmin,abs(a[i][k]-a[j][k]));
dis[i][j]=mmin;
mmin=INF;
for(int k=1;k<m;k++) mmin=min(mmin,abs(a[i][k+1]-a[j][k]));//i是首行
tr[i][j]=mmin;
}
}
int k=0;
for(int i=1;i<=n;i++)
{
dp[1<<(i-1)][i][i]=tr[i][i];//
}
for(int i=1;i<(1<<n);i++)
{
int cnt=0;//记录这个状态选了几行
for(int s=1;s<=n;s++) if(i & (1<<(s-1))) cnt++;
for(int st=1;st<=n;st++)
{
for(int ed=1;ed<=n;ed++)
{
if(dp[i][st][ed] != -1) //合法状态
{
for(int to=1;to<=n;to++) //下一行选啥
{
if(i & (1<<(to-1))) continue;
int sta=i | (1<<(to-1));
if(cnt==1) //此时st==ed 不能拿dp[i][st][ed]去更新 因为这是首末行的解
dp[sta][st][to]=max(dp[sta][st][to],dis[ed][to]);
else dp[sta][st][to]=max(dp[sta][st][to],min(dp[i][st][ed],dis[ed][to]));
}
}
}
}
}
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
dp[(1<<n)-1][i][j]=min(dp[(1<<n)-1][i][j],tr[i][j]);//考虑首末行
k=max(k,dp[(1<<n)-1][i][j]);
}
}
cout<<k;
return 0;
}