子集DP

---

1.子集DP什么鬼

今天第一次见到这个名字，哇，新鲜算法。。。

子集dp 英文名:subset DP

第一次见，还不能总结，看例题吧。

2.例题

原题：Codeforces Round #585 (Div. 2) E. Marbles
简洁中文版：H-树叶

分析（以中文版题面为依据）

（以下是我自己的中二分析）
对于每片叶子，只关心别的叶子是否与它同色，而不在乎具体的颜色。
而20这个数太小了，一定有猫腻，莫非要状压dp？可是这个dp如何表示呢？状压只有0和1，怎么也表示不了啊。

（以下是正经题解版）

• $20$这个数太小了，因此可以用指数级的算法。
• 我们先求一个数组$exc[i][j]$，意为忽略掉其他颜色，将颜色$i$和$j$摆成 “$i$在前，$j$在后” 的位置关系时需要的交换次数（$exc$即exchange交换）。
• 之后状压$s$，代表某些颜色已经各自连在一起的状态（忽略其他颜色）；$f[s]$为到达状态$s$时至少需要的交换数（忽略其他颜色）。这样，便可由$f[s]$转移到$f[s|1<<i],0<=i<20$。
  $f[s|1<<i]=f[s]+{\sum }_{s包含的颜色j}exc[i][j]$
• O(2^20)

code

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
using namespace std;
typedef long long LL;
const int N=4e5+10;
const LL INF=1e15;

int n;
vector<int> v[20];
LL exc[20][20], f[1<<20];

int main()
{
    ///input
    scanf("%d", &n);
    for(int i=1; i<=n; ++i)
    {
        int x;
        scanf("%d", &x);
        v[x-1].push_back(i);
    }

    ///exc j,i->i,j
    for(int i=0; i<20; ++i)
    {
        if(!v[i].size()) continue;
        int szi=v[i].size();
        for(int j=0; j<20; ++j)
        {
            if(i==j||!v[j].size()) continue;
            LL sum=0; int wj=-1, szj=v[j].size();
            for(int wi=0; wi<szi; ++wi)
            {
                while(wj+1<szj&&v[j][wj+1]<v[i][wi]) wj++;
                sum+=wj+1;
            }
            exc[i][j]=sum;
        }
    }

    ///f
    int stu=1<<20;
    fill(f,f+stu,INF); f[0]=0;
    for(int i=0; i<stu; ++i)
    {

        bool ok=true;
        for(int j=0; j<20; ++j)
            if((i&(1<<j))&&!v[j].size()) ok=false;
        if(!ok) continue;
        for(int j=0; j<20; ++j)
        {
            if((i&(1<<j))||!v[j].size()) continue;
            LL sum=0;
            for(int k=0; k<20; ++k)
            {
                if(!(i&(1<<k))) continue;
                sum+=exc[j][k];
            }
            f[i|(1<<j)]=min(f[i|(1<<j)],f[i]+sum);
        }
    }
    ///output
    int ans=0;
    for(int i=0; i<20; ++i)
        if(v[i].size()) ans|=1<<i;
    cout<<f[ans]<<endl;

    return 0;
}

2020.1.31