动规-Subset Sums

题目如下：

对于从1到N的连续整数集合，能划分成两个子集合，且保证每个集合的数字和是相等的。

举个例子，如果N=3，对于{1，2，3}能划分成两个子集合，他们每个的所有数字和是相等的：

{3}and{1,2}

这是唯一一种分发（交换集合位置被认为是同一种划分方案，因此不会增加划分方案总数）

如果N=7，有四种方法能划分集合{1，2，3，4，5，6，7}，每一种分发的子集合各数字和是相等的:

{1,6,7}and {2,3,4,5} {注 1+6+7=2+3+4+5}

{2,5,7}and {1,3,4,6}

{3,4,7}and {1,2,5,6}

{1,2,4,7}and {3,5,6}

给出N，你的程序应该输出划分方案总数，如果不存在这样的划分方案，则输出0。程序不能预存结果直接输出。

PROGRAMNAME: subset

INPUTFORMAT

输入文件只有一行，且只有一个整数N

SAMPLEINPUT (file subset . in)

7

OUTPUTFORMAT

输出划分方案总数，如果不存在则输出0。

SAMPLEOUTPUT (file subset.out)

4

 

题目类型：动态规划

思路：F[i][j]表示用1到i中的数拼出和为j的方案数。方程为F[i][j]=F[i-1][j]+F[i-1][j-i]，边界状态为F[i][0]=1。

对于每个数字i，只有两种情况，用和不用。不用的方案数即为F[i-1][j]，用的方案数为F[i-1][j-i]。

由于题目要求两个集合个数相等，所以从1到n的和必须为偶数，即n%4==2或3  且每个集合的和应为(n+1)*n/4。由于交换集合位置被认为是同一种划分方案，所以最终答案应为f[n][(n+1)*n/4]/2。

 

 

#include <fstream>

 

#define I_F "subset.in"

#define O_F "subset.out"

#define MAXn (39+1)

 

short n;

long long ans;

 

void Input();

void Dynap();

void Output();

 

int main()

{

    Input();

    if ((n%4==0)||(n%4==3))

        Dynap();

    Output();

    return 0;

}

 

void Input()

{

    std::ifstream fin(I_F);

    fin>>n;

    fin.close();

}

 

void Dynap()

{

    long long f[MAXn][(1+MAXn)*MAXn/2]={{0}};

    for (short i=0; i<=n; f[i++][0]=1);

    f[1][1]=1;

    for (short i=2; i<=n; i++)

        for (int j=0; j<=(i+1)*i/2; j++)

            f[i][j]=f[i-1][j]+f[i-1][j-i];

    ans=f[n][(1+n)*n/4]/2;

}

 

void Output()

{

    std::ofstream fout(O_F);

    fout<<ans<<std::endl;

    fout.close();

}