ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

P2258 [NOIP2014 普及组] 子矩阵

2021-05-04 09:06:12  阅读:188  来源: 互联网

标签:NOIP2014 20 P2258 int memset 矩阵 dfs sizeof include


题目

题目

思路

暴力枚举每一列情况,然后设 f x , y f_{x,y} fx,y​为前y行选x个的最优解, d x d_x dx​为x列所需代价 e x , y e_{x,y} ex,y​为x,y2行相邻的代价,则有:
f i , j = m i n ( f i − 1 , k + d j + e k , j ) ( 1 < = i < = c , i < = j < = m , 0 < = k < j ) f_{i,j}=min(f_{i-1,k}+d_j+e_{k,j})(1<=i<=c,i<=j<=m,\color{red}{0<=k<j}) fi,j​=min(fi−1,k​+dj​+ek,j​)(1<=i<=c,i<=j<=m,0<=k<j)
code:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring> 
using namespace std;
int n,m,r,c,a[20][20],b[20],d[20][20],e[20][20],f[20],ans=1<<30;
void dp()
{
	memset(b,0,sizeof(b));
	memset(e,0,sizeof(e));
	memset(d,127,sizeof(d));
	for (int i=1;i<=m;i++) for (int j=1;j<r;j++) b[i]+=abs(a[f[j]][i]-a[f[j+1]][i]);
	for (int i=1;i<=m;i++) for (int j=i+1;j<=m;j++) for (int k=1;k<=r;k++) e[i][j]+=abs(a[f[k]][i]-a[f[k]][j]);
	d[0][0]=0;
	for (int i=1;i<=c;i++)
	{
		for (int j=i;j<=m;j++)
		{
			for (int k=0;k<j;k++) d[i][j]=min(d[i][j],d[i-1][k]+b[j]+e[k][j]);
			if (i==c) ans=min(ans,d[c][j]);
		}
	}
	return;
}
void dfs(int x,int y)
{
	if (y>r)
	{
		dp();
		return;
	}
	if (x>n) return;
	dfs(x+1,y);
	f[y]=x;
	dfs(x+1,y+1);
	return;
}
signed main()
{
	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,1);
	cout<<ans;
	return 0;
}

标签:NOIP2014,20,P2258,int,memset,矩阵,dfs,sizeof,include
来源: https://blog.csdn.net/weixin_49843717/article/details/116393745

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有