luoguP1606 [USACO07FEB]荷叶塘Lilypad Pond--最短路计数

题目大意：添加最少石头的方案数。

一种朴素的想法是是每个节点向外建边，是荷花建0边，添加建1边，然后跑最短路计数，这样只能得部分分。为什么呢？比如可以多个荷花节点向同一个水节点建边，实际是一种放置石头得方案，但这样就会统计成多种。
因此正确得建图方法是一个荷花联通块，添加一个石头可以到达哪里。，然后再统计最短路条数就是不同得方案数了。

建图抄了dfs的，又写了个bfs的。原理是一样的。

参考代码：

#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long LL;
const int N=31,inf=0x3f3f3f3f;
const int fx[8][2]= {2,-1,2,1,-2,-1,-2,1,1,2,-1,2,1,-2,-1,-2};

struct node {
	int to,next;
} e[500010];
int cnt,head[N*N];
void add(int u,int v) {
	e[++cnt].to=v;
	e[cnt].next=head[u];
	head[u]=cnt;
}
int mp[N][N],dis[N*N],p[N][N],inq[N*N],n,m;
LL f[N*N];
bool vis[N][N];
void dfs(int po,int x,int y) {//dfs建边。 
	vis[x][y]=1;
	for(int i=0; i<8; i++) {
		int dx=x+fx[i][0],dy=y+fx[i][1];
		if(dx<1 || dy<1 || dx>n || dy>m || vis[dx][dy])continue;
		if(mp[dx][dy]==1)dfs(po,dx,dy);//这个点是花，直接跳过去。
		else vis[dx][dy]=1,add(po,p[dx][dy]);//这个点不是花，加一个花
	}
}
void bfs(int po,int x,int y) { //bfs建边，处理荷叶联通块，荷叶直接进队列，向水或终点建边。 
	queue<int>qx,qy;
	qx.push(x),qy.push(y);
	vis[x][y]=1;//进过队列不再进队列 
	while(!qx.empty()) {
		int x=qx.front(),y=qy.front();
		qx.pop(),qy.pop();
		for(int i=0; i<8; i++) {
			int dx=x+fx[i][0],dy=y+fx[i][1];
			if(dx<1 || dy<1 || dx>n || dy>m || vis[dx][dy])continue;
			vis[dx][dy]=1;
			if(mp[dx][dy]==1)qx.push(dx),qy.push(dy);//
			else if(mp[dx][dy]!=2)add(po,p[dx][dy]);
		}
	}
}
queue<int>q;
int main() {
	int i,j,x,y,k,st,ed;
	cin>>n>>m;
	for(i=1; i<=n; i++)
		for(j=1; j<=m; j++)p[i][j]=(i-1)*m+j;
	for(i=1; i<=n; i++)
		for(j=1; j<=m; j++) {
			cin>>mp[i][j];
			if(mp[i][j]==3)st=p[i][j];
			if(mp[i][j]==4)ed=p[i][j];

		}
	for(i=1; i<=n; i++)
		for(j=1; j<=m; j++)
			if(mp[i][j]!=2) {
				memset(vis,0,sizeof(vis));
				bfs(p[i][j],i,j);
			}
	memset(dis,127,sizeof(dis));
	memset(inq,0,sizeof(inq));
	dis[st]=0,inq[st]=1;
	f[st]=1;
	q.push(st);
	while(!q.empty()) {
		int v,u=q.front();
		q.pop();
		inq[u]=0;
		for(int i=head[u]; i; i=e[i].next) {
			v=e[i].to;
			if(dis[v]>dis[u]+1) {
				dis[v]=dis[u]+1,f[v]=f[u];
				if(!inq[v]) {
					inq[v]=1;
					q.push(v);
				}
			} else if(dis[v]==dis[u]+1)
				f[v]+=f[u];
		}
	}
	if(dis[ed]<inf)printf("%d\n%lld\n",dis[ed]-1,f[ed]);
	else printf("-1\n");
	return 0;
}