POJ3271/BZOJ 1698: [Usaco2007 Feb]Lilypad Pond 荷叶池塘

这道题“真的”好玩 手动笑cry

一开始 你可能会想到直接bfs 但是 因为原本荷叶的存在 你会很烦

第一问还好说 统计方案数就GG

然后 你可能会想到把这些荷叶缩成连通块。。。

可是！！！！ 你又会发现有问题

比如这种

0 0 y 0 0

1 0 0 0 1

0 0 x 0 0

x和y都是清水 从x去到y 可以经过左边的联通块去 可以经过右边的联通块去 答案有两种但其实就只有x，y放荷叶一种

是不是很TM爽！

最后 我来告诉你这个世界的真相

因为它是统计放荷叶的方案数啊

我们可以去避免遇到那些原本存在的嘛  我们就预处理出每个起点以及水位置的花费1个荷叶能到的地方 然后建边 跑spfa就可以了~

为什么呢...因为这样就肯定没有环了 而且上面那种情况其实就是从x直接到y 就不怕算多情况了

觉得这个代码还是挺好看的（毕竟挺快的）

#include<cmath>#include<queue>#include<cstdio>#include<cstring>#include<cstdlib>#include<algorithm>#define me(a,x) memset(a,x,sizeof a)using namespace std;typedef long long LL;const int N=31,inf=1e9;const int dx[8]={2,2,-2,-2,1,-1,1,-1};const int dy[8]={-1,1,-1,1,2,2,-2,-2};inline int read(){    int x=0,f=1; char ch=getchar();    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}    return x*f;}struct node{int y,next;}a[500010]; int len,first[N*N];void ins(int x,int y){    a[++len].y=y,a[len].next=first[x],first[x]=len;}int mp[N][N],d[N*N],p[N][N],n,m; LL f[N*N];bool b[N][N];void dfs(int o,int x,int y){    b[x][y]=1;    for(int k=0;k<8;k++)    {        int xx=x+dx[k],yy=y+dy[k];        if(xx<1 || yy<1 || xx>n || yy>m || b[xx][yy])continue;        if(mp[xx][yy]==1)dfs(o,xx,yy);        else b[xx][yy]=1,ins(o,p[xx][yy]);    }}queue<int>q; bool v[N*N];int main(){    int i,j,x,y,k,st,ed;    n=read(),m=read();    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++)    {        mp[i][j]=read();        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]==0||mp[i][j]==3)me(b,0),dfs(p[i][j],i,j);    me(d,63); d[st]=0,v[st]=f[st]=1; q.push(st);    while(!q.empty())    {        x=q.front(); q.pop(); v[x]=0;        for(k=first[x];k;k=a[k].next)        {            y=a[k].y;            if(d[y]>d[x]+1)            {                d[y]=d[x]+1,f[y]=f[x];                if(!v[y]){v[y]=1; q.push(y);}            }            else if(d[y]==d[x]+1)                f[y]+=f[x];        }    }    if(d[ed]<inf)printf("%d\n%lld\n",d[ed]-1,f[ed]);    else printf("-1\n");    return 0;}