题目链接:

https://www.lydsy.com/JudgeOnline/problem.php?id=1084

题目大意:

这里有一个n*m的矩阵,请你选出其中k个子矩阵,使得这个k个子矩阵分值之和最大。注意:选出的k个子矩阵不能相互重叠。n<=100 m<=2

思路:

m=1时,就是数组选出k个连续子段和最大。

dp[i][j]表示前i个数中已经取了j个连续子段和的最优解。

dp[i][j] = dp[i-1][j] 不取这个数

dp[i][j] = dp[start-1][j-1] + s[start]+...+s[i] 取这个数,并且从start到i作为第j个连续段(求区间和直接用前缀和求差代替)

m=2时

dp[i][j][k]表示第一列前i个数 第2列前j个数,已经取了k个子矩阵的最优解

dp[i][j][k] = max(dp[i-1][j][k], dp[i][j-1][k]) 不取这个数

dp[i][j][k] = dp[start-1][j][k-1] + s[start][1]+...+s[i][1] 第一列从start到i取出来作为第k个子矩阵

dp[i][j][k] = dp[i][start-1][k-1] + s[start][2]+...+s[i][2] 第二列从start到i取出来作为第k个子矩阵

dp[i][j][k] = dp[start-1][start-1][k-1] + s[start][1]+...+s[i][1]+s[start][2]+...+s[j][2]当且仅当i==j 两列同时取。

同样的,区间求和用前缀和快速求出。

时间复杂度为O(n^3*k)

 #include<bits/stdc++.h>
#define IOS ios::sync_with_stdio(false);//不可再使用scanf printf
#define Max(a, b) ((a) > (b) ? (a) : (b))//禁用于函数,会超时
#define Min(a, b) ((a) < (b) ? (a) : (b))
#define Mem(a) memset(a, 0, sizeof(a))
#define Dis(x, y, x1, y1) ((x - x1) * (x - x1) + (y - y1) * (y - y1))
#define MID(l, r) ((l) + ((r) - (l)) / 2)
#define lson ((o)<<1)
#define rson ((o)<<1|1)
#define Accepted 0
#pragma comment(linker, "/STACK:102400000,102400000")//栈外挂
using namespace std;
inline int read()
{
int x=,f=;char ch=getchar();
while (ch<''||ch>''){if (ch=='-') f=-;ch=getchar();}
while (ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
return x*f;
} typedef long long ll;
const int maxn = + ;
const int MOD = ;//const引用更快,宏定义也更快
const int INF = 1e9 + ;
const double eps = 1e-;
int dp1[maxn][];//m = 1 dp1[i][j]表示前i个数中取了j个连续段的最大值
int dp2[maxn][maxn][];//m = 2 dp2[i][j][k]表示第一列前i个数第二列前j个数取了k个连续段的最大值
int s[maxn];
int sum[maxn][];
int main()
{
int n, m, K;
scanf("%d%d%d", &n, &m, &K);
if(m == )
{
for(int i = ; i <= n; i++)scanf("%d", &s[i]), s[i] += s[i - ];//直接记录前缀和
for(int i = ; i <= n; i++)
for(int j = ; j <= K; j++)dp1[i][j] = -INF;//dp[i][0]均为初始化的0
for(int i = ; i <= n; i++)
{
for(int j = ; j <= K; j++)
{
dp1[i][j] = dp1[i - ][j];//不取这个数字
for(int start = ; start <= i; start++)//从start开始一直取到i作为第j段
{
dp1[i][j] = max(dp1[i][j], dp1[start - ][j - ] + s[i] - s[start - ]);
}
}
}
printf("%d\n", dp1[n][K]);
}
else
{
for(int i = ; i <= n; i++)
{
for(int j = ; j <= m; j++)
{
scanf("%d", &sum[i][j]);
sum[i][j] += sum[i - ][j];//记录每一列的前缀和
}
}
for(int i = ; i <= n; i++)for(int j = ; j <= n; j++)for(int k = ; k <= K; k++)dp2[i][j][k] = -INF;
for(int i = ; i <= n; i++)
for(int j = ; j <= n; j++)for(int k = ; k <= K; k++)
{
dp2[i][j][k] = max(dp2[i - ][j][k], dp2[i][j - ][k]);//不取这个数字
for(int start = ; start <= i; start++)//第一列从start开始一直取到i作为第k个矩阵
dp2[i][j][k] = max(dp2[i][j][k], dp2[start - ][j][k - ] + sum[i][] - sum[start - ][]);
for(int start = ; start <= j; start++)//第二列从start开始一直取到j作为第k个矩阵
dp2[i][j][k] = max(dp2[i][j][k], dp2[i][start - ][k - ] + sum[j][] - sum[start - ][]);
if(i == j)//两列从start开始一直取到i作为第k个矩阵
{
for(int start = ; start <= i; start++)
dp2[i][j][k] = max(dp2[i][j][k], dp2[start - ][start - ][k - ] + sum[i][] + sum[i][] - sum[start - ][] - sum[start - ][]);
}
}
printf("%d\n", dp2[n][n][K]);
}
return Accepted;
}

BZOJ 1084 最大子矩阵 dp的更多相关文章

  1. BZOJ 1084 最大子矩阵

    http://www.lydsy.com/JudgeOnline/problem.php?id=1084 思路:分m=1和m=2操作 #include<algorithm> #includ ...

  2. [SCOI2005][BZOJ 1084]最大子矩阵

    Description 这里有一个n*m的矩阵,请你选出其中k个子矩阵,使得这个k个子矩阵分值之和最大.注意:选出的k个子矩阵不能相互重叠. Input 第一行为n,m,k(1≤n≤100,1≤m≤2 ...

  3. [BZOJ 1084] [SCOI2005] 最大子矩阵 【DP】

    题目链接:BZOJ - 1084 题目分析 我看的是神犇BLADEVIL的题解. 1)对于 m = 1 的情况, 首先可能不取 Map[i][1],先 f[i][k] = f[i - 1][k];   ...

  4. 【SCOI2005】 最大子矩阵 BZOJ 1084

    Description 这里有一个n*m的矩阵,请你选出其中k个子矩阵,使得这个k个子矩阵分值之和最大.注意:选出的k个子矩阵不能相互重叠. Input 第一行为n,m,k(1≤n≤100,1≤m≤2 ...

  5. BZOJ 1084: [SCOI2005]最大子矩阵 DP

    1084: [SCOI2005]最大子矩阵 题目连接: http://www.lydsy.com/JudgeOnline/problem.php?id=1084 Description 这里有一个n* ...

  6. 【BZOJ 1084】 1084: [SCOI2005]最大子矩阵 (DP)

    1084: [SCOI2005]最大子矩阵 Description 这里有一个n*m的矩阵,请你选出其中k个子矩阵,使得这个k个子矩阵分值之和最大.注意:选出的k个子矩阵不能相互重叠. Input 第 ...

  7. 【BZOJ 1084】 [SCOI2005]最大子矩阵(DP)

    题链 http://www.lydsy.com/JudgeOnline/problem.php?id=1084 Description 这里有一个n*m的矩阵,请你选出其中k个子矩阵,使得这个k个子矩 ...

  8. bzoj 1084: [SCOI2005]最大子矩阵【dp】

    分情况讨论,m=1的时候比较简单,设f[i][j]为到i选了j个矩形,前缀和转移一下就行了 m=2,设f[i][j][k]为1行前i个,2行前j个,一共选了k个,i!=j的时候各自转移同m=1,否则转 ...

  9. BZOJ 1084 (SCOI 2005) 最大子矩阵

    1084: [SCOI2005]最大子矩阵 Time Limit: 10 Sec Memory Limit: 162 MB Submit: 3560 Solved: 1779 [Submit][Sta ...

随机推荐

  1. JAVA-2NIO之Channel

    注意:转载自并发编程网 – ifeve.com本文链接地址: Java NIO系列教程(二) Channel Channel Java NIO的通道类似流,但又有些不同: 既可以从通道中读取数据,又可 ...

  2. 说一说HTTP

    什么是URI和URL URI用字符串标示某一互联网资源,而URL表示资源的地点.可见URL是URI的子集. URI要使用涵盖全部必要信息的URI.绝对URL以及相对URL.相对URL是指从浏览器中基本 ...

  3. RxJS库

    介绍 RxJS是一个异步编程的库,同时它通过observable序列来实现基于事件的编程.它提供了一个核心的类型:Observable,几个辅助类型(Observer,Schedulers,Subje ...

  4. io多大合适

    MySQL的cpu iowait值,<=1/(number of cores).例如,如果是8核的,io wait 应<=12%. "Monitor DB CPU IO wait ...

  5. Web开发 学习积累20161018

    项目 一.项目做的是什么 业务逻辑 -> 增删改查 二.什么是面向对象编程,它有哪些好处 oop:object oriented programming <>核心思想:使用人类思考问 ...

  6. 系统分析与设计 homework2

    1. 简述瀑布模型.增量模型.螺旋模型(含原型方法)的优缺点. 瀑布模型 优点: 降低了软件开发的复杂度,提高软件开发过程中的透明性,提高软件开发的可管理性. 为项目提供了按阶段划分的检查点. 当前一 ...

  7. HDU-2046 骨牌铺方格【递推】

    http://acm.hdu.edu.cn/showproblem.php?pid=2046 和前面的一样,a[i] = a[i-1] + a[i-2] #include<iostream> ...

  8. js在ie6下的一个bug—未结束标签的错误

    在IE6下,如果在body标签没结束前,用代码获取body对象就会出现错误.如: <html> <head> <script type="text/javasc ...

  9. Linux常见操作指令(转载)

    Linux,免费开源,多用户多任务系统.基于Linux有多个版本的衍生.RedHat.Ubuntu.Debian 安装VMware或VirtualBox虚拟机.具体安装步骤,找百度. 再安装Ubunt ...

  10. 关于 Blog 修改

    关于 Blog 修改 本 Blog 使用的是 WordPress,每次升级 WordPress 都需要修改文件,以修正一些问题,因此做个总记录,便于自己修改. 解决 WordPress 无法打开中文链 ...