P1879 Corn Fields G （状压dp

stay fool

已于 2023-10-05 17:34:31 修改

阅读量119

点赞数

分类专栏： dp 文章标签：算法数据结构

于 2023-07-20 13:40:20 首次发布

本文链接：https://blog.csdn.net/m0_75125820/article/details/131828497

版权

dp 专栏收录该内容

89 篇文章 0 订阅

订阅专栏

该代码示例是用C++编写的，主要涉及位操作和动态规划。程序首先定义了一些数据类型别名，然后处理输入，构建状态数组。通过检查相邻位是否都是0来确定状态的有效性。接下来，通过ff函数计算两个状态的最小公共祖先，并在grass_st函数中判断草丛（grass）和状态（st）是否兼容。最后，使用动态规划方法（f[i][x]）计算满足条件的状态数量。

摘要由CSDN通过智能技术生成

#include <bits/stdc++.h>
using namespace std;
using VI = vector<int>;
using PII = pair<int,int>;
using ll = long long;
using ull = unsigned long long;
int n,m;
int f[20][1<<13];
int mp[110];
int cnt[1<<13];
const int mod = 100000000;
bool check(int x){
    //当前状态没有相邻的1
    if(!(x & x>>1)){
        return true;
    }
    return false;

}
int ff(int x,int y){

    while(x && y){
        if((x & 1) == 1 && (y & 1) == 1) return 0;
        x = x>>1;
        y = y>>1;
    }
    return 1;
}
int grass_st(int grass ,int st){
    while(grass || st){
        if((grass & 1) ==0 && (st & 1) == 1) return 0;
        grass = grass >>1;
        st = st>>1;
    }
    return 1;

}
int main(){
    cin>>m>>n;
    for(int i=1;i<=m;i++){
        int st = 0;
        for(int j=0;j<n;j++){
            int x;
            cin>>x;
            if(x) st |= (1<<j);
        }
        mp[i] = st;
    }
    mp[m+1] = 0;
    VI state;
    for(int i=0;i<=(1<<n)-1;i++){
        if(check((i))) state.push_back(i);
    }
    VI head[1<<13];
    for(auto x :state){
        for(auto y : state){
            if(ff(x,y)) {
                head[x].push_back(y);
                //cout << x << " " << y << endl;
            }
        }
    }
    VI grass[1<<13];
    for(int i=1;i<=m+1;i++){
        for(auto y : state){
            if(grass_st(mp[i],y)){
                grass[i].push_back(y);
                //cout<<y<<" ";
            }
        }
        //cout<<endl;
    }


    memset(f,0,sizeof f);
    f[0][0] = 1;

    for(int i=1;i<=m+1;i++){
        for(auto x : grass[i])
        {
            {
                for(auto y : head[x])
                {
                    f[i][x] = (f[i-1][y] + f[i][x]) % mod;
                }
            }
        }
    }

    cout<<f[m+1][0] % mod;
}

感觉思维上好像难度不大，主要就是这个预处理合理状态得仔细想想