Codeforces Round #531 F. Elongated Matrix （状压dp）

传送门

题意：给定了一个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;
}