1. 程式人生 > >[luoguP3694] 邦邦的大合唱站隊/簽到題(狀壓DP)

[luoguP3694] 邦邦的大合唱站隊/簽到題(狀壓DP)

pac sizeof size 前綴和 for ensure -m open https

傳送門

來自kkk的題解:

70分做法:枚舉每個學校順序,暴力。

100分:狀壓dp。從隊列頭到尾DP,

狀態:f[i]表示i狀態下最小的出列(不一致)的個數。

比如f[1101]表示從頭到位為1/3/4樂隊的偶像的最小出列個數。

f[i]=min(f[i\ xor\ 2^j]+num[j]-(sum[length][j]-sum[length-num[j]][j]));f[i]=min(f[i xor 2?j??]+num[j](sum[length][j]sum[lengthnum[j]][j]));

j表示團隊編號,sum表示某種團隊的前綴和,length表示到此已經排到的長度。

#include <cstdio>
#include <cstring>
#include <iostream>
#define M 21
#define N 100001
#define min(x, y) ((x) < (y) ? (x) : (y))

int n, m, S, cnt;
int sum[N][M], num[M], f[1 << M];

inline int read()
{
	int x = 0, f = 1;
	char ch = getchar();
	for(; !isdigit(ch); ch = getchar()) if(ch == ‘-‘) f = -1;
	for(; isdigit(ch); ch = getchar()) x = (x << 1) + (x << 3) + ch - ‘0‘;
	return x * f;
}

int main()
{
	int i, j, x;
	n = read();
	m = read();
	memset(f, 127, sizeof(f));
	for(i = 1; i <= n; i++)
	{
		x = read();
		num[x]++;
		for(j = 1; j <= m; j++)
			sum[i][j] = sum[i - 1][j] + (x == j);
	}
	f[0] = 0;
	for(i = 1; i < (1 << m); i++)
	{
		S = 0;
		for(j = 1; j <= m; j++)
			if(i & (1 << j - 1))
				S += num[j];
		for(j = 1; j <= m; j++)
			if(i & (1 << j - 1))
			{
				cnt = num[j] - sum[S][j] + sum[S - num[j]][j];				
				f[i] = min(f[i], f[i ^ (1 << j - 1)] + cnt);
			}
	}
	printf("%d\n", f[(1 << m) - 1]);
	return 0;
}

  

[luoguP3694] 邦邦的大合唱站隊/簽到題(狀壓DP)