[SDOI2005]动物园 匹配

一开始读错题以为动物不能停，，如果不能停的话就需要黑白染色原理、

这个题只需要广搜最短路判断后直接dinic就可以了

注：为保障程序兼容性。尽量用传统写法，如最短距离数组初始化为INF，而不是-1，不然判断会想当然出错

码：

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
int tot=-1,i,j,n,m,yuan[100005],ky[105][105],xx[100005],bh[100005],gs[100005],yy[100005],D[100005],xs,ys,zhong[100005],xia[100005],v[100005],hou[100005],d[105][105],S,T,p,r,x[100005],y[100005],t[100005],daan[105];
char tu[105][105];
queue<int>q,qx,qy;
void jian(int a,int b,int c)
{
	++tot,hou[tot]=yuan[a],yuan[a]=tot,zhong[tot]=b,v[tot]=c;
}
void jia(int a,int b,int c)
{
	jian(a,b,c);jian(b,a,0);
}
void bfs()
{
	memset(d,-1,sizeof(d));
   qx.push(xs);
   qy.push(ys);
   d[xs][ys]=0;
   while(!qx.empty())	
	{
		int stx=qx.front(),sty=qy.front();
		qx.pop();qy.pop();
	if(sty>1&&tu[stx][sty-1]!='*'&&d[stx][sty-1]==-1)
	{
	d[stx][sty-1]=d[stx][sty]+1;	
	qx.push(stx);qy.push(sty-1);
	}	
	if(stx>1&&tu[stx-1][sty]!='*'&&d[stx-1][sty]==-1)
	{
	d[stx-1][sty]=d[stx][sty]+1;	
	qx.push(stx-1);qy.push(sty);
	}			
	if(sty<n&&tu[stx][sty+1]!='*'&&d[stx][sty+1]==-1)
	{
	d[stx][sty+1]=d[stx][sty]+1;	
	qx.push(stx);qy.push(sty+1);
	}	
	if(stx<n&&tu[stx+1][sty]!='*'&&d[stx+1][sty]==-1)
	{
	d[stx+1][sty]=d[stx][sty]+1;	
	qx.push(stx+1);qy.push(sty);
	}			
    }
}
bool bfs2()
{
	for(i=1;i<=T;i++)xia[i]=yuan[i],D[i]=99999999;
	q.push(S);
	D[S]=1;
	while(!q.empty())
	{
		int st=q.front();
		q.pop();
		for(i=xia[st];i!=-1;i=hou[i])
		{
			if(v[i]==0)continue;
			int nd=zhong[i];
			if(D[nd]!=99999999)continue;
			D[nd]=D[st]+1;
			q.push(nd);		
		}		
	}
if(D[T]==99999999)	return 0;
else return 1;
}
int dfs(int o,int limit)
{
	if(!limit||o==T)return limit;
	int i,f=0,flow=0,nd;
	for(i=xia[o];i!=-1;i=hou[i])
	{
		xia[o]=i;
		 nd=zhong[i];
		 if(D[nd]==D[o]+1&&(f=dfs(nd,min(limit,v[i]))))
		 {
		 	flow+=f;
		 	limit-=f;
		 	v[i]-=f;
		 	v[i^1]+=f;
			 if(!limit)break; 	
		 }	
	}	
	return flow;
}
void dinic()
{
	while(bfs2())
	{
		dfs(S,1000000);		
	}	
}
int jdz(int o)
{
	if(o<0)return -o;
	else return o;
}
int main()
{
	memset(yuan,-1,sizeof(yuan));
	scanf("%d",&n);
	for(i=1;i<=n;i++)
	for(j=1;j<=n;j++)
	{
		scanf("%c",&tu[i][j]);
		while(tu[i][j]!='*'&&tu[i][j]!='.')scanf("%c",&tu[i][j]);	
	}
	scanf("%d",&p);
	for(i=1;i<=p;i++)
	{
		scanf("%d%d",&xx[i],&yy[i]);
	}
		for(i=1;i<=p;i++)
	{
		scanf("%d",&v[i]);		
	}
scanf("%d",&r);
for(i=1;i<=r;i++)
{
	scanf("%d%d%d%d",&t[i],&x[i],&y[i],&bh[i]);
	gs[bh[i]]++;
}
for(i=1;i<=p;i++)	
{
	xs=xx[i];ys=yy[i];
bfs();	
	for(j=1;j<=r;j++)
	{
		int jl=v[bh[j]]*t[j];
		if(d[x[j]][y[j]]<=jl&&d[x[j]][y[j]]!=-1)
	{
	ky[i][bh[j]]++;		
	}	
	}
}
	S=p*2+1;T=p*2+2;
	for(i=1;i<=p;i++)
	{
		jia(S,i,1);
		jia(p+i,T,1);
	for(j=1;j<=p;j++)	
	{
	if(ky[i][j]==gs[j])	
	jia(i,p+j,1);	
	}	
	}
	dinic();
		for(i=1;i<=p;i++)
		{
			for(j=xia[i];j!=-1;j=hou[j])
			{
				int nd=zhong[j];
				if(nd==S||v[j])continue;
			   daan[nd-p]=i;						
			}			
		}
		for(i=1;i<=p;i++)
		{
			printf("%d %d %d\n",i,xx[daan[i]],yy[daan[i]]);			
		}	
}