ch4(DP)——区间DP、状态压缩DP

1.自然语言描述
区间DP：一般是给定一个序列，然后进行区间内元素合并，问合并整个序列的最大价值/最小耗费是多少。区间DP问题有着明显的特征，问题的思考方式也有着确定的模板。

状态压缩DP：一般是将问题中的状态，用一个二进制数来表示，将状态的枚举转变为位运算；如果用一般的暴力搜索，复杂度将会是出奇的高；这个时候就考虑能不能将每个状态压缩为一个二进制数，然后通过问题本身的信息来降低复杂度。

2.代码描述
题目：Acwing.282 石子合并 题目链接

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

const int MAXN=310;

int n;
int s[MAXN];
int f[MAXN][MAXN];//f[i][j]:区间[i,j]内所有石子堆合并成一个堆所需要的最小代价

int main(void)
{
    cin>>n;
    
    for(int i=1;i<=n;i++)
        cin>>s[i];
        
    for(int i=1;i<=n;i++)
        s[i]+=s[i-1];
        
    for(int len=2;len<=n;len++)//len表示区间长度,从小到大枚举区间长度使得小区间内的f先被算出来,后续大区间计算f要用
        for(int i=1;i+len-1<=n;i++){//枚举起点
            int l=i,r=i+len-1;
            f[l][r]=1e9;//初始化f以一个很大的数,因为要取最小值
            for(int j=l;j<r;j++)//枚举分界点
                f[l][r]=min(f[l][r],f[i][j]+f[j+1][r]+s[r]-s[l-1]);//状态转移方程
        }
        
    cout<<f[1][n]<<endl;
    
    return 0;
}

题目：Acwing.291 蒙德里安的梦想题目链接

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
typedef long long LL;
const int MAXN=12,MAXM=1<<MAXN;

//一列一列考虑,从1~m列,将每一列的状态抽象为一个n位二进制数,只需要考虑横向长方形填满的情况,然后纵向长方形的填充位置也就确定了
LL f[MAXN][MAXM];//f(i,j):表示第i列中已经被i-1列放的横向长方形占了的j状态情况下放置横向长方形的方案数
                    //j为n位二进制数,用j代表第i列放置横向长方形的方案
                    //设前一列状态f(i-1,k),当j&k==0时说明枚举的j和k不冲突,也就是说没有在第i列已经被前一列的长方形占了的地方放置长方形
                    //j|k的结果二进制数中,只有0的个数为偶数,才能放得下纵向长方形,所以必须预处理出数据范围内所有状态的含0奇偶性,st数组应运而生
int n,m;
bool st[MAXM];//st[i]:表示某行状态i二进制下的连续0是否是偶数个,是则为true
//只有连续0个数为偶数的状态才可取

int main(void)
{
    while(cin>>n>>m,n||m){
        
        memset(f,0,sizeof(f));
        
        //预处理所有状态的含0奇偶状态
        for(int i=0;i<1<<n;i++){
            st[i]=true;
            int cnt=0;//cnt临时记录连续个0的个数
            for(int j=0;j<n;j++)
                if(i>>j&1){//枚举n位,若当前二进制位是1
                    if(cnt&1)//若发现连续0数量为奇数,则说明这个状态不可取
                        st[i]=false;
                    cnt=0;//清空计数
                }else
                    cnt++;
            if(cnt&1)
                st[i]=false;
        }
        
        f[0][0]=1;//规定0列不放横向长方形
        for(int i=1;i<=m;i++)
            for(int j=0;j<1<<n;j++)
                for(int k=0;k<1<<n;k++)
                    if((j&k)==0&&(st[j|k]))
                        f[i][j]+=f[i-1][k];
                        
        cout<<f[m][0]<<endl;//最后一列也就是第m列不能放任何横向长方形
        
    }
    return 0;
}

题目：Acwing.91 最短Hamilton路径题目链接

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

const int MAXN=21,MAXM=1<<MAXN;

int n,f[MAXM][MAXN];//f(i,j):表示从0~j经过的点状态为i的最小路径,i二进制的每一位表示该位对应的点是否经过
int w[MAXN][MAXN];

int main(void)
{
    cin>>n;
    
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            cin>>w[i][j];
            
    memset(f,0x3f,sizeof(f));//初始化
    f[1][0]=0;//开始位于0点,状态i的第0位为1,最后点为0,耗费为0
    for(int i=0;i<1<<n;i++)
        for(int j=0;j<n;j++)
            if(i>>j&1)
                for(int k=0;k<n;k++)
                    if((i-(1<<j))>>k&1)
                        f[i][j]=min(f[i][j],f[i-(1<<j)][k]+w[k][j]);
                        
    cout<<f[(1<<n)-1][n-1]<<endl;//最终所有点都被选,即状态i的所有位全为1,并且最后抵达点j为n-1的状态为最小哈密顿路径
    
    return 0;
}