并查集——vijos1944琵琶湖

vijos1944

下午除了订正只做了这题啊。。

这题我感觉我的xjb算法已经有点用处了……虽然我不知道别人怎么搞的，但是我download觉得我的做法很玄

首先，并查集并不支持断开，所以我们不可能去顺序模拟题意

因此倒序

在Ti时会有哪些点浮出水面

那首先我们二分出点(x,y)在什么时刻浮出水面，然后把压缩后的点向这一个时刻连一条边……

是不是很奇怪……而且有可能会有自边是吧。不管他，反正当森林看

这个用倒挂才可以存下……我本来不知道怎么存的下，因为第一反应显然是邻接表啊。然后想，如果邻接表可以存的话那边表应该也可以存

于是边表

暴力的做法是直接出来一个点，上下左右跑四次并查集，然后n*m判一遍。……可能和直接BFS差不多？

然后一个点浮出水面后，首先对答案的贡献是1，然后去考虑他把上下左右连接起来对答案的贡献，一次有效对答案的贡献是-1，然后直接用并查集维护就可以了。

#include <iostream>#include <cstdio>#include <cmath>#include <algorithm>#include <queue>#include <string>#include <map>#include <cstring>#include <ctime>#include <vector>#define inf 1e9#define ll long long#define For(i,j,k) for(int i=j;i<=k;i++)#define Dow(i,j,k) for(int i=k;i>=j;i--)using namespace std;int a[1001][1001],que[100001],n,m,T,fa[1000001],poi[1000001],f[1000001],nxt[1000001],cnt,ANS[100001],size[1000001],ans;int dx[]={0,1,0,-1,0};int dy[]={0,0,1,0,-1};inline int c(int x,int y){return (x-1)*m+y;}inline void add(int x,int y){poi[++cnt]=y;nxt[cnt]=f[x];f[x]=cnt;}inline int get(int x){return fa[x]==x?x:fa[x]=get(fa[x]);}inline void merge(int x,int y){int tx=get(x),ty=get(y);fa[tx]=ty;if(tx!=ty)ans-=1;}inline void doit(int x){for(int i=f[x];i;i=nxt[i]){int ty=(poi[i]-1)%m+1,tx=(poi[i]-1)/m+1;fa[poi[i]]=poi[i];size[poi[i]]=1;ans++;For(dir,1,4){int px=tx+dx[dir],py=ty+dy[dir];if(px>=1&&py>=1&&px<=n&&py<=m)if(fa[c(px,py)])merge(c(tx,ty),c(px,py));}}ANS[x]=ans;}inline void get(int x,int y){int l=1,r=T,t=0;while(l<=r){int mid=(l+r)/2;if(que[mid]<a[x][y]) t=max(t,mid),l=mid+1;else r=mid-1;}add(t,c(x,y));}int main(){scanf("%d%d",&n,&m);For(i,1,n)For(j,1,m)scanf("%d",&a[i][j]);scanf("%d",&T);For(i,1,T)scanf("%d",&que[i]);For(i,1,n)For(j,1,m)get(i,j);Dow(i,1,T)doit(i);For(i,1,T)printf("%d ",ANS[i]);}