剛開始想這道題想了很久,然後就有個思路,沒立即敲。
然後敲的時候將起跳後忍受能力減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;
}