P3694 邦邦的大合唱站队（集合类状压dp）

 邦邦的大合唱站队 - 洛谷

关键：由于是线性安排的，通过枚举所有M的选择情况+最后一个不同点，可以达到枚举M的所有摆放顺序

分析：由于M的范围是1~20,所以想到用状压dp来做

由于按照出队后，出队的偶像可以归队到任意的空位中，且我们需要同一乐队的偶像连续站在一起

• 性质1：一队偶像安排到的位置可以是任意的,在他之前和之后的偶像都可以是任意的

因为位置可能是任意的，那么我们可以先从头开始安排，一直安排到尾。

因为我们是使用状压dp来做的，通过枚举所有  2^j   （0<=j<=M），就可以达到枚举排队的顺序的效果（例如 1011，可以是 134，143，341...，通过不断枚举最后一个分割点可以达到左边的所有状态  参考：最短Hamilton路径）

• 性质2：通过枚举最后一次不同点，可以枚举M个团队所有的安排顺序

因为要排一队团队的时候，由于我们已经确定好了他的位置，设队伍长度为 num[j],当前已经排到的队伍总长度(包括队伍j)为lenth。那么如图当前要出队的人是  框内所有不是该团队的偶像 +  框外所有是该团队的偶像 

由于当前框内的偶像可能是其他团内的偶像，所以可能重复计算，那么这里只计算 框内所有不是该团队的偶像 ，由于从头安排到尾，所以不在框内的该团队的偶像会在其他团队计算时被计算到。

• 性质3：排列当前团队时  只计算  不是该团队偶像的个数

因为框内实际上是一个区间，那么就需要O(1)求区间和，想到前缀和

• 性质4：前缀和优化，一个区间内所有该偶像的数量=s[r][i]-s[l][i]

下面是dp的分析图

 状态表示：当前所有团队的安排状态是i，且最后一个摆好的团队是j的所有集合

属性：min

初始化：0x3f                     边界：dp[0][0] =0 全部都没有摆

状态计算/子集划分：根据上一个子集中最后一个到达的点作为划分依据，可以分为不重不漏的多个子集。

• dp[i][j]=min(dp[i][j],dp[i^(1<<j)][k]+num[j]-(s[lenth][j]-s[lenth-num[j][j])

由于会mle，所以可以进行空间优化，只要把第二维去掉就可以了

• dp[i]=min(dp[i],dp[i^(1<<j)]+num[j]-(s[lenth]-s[lenth-num[j]])

#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
const int N = 1e5 + 10, M = 21;
int num[N];
int s[N][M];
int dp[1 << M];
int main()
{
	int n, m;
	cin >> n >> m;
	for (int i = 1;i <= n;i++)
	{
		int a;
		cin >> a;a--;
		s[i][a]++;
		num[a]++;
		for (int j = 0;j < m;j++)
			s[i][j] += s[i - 1][j];
	}

	memset(dp, 0x3f, sizeof dp);
	dp[0] = 0;

	for(int i=1;i<1<<m;i++)
		for (int j = 0;j < m;j++)
			if (i >> j & 1)
			{
				int lenth = 0;
				
				for (int k = 0;k < m;k++)
					if (i >> k & 1)
						lenth += num[k];
				//cout << lenth << ' ';
						dp[i] = min(dp[i], dp[i ^ (1 << j)] + num[j] - (s[lenth][j] - s[lenth - num[j]][j]));
				
			}
	int res = 0x3f3f3f;
	for (int i = 0;i < m;i++)
		res = min(res, dp[(1 << m) - 1]);
	cout << res;
}