状压dp基础题目

话痨fw的碎碎念

一直觉得dp是玄学
原来状压dp才是 真的玄学
fw落泪

状压dp个人理解，就是把某些状态用01串(看作二进制)表示出来，再将其压缩为十进制表示,提到01就少不了神奇的位运算 了。

位运算 来自大佬的博客
状压dp 依旧来自大佬.

从上面大佬的博客扣来的灰常有用的位运算神奇操作

题目

P1879 [USACO06NOV]Corn Fields G

大佬的题解.
思路:
预处理:
(1)将每行是否可以种草转化为十进制.
(2)找出只考虑每行时可以种草的十进制表示
遍历每行及每种可以种草的状态，当符合要求时找到上一行与之不矛盾的状态进行累加。

题目链接.

#include<bits/stdc++.h>
#define ll long long
using namespace std;

ll mod=100000000;
ll n,m;
int mp[15];
ll dp[15][1<<13];
int ad[2<<13],cnt;

void init()
{
    for(int i=0;i<(1<<m);i++)
    {
        if((!(i&(i<<1)))&&(!(i&(i>>1))))ad[cnt++]=i;
    }
}

int main()
{
    ios::sync_with_stdio(0);
    //cin.tie(0);cout.tie(0);

    cin>>n>>m;
    init();
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            int x;cin>>x;
            mp[i]=(mp[i]<<1)+x;
        }
    }

    dp[0][0]=1;
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<cnt;j++)
        {
            if((mp[i]&ad[j])==ad[j])
            {
                for(int k=0;k<cnt;k++)
                {
                    if((ad[j]&ad[k])==0)
                    {
                        dp[i][j]+=dp[i-1][k];
                        dp[i][j]%=mod;
                    }
                }
            }

        }
    }
    ll ans=0;
    for(int i=0;i<cnt;i++)
    {
        ans=(ans+dp[n][i])%mod;
    }
    cout<<ans<<"\n";
    return 0;
}

Mondriaan’s Dream

题目链接.
思路:将每行二进制表示出后压缩为十进制,根据本行的放置对下行的影响依次递推。

#include<iostream>
#include<cstring>
#define ll long long
using namespace std;

ll dp[15][1<<12];
int n,m;

void dfs(int i,int j,int now,int ne)
{
    if(j==m)
    {
        dp[i+1][ne]+=dp[i][now];
        return ;
    }
    
    //该位置已被上行填充
    if(((1<<j)&now)>0)
        dfs(i,j+1,now,ne);
    //该位置未被上行填充
    if(((1<<j)&now)==0)
        dfs(i,j+1,now,ne|(1<<j));
     //该位置及其下一个位置均未未被上行填充
    if(((j+1)<m)&&(((1<<j)&now)==0)&&(((1<<(j+1))&now)==0))
        dfs(i,j+2,now,ne);

}

int main()
{
    while(cin>>n>>m&&(n||m))
    {
        memset(dp,0,sizeof dp);
        dp[1][0]=1;
        for(int i=1;i<=n;i++)
        {
            for(int j=0;j<(1<<m);j++)
            {
                if(dp[i][j])dfs(i,0,j,0);
            }
        }
        cout<<dp[n+1][0]<<"\n";
    }
    return 0;
}

三进制状压 HDU - 3001

题目链接.

#include<bits/stdc++.h>
#define ll long long
#define INF 0x3f3f3f3f
using namespace std;

int n,m;
int mp[12][12],dp[80000][12];
int a[12],b[80000][12];

void init()
{
    a[1]=1;
    for(int i=2;i<=11;i++)a[i]=a[i-1]*3;
    for(int i=0;i<a[11];i++)
    {
        int x=i;
        for(int j=1;j<=10;j++)
        {
            b[i][j]=x%3;
            x=x/3;
        }
    }
}

int main()
{
    ios::sync_with_stdio(0);
    //cin.tie(0);cout.tie(0);
    init();

    while(cin>>n>>m)
    {
        memset(dp,INF,sizeof dp);
        memset(mp,INF,sizeof mp);
        while(m--)
        {
            int a,b,c;cin>>a>>b>>c;
            if(c<mp[a][b])
                mp[a][b]=mp[b][a]=c;
        }
        for(int i=1;i<=n;i++)dp[a[i]][i]=0;

        int ans=INF;
        for(int i=0;i<a[n+1];i++)
        {
            int flag=1;
            for(int j=1;j<=n;j++)
            {
                if(b[i][j]==0)flag=0;
                if(dp[i][j]==INF)continue;
                for(int k=1;k<=n;k++)
                {
                    if(j!=k&&mp[j][k]<INF&&b[i][k]<2)
                    {
                        dp[i+a[k]][k]=min(dp[i+a[k]][k],dp[i][j]+mp[j][k]);
                    }
                }
            }
            if(flag)
            {
                for(int j=1;j<=n;j++)
                    ans=min(ans,dp[i][j]);
            }
        }
        if(ans==INF)cout<<"-1\n";
        else cout<<ans<<"\n";
    }

    return 0;
}