这周开始新的专题开启了,前几天还在看博客,得先了解状压DP到底解决的是什么问题,主要还是看对于自己来说新型的题目,例如H题,题目的意思是研究人员发现有N个原子两两组合会发生反应产生能量,且其中一个会被和谐掉,给出一个矩阵map,其中map[i][j]表示原子i与j发生反应且原子j被和谐掉所释放的能量,问给定的n个原子反应最多能产生多少能量。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int n;
int dp[1<<11];
int power[11][11];
int main(){
while(cin>>n){
if(n==0)break;
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
cin>>power[i][j];
memset(dp,0,sizeof(dp));
for(int s=(1<<n)-1;s>=0;s--){
for(int i=0;i<n;i++)
if(s&(1<<i)){
for(int j=0;j<n;j++){
if(s&(1<<j)&&i!=j)
dp[s^(1<<j)]=max(dp[s^(1<<j)],dp[s]+power[i][j]);
}
}
}
int maxx=0;
for(int i=0;i<(1<<n)-1;i++)
maxx=max(maxx,dp[i]);
cout<<maxx<<endl;
}
}
还有A,这道题的大概意思是给出一个n*m的棋盘,问用1*2的多米诺骨牌将其完全覆盖有多少种不同的方案。
分析:经典的状态压缩动态规划。n和m最大都只有11,所以可以采用状态压缩,即用一个二进制数表示每一行的状态。定义f[i,j]表示前i-1行已经完全覆盖,第i行的状态为j的时候的方案数。则f[i,j]=∑f[i-1,k]。k能通过放置骨牌使第i-1行被完全覆盖,且第i行状态为j。这个状态转移可以用dfs来实现。在开始动态规划之前,要先用dfs预处理出第一行可能的状态。对于第一行,只有两种可能操作:
1、不放。-->dfs(x+1,now shl 1);
2、横放。-->dfs(x+2,now shl 2 or 3);
当x等于m+1时,说明得到了一种可行状态,inc(f[1,now]);当x大于m时,退出dfs。
状态转移时的dfs过程与此类似,只是多了一种竖放操作而已,运用位运算实现 。