bzoj2064 分裂

$题目传送门bzoj2064$

---

题目

$2064: 分裂$

$Time Limit: 10 Sec$
$Memory Limit: 64 MB$

$Description$

和久必分，分久必和。。。 题目描述： 中国历史上上分分和和次数非常多。。通读中国历史的WJMZBMR表示毫无压力。同时经常搞OI的他把这个变成了一个数学模型。 假设中国的国土总和是不变的。 每个国家都可以用他的国土面积代替， 又两种可能，一种是两个国家合并为一个，那么新国家的面积为两者之和。 一种是一个国家分裂为两个，那么两个新国家的面积之和为原国家的面积。 WJMZBMR现在知道了很遥远的过去中国的状态，又知道了中国现在的状态，想知道至少要几次操作（分裂和合并各算一次操作），能让中国从当时状态到达现在的状态。

$Input$

第一行一个数$n_{1}$，表示当时的块数，接下来$n_{1}$个数分别表示各块的面积。 第二行一个数$n_{2}$，表示现在的块，接下来$n_{2}$个数分别表示各块的面积。

$Output$

一行一个数表示最小次数。

$Sample Input$

1 6
3 1 2 3

$Sample Output$

2

$HINT$

$对于30\%的数据，n_{1},n_{2}\leq6$
$对于100\%的数据，n_{1},n_{2}\leq10，每个数\leq50$

---

题解

$([A]为bool表达式，A成立时返回1，A不成立时返回0。博主在这声明，下文可能用到。)$
状压DP！
状压DP！！
状压DP！！！
因为$n_{1},n_{2}\leq10$，所以可以想到时空复杂度约为$O(2^{n})$。$(n=n_{1}+n_{2})$
这样子就可以想到状压DP。
关键是怎么压缩呢？
我想一想，想一想，想一想……想到对于一个过去状态的子集和一个现在状态的子集，如果两个子集元素之和相等，那么可以先把过去全部$m_{1}$个状态合并为一个，再把一个分裂为现在全部$m_{2}$个状态，这样的代价是$(m_{1}-1)+(m_{2}-1)=m_{1}+m_{2}-2$。而$\sum_{i=1}^{m}(n_{1}+n_{2}-2)=n-2m$。因为$n$是定值，所以$m$越大越好。
但是我们如何判断一个状态是否合法呢？我们可以发现$s_{1}-s_{2}=0$。那么我们只需要把现在的状态变为负数，子集元素的和为$0$时状态合法。
接着我们需要求出最大的$m$。我们用$f(i)$来表示状态为$i$时的$m$。首先我们可以发现$f(2^{i})=0，f(S）=max\{f(T)+f(S-T)\}(T\subset S)$。那么我们就发现了$O(3^{n})$的真·枚举子集的算法。
但是容易发现无法AC……
接下来我们要优化算法。我们可以发现$f(S)\geq f(T)(T\subset S)$，接着可以推出$f(S)=max\{f(T)+[Sum(S-T)=0]\}(T\subset S)$，进而推出$f(S)=max\{f(T)\}+[Sum(S)=0](T\subset S)$，进而推出$f(S)=max\{f(T)\}+[Sum(S)=0](T=S-2^{i})$(因为这样的$f(T)$包含了所有的较优解)。
这样，我们就成功的写出了$O(n·2^{n})$的假·枚举子集的算法。

---

总结

对于这一类的状压DP，我们需要处理的是将最终状态变为多个合法状态的最大状态个数。当一个状态的DP值大于等于所有子集的DP值，并且$f(2^{i})=0$时，就可以将$O(3^{n})$变为$O(n·2^{n})$的复杂度。但是这种方法慎用，否则可能会有惊喜。

---

标程

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int N=1050005;
int s[N], f[N];
int n, m, t;
int main()
{
    scanf("%d", &n);
    for (int i=1; i<=n; i++)
        scanf("%d", &s[1<<(i-1)]);
    scanf("%d", &m);
    for (int i=n+1; i<=n+m; i++)
    {
        scanf("%d", &s[1<<(i-1)]);
        s[1<<(i-1)]=-s[1<<(i-1)];
    }
    n+=m;
    t=(1<<n)-1;
    for (int i=1; i<=t; i++)
    {
        int tt=i&(-i);
        s[i]=s[tt]+s[i^tt];
        for (int j=1; j<=n; j++)
        {
            tt=(i&(1<<(j-1)));
            if (tt) f[i]=max(f[i], f[i^tt]);
        }
        if (!s[i]) f[i]++;
    }
    printf("%d\n", n-(f[t]<<1));
    return 0;
}

---

友情链接

$bzoj2064题解报告$
（大神比较喜欢绿色）
$lvzelong2014$大神解题报告
若有疑惑，问大神去。

---