学步园

题意：价值为1，2，3，4，5，6的石子分别有n1，n2，n3，n4，n5，n6个，问能否把石子分成价值相等的两份（石子总个数 <= 20000）。

题目链接：http://poj.org/problem?id=1014

——>>转换为二进制拆分+01背包。。

二进制拆分：把第i种石子拆分成ni份，对于每一份，可取可不取，可用01背包，但这样时间开销较大。有这样一个事实：对于一个正整数n，小于等于n的任何一个正整数都可以用1, 2, 4, 8, 16, ..., 2^k, n-(2^(k+1)-1)中的几个数的和表示（想想n的二进制表示，小于等于n的正整数中，对于 <=2^k 的正整数，二进制位对应就是；对于 > 2^k 且 <= n 的正整数x，设前面的1 + 2 + 4 + 8 + 16 + ... + 2^k = s1, n-(2^(k+1)-1) =
s2，则n = s1 + s2，则x = s1 + s2 - y，y是s1中某些数的和。），这样同样可以用01背包且时间开销减小了。。。

——>>将同一种石子重新分堆，分别分成2^0个，2^1个，2^2个，2^3个……

状态：设dp[i][j]为从前i堆石子取石子组成价值不超过 j 的最大价值是多少。

状态转移方程：dp[i][j] = max(dp[i - 1][j], dp[i - 1][j - nVal[i]] + nVal[i]);

使用滚动数组优化为一维。。

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

const int maxn = 20000 + 10;

int w[maxn], cnt, d[maxn*6];

void init() {
    cnt = 0;
}

void binary_split(int *a, int n) {
    for(int i = 1; i <= n; i++) {
        int sum = 0;
        for(int j = 1; sum+j <= a[i]; j <<= 1) {
            w[++cnt] = i * j;
            sum += j;
        }
        if(sum < a[i]) w[++cnt] = i * (a[i] - sum);
    }
}

bool dp(int sum) {
    memset(d, 0, sizeof(d));
    for(int i = 1; i <= cnt; i++) {
        for(int j = sum; j >= 0; j--)
            if(j >= w[i]) d[j] = max(d[j], d[j-w[i]] + w[i]);
    }
    return d[sum] == sum;
}

int main()
{
    int a[7], kase = 0;
    while(scanf("%d%d%d%d%d%d", a+1, a+2, a+3, a+4, a+5, a+6) == 6) {
        if(!a[1] && !a[2] && !a[3] && !a[4] && !a[5] && !a[6]) return 0;
        printf("Collection #%d:\n", ++kase);
        int sum = 0;
        for(int i = 1; i <= 6; i++) sum += i * a[i];
        if(sum&1) {
            puts("Can't be divided.\n");
        }
        else {
            init();
            binary_split(a, 6);
            dp(sum/2) ? puts("Can be divided.\n") : puts("Can't be divided.\n");
        }
    }
    return 0;
}

2014.9.24

#include <cstdio>
#include <cstring>
#include <algorithm>

using std::max;

const int MAXN = 6;
const int MAXV = 100;
const int MAXDP = 210000 + 10;

int nCnt[MAXN + 1];
int nVal[MAXV], nNewCnt;
int dp[MAXDP];

void Split()
{
    nNewCnt = 0;
    for (int i = 1; i <= MAXN; ++i)
    {
        for (int j = 1; nCnt[i] - j >= 0; nCnt[i] -= j, j <<= 1)
        {
            nVal[++nNewCnt] = i * j;
        }
        if (nCnt[i]) nVal[++nNewCnt] = i * nCnt[i];
    }
}

bool ZeroOnePack(int target)
{
    memset(dp, 0, sizeof(dp));

    for (int i = 1; i <= nNewCnt; ++i)
    {
        for (int j = target; j >= 0; --j)
        {
            if (j >= nVal[i])
            {
                dp[j] = max(dp[j], dp[j - nVal[i]] + nVal[i]);
            }
        }
    }

    return dp[target] == target;
}

int main()
{
    int nCase = 0;

    while (scanf("%d%d%d%d%d%d", nCnt + 1, nCnt + 2, nCnt + 3, nCnt + 4, nCnt + 5, nCnt + 6) == 6)
    {
        if (!nCnt[1] && !nCnt[2] && !nCnt[3] && !nCnt[4] && !nCnt[5] && !nCnt[6])
        {
            break;
        }

        int nSumValue = 0;
        for (int i = 1; i <= MAXN; ++i)
        {
            nSumValue += i * nCnt[i];
        }

        printf("Collection #%d:\n", ++nCase);
        if (nSumValue & 1)
        {
            puts("Can't be divided.");
        }
        else
        {
            Split();
            ZeroOnePack(nSumValue >> 1) ? puts("Can be divided.") : puts("Can't be divided.");
        }
        puts("");
    }

    return 0;
}