poj 1014 Dividing-动态规划

题目大意:有六种珠子价值分别为1,2,...,6，给定每种珠子的总数，问是否能将这些珠子分解成价值相等的两堆

因为珠子的数目很大所以直接DP肯定会超时，所以要设法减小数据的规模，刚开始听人讲要取模做，觉得好像应该是这么回事的样子，但是自己又不能证明这种方法是对的，然后敲了个代码，WA，把人家取模的例程submit，WA.....最后看了poj上某大神的反例才发现取模的做法是错的....

Posted by dashashi at 2009-12-07 15:42:13 on Problem 1014

---

假设对K取模。如果K是偶数，那么对于数据K+1 0 1 0 0 0，取模之后是无法分割的，但是实际上是可以的。答案是K/2+2。

如果K是奇数(显然是不行的……也给个简单的反例)，那么对于K+1 1 0 0 0 0，取模之后无法分割，答案也是K/2+2，(不过这里的K/2是把多出来的0.5舍去了的)
是测试数据太弱了

重新思考了之后想出了一种方法，不保证是对的...AC了但是也不排除数据太弱的可能性，所以大家如果发现了什么问题就直接说吧...

对于每一种珠子的个数组合，如果它能分解价值相等的两堆的话不一定只有一种分解方式，如果某一个价值的珠子在一种分解方式中被分解成m1和m2，不妨设m1>m2，那么这种价值的珠子总数减少2 * m2个，仍然能被分解，m1 - m2 和 0，如果可以证明分解之后m1 - m2 < 常数l的话，数据规模就可以减小了

先假设减去公共个数之后，可分解成

0 k2 k3 0 0 k6

k1 0 0 k4 k5 0 (每个k在上下出现的位置是随机的)

以k6为例，在这种情况下应有6 * k6 <= 1 * k1 + 4 * k4 + 5 * k5，即k6的大小是受下派总价值控制的

若k1的个数大于6个，则可以将这6个移动到上排，将一个6号珠子移动下来，仍保持两边价值相等，这样再减去公共的数目，应有k1 < 6

同理(若这些k都存在于下排且k6足够大）k2 < 3, k3 < 2, k4 < 3, k5 < 6，存在这样的一种分解方式

在这种分解方式下k6取最大值应当是每一个ki都处在下排，但若1,5 2,4同时存在，它们相加也可消去一个6，故在最大的情况下

6 * k6 <= 3 * k3 + 4 * k4 + 5 * k5 <= 3 * 1 + 4 * 2 + 5 * 5 = 36， 即k6 <= 6

这样我们证明了（或者说，说明了）如果能分解那么一定存在一种分解方式使两堆中6号珠子个数的差<=6

同理有k1, k2, k3, k4, k5 < 8（没严格去计算，这个数也肯定不是上确界）

这样对于每一种珠子，可以将多余8个的部分成对的刨去，问题就简化了

（做完之后查了一下，发现这道题还可以用dfs和二进制优化的方式做...前者没敢用...总觉得会TLE...后者学习了，感觉好厉害的样子...）

下面是代码(之前把代码贴错了...囧...)

#include <iostream>
#include <cstring>
using namespace std;

int f[7][150];
int main()
{
    int n[7], i, j, sum, cnt = 0, k;
    for (;;)
    {
        sum = 0;
        ++cnt;
        for (i = 1; i < 7; ++i)
        {
            cin >> n[i];
            if (n[i] > 7)
            {
                n[i] = n[i] % 2 + 6;
            }
            sum += n[i] * i;
        }
        if (0 == sum) break;
        memset(f, 0, sizeof(f));
        f[0][0] = 1;
        cout << "Collection #" << cnt << ":" << endl;
        if (sum % 2)
        {
            cout << "Can't be divided." << endl;
            cout << endl;
            continue;
        }
        for (i = 1; i < 7; ++i)
        {
            for (j = 0; j < 150; ++j)
            {
                if (f[i - 1][j])
                {
                    for (k = 0; k <= n[i]; ++k)
                    {
                        f[i][j + k * i] = 1;
                    }
                }
            }
        }
        if (f[6][sum / 2])
        {
            cout << "Can be divided." << endl;
        }
        else
        {
            cout << "Can't be divided." << endl;
        }
        cout << endl;
    }
    return 0;
}