有一种神奇的算法叫做floodfill
就是一个n*m的矩阵,a[i][j]为当前高度,我们可以任选一个点倒水,开始bfs,如果要搜的点没有被搜到过,并且高度小于当前的点,我们就把这个点加入队列中
而这道题我们可以用floodfill做
先假设有解,也就是说假设干旱区里的所有城区都能建水利措施。
对于第一行湖泊,显然我们bfs往下floodfill覆盖的最后一行肯定是个区间,那么对于第一行两个点floodfill后覆盖的最后一行的两个区间,如果左边端点floodfill的右边一部分与右边端点floodfill的左边一部分有相交部分,那么这两个floodfill最后一行所覆盖的区间的左端点肯定是相同的,但是我们选择建蓄水厂肯定是选择右边端点,因为右边端点的floodfill能够多覆盖一部分。
所以我们从左往右枚举最后一行的端点,从最后一行的端点先倒着做一遍floodfill,所覆盖的第一行的区间中选择最右边的端点,再正着做一遍floodfill,直接标记floodfill中扫到的最后一行的区间里的端点,再把这个区间的右端点的右边一个点进行倒着做一遍floodfill......就这样循环做下去直到(n,m)被标记
那么怎么判断无解呢?显然,我们把第一行所有端点都做一遍正着的floodfill,也就是说能建蓄水厂的都建起来。再统计最后一行的没被标记的点数,如果有被标记的,就按照无解的输出形式输出,否则就是有解,按照有解的办法解决问题
可能解释有些不清,可以线下找我)
那么就是我们的代码,其中booms函数是正着floodfill,boom函数是倒着floodfill
Code:
#include<iostream>
#include<cstdio>
using namespace std;
int dx[]={0,0,1,-1};
int dy[]={1,-1,0,0};
int n,m,a[10001][10001];
bool vis[10001][10001];
bool viss[10001][10001];
int qx[400001],qy[400001];
int ans;
int boom(int p)
{
int tail=1,head=1;
vis[n][p]=1;
qx[1]=n;
qy[1]=p;
while(head<=tail)
{
int x=qx[head],y=qy[head];
head++;
for(int i=0;i<4;i++)
{
int tx=x+dx[i],ty=y+dy[i];
if(tx>=1&&ty>=1&&tx<=n&&ty<=m&&!vis[tx][ty]&&a[x][y]<a[tx][ty])
{
qx[++tail]=tx;
qy[tail]=ty;
vis[tx][ty]=1;
}
}
}
for(int i=m;i>=1;i--)
if(vis[1][i])
return i;
return 0;
}
void booms(int p)
{
int tail=1,head=1;
qx[1]=1;
qy[1]=p;
viss[1][p]=1;
while(head<=tail)
{
int x=qx[head],y=qy[head];
head++;
for(int i=0;i<4;i++)
{
int tx=x+dx[i],ty=y+dy[i];
if(tx>=1&&ty>=1&&tx<=n&&ty<=m&&!viss[tx][ty]&&a[x][y]>a[tx][ty])
{
qx[++tail]=tx;
qy[tail]=ty;
viss[tx][ty]=1;
}
}
}
}
inline int read()
{
int r=0,w=1;
char ch=getchar();
while(ch<'0'||ch>'9'){
if(ch=='-') w=~w+1;
ch=getchar();
}
while(ch>='0'&&ch<='9'){
r=(r<<3)+(r<<1)+ch-(1<<4)-(1<<5);
ch=getchar();
}
return r*w;
}
int main()
{
//freopen("P1514_5.in","r",stdin);
n=read();m=read();
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
a[i][j]=read();
for(int i=1;i<=m;i++)
booms(i);
ans=0;
for(int i=1;i<=m;i++)
if(!viss[n][i])
ans++;
if(ans){
cout<<0<<endl<<ans;
return 0;
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
vis[i][j]=viss[i][j]=0;
int right=1;
while(!viss[n][m]&&!vis[n][m])
{
int x=boom(right);
booms(x);
for(int i=m;i>=1;i--)
if(!viss[n][i])
right=i;
ans++;
}
cout<<1<<endl<<ans;
return 0;
}