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)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • 域名怎么与主机空间绑定的_域名绑定虚拟主机

    域名怎么与主机空间绑定的_域名绑定虚拟主机域名怎么与主机空间绑定2008-07-2623:41由于各种原因,我们有时候需要在一个IP地址上建立多个web站点,在IIS中,我们可能通过简单的设置达到这个目标。  在IIS中,每个Web站点都具有唯一的、由三个部分组成的标识,用来接收和响应请求:  1、IP地址  2、端口号  3、主机头名。  在IIS中,在一个IP地址上建立多个独

    2022年10月15日
    0
  • php 字符串转换时间_php 字符时间如何转换「建议收藏」

    php 字符串转换时间_php 字符时间如何转换「建议收藏」php字符时间转换的方法:1、通过php中的“strtotime()”函数将任何英文文本的日期时间描述解析为时间戳;2、使用php中的“mktime()”函数从日期取得时间戳即可。本文操作环境:windows7系统、PHP5.6版,DELLG3电脑。php字符串转时间戳PHP提供了函数可以方便的将各种形式的日期转换为时间戳,该类函数主要是:strtotime():将任何英文文本的日期时间描述解…

    2022年6月2日
    29
  • linux下安装tomcat8

    linux下安装tomcat8(一)安装JDK环境64位JDK包:http://download.oracle.com/otn-pub/java/jdk/8u45-b14/jdk-8u45-linux-x64.tar.gz1.用 SecureCRT链接上linux,用命令直接下载 wget http://download.oracle.com/otn-pub/java/jdk/8u45-b14/jdk-8u45

    2022年6月2日
    39
  • HTTP、HTTPS等常用的默认端口号

    HTTP、HTTPS等常用的默认端口号最近部署一个接口访问对方域名,https开头Url,端口配成80,犯了个低级错误。重新学习一下。端口号标识了一个主机上进行通信的不同的应用程序。1,HTTP服务器,默认端口号为80/tcp(木马Executor开放此端口)2,HTTPS(securelytransferringwebpages)服务器,默认端口号为443/tcp  443/udp3,HTTP协议代理服务器常用…

    2022年6月29日
    22
  • 1024,节日快乐!

    1024,节日快乐!每个你不了解的新世界,都值得说一声“Hello,World!”,1024,祝大家节日快乐!少熬夜,多赚钱,护好发!

    2022年6月11日
    34
  • Java中executeUpdate和executeQuery的区别

    Java中executeUpdate和executeQuery的区别executeQueryexecuteQuery()方法是在PreparedStatement对象中执行SQL查询,并返回该查询生成的ResultSet对象,executeUpdate该方法是在PreparedStatement对象中执行SQL语句,该语句必须是一个SQL数据操作语言,例如INSERT,UPDATE,DELETE等语句,或者是没有返回内容的SQL语句….

    2022年10月20日
    0

发表回复

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

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