这是一个非常神奇的题;
感觉像一个模拟搜索;
但是竟然可以用网络流来解决;
直接粘题解把:
如果不能走通的话,必然说明能够从右上角(图外面)沿雷“跳” ,一直可以“跳”左
下角(图外面) ,因此建好图之后求一个最小割就可以得到结果了。但是关键在于:1.哪些
雷之间可以相互“跳” ?2.哪些雷可以从右上角“跳”过去,哪些雷可以“跳”到左下角?
第二个问题很好办,如果地雷的范围能接触到最上或者最右的格子,就可以从右上角跳
到这个雷上,如果地雷的范围能接触到最下或者最右的格子,就可以跳到左下角。
第一个问题需要分类讨论一下,如果两个雷在同一水平线或者竖直线上,当两个雷的距
离不超过 2*K+1 时可以认为两个雷是连通的,但如果不在同一水平线或者竖直线上时,当
两个雷的距离不超过 2*K+2 时可以认为两个雷是连通的。
注意雷要拆点,容量为 1,其他边容量为 INF
代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define maxn 1550
#define maxm 40009
#define inf 9999999
using namespace std; struct edge
{
int from,to,cap,flow;
edge() {}
edge(int from,int to,int cap,int flow)
:from(from),to(to),cap(cap),flow(flow) {}
}; struct dinic
{
int n,m,s,t;
edge edges[maxm*];
int head[maxn];
int next[maxm*];
bool inq[maxn];
int d[maxn];
int cur[maxn];
void init(int n)
{
this->n=n;
m=;
memset(head,-,sizeof(head[])*(n+));
} void addedge(int from,int to,int cap)
{
next[m]=head[from];
edges[m]=edge(from,to,cap,);
head[from]=m++;
next[m]=head[to];
edges[m]=edge(to,from,,);
head[to]=m++;
}
bool bfs()
{
memset(inq,,sizeof(inq[])*(n+));
queue<int>q;
q.push(s);
d[s]=;
inq[s]=;
while(!q.empty())
{
int u=q.front();
q.pop();
for(int i=head[u]; i!=-; i=next[i])
{
edge& e=edges[i];
int v=e.to;
if(!inq[v]&&e.cap>e.flow)
{
inq[v]=;
d[v]=d[u]+;
q.push(v);
if(v==t)return ;
}
}
}
return ;
}
int dfs(int u,int a)
{
if(u==t||a==)return a;
int flow=,f;
for(int&i=cur[u]; i!=-; i=next[i])
{
edge& e=edges[i];
int v=e.to;
if(d[u]+==d[v]&&(f=dfs(v,min(a,e.cap-e.flow)))>)
{
e.flow+=f;
edges[i^].flow-=f;
flow+=f;
a-=f;
if(a==)break;
}
}
return flow;
}
int maxflow(int s,int t)
{
this->s=s;
this->t=t;
int flow=;
while(bfs())
{
memcpy(cur,head,sizeof(head[])*(n+));
flow+=dfs(s,inf);
}
return flow;
}
}; dinic solve; struct node
{
int x,y;
} no[maxn];
char s[];
int main()
{
// freopen("test0.in","r",stdin);
int t;
int n,m,d;
scanf("%d",&t);
while(t--)
{
int cnt=;
scanf("%d%d%d",&n,&m,&d);
for(int i=; i<=n; i++)
{
scanf("%s",s+);
for(int j=; j<=m; j++)
{
if(s[j]=='*')
{
no[cnt].x=i;
no[cnt++].y=j;
}
}
}
int sr=;
int tr=*cnt+;
solve.init(*cnt+);
for(int i=; i<cnt; i++)
{
solve.addedge(*i,*i+,);
if((no[i].x-)<=d||((m-no[i].y)<=d))
solve.addedge(sr,*i,inf-);
if((n-no[i].x)<=d||(no[i].y-)<=d)
solve.addedge(*i+,tr,inf-);
}
for(int i=; i<cnt; i++)
for(int j=; j<cnt; j++)
{
int dis=abs(no[i].x-no[j].x)+abs(no[i].y-no[j].y);
if(no[i].x==no[j].x||no[i].y==no[j].y)
{
if(dis<=*d+)
solve.addedge(*i+,*j,inf-);
}
else
{
if(dis<=*d+)
solve.addedge(*i+,*j,inf-);
}
}
printf("%d\n",solve.maxflow(sr,tr));
}
return ;
}