最大子列和,及其扩展最大子矩阵和,都是比较经典的dp,把这两题写在一起,共大家参考学习。
~~~~
lz弱菜啊,到现在还只能写这么水的DP...orz。
~~~~
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1003
http://acm.hdu.edu.cn/showproblem.php?pid=1081
题意应该比较好理解,不解释了。
令dp[i]表示以第i个字符开始的的最大子列和。
故状态转移就是:dp[i]=max(dp[i+1]+a[i],a[i]);
1003:
#include<cstdio>
#include<iostream>
#include<cstring>
#define INF 1111
#define N 111111
using namespace std;
int a[N],dp[N];
int main()
{
int t,T=0;
scanf("%d",&t);
while(t--)
{
int n;
scanf("%d",&n);
memset(dp,0,sizeof(dp));
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
dp[n]=a[n];
if(T) puts("");
printf("Case %d:\n",++T);
int ans=dp[n],l=n,r=n;
for(int i=n-1;i>=1;i--)
{
dp[i]=max(dp[i+1]+a[i],a[i]); //状态转移
if(ans<=dp[i])
{
ans=dp[i];
l=i; //记录最大子序列区间的左端点
}
}
int s=0;
for(int i=l;i<=n;i++) //根据左端点找右端点
{
s += a[i];
if(s == ans)
{
r=i;
break;
}
}
printf("%d %d %d\n",ans,l,r);
}
}
还有一种写法:
但其实第一种思想更重要。
#include<cstdio> //这里没有输出答案的区间
int main()
{
int n,test,ans,t,a,i;
scanf("%d",&test);
while(test--)
{
scanf("%d",&n);
scanf("%d",&a);
ans=t=a; //t记录当前保存的最大字段和
for(i=1;i<n;i++)
{
scanf("%d",&a);
if(t<0) t=a;
else t=t+a;
if(ans<t) ans=t;
}
printf("%d\n",ans);
}
return 0;
}
~~~~
HDU 1081:
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 111
#define INF 0x3f3f3f3f
using namespace std;
int n;
int g[N][N],sub[N],f[N];
int dp()
{
int Max=-INF;
f[n-1]=sub[n-1];
for(int i=n-2;i>=0;i--)
{
f[i]=max(f[i+1]+sub[i],sub[i]);
Max=max(Max,f[i]);
}
return Max;
}
int qyj()
{
int ans=-INF;
for(int i=0;i<n;i++) //枚举行
{
memset(sub,0,sizeof(sub)); //~~
for(int j=i;j<n;j++) //以该行为边界扩展的的所有子矩阵
{
for(int k=0;k<n;k++)
{
sub[k] += g[j][k]; //子矩阵
}
int sum=dp(); //~~
ans=max(ans,sum);
}
}
return ans;
}
int main()
{
while(scanf("%d",&n)==1)
{
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
scanf("%d",&g[i][j]);
int ans=qyj();
printf("%d\n",ans);
}
return 0;
}
1081不理解的可以手写画画得到子矩阵的那段代码、