BZOJ 1231: [Usaco2008 Nov]mixup2 混乱的奶牛

状压dp入门题。

写状压dp的时候，其实有很多小优化技巧，最方便的也是最没啥用的一个（不过也有点用，上次就有一次因为没优化被卡常），就是数组定义的时候尽量定义为0—n-1，这样的话就不用1<<(i-1)了，直接1<<i即可。

如果把加了这个小优化的话，bzoj快20ms，洛谷快9ms。

还有一个更没啥用的小优化叫做寻址优化，就是f数组，把小的N放在大的1<<N前面，听大佬们说能快一点儿。

#include <bits/stdc++.h>
#define int long long
using namespace std;
const int N=16;
int n,k,max_statue,ans;
int a[N],f[N][1<<N];
//f[p][i]表示：已选的奶牛状态为p，最后一头奶牛为i的方案数 
//状态是显然的，因为转移过程和最后一头奶牛有关所以必须记录一个最后一头奶牛是什么的数组
//转移过程也是很显然的
//初始化想了一会。 
signed main(){
	scanf("%lld%lld",&n,&k);
	for (register int i=0; i<n; ++i) scanf("%lld",&a[i]); 
	
	max_statue=(1<<n)-1;
	for (register int i=0; i<n; ++i) f[i][1<<i]=1;
	for (register int p=1; p<max_statue; ++p)
	{
		for (register int i=0; i<n; ++i)
		if ((1<<i)&p)
		{
			for (register int j=0; j<n; ++j)
			if (abs(a[i]-a[j])>k)
			{
				if (!((1<<j)&p))
				{
					f[j][p|(1<<j)]+=f[i][p];
				}
			}
		}
	}
	
	for (register int i=0; i<n; ++i) ans+=f[i][max_statue];
	printf("%lld\n",ans);
return 0;
}