luogu P3694 邦邦的大合唱站队

题面传送门
挺入门的状压\(dp\)题。
先枚举状态，然后枚举最后一个合并的是哪一个，状态转移方程:
\(f_{i}=min(f_{i-2^j}+\sum\limits_{k=1}^{n}{1[a_k=j]}-\sum\limits_{z=\sum\limits_{s=1}^{m}(\sum\limits_{s=1}^{n}1[a_s=s])[i\&2^{s-1}]-\sum\limits_{s=1}^{n}1[a_s=j]}^{\sum\limits_{s=1}^{m}\sum\limits_{s=1}^{n}1[a_s=s])[i\&2^{s-1}]}{1[a_z=j]})\)
代码实现:

#include<cstdio>
#define min(a,b) ((a)<(b)?(a):(b))
using namespace std;
int n,m,a[100039],q[39][100039],x,y,ans,tot,pus,dp[1100039];
int main(){
	register int i,j;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++) scanf("%d",&a[i]),q[a[i]][i]++;
	for(i=1;i<=m;i++){
		for(j=1;j<=n;j++) q[i][j]+=q[i][j-1];
	}
	for(i=1;i<=(1<<m)-1;i++){
		ans=0;
		for(j=1;j<=m;j++) if(i&(1<<j-1)) ans+=q[j][n];
		dp[i]=1e9;
		for(j=1;j<=m;j++){
			if(i&(1<<j-1)){
				dp[i]=min(dp[i],dp[i-(1<<j-1)]+q[j][n]-q[j][ans]+q[j][ans-q[j][n]]);
			}
		}
	}
	printf("%d\n",dp[(1<<m)-1]);
}