P2258 [NOIP2014 普及组] 子矩阵 day16

你确定不把这道题升紫？

还是说我以前做的蓝都太水了。。。

暴力只有五十分

我们思考枚举行 然后贪心

好的显然没有任何贪心规律 爆搜肯定是会T的

那怎么办呢？？？

别急！！！我知道你很急！！！因为我也没想到啊！！！

我们尝试推出一个状态转移方程

因为行确定所以每个纵向分数也确定了！

我们常规的dp[i][j]一般表示的都是你前i行选择了j个的 max or min or cnt

但是这题不行啊 因为你不确定j前面选择了谁所以就很难写啊！

那我们只好调整一下dp[i][j]用i来表示你前面i个里面 并且选了i这个数！！！！

一共选了j个来推 就非常好推了

我们可以先预处理出 up[i]表示 第i列的纵向分数和总和

lft[i][j]表示第i个列到第j列的横向值 （可以模拟一下啊 一个子矩阵的sum就是up+lft）

然后我们就可以来推dp[i][j]=枚举前面i个点dp[k][j-1]+up[i]+lft[k][i]

这样子看起来还是非常暴力啊！

实际上的时间复杂度只有O（C2n*n^3）因为我们dp的时候相当于加了格外多的剪枝 把不够优秀的开局全部筛掉了

#include <bits/stdc++.h>
using namespace std;
const int N = 18;
const int M = 1e4+10;
const int mod = 1e9+7;
#define int long long
#define endl '\n'
#define Endl '\n'
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int a[N][N],ch[N],n,m,r,c,ans=2e9,dp[N][N],up[N],lft[N][N],cnt;
int calc(){
    memset(dp,0x3f3f,sizeof dp);
    int res=inf;
    for(int i=1;i<=m;i++){
        dp[i][1]=up[i];
        for(int j=2;j<=c;j++){
            for(int k=1;k<i;k++){
                dp[i][j]=min(dp[k][j-1]+up[i]+lft[k][i],dp[i][j]);
            }
        }
        res=min(res,dp[i][c]);
    }
    return res;
}
void pre(){
    memset(up,0,sizeof up);
    memset(lft,0,sizeof lft);
    for(int i=1;i<=m;i++)
        for(int j=2;j<=r;j++)
            up[i]+=abs(a[ch[j]][i]-a[ch[j-1]][i]);
    for(int i=1;i<=r;i++)
        for(int j=1;j<=m;j++)
            for(int k=j+1;k<=m;k++)
                lft[j][k]+=abs(a[ch[i]][j]-a[ch[i]][k]);
}
void dfs(int step){
    if(cnt==r){
        pre();
        ans=min(ans,calc());
        return;
    }
    if(step==n+1)return;
    for(int i=ch[step-1]+1;i<=n;i++){
        ch[step]=i;
        cnt++;
        dfs(step+1);
        ch[step]=0;
        cnt--;
    }
}
signed main(){
    fast
    cin>>n>>m>>r>>c;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            cin>>a[i][j];
        }
    }
    dfs(1);
    cout<<ans<<endl;
    return 0^0;
}

想了一下普遍的dp想法好像还是可以 我发了再去写下试试看