hdu5691状压dp

传送门

• 状态转移 dp[s|(1<<(j-1))][j]=max(dp[s|(1<<(j-1))][j],dp[s][i]+ai*aj)
• 有虚假态

#include<bits/stdc++.h>
using namespace std;
#define inf 0x80000000
int dp[1<<20][20];//dp[i][j],i为最大1<<n的二进制数，表示这n个数的当前选择方案，j表示第j个数为结尾 
int a[20],p[20];//a[i]标记数值，p[i]标记位置 
int main()
{
	int t,n,s,ans;
	int cas=1;
	scanf("%d",&t);
	while(t--){
		scanf("%d",&n);
		for(int i=0;i<n;i++) scanf("%d%d",&a[i],&p[i]);
		for(int i=0;i<(1<<n);i++)
			for(int j=0;j<n;j++)
				dp[i][j]=inf;
		for(int i=0;i<n;i++)
			if(p[i]==0||p[i]==-1) dp[1<<i][i]=0;
		for(int i=0;i<(1<<n);i++)
			for(int j=0;j<n;j++)
				if(dp[i][j]!=-inf)
					for(int k=0;k<n;k++)
						if(((1<<k)&i)==0&&(p[k]==__builtin_popcount(i)||p[k]==-1))
							dp[i|(1<<k)][k]=max(dp[i|(1<<k)][k],dp[i][j]+a[j]*a[k]);
		ans=inf;
		for(int i=0;i<n;i++) ans=max(ans,dp[(1<<n)-1][i]);
		printf("Case #%d:\n%d\n",cas++,ans);
	}
	return 0;
}

• 无虚假态

#include<bits/stdc++.h>
using namespace std;
const int inf=0x80000000; 
const int maxn=(1<<16)+5;
int n,a[20],p[20];
int h[20];//位置0，1，2……上预存放的数的下标
int c[maxn];//a的下标与二进制位的映射，例c[0001]=1,c[0010]=2,c[0100]=3 
int dp[maxn][20];
//自定义函数 
int bit(int i){
	int cnt=0;
	while(i) i&=i-1,cnt++;
	return cnt;
}//内建函数__builtin_popcount(s); 
int main()
{
	int t,ans,S;
	int cas=1;
	for(int i=0;i<16;i++) c[1<<i]=i+1;//初始化c
	scanf("%d",&t);
	while(t--){
		scanf("%d",&n);
		S=(1<<n)-1;
		memset(h,-1,sizeof(h));
		for(int i=1;i<=n;i++){
			scanf("%d%d",&a[i],&p[i]);
			if(~p[i]) h[p[i]]=i;
		}
		for(int i=0;i<maxn;i++)
			for(int j=0;j<20;j++)
				dp[i][j]=inf;
		if(~h[0]) dp[1<<(h[0]-1)][h[0]]=0;//0位置有预存放的数 
		else{
			for(int i=1;i<=n;i++) if(p[i]==-1) dp[1<<(i-1)][i]=0;
		}
		for(int nxt,i=1;i<S;i++){
			nxt=bit(i);
			for(int lst,j=i;j>0;j&=j-1){//枚举尾元素0101010 & (0101010-1) 
				lst=-j&j;//最低位1,0110->0010
				if(dp[i][c[lst]]==inf) continue;//以元素a[c[lst]]结尾的状态未出现过
				if(~h[nxt]){
					dp[i|(1<<(h[nxt]-1))][h[nxt]]=max(dp[i|(1<<(h[nxt]-1))][h[nxt]],dp[i][c[lst]]+a[c[lst]]*a[h[nxt]]);
					continue;
				} 
				for(int k=1;k<=n;k++){//下一个位置没有预排的数 
					if(1<<(k-1)&i) continue;//a[k]已经排好
					if(~p[k]&&p[k]!=nxt) continue;//a[k]预排到其他位置 
					dp[i|(1<<(k-1))][k]=max(dp[i|(1<<(k-1))][k],dp[i][c[lst]]+a[c[lst]]*a[k]);
				}
			}
		}
		ans=inf;
		for(int i=1;i<=n;i++) ans=max(ans,dp[S][i]);
		printf("Case #%d:\n%d\n",cas++,ans);
	} 
	return 0;
}