【每日一题】德玛西亚万岁

德玛西亚万岁

https://ac.nowcoder.com/acm/problem/15034

solution

状压dp。

用数字来表示一行方案,对于一个状态x,如果x的第i位(二进制下)为1,表示这一行的第i个位置安排士兵。否则表示不安排士兵。

先把每一行可能的状态都搜出来,也就是在满足相邻两个二进制位不都为1的条件下,搜出所有的方案。可以发现一行的方案数最多不超过400。

然后将两个可以转移的状态之间两两连边。x可以转移到y表示的意思是x这个状态可以和y这个状态相邻。判断两个状态之间能不能转移显然就是看x和y是不是有某个二进制位同时为1,即只有时,x可以转移到y。

表示前i行,第行的状态是y的方案数,注意到,第i行的状态可以是y的前提是y这个状态与所给地形之间不存在某个二进制位置i满足(b表示第i行的地形状态)。。

考虑转移就从所有的转移过来就行了。转移的条件是j这个状态与k这个状态之间有边。

code

/*
* @Author: wxyww
* @Date:   2020-06-02 19:42:43
* @Last Modified time: 2020-06-02 19:54:25
*/
#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include<ctime>
#include<cmath>
using namespace std;
typedef long long ll;
const int N = 1010,mod = 100000000;
ll read() {
    ll x = 0,f = 1;char c = getchar();
    while(c < '0' || c > '9') {
        if(c == '-') f = -1; c = getchar();
    }
    while(c >= '0' && c <= '9') {
        x = x * 10 + c - '0'; c = getchar();
    }
    return x * f;
}
int a[N],ma[N][N],f[20][N],n,m,b[20];

void solve() {    
    int cnt = 0;
    for(int i = 0;i < (1 << m);++i) {
        int flag = 0;
        for(int j = 1;j < m;++j) {
            if((i >> j & 1) && (i >> (j - 1) & 1)) {flag = 1;break;}
        }
        if(!flag) a[++cnt] = i;
    }

    memset(b,0,sizeof(b));
    for(int i = 1;i <= n;++i) {
        for(int j = 0;j < m;++j) {
            int x = read() ^ 1;
            b[i] |= x << j;
        }
    }

    memset(ma,0,sizeof(ma));
    for(int i = 1;i <= cnt;++i) {
        for(int j = i;j <= cnt;++j) {
            if(!(a[i] & a[j])) ma[i][j] = ma[j][i] = 1;
        }
    }
    memset(f,0,sizeof(f));
    f[0][1] = 1;
    // cout<<ma[1][2]<<endl;
    for(int i = 1;i <= n;++i) {
        for(int j = 1;j <= cnt;++j) {
            if(a[j] & b[i]) continue;
            // if(i == 1 && j == 1) puts("!!!");
            for(int k = 1;k <= cnt;++k) {
                if(ma[k][j]) {
                    // if(k == 1 && i == 1) puts("!!!");
                    f[i][j] += f[i - 1][k],f[i][j] %= mod;
                }
            }
        }
    }
    ll ans = 0;
    for(int i = 1;i <= cnt;++i) 
        ans += f[n][i],ans %= mod;
    cout<<ans<<endl;
}

int main() {
    while(~scanf("%d%d",&n,&m))
    solve();
    return 0;
}
全部评论

相关推荐

昨天 17:22
已编辑
西安交通大学 Java
华为 昇腾 ai软件开发 薪资20k x (14-16),职级13A,5%公积金,c/cpp
BLOOMING7:闭眼滴滴,华子给的又少又累
点赞 评论 收藏
分享
oppo 应用软开 22*15+0.5*12
拿到了ssp完美:真的坎坷,但是你至少拿到这么多offer了!
点赞 评论 收藏
分享
沉淀一会:**圣经 1.同学你面试评价不错,概率很大,请耐心等待;2.你的排名比较靠前,不要担心,耐心等待;3.问题不大,正在审批,不要着急签其他公司,等等我们!4.预计9月中下旬,安心过节;5.下周会有结果,请耐心等待下;6.可能国庆节前后,一有结果我马上通知你;7.预计10月中旬,再坚持一下;8.正在走流程,就这两天了;9.同学,结果我也不知道,你如果查到了也告诉我一声;10.同学你出线不明朗,建议签其他公司保底!11.同学你找了哪些公司,我也在找工作。
点赞 评论 收藏
分享
评论
2
收藏
分享
牛客网
牛客企业服务