hdu1024—最大m段不重合子序列和:策略优化

hdu1024—最大m段不重合子序列和:策略优化

http://acm.hdu.edu.cn/showproblem.php?pid=1024(到这里提交)

Max Sum Plus Plus

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 4977    Accepted Submission(s): 1627

Problem Description
Now I think you have got an AC in Ignatius.L’s “Max Sum” problem. To be a brave ACMer, we always challenge ourselves to more difficult problems. Now you are faced with a more difficult problem.

Given a consecutive number sequence
S1, S2, S3, S4 … Sx, … Sn(
1 ≤ x ≤ n ≤ 1,000,000, -32768 ≤ Sx ≤ 32767). We define a function sum(i, j) =
Si + … + Sj (
1 ≤ i ≤ j ≤ n).

Now given an integer m (m > 0), your task is to find m pairs of i and j which make
sum(i1, j1) + sum(i2, j2) + sum(i3, j3) + … + sum(im, jm) maximal (
ix ≤ iy ≤ jx or ix ≤ jy ≤ jx is not allowed).

But I`m lazy, I don’t want to write a special-judge module, so you don’t have to output m pairs of i and j, just output the maximal summation of sum(i
x, j
x)(
1 ≤ x ≤ m) instead. ^_^

 

 

Input
Each test case will begin with two integers m and n, followed by n integers S
1, S
2, S
3 … S
n.

Process to the end of file.
 

 

Output
Output the maximal summation described above in one line.
Sample Input

1 3 1 2 32 6 -1 4 -2 3 -2 3

Sample Output

68

Hint

Huge input, scanf and dynamic programming is recommended.

 最终程序代码:

 

ContractedBlock.gif
ExpandedBlockStart.gif
代码


1
#include
<
stdio.h
>


2
 
int
m,n,a[
1000001
],ans;

3
 
int
f[
1000001
],sum[
1000001
],b[
2
][
1000001
];

4
 
int
main(){


5

int
i,j,k,i1,i2;

6

while
(scanf(

%d%d

,
&
m,
&
n)
!=
EOF)

7
{


8
sum[
0
]
=
0
;

9

for
(i
=
1
;i
<=
n;i
++
)

10
{


11
scanf(

%d

,
&
a[i]);

12
sum[i]
=
sum[i

1
]
+
a[i];

13
}

14


15

for
(j
=
0
;j
<=
n;j
++
)b[0][j]=0;————————前0个数分成j段的最大和为0

16
for ( j = 0; j <=n; j++) b[1][j]= 2000000000; //已改动


17
i1
=
1
;ans
=-
2000000000
;

18


19

for
(i
=
1
;i
<=
m;i
++
)

20
{ b[i1][j]每一次都应该初始化,因为空间是循环利用的

21

b[i1][i-1]
=


2000000000
; //已改动

22

for
(j
=
i;j
<=
n;j
++
)——————从数字数不少于分段数的地方开始

23
{


24

if
(i
==
j)f[j]
=
sum[j];

25

else


26
{

27
f[j]
=
f[j

1
]
+
a[j];

28

if
(f[j]
<
b[
1

i1][j

1
]
+
a[j])f[j]
=
b[
1

i1][j

1
]
+
a[j];

29
}

30
b[i1][j]
=
b[i1][j

1
];

31

if
(b[i1][j]
<
f[j])b[i1][j]
=
f[j];

32
}

33
i1
=
1

i1;—————i1总是0-1交替变化

34
}

35

for
(i
=
m;i
<=
n;i
++)——————此时f[i]中存储前i个数分成m段得到的最大和


36

if
(ans
<
f[i])ans
=
f[i];

37
printf(

%d\n

,ans);

38
}

39

return

0
;

40
}

 

 

此题难在两点:1.找到状态转移方程;2.优化程序。

状态转移方程:f[i][j]=max{f[i-1][j],max{f[i-1][k]}}+a[j];         (i-1=<k<=j-1)

f[i][j]代表将前j个数分成i段可以得到的最大和,并且这i段中必须包括a[j].它可以由两个状态转移到:a[j]与其他数一起包含

到第j段中,所得到的最大和(f[i-1][j]+a[j]);a[j]独立成段,与前i-1个数分成的j-1段得到的最大和(max{f[i-1][k]}+a[j])。

 

由于问题所涉及的规模实在太大,所以必须用滚动数组存储状态。同时设置sum[]数组,其中存储当i=j即分段数等于数字数时的

状态值,也就是所有j个数的和。主要代码见下:

ContractedBlock.gif
ExpandedBlockStart.gif
代码


1
i1
=
1
;

2
ans
=-
(
1
<<
30
);

3

for
(i
=
1
;i
<=
m;i
++
)

4
{

5

for
(j
=
i;j
<=
n;j
++
)

6
{


7

if
(i
==
j)f[i1][j]
=
sum[j];

8

else


9
{

10
f[i1][j]
=
f[i1][j

1
]
+
a[j];

11

for(k=i-1;k<=j-1;k++
)

12

if(f[i1][j]<f[1-i1][k]+a[j])f[i1][j]=f[1-i1][k]+
a[j];

13
}

14
}

15
i1
=
1

i1;

16
}

17
j
=
m
%
2
;

18

for
(i
=
m;i
<=
n;i
++
)

19

if
(ans
<
f[j][i])ans
=
f[j][i];

20
printf(

%d\n

,ans);

红色代码区就是上述转移状态的第二个状态。这个代码的效率不高,很大程度上是因为这个循环。优化办法如下:

设置一个数组b[i][j],它是滚动式的,用来存储前j个数划为i段得到的最大和,但不一定包括a[j]。这样转移方程就可以变

为这样: f[i][j]=max{f[i-1][j],b[i-1][j-1]}+a[j];

             b[i][j]=max{b[i][j-1],f[i][j]}

这样的话,在程序运行是就包含了两个转移方程,b[i][j]的值由两种状态转移可得:前j个数分为i段选a[j]得到的最大和(f[i][j]);

不选a[j]得到的最大和,即前j-1个数分为i段得到的最大和(b[i][j-1])。此时f[i][j]还可以省去一维,变为f[j].最终优

化后的程序见最上部,还有若干细节问题在程序旁会有注释。

 

转载于:https://www.cnblogs.com/aiyite826/archive/2010/07/24/1784449.html

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/110749.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • Oscium双频带频谱分析仪「建议收藏」

    Oscium双频带频谱分析仪「建议收藏」新一代Wipry5X为频谱分析仪WipryCombo的升级版本,主要针对2.4&5GHz两个频段的wifi信号进行测试,同时兼容苹果iOS和Android双系统;具有双带频谱分析功能,智能的触摸屏操作,简便快捷,口袋型轻巧的设计,对数据记录分析,随时发送E-mail功能,让你的测试测量变得随时随地,不再受时间、空间、地理位置的影响。WiPry5X性能指标

    2022年8月11日
    9
  • visio2013密钥激活_visio2013不激活可以用吗

    visio2013密钥激活_visio2013不激活可以用吗  2NYF6-QG2CY-9F8XC-GWMBW-29VV8  FJ2N7-W8TXC-JB8KB-DCQ7Q-7T7V3  VXX6C-DN3HQ-3CRXG-RF4KT-YG7V3  B3C7Q-D6NH2-2VRFW-HHWDG-FVQB6  TCWJK-N6GFH-82BP9-HV7YQ-T6KMQ  KD8CP-DN968-RGQM6-PMFMK-M96YQ  …

    2025年8月28日
    7
  • vsftp日志查看_vsftp日志配置

    vsftp日志查看_vsftp日志配置1、开始vsftp记录日志。修改/etc/vsftpd/vsftpd.conf如下:xferlog_enable=YESxferlog_std_format=YESxferlog_file=/var/log/xferlogFTP服务器的日志设置,可以通过修改主配置文件/etc/vsftpd.conf实现。主配置文件中与日志设置有关的选项包括xferlog_enable、xferlog_file…

    2025年6月13日
    8
  • map怎么转json对象_json怎么获取map

    map怎么转json对象_json怎么获取map如何把JSON对象转为map对象呢?JSON对象保存在大括号内。就像在JavaScript中,对象可以保存多个键/值对。Map对象保存键/值对,是键/值对的集合。任何值(对象或者原始值)都可以作为一个键或一个值。Object结构提供了“字符串—值”的对应,Map结构提供了“值—值”的对应。javascript将JSON对象转为map对象可以利用阿里巴巴封装的FastJSON来转换。有多种…

    2022年8月23日
    7
  • git 使用总结

    git 使用总结

    2021年11月7日
    87
  • android10锁屏时钟样式,三星s10息屏时钟[通俗易懂]

    android10锁屏时钟样式,三星s10息屏时钟[通俗易懂]三星s10息屏时钟设置软件相当的方便强大,是不是设置息屏时钟太麻烦了,或者横竖屏总是错乱不尽人意?不用再在网上找各种攻略啦,这款软件可以直接帮助你修改息屏时钟,一键操作设置即可!同时这款软件中还有这更多的炫酷功能,图标壁纸时钟主题等等自定义和制作功能相当的强大!喜欢的朋友们不要错过哦!赶快下载三星s10息屏时钟开始体验吧!三星s10息屏时钟介绍三星s10息屏时钟软件app是一款提取版三星s10息屏…

    2022年9月29日
    5

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号