36 状态压缩DP 树形DP 记忆化搜索DP 数位统计DP

1 状态压缩DP

1.1 蒙德里安的梦想

/*
学acwing的算法基础课学来的，喜欢的话多多支持呀。
*/
//状态压缩DP 蒙德里安的梦想  N×M的棋盘分割成若干个1×2的的长方形，有多少种方案
#include<bits/stdc++.h>
using namespace std;
const int N=12,M=1<<N;
int n,m;
long long f[N][M];//第一维表示列,第二维表示所有可能的状态
bool st[M];//存储每种状态是否有奇数个连续的0，如果奇数个0是无效状态，如果是偶数个零置为true
int main()
{
   while(cin>>n>>m,n||m)
    {
        memset(f,0,sizeof f);//全部初始化为0，因为是连续读入，这里是一个清空操作
        for(int i=0;i<1<<n;i++)//对于每种状态，先预处理每列不能有奇数个连续的0
        {
            st[i]=true;// 某种状态没有奇数个连续的0则标记为true
            int cnt=0;//记录连续的0的个数
            for(int j=0;j<n;j++)//遍历这一列，从上到下
                if(i>>j&1)//i>>j位运算，表示i（i在此处是一种状态）的二进制数的第j位； &1为判断该位是否为1，如果为1进入if
                {
                   if(cnt&1) st[i]=false;//这一位为1，看前面连续的0的个数，如果是奇数（cnt &1为真）则该状态不合法
                    cnt=0;//然后清0
                }
                else cnt++;//否则的话该位还是0，则统计连续0的计数器++
        if(cnt&1) st[i]=false;//假如最下面的那一段判断一下连续0的个数
        }
    }
    f[0][0]=1;// 这里需要回忆状态表示的定义，按定义这里是：前第-1列都摆好，且从-1列到第0列伸出来的状态为0的方案数
    for(int i=1;i<=m;i++) //遍历每一列:第i列合法范围是(0~m-1列)
        for(int j=0;j<1<<n;j++)//遍历当前列（第i列）所有状态j
          for(int k=0;k<1<<n;k++)  // 遍历第i-1列的状态k，如果“真正”可行，就转移
            if(!(j&k)&&st[j|k])
               f[i][j]+=f[i-1][k];  // 当前列的方案数就等于之前的第i-1列所有状态k的累加
    cout<<f[m][0]<<endl; //f[m][0]表示 前m-1列都处理完，并且第m-1列没有伸出来的所有方案数,也就是总方案数
    return 0;
}
/*到达胜利之前无法回头*/

1.2 最短Hamilton路径

/*
学acwing的算法基础课学来的，喜欢的话多多支持呀。
*/
//状态压缩DP 最短Hamilton路径 给定一张n个点的带权无向图，点从0∼n−1标号求起点0到终点n−1的最短Hamilton路径
//Hamilton路径的定义是从0到n−1不重不漏地经过每个点恰好一次
#include<bits/stdc++.h>
using namespace std;
const int N=20,M=1<<N;
int n;
int w[N][N];//存i到j点的距离
int f[M][N];//所以从0走到j，走过的所有点是i的所以路径
int main()
{
   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这个1个点，从0走到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;//输出结果
    return 0;
}
/*到达胜利之前无法回头*/

2 树形DP

/*
学acwing的算法基础课学来的，喜欢的话多多支持呀。
*/
//树形DP 没有上司的舞会
/*Ural 大学有 N 名职员，编号为 1∼N。
他们的关系就像一棵以校长为根的树，父节点就是子节点的直接上司。
每个职员有一个快乐指数，用整数 Hi 给出，其中 1≤i≤N。
现在要召开一场周年庆宴会，不过，没有职员愿意和直接上司一起参会。
在满足这个条件的前提下，主办方希望邀请一部分职员参会，使得所有参会职员的快乐指数总和最大，求这个最大值
*/ 
#include<bits/stdc++.h>
using namespace std;
const int N=6010;
int n;
int happy[N];//存的是每个点的开心度
int h[N],e[N],ne[N],idx;//领接表的存储方式
int f[N][2];//存状态
bool has_father[N];//存上司
void add(int a,int b)//链接a->b
{
    e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void dfs(int u)
{
    f[u][1]=happy[u];//选这个点
    for(int i=h[u];i!=-1;i=ne[i])//遍历这个点的所有
    {
        int j=e[i];
        dfs(j);//在遍历上司的上司
        //状态计算得来的
        f[u][0]+=max(f[j][0],f[j][1]);
        f[u][1]+=f[j][0];
    }
}
int main()
{
   cin>>n;
   for(int i=1;i<=n;i++) cin>>happy[i];
   memset(h,-1,sizeof h);//初始化为空结点
   for(int i=0;i<n-1;i++)
    {
        int a,b;
        cin>>a>>b;
        has_father[a]=true;//标记父节点
        add(b,a);
    }
    int root=1;
    while(has_father[root]) root++;//找一个没有下司的点
    dfs(root);
    printf("%d\n",max(f[root][0],f[root][1]));//输出选或者不选的最大值
    return 0;
}
/*到达胜利之前无法回头*/

3 记忆化搜索DP 

/*
学acwing的算法基础课学来的，喜欢的话多多支持呀。
*/
//记忆化搜索DP 滑雪
#include<bits/stdc++.h>
using namespace std;
const int N=310;
int n,m;
int h[N][N];//存图
int f[N][N];//存状态
int dx[4]={-1,0,1,0},dy[4]={0,1,0,-1};//定义上下左右
int dp(int x,int y)
{
    int &v=f[x][y];//这里取f[x][y]的地址，相当于v=f[x][y]
    if(v!=-1) return v;//假如被求过则直接返回，这里则为记忆化搜索
    v=1;//令f[x][y]=1
    for(int i=0;i<4;i++)//枚举一下每一步
    {
        int a=x+dx[i],b=y+dy[i];//下一步方位
        if(a>=1&&a<=n&&b>=1&&b<=m&&h[a][b]<h[x][y])//假如合法并且满足条件
            v=max(v,dp(a,b)+1);//则走下一步，每次取最大
    }
    return v;//返回f[x][y],也就是每个点的最大步
}
int main()
{
   cin>>n>>m;
   for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
         cin>>h[i][j];
   memset(f,-1,sizeof f);//初始化为-1
   int res=0;//用来记录答案
   for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
      res=max(res,dp(i,j));//求一遍每个点的值，再取最大
   cout<<res<<endl;
    return 0;
}
/*到达胜利之前无法回头*/

4 数位统计DP

 

/*
学acwing的算法基础课学来的，喜欢的话多多支持呀。
*/
//数位统计DP 计数问题
#include<bits/stdc++.h>
using namespace std;
int get(vector<int>num,int l,int r)//用来获取l到r之间位数的值
{
    int res=0;
    for(int i=l;i>=r;i--) res=res*10+num[i];
    return res;
}
int count(int n,int x)
{
    if(!n) return 0;
    vector<int>num;
    while(n)//获取n的每一位
    {
        num.push_back(n%10);
        n/=10;
    }
    n=num.size();
    int res=0;
    for(int i=n-1-!x;i>=0;i--)
    {
        if(i<n-1)//假如第i位小于n-1
        {
            res+=get(num,n-1,i+1)*pow(10,i);//则第一种情况
            if(!x) res-=pow(10,i);//假如只是0的话，就减多一次
        }
        if (num[i]==x) res+=get(num,i-1,0)+1;//第二的第二种情况
        else if (num[i]>x)res+=pow(10,i);//第二的第三种情况
    }
    return res;//返回最终答案
}
int main()
{
    int a,b;
    while(cin>>a>>b,a)
    {
        if (a>b) swap(a,b);
        for (int i=0;i<=9;i++)//枚举每一位
        cout<<count(b,i)-count(a-1,i)<<' ';//前缀和的方式求出
        puts("");
    }
    return 0;
}

/*到达胜利之前无法回头*/