刚开始想这道题想了很久,然后就有个思路,没立即敲。
然后敲的时候将起跳后忍受能力减1,记为降落时忍受能力减1。
增加一个源点n*m * 2和汇点n *m *2+1。
首先,对于每个点node,将其分为两个点node1和node2,node1连向node2的权值就等于这个点的承受能力(当然承受能力为0就不连了,因为没有柱子),然后如果可以从i跳到j,就将i2连向j1,权值无穷大,这里可以是400。对于每一个可以直接跳到外面的点,将其(这时连的是node2)连向汇点,权值为无穷。
对于每一个有Lizard的点,从源点连向它(这时连的是node1),权值为1。(如果题意是起跳后忍受能力减1,则应该是从源点连向node2)。
然后吐槽一下,wa了几发,原来是欧拉距离,不是哈曼顿距离。。。。。。
还有输出格式。
#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <stdlib.h>
#include <stack>
#include <vector>
#include <string.h>
#include <queue>
#include <math.h>
#define msc(X) memset(X,-1,sizeof(X))
#define ms(X) memset(X,0,sizeof(X))
typedef long long LL;
using namespace std;
const int INF=0x3f3f3f3f;
const int MAXN=805;
const int MAXM=32222;
struct _Edge
{
int to,next,cap,flow;
}edge[MAXM];
int tol,hd[MAXN];
void inti(void)
{
tol=0;
msc(hd);
}
void addedge(int u,int v,int w,int rw=0)
{
edge[tol].to=v;edge[tol].cap=w;edge[tol].flow=0;
edge[tol].next=hd[u];hd[u]=tol++;
edge[tol].to=u;edge[tol].cap=rw;edge[tol].flow=0;
edge[tol].next=hd[v];hd[v]=tol++;
}
int Q[MAXN];
int dep[MAXN],cur[MAXN],sta[MAXN];
bool bfs(int s ,int t,int n)
{
int front=0,tail=0;
msc(dep);
dep[s]=0;
Q[tail++]=s;
while(front<tail){
int u=Q[front++];
for(int i=hd[u];i!=-1;i=edge[i].next)
{
int v=edge[i].to;
if(edge[i].cap>edge[i].flow&&dep[v]==-1){
dep[v]=dep[u]+1;
if(v==t) return true;
Q[tail++]=v;
}
}
}
return false;
}
int dinic(int s,int t,int n)
{
int maxflow=0;
while(bfs(s,t,n)){
for(int i=0;i<n;i++) cur[i]=hd[i];
int u=s,tail=0;
while(cur[s]!=-1){
if(u==t){
int tp=INF;
for(int i=tail-1;i>=0;i--)
tp=min(tp,edge[sta[i]].cap-edge[sta[i]].flow);
maxflow+=tp;
for(int i=tail-1;i>=0;i--)
{
edge[sta[i]].flow+=tp;
edge[sta[i]^1].flow-=tp;
if(edge[sta[i]].cap-edge[sta[i]].flow==0)
tail=i;
}
u=edge[sta[tail]^1].to;
}
else if(cur[u]!=-1&&edge[cur[u]].cap>edge[cur[u]].flow
&&dep[u]+1==dep[edge[cur[u]].to]){
sta[tail++]=cur[u];
u=edge[cur[u]].to;
}
else {
while(u!=s&&cur[u]==-1)
u=edge[sta[--tail]^1].to;
cur[u]=edge[cur[u]].next;
}
}
}
return maxflow;
}
char maze[22][22];
#define mpow(X) ((X)*(X))
inline int dis(int i,int j,int ni,int nj)
{return (int)ceil(sqrt(mpow(i-ni)+mpow(j-nj)));}
int main(int argc, char const *argv[])
{
int t,ti=0;
scanf("%d",&t);
while(++ti<=t){
int n,d;
scanf("%d %d",&n,&d);
inti();
for(int i=0;i<n;i++)
scanf("%s",maze[i]);
int m=strlen(maze[0]);
for(int i=0;i<n;i++)
for(int j=0;j<m;j++)
if(maze[i][j]!='0'){
for(int iii=0;iii<n;iii++)
for(int jjj=0;jjj<m;jjj++)
if(maze[iii][jjj]!='0'&&!(i==iii&&j==jjj))
{
if(dis(i,j,iii,jjj)<=d)
addedge((i*m+j)<<1|1,(iii*m+jjj)<<1,400,0);//最多400只
}
if(i<d||j<d||i+d>=n||j+d>=m)
addedge((i*m+j)<<1|1,n*m*2+1,400,0);
addedge((i*m+j)<<1,(i*m+j)<<1|1,maze[i][j]-'0',0);
}
for(int i=0;i<n;i++)
scanf("%s",maze[i]);
int cnt=0;
for(int i=0;i<n;i++)
for(int j=0;j<m;j++)
if(maze[i][j]=='L')
addedge(n*m*2,(i*m+j)<<1,1,0),cnt++;
int res=dinic(n*m*2,n*m*2+1,n*m*2+2);
printf("Case #%d: ",ti);
if(res==cnt) puts("no lizard was left behind.");
else if(res+1==cnt) puts("1 lizard was left behind.");
else printf("%d lizards were left behind.\n",cnt-res);
}
return 0;
}