暴搜 - 状态压缩 - Scaring the Birds - HDU - 4462

暴搜 - 状态压缩 - Scaring the Birds - HDU - 4462

题意：

$$\begin{gathered} 多组测试用例。 \\ 每组包括方阵的阶n，允许摆放稻草人的m个位置(x_i,y_i),以及每个位置能够覆盖的范围r_i,i\in [1,m]。 \end{gathered}$$

$要求至少摆放多少个稻草人，才能覆盖整个n\times n的稻田。$

$这里的覆盖范围是曼哈顿距离(两点间横纵坐标差的和)。$

Sample Input：
4
2
2 2 3 3
1 3
4
2
2 2 3 3
1 4
0
Sample Output：
-1
1

数据范围：

$$\begin{gathered} 2<=n<=50，0<=m<=10，1<=x_i,y_i<=n，r_i<=2n。 \\ Timelimit：1000ms，Memorylimit：32768kB \end{gathered}$$

---

分析：

$共有m个位置可选择，可用一个长度为m的二进制数来表示当前方案。$

$$\begin{gathered} 共有2^m种方案，每种方案在最坏、最暴力的情况会全图搜索10次， \\ 计算次数上限大约为2^{10}\times 50\times 50\times 10\approx 2.5\times 10^7。 \end{gathered}$$

$尽管是最暴力的方式枚举，仍然可以在1s内跑完。事实上测试不到50ms就跑完了。$

具体落实：

$$\begin{gathered} ①、从[0,2^m-1]，每次枚举一种状态，对每种状态用一个vis数组标记哪些点能够被覆盖到， \\ 这里需要注意，允许摆放稻草人的位置可以不必覆盖。 \end{gathered}$$

$$\begin{gathered} ②、枚举当前状态上的1，从这个状态所表示的j点的位置开始暴力标记一个以(x_j,y_j)为中心，2r_j为边长的方阵， \\ 只要这个方阵内的点在半径覆盖范围内，就将其标记。即方阵内满足|x-x_j|+|y-y_j|<=r_j的所有点。 \end{gathered}$$

$③、统计当前状态标记了方阵上多少个点，若所有点均被覆盖了，就更新一下答案。$

代码：

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<vector>

#define P pair<int,int>
#define x first
#define y second
#define ll long long
#define inf 0x3f3f3f3f

using namespace std;

const int N=55,M=1<<10;

int n,m,r[N];
P v[15];
bool vis[N][N];

void dfs(int x,int y,int r)
{
    for(int i=max(1,x-r);i<=min(x+r,n);i++)
        for(int j=max(1,y-r);j<=min(y+r,n);j++)
            if(abs(i-x)+abs(j-y)<=r)
                 vis[i][j]=true;
}

int main()
{
    while(~scanf("%d",&n),n)
    {
        scanf("%d",&m);
        for(int i=0;i<m;i++) scanf("%d%d",&v[i].x,&v[i].y);
        for(int i=0;i<m;i++) scanf("%d",&r[i]);

        int ans=inf;
        for(int i=0;i<M;i++)
        {
            memset(vis,false,sizeof vis);
            for(int k=0;k<m;k++) vis[v[k].x][v[k].y]=true;
            int cnt=0;
            for(int j=0;j<m;j++)
                if((i>>j)&1)
                {
                    dfs(v[j].x,v[j].y,r[j]);
                    cnt++;
                }

            int tmp=0;
            for(int k=1;k<=n;k++)
                for(int j=1;j<=n;j++)
                    if(vis[k][j])
                        tmp++;
            if(tmp==n*n) ans=min(ans,cnt);
        }

        if(ans!=inf) printf("%d\n",ans);
        else puts("-1");
    }

    return 0;
}

---