题意：

有一个n∗nn*nn∗n的网格图，有若干位置有障碍，有qqq个询问，每次询问给出两个点坐标x1,y1,x2,y2x1,y1,x2,y2x1,y1,x2,y2,询问把一个正方形箱子从x1,x2x1,x2x1,x2推到y1,y2y1,y2y1,y2的时候，箱子最大能是多大。

n≤1e3,q≤3e5n\le1e3,q\le3e5n≤1e3,q≤3e5

思路：

首先我们可以预处理出来每个点能放的最大的正方形边长是多少，即为数组aaa，这个可以通过二分+前缀和或者多源bfsbfsbfs来处理出来。

考虑用我们处理的东西来干点什么。

二分+bfs?+bfs?+bfs?

这个思路显然是可以的，但是复杂度不允许，而且不是很好优化，需要另寻思路。

我们从bfsbfsbfs得到启发，起点和终点必须联通？问题就是从起点到终点走过的路径中aaa的最小值，最大化最小值，这个不就是最大生成树的最小边吗？所以我们将这个图建成克鲁斯卡尔重构树，按照最大生成树来建，两点之间的最大的最小值就是两点的lcalcalca，这样问题就顺利解决啦~

// 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 克鲁斯卡尔重构树 + 建图相关推荐

Gym - 101173H Hangar Hurdles(bfs+克鲁斯卡尔重构树)
题目链接:点击查看题目大意:给出一个 n∗nn*nn∗n 的矩阵,有些位置存在障碍物,现在有 qqq 次询问,每次询问给出两个点 st=(x1,y1),ed=(x2,y2)st=(x1,y1),ed ...
Central Europe Regional Contest 2019 J. Saba1000kg （并查集+根号讨论）
链接:https://ac.nowcoder.com/acm/contest/7817/I 来源:牛客网时间限制:C/C++ 2秒,其他语言4秒空间限制:C/C++ 262144K,其他语言524 ...
Central Europe Regional Contest 2012 Problem c: Chemist’s vows
字符串处理的题目: 学习了一下string类的一些用法: 这个代码花的时间很长,其实可以更加优化: 代码: 1 #include<iostream> 2 #include<strin ...
2015-2016 ACM-ICPC, Central Europe Regional Contest (CERC 15) F.Frightful Formula（大数NTT）
题目链接:http://codeforces.com/gym/101480/attachments mark一下O(1)快速乘,啥原理啊?还有这个长度计算,貌似也得好好学学啊? 代码: #includ ...
ICPC Central Europe Regional Contest 2019 K. K==S（AC自动机+矩阵快速幂）
Progressive hard octave rock tunes (so-called "phorts") are written using a specifific mus ...
2021 ICPC Southeastern Europe Regional Contest Werewolves（树上背包）
2021 ICPC Southeastern Europe Regional Contest Werewolves(树上背包) 链接题意:给出一个n个节点的树(n≤3000n\le3000n≤300 ...
ICPC Central Russia Regional Contest (CRRC 19)
ICPC Central Russia Regional Contest (CRRC 19) 文章目录 ICPC Central Russia Regional Contest (CRRC 19) A ...
2021 ICPC Southeastern Europe Regional Contest（更新至六题）
2021 ICPC Southeastern Europe Regional Contest A题签到 A. King of String Comparison 题意:给两个字符串,找出有多少对(l, ...
2015-2016 ACM-ICPC Southwestern Europe Regional Contest (SWERC 15)
C. Canvas Painting 合并果子. E. Wooden Signs \(dp(i,l,r)\)表示第\(i\)块木板的长度区间为\([l,r]\)的方案数,根据题意,\(l\)或\(r\ ...

2016-2017 Central Europe Regional Contest Hangar Hurdles 克鲁斯卡尔重构树 + 建图

文章目录

题意：

思路：

2016-2017 Central Europe Regional Contest Hangar Hurdles 克鲁斯卡尔重构树 + 建图相关推荐

最新文章

热门文章