傳送門

文章目錄

• 題意：
• 思路：

題意：

有一個$n?n$的網格圖，有若干位置有障礙，有$q$個詢問，每次詢問給出兩個點坐標$x1,y1,x2,y2$,詢問把一個正方形箱子從$x1,x2$推到$y1,y2$的時候，箱子最大能是多大。

$n\le 1e3,q\le 3e5$

思路：

首先我們可以預處理出來每個點能放的最大的正方形邊長是多少，即為數組$a$，這個可以通過二分+前綴和或者多源$bfs$來處理出來。

考慮用我們處理的東西來干點什么。

二分$+bfs?$

這個思路顯然是可以的，但是復雜度不允許，而且不是很好優化，需要另尋思路。

我們從$bfs$得到啟發，起點和終點必須聯通？問題就是從起點到終點走過的路徑中$a$的最小值，最大化最小值，這個不就是最大生成樹的最小邊嗎？所以我們將這個圖建成克魯斯卡爾重構樹，按照最大生成樹來建，兩點之間的最大的最小值就是兩點的$lca$，這樣問題就順利解決啦~

// Problem: G - Hangar Hurdles // Contest: Virtual Judge - Camp 2 // URL: https://vjudge.net/contest/454780#problem/G // Memory Limit: 524 MB // Time Limit: 4000 ms // // Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math") //#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native") //#pragma GCC optimize(2) #include<cstdio> #include<iostream> #include<string> #include<cstring> #include<map> #include<cmath> #include<cctype> #include<vector> #include<set> #include<queue> #include<algorithm> #include<sstream> #include<ctime> #include<cstdlib> #include<random> #include<cassert> #define X first #define Y second #define L (u<<1) #define R (u<<1|1) #define pb push_back #define mk make_pair #define Mid ((tr[u].l+tr[u].r)>>1) #define Len(u) (tr[u].r-tr[u].l+1) #define random(a,b) ((a)+rand()%((b)-(a)+1)) #define db puts("---") using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); } //void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); } //void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL; typedef unsigned long long ULL; typedef pair<int,int> PII;const int N=1010,mod=1e9+7,INF=0x3f3f3f3f; const double eps=1e-6;int n; int p[N*N*2],weigth[N*N*2]; int depth[N*N*2]; int fa[N*N*2][22]; int f[N][N],a[N][N]; char s[N][N]; vector<int>v[N*N*2]; struct Edge {int a,b,w;bool operator < (const Edge &W) const {return w>W.w;} }edge[N*N*2];int find(int x) {return x==p[x]? x:p[x]=find(p[x]); }int get(int x1,int y1,int x2,int y2) {// if(x1>x2) swap(x1,x2),swap(y1,y2);// if(y1>y2) swap(y1,y2);return f[x2][y2]-f[x1-1][y2]-f[x2][y1-1]+f[x1-1][y1-1]; }bool check(int x,int y,int mid) {int x1=x-mid,y1=y-mid;int x2=x+mid,y2=y+mid;if(x1<1||x2>n||y1<1||y2>n) return false;if(get(x1,y1,x2,y2)==(x2-x1+1)*(y2-y1+1)) return true;return false; } int get(int x,int y) {return (x-1)*n+y; }void dfs(int u,int f) {fa[u][0]=f; depth[u]=depth[f]+1;for(int i=1;i<=20;i++) fa[u][i]=fa[fa[u][i-1]][i-1];for(auto x:v[u]) if(x!=f) dfs(x,u); }int lca(int a,int b) {if(depth[a]<depth[b]) swap(a,b);for(int i=20;i>=0;i--) if(depth[fa[a][i]]>=depth[b]) a=fa[a][i];if(a==b) return a;for(int i=20;i>=0;i--) if(fa[a][i]!=fa[b][i]) a=fa[a][i],b=fa[b][i];return fa[a][0]; }int main() { // ios::sync_with_stdio(false); // cin.tie(0);scanf("%d",&n);for(int i=1;i<=n;i++) scanf("%s",s[i]+1);for(int i=1;i<=n;i++) {for(int j=1;j<=n;j++) {f[i][j]+=f[i-1][j]; f[i][j]+=f[i][j-1];f[i][j]-=f[i-1][j-1]; f[i][j]+=s[i][j]=='.';}}for(int i=1;i<=n;i++) {for(int j=1;j<=n;j++) {a[i][j]=0;if(s[i][j]=='#') continue;int l=0,r=n,ans=-1;while(l<=r) {int mid=(l+r)>>1;if(check(i,j,mid)) l=mid+1,ans=mid;else r=mid-1;}a[i][j]=1+ans*2;weigth[get(i,j)]=1+ans*2;}}int tot=0;for(int i=1;i<=n;i++) {for(int j=1;j<=n;j++) {// if(s[i][j]=='#') continue;if(i+1<=n) {edge[++tot]={get(i,j),get(i+1,j),min(a[i][j],a[i+1][j])};// edge[++tot]={get(i,j),get(i+1,j),a[i+1][j]};} if(j+1<=n) {edge[++tot]={get(i,j),get(i,j+1),min(a[i][j],a[i][j+1])};// edge[++tot]={get(i,j),get(i,j+1),a[i][j+1]};}}}for(int i=1;i<N;i++) p[i]=i;sort(edge+1,edge+1+tot);int m=tot; tot=n*n;for(int i=1;i<=m;i++){int a=edge[i].a,b=edge[i].b,w=edge[i].w;int pa=find(a),pb=find(b);if(pa==pb) continue;tot++; p[pa]=tot; p[pb]=tot;v[tot].pb(pa); v[tot].pb(pb);weigth[tot]=w;}// weigth[0]=INF;dfs(tot,0);int q; scanf("%d",&q);while(q--) {int x1,y1,x2,y2;scanf("%d%d%d%d",&x1,&y1,&x2,&y2);int a=get(x1,y1),b=get(x2,y2);// if(find(a)!=find(b)) {// puts("0");// continue;// }printf("%d\n",weigth[lca(a,b)]);}return 0; } /**/

總結

以上是生活随笔為你收集整理的2016-2017 Central Europe Regional Contest Hangar Hurdles 克鲁斯卡尔重构树 + 建图的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。