noip2015普及组第四题

来源：互联网发布：用户网络数据监控编辑：程序博客网时间：2024/05/21 21:43

不会写全动归，只能写半动归半暴搜。主要DFS枚举行数，dp目前矩阵，类似求最大和，预处理好每一列自身的值以及各行间的值就可以了。

#include<cstdio>#include<cstring>#include<cstdlib>#include<algorithm>#define maxn 214748364using namespace std;int n, m, r, c, ans=maxn;int v[18], a[18][18], s1[18], s2[18][18], dp[18][18];int dpp(){memset(s1,0,sizeof(s1));memset(s2,0,sizeof(s2));memset(dp,0,sizeof(dp));for(int i=1;i<=m;i++){    for(int j=1;j<v[0];j++) s1[i]+=abs(a[v[j]][i]-a[v[j+1]][i]);    for(int j=i+1;j<=m;j++)    for(int k=1;k<=v[0];k++) s2[i][j]+=abs(a[v[k]][i]-a[v[k]][j]);}for(int i=1;i<=m;i++)for(int j=1;j<=m;j++) dp[i][j]=maxn;for(int i=1;i<=m;i++) dp[1][i]=s1[i];for(int i=2;i<=c;i++)for(int j=i;j<=m;j++)for(int k=j-1;k>0;k--)dp[i][j]=min(dp[i][j],dp[i-1][k]+s1[j]+s2[k][j]);int ans=maxn;for(int i=1;i<=m;i++) ans=min(ans,dp[c][i]);return ans;}void dfs(int pos, int z){if(z==r){ans=min(ans,dpp());return;}for(int i=pos+1;i<=m;i++){v[++v[0]]=i;dfs(i,z+1);v[0]--;}}int main(){scanf("%d %d %d %d", &n, &m, &r, &c);for(int i=1;i<=n;i++)for(int j=1;j<=m;j++) scanf("%d", &a[i][j]);dfs(0,0);printf("%d\n", ans);return 0;}

1 0