计算拆分方案.docx

【题6】计算拆分方案

任何大于1的自然数n都可以拆分成若干个小于n的自然数之和。输入n，计算和输出n的不同拆分方案。

题解

我们用数组a[1..100]存储已完成的拆分，所有自然数按照递增顺序排列：n=a[1]+a[2]+…a[k](a[1]

≤a[2]≤…≤a[k])。这个多项式能否延长，关键看a[k]能否继续拆分成m和a[k]-m。按照递增要求，当

a[k-1]≤m≤?a[k]?

时，多项式可延长至 n=a[1]+a[2]+…a[k]’+a[k+1]’(a[1]≤a[2]≤…≤a[k]’≤

2? ?

2

? ?

a[k+1]’（a[k]’=m，a[k+1]’=a[k]-m）。设

⑴状态kx

kx为待拆分的元素序号，即多项式中最后一个自然数的序号。由于需要输出每一个拆分方案，因此当前的拆分方案n=a[1]+a[2]+…a[kx]不存在什么边界条件。

⑵搜索范围

设a[kx]拆分成m和a[kx]-m。我们将m作为算符。按照递增要求，a[kx-1]≤m≤?a[kx]?

。只要m在这

2? ?

2

? ?

个范围内，a[kx]即可拆分成a[kx]’=m和a[kx+1]=a[kx]-m，该搜索范围隐含了约束条件。我们用递归程序sum(kx)描述拆分方案的计算过程：

proceduresum(kx:integer)； {从a[kx]出发，递归计算拆分方案}var

k,l,m:integer；begin

s←s+1；

{累计拆分方案数}

write(’sumN0’,s:3,’:’,n,’=’)；

{输出当前的拆分方案}

fork←1tokx-1dowrite(a[k],’+’)；

writeln(a[kx])；

k←kx；l←a[k]；

form←a[k-1]toldiv2do

{

拆分a[kx]=a[k]’+a[k+1]’}

{枚举a[k]’的所有可能值}

begin

a[k]←m；a[k+1]←l-m；

sum(k+1)； {递归拆分a[kx+1]}

end；{for}end；{sum}

调用sum过程的前提是a[1]和a[2]已知。我们按照递增要求，枚举a[1]和a[2]的所有可能值：1≤

a[1]≤i,a[2]=n-a[1](1≤i≤?n?

?2?

),并分别调用sum(2)，便可计算出所有拆分方案：

? ?

ifn=2then

fori←1tondiv2do {枚举a[1]和a[2]的所有可能情况}begin

a[1]←i；a[2]←n-a[1]；

sum(2)； {从当前的a[2]出发进行拆分}

end；{for}