计蒜客蓝桥杯模拟赛五J. 程序设计：放置守卫

在一张 n 行 m 列的方格地图上放置一些守卫，每个守卫能守护上、左、右三个方向上相邻的方格和自己所在的方格。如下图，红色的方格放置守卫，绿色的方格为该守卫守护的区域。

现在要求在地图上放置若干个守卫，让每个方格至少被一个守卫守护（可以同时被多个守卫守护），但是有些方格上不能放置守卫（这个方格也需要被守护），求出最少需要多少个守卫才能满足条件。

输入格式

第一行输入两个整数 n, m。

接下来输入一个 n×m 的矩阵。矩阵中元素为 0 表示该位置不能放置守卫，为 1 表示该位置能放置守卫。元素之间用空格隔开。

数据约定： 所有数据保证一定有一种方案满足条件。

对于 20% 的数据： 1≤n,m≤5；

对于 50% 的数据： 1≤n≤20,1≤m≤10；

对于100% 的数据： 1≤n∗m≤1000,1≤m≤15。

输出格式

输出最少需要放置的守卫数量。

没有打这个。。不过学校群有人发出来了。做了下，没得交，大概思路应该是对的

做法：

我们做一个dp，假设格子移动到了图中标记为1的地方，我们定义dp[now][sta]为当移动到这个格子，而且(4,3,2,1,0)的状态为sta的二进制表示（顺序也是一样的，1表示已放置，0表示未放置），并且这个分界线以上的全部被放置了的最小放置守卫数，这个分界线是这样的，当转移到（i，j）时，分界线就是(i,0)(i,1)...(i,j)(i-1,j)(i-1,j+1)..(i-1,m-1)转移比较简单，初始化第一个格子的时候，我们知道第一个格子是没有格子的，但是我们把那里当成已被守卫的格子就好了所以就是dp[now][(1 << m+1) - 2] = 0,其他定义为-1。

举个例子吧,下面的数字外面一层表示行列数，里面的表示分界线

0 1 2 3

0  _ _ 3 4

1  0 1 2 _

图1

我们知道我们当前是转移到(1,2)这个位置的，当这个位置可以放守卫，我们选择放一个守卫之后

0 1 2 3

0  _ _ _ 4

1  0 1 2 3

图2

现在转移到下一个格子了，但是刚刚在(1,2)放了守卫，所以(1,1)(1,2)(1,3)(0,2)这四个位置的状态都要变成了1，所以我们的sta需要或一个(1 << (j-1))(1 << (j))(1 << (j+1))如果j-1小于0就不用了，(0,2)这个位置已经到了分界线外。所以如果当前格子上面是没有守卫的，就必须放一个。因为我们定义的dp是分界线上必须全部被守卫。

如果当前格子有守卫，我们可以选择不放，如果不放的话我们需要把sta^(1 << (j + 1))，因为我们可以知道图1的分界线3转移到下一个的图2时，这个位置原本是没有守卫的，所以要置为0。

这上面的转移是一种情况，但是当转移需要换行的时候，转移的式子会有所改变，不过手推很容易推出来。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
#include <time.h>
#include <string>
#include <stack>
#include <set>
#include <map>
#include <iostream>
#include <bitset>
#include <algorithm>
using namespace std;
#define MP make_pair
#define rep(i,n) for(int i = 0; i < (n); i++)
#define per(i,n) for(int i = (n-1); i >= 0; i--)
#define rep1(i,n) for(int i = 1; i <= (n); i++)
#define per1(i,n) for(int i = (n); i > 0; i--)

#define PB push_back
#define mst(a,b) memset((a),(b),sizeof(a))
typedef long long LL;
typedef unsigned long long uLL;
typedef pair<int, int> Pii;
typedef vector<int> Vi;
typedef vector<Pii> Vii;
const int inf = 0x3f3f3f3f;
const LL INF = (1uLL << ) - ;
const LL mod = ;
const int N = 1e5 + ;
const double Pi = acos(-1.0);
const int maxn = ( << ) + ;
const uLL Hashmod = ;
int dp[][maxn];
int maz[][];
int n, m;
inline bool up(int sta, int pos) {
    return (sta & ( << (pos + ))) != ;
}
inline void uMin(int &A, int B) {
    if(B == -)return;
    if(A == -)A = B;
    else A = min(A, B);
}
inline int Set(int sta, int pos) {
    if(pos -  >= )sta |= ( << (pos - ));
    sta |= ( << pos);
    sta |= ( << (pos + ));
    if(pos == m - ) {
        sta <<= ;
        sta %= ( << (m + ));
    }
    return sta;
}
inline int noSet(int sta,int pos){
    if(pos != m - )sta ^= ( << (pos + ));
    else sta <<= ,sta %= ( << (m + ));
    return sta;
}
//inline void show(int x) {
//    putchar('0' + x % 2);
//    if(x > 1)show(x / 2);
//}
int main() {
#ifdef local
    freopen("input", "r", stdin);
    //freopen("w","w",stdout);
#endif
    ios::sync_with_stdio();
    cin.tie();
    cin >> n >> m;
    rep(i, n)rep(j, m)cin >> maz[i][j];
    int nx = , now = ;
    mst(dp[nx], -);
    dp[nx][( << (m + )) - ] = ;
    rep(i, n) {
        rep(j, m) {
            swap(nx, now);
            mst(dp[nx], -);
            rep(sta, ( << (m + ))) {
                if(dp[now][sta] != -) {
                    if(up(sta, j)) uMin(dp[nx][noSet(sta,j)], dp[now][sta]);
                    if(maz[i][j]) uMin(dp[nx][Set(sta, j)], dp[now][sta] + );
                }
            }
        }
    }
    int ans = dp[now][( << (m + )) - ];
    uMin(ans,dp[nx][( << (m + )) - ]);
    cout << ans << endl;
}

...因为并没有找到可以提交的地方，所以代码也只是随便测了几组，不过思路好像没什么问题，如果找到代码错误的话可以留言或者私信交流一下