gym102769B. Bounding Wall ccpc2020秦皇岛

最新推荐文章于 2022-05-11 20:55:35 发布
最新推荐文章于 2022-05-11 20:55:35 发布
阅读量224 收藏
点赞数
分类专栏： 树状数组
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/liufengwei1/article/details/109394355
版权

https://codeforces.com/gym/102769/problem/B

写完到调试到对拍出来前前后后用了6个小时，写这种码农题还是功力不够。。。现场感觉是不可能写出来的，现场也只有pku和渡渡鸟过了这题

思路挺简单的，就枚举这个点在分别的左右上下边界上

大致讲一下在左边界上怎么求最大的矩形

首先l[i][j]表示i能连续最左边到哪，也就是最小的j，r[i][j]为最优边，u[i][j]为最上边也就是最小的x，d[i][j]最下边

然后先求出上边最高可以到哪，由于(x,y)在左边界，我们就去枚举右边界的j，然后首先a[x][j]必须是有的，然后对于他能拓展到的最高的ua[x][j]，x最小可以是多少，这个东西可以在树状数组上二分完成，树状数组存长度为下标，前缀和为有多少条a(i,y)可以拓展到j的，也就是需要预处理(x,y)向上走也就是(i<x),y)可以向右拓展的最远距离，先把他们全部加入树状数组中，再把他们加入向右最远距离的消除vector中

然后枚举右边界j增加的时候，先树状数组上二分找到能够到当前j中，且左边界和右边界向上都能够着的区间中，最大长度是多少，然后再看有哪些这里过后就够不到j+1的(i<x,y)向右拓展的长度，吧他们从树状数组中删掉。

这样我们就求出了(x,y)为左边界，且对于每个j能向上最多到多远，然后再求一遍向下最多能到多远，这就求出了左边界的最大矩阵。然后类似地求上下左右就行了

我吐了，应该有简单一点的写法但我不会写，比如四个方向都抽象一下放到一个函数里面，但感觉不会啊 

#include<bits/stdc++.h>
using namespace std;

const int maxl=1010;

int n,m,q,len;
int b[2][maxl],c[2][maxl],bit[2][maxl];
int a[maxl][maxl];
int la[maxl][maxl],ra[maxl][maxl],ua[maxl][maxl],da[maxl][maxl];
char s[maxl];
vector <int> out[2][maxl];

inline void prework()
{
	scanf("%d%d%d",&n,&m,&q);
	len=max(n,m);
	for(int i=1;i<=n;i++)
	{
		scanf("%s",s+1);
		for(int j=1;j<=m;j++)
			a[i][j]=(s[j]=='#')?1:0;
	}
	for(int i=1;i<=n+1;i++) a[i][m+1]=0;
	for(int i=1;i<=m+1;i++) a[n+1][i]=0;
	for(int i=1;i<=n;i++)
	{	
		for(int j=1;j<=m;j++)
		if(a[i][j]==0)
			la[i][j]=0;
		else if(j-1>0 && a[i][j-1])
			la[i][j]=la[i][j-1];
		else
			la[i][j]=j;
		for(int j=m;j>=1;j--)
		if(a[i][j]==0)
			ra[i][j]=0;
		else if(j+1<=m && a[i][j+1])
			ra[i][j]=ra[i][j+1];
		else
			ra[i][j]=j;
	}
	for(int j=1;j<=m;j++)
	{
		for(int i=1;i<=n;i++)
		if(a[i][j]==0)
			ua[i][j]=0;
		else if(i-1>0 && a[i-1][j])
			ua[i][j]=ua[i-1][j];
		else
			ua[i][j]=i;
		for(int i=n;i>=1;i--)
		if(a[i][j]==0)
			da[i][j]=0;
		else if(i+1<=n && a[i+1][j])
			da[i][j]=da[i+1][j];
		else
			da[i][j]=i;
	}
}

inline void chg(int x,int y)
{
	if(a[x][y]==1)
	{
		ra[x][y]=(y+1<=m && a[x][y+1])?ra[x][y+1]:y;
		for(int j=y-1;j>=1 && a[x][j];j--)
		if(a[x][j])
			ra[x][j]=ra[x][j+1];
		la[x][y]=(y-1>=1 && a[x][y-1])?la[x][y-1]:y;
		for(int j=y+1;j<=m && a[x][j];j++)
			la[x][j]=la[x][j-1];
		ua[x][y]=(x-1>=1 && a[x-1][y])?ua[x-1][y]:x;
		for(int i=x+1;i<=n && a[i][y];i++)
			ua[i][y]=ua[i-1][y];
		da[x][y]=(x+1<=n && a[x+1][y])?da[x+1][y]:x;
		for(int i=x-1;i>=1 && a[i][y];i--)
			da[i][y]=da[i+1][y];
	}
	else
	{
		ua[x][y]=da[x][y]=la[x][y]=ra[x][y]=0;
		if(y-1>0 && a[x][y-1])
			for(int j=y-1;j>=1 && a[x][j];j--)
				ra[x][j]=y-1;
		if(y+1<=m && a[x][y+1])
			for(int j=y+1;j<=m && a[x][j];j++)
				la[x][j]=y+1;
		if(x+1<=n && a[x+1][y]);
			for(int i=x+1;i<=n && a[i][y];i++)
				ua[i][y]=x+1;
		if(x-1>0 && a[x-1][y])
			for(int i=x-1;i>=1 && a[i][y];i--)
				da[i][y]=x-1;
	}
}

inline void add(int i,int x,int bit[])
{
	while(i<=len)
	{
		bit[i]+=x;
		i+=i&-i;
	}
}
inline int sum(int i,int bit[])
{
	int ret=0;
	while(i)
	{
		ret+=bit[i];
		i-=i&-i;
	}
	return ret;
}
inline int qry(int d,int bit[])
{
	int ret=0,sum=0;
	for(int i=10;i>=0;i--)
	if(ret+(1<<i)<=len && sum+bit[ret+(1<<i)]<d)
		ret+=1<<i,sum+=bit[ret];
	return ret+1;
}

inline void init(int a[2][maxl])
{
	for(int i=1;i<=len;i++)
	{
		bit[0][i]=bit[1][i]=0;
		out[0][i].clear();out[1][i].clear();
		a[0][i]=a[1][i]=0;
	}
}

inline void mainwork()
{
	int op,x,y,l1,l2,ans,d;
	for(int cas=1;cas<=q;cas++)
	{
		scanf("%d%d%d",&op,&x,&y);
		if(op==1)
		{
			a[x][y]^=1;
			chg(x,y);
		}
		else
		{
			if(!a[x][y])
			{
				puts("0");
				continue;
			}
			ans=1;
			init(b);
			l1=b[0][0]=b[1][0]=x-ua[x][y]+1;//up down left 0/right 1
			for(int i=x;i>=1 && a[i][y];i--)
			{
				if(y-la[i][y]>0)
					add(x-i+1,1,bit[0]),out[0][y-la[i][y]].push_back(x-i+1);	
				if(ra[i][y]-y>0)
					add(x-i+1,1,bit[1]),out[1][ra[i][y]-y].push_back(x-i+1);
			}
			for(int j=y-1;j>=1;j--)
			{
				if(a[x][j])
				{
					l2=x-ua[x][j]+1;
					d=sum(min(l1,l2),bit[0]);
					b[0][y-j]=(d==0)?0:qry(d,bit[0]);
				}
				for(int id:out[0][y-j])
					add(id,-1,bit[0]);
			}
			for(int j=y+1;j<=m;j++)
			{
				if(a[x][j])
				{	
					l2=x-ua[x][j]+1;
					d=sum(min(l1,l2),bit[1]);
					b[1][j-y]=(d==0)?0:qry(d,bit[1]);
				}
				for(int id:out[1][j-y])
					add(id,-1,bit[1]);
			}
			init(c);
			l1=c[0][0]=c[1][0]=da[x][y]-x+1;
			for(int i=x;i<=n && a[i][y];i++)
			{	
				if(y-la[i][y]>0)
					add(i-x+1,1,bit[0]),out[0][y-la[i][y]].push_back(i-x+1);
				if(ra[i][y]-y>0)
					add(i-x+1,1,bit[1]),out[1][ra[i][y]-y].push_back(i-x+1);
			}
			for(int j=y-1;j>=1;j--)
			{
				if(a[x][j])
				{
					l2=da[x][j]-x+1;
					d=sum(min(l1,l2),bit[0]);
					c[0][y-j]=(d==0)?0:qry(d,bit[0]);
				}
				for(int id:out[0][y-j])
					add(id,-1,bit[0]);
			}
			for(int j=y+1;j<=m;j++)
			{
				if(a[x][j])
				{
					l2=da[x][j]-x+1;
					d=sum(min(l1,l2),bit[1]);
					c[1][j-y]=(d==0)?0:qry(d,bit[1]);
				}
				for(int id:out[1][j-y])
					add(id,-1,bit[1]);
			}
			for(int i=0;i<=len;i++)
			{
				if(b[0][i] && c[0][i])
					ans=max(ans,(b[0][i]+c[0][i]-1)*(i+1));
				if(b[1][i] && c[1][i])
					ans=max(ans,(b[1][i]+c[1][i]-1)*(i+1));
			}
			init(b);
			l1=b[0][0]=b[1][0]=y-la[x][y]+1;  // left right up 0/down 1
			for(int j=y;j>=1 && a[x][j];j--)
			{
				if(x-ua[x][j]>0)
					add(y-j+1,1,bit[0]),out[0][x-ua[x][j]].push_back(y-j+1);
				if(da[x][j]-x>0)
					add(y-j+1,1,bit[1]),out[1][da[x][j]-x].push_back(y-j+1);
			}
			for(int i=x-1;i>=1;i--)
			{
				if(a[i][y])
				{
					l2=y-la[i][y]+1;
					d=sum(min(l1,l2),bit[0]);
					b[0][x-i]=(d==0)?0:qry(d,bit[0]);
				}
				for(int id:out[0][x-i])
					add(id,-1,bit[0]);
			}
			for(int i=x+1;i<=n;i++)
			{
				if(a[i][y])
				{
					l2=y-la[i][y]+1;
					d=sum(min(l1,l2),bit[1]);
					b[1][i-x]=(d==0)?0:qry(d,bit[1]);
				}
				for(int id:out[1][i-x])
					add(id,-1,bit[1]);
			}
			init(c);
			l1=c[0][0]=c[1][0]=ra[x][y]-y+1;
			for(int j=y;j<=m && a[x][j];j++)
			{
				if(x-ua[x][j]>0)
					add(j-y+1,1,bit[0]),out[0][x-ua[x][j]].push_back(j-y+1);
				if(da[x][j]-x>0)
					add(j-y+1,1,bit[1]),out[1][da[x][j]-x].push_back(j-y+1);
			}
			for(int i=x-1;i>=1;i--)
			{
				if(a[i][y])
				{
					l2=ra[i][y]-y+1;
					d=sum(min(l1,l2),bit[0]);
					c[0][x-i]=(d==0)?0:qry(d,bit[0]);
				}
				for(int id:out[0][x-i])
					add(id,-1,bit[0]);
			}
			for(int i=x+1;i<=n;i++)
			{
				if(a[i][y])
				{
					l2=ra[i][y]-y+1;
					d=sum(min(l1,l2),bit[1]);
					c[1][i-x]=(d==0)?0:qry(d,bit[1]);
				}
				for(int id:out[1][i-x])
					add(id,-1,bit[1]);
			}
			for(int i=0;i<=len;i++)
			{
				if(b[0][i] && c[0][i])
					ans=max(ans,(b[0][i]+c[0][i]-1)*(i+1));
				if(b[1][i] && c[1][i])
					ans=max(ans,(b[1][i]+c[1][i]-1)*(i+1));
			}
			printf("%d\n",ans);
		}
	}
}

int main()
{
	int t;
	scanf("%d",&t);
	for(int i=1;i<=t;i++)
	{
		prework();
		printf("Case #%d:\n",i);
		mainwork();
	}
	return 0;
}

 

 

二分抄代码
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
gym-0.13.1.tar.gz
07-29
env = gym.make('CartPole-v0')`。 四、强化学习的基础知识 强化学习是一种机器学习方法，其中智能体通过与环境的交互学习最优策略。gym库中的环境就是这种互动的平台。强化学习的关键概念包括： - **状态...
gym.github.io
05-11
"GYM.github.io-main"的目录结构可能反映了这样的组织方式，其中可能有单独的文件夹存放图片、样式表和脚本文件。 总之，"GYM.github.io"提供了一个学习和实践HTML的平台，特别是对于初学者来说，可以通过查看和...
2020 CCPC 秦皇岛 B.Bounding Wall
Wings
04-29 143
我搬运我自己应该算原创吧 题目 题意 题解 题意 TTT组测试. 给出n×mn \times mn×m的010101矩阵MMM. qqq次操作. 1 x y 将矩阵中(x,y)(x, y)(x,y)取反 2 x y 表示查询 查询如下: 找一个面积最大的矩形边框(平行与坐标轴), 满足: 边框上的点都是000 (x,y)(x, y)(x,y)在边框上 1≤T,n,m,q≤1031 \le T, n, m, q \le 10^31≤T,n,m,q≤103 题解 先来看询问. 首先很容易发现,
2020秦皇岛CCPC 7-2 Bounding Wall(线段树)
0ng的博客
10-24 163
维护每个点最左边/右边/上边的’.'位置 每次查询建立线段树然后枚举高度查找两边长度即可 #include<bits/stdc++.h> using namespace std; #define memarray(a,b) memset(a,b,sizeof(a)) #define pb push_back const long long mod=1e9+7; const int MAXN=1e3+10; const int N=3e2+10; const int INF=0x3f; int n
2020秦皇岛CCPC 7-2 Bounding Wall（树状数组+二分，并查集）
tomjobs的博客
10-23 290
题意： 一个n*m的方格，每个格子可以是’#‘或者’.‘ 。每次操作可以改变一个格子的状态，或者询问经过(x,y)的最大方框面积，且这个方框不能经过’.’ 。 思路： 思路参考自大佬 每次询问一个点(x,y)(x,y)(x,y)，我们只需要考虑以这个点为矩形下边的情况，然后枚举上面那条边。其他的情况只需要将方格旋转90°/180°/270°即可。 我们维护每个点最左边/右边/上边的’.'位置。那么每次修改点，只需要单独维护出这个点左边/右边/上边的位置，复杂度O(n+m)O(n+m)O(n+m)。 你想想
2020 China Collegiate Programming Contest Qinhuangdao Site B. Bounding Wall（思维+并查集）
LSD20164388的博客
10-28 429
传送门 B. Bounding Wall time limit per test4.0 s memory limit per test512 megabytes inputstandard input outputstandard output Alex is a professional computer game player. These days, Alex is playing a war strategy game. His land is a rectangular grid with n r
2020 CCPC秦皇岛 D - Exam Results Gym - 102769E 线段树+离散化
qq_53680405的博客
10-19 341
链接https://codeforces.com/gym/102769/standings 题目大意 有n个学生 每个学生有两个成绩 你可以选择一个为这个学生的最终成绩 设这n个学生的最终成绩中最大的为MAX 问你如何选择 让数量最多的学生最终成绩不低于MAX的百分之P 题目思路 每个学生的两个成绩 ai 和 bi ai >= bi 我们考虑对每个学生划一个区间 代表着让这个学生满足题意的max区间 例如某学生ai = 3 bi =6 P=60 那么对于该学生的满足题意的区间.
2020年CCPC秦皇岛站部分题解
Mr.RainsdRop的胡思乱想
11-25 1136
前言 tmd，哪个天杀的狗东西想出来的这么一个装饰品！！！ 我刷了半小时！！！ 就一只！！！ 还让我给弄成这么个破玩意！！！ 太上头了，得来一道线段树消消气。 Fixing Banners(DFS) 比赛链接：https://codeforces.com/gym/102769/problem/F、 题目大意 现在有666个字符串，你需要从每一个字符串中各取出一个字符，最后拼成"harbin"(顺序无所谓，只要最后得到的666个字符是这666个即可)。 问是否可以完成任务，可以输出"YesYesYes"，
GYM-PROJECT.rar_gym_gym\.project_project management
09-21
在IT行业中，项目管理是一项至关重要的技能，尤其在健身房（GYM）管理这样的领域中，高效、有序的项目管理能够确保健身房的运营顺畅，会员服务优质。本文将深入探讨"GYM-PROJECT.rar"这个压缩包文件所代表的健身房...
Python库 | gym-0.5.2.tar.gz
03-07
安装完成后，可以使用`import gym`来引入库，并通过`gym.make()`函数创建特定的环境，例如经典的CartPole平衡任务，只需输入`env = gym.make('CartPole-v0')`。 Gym库包含了多个内置的环境，这些环境分为两类：无尽...
PyPI 官网下载 | gym-0.4.10.tar.gz
02-01
《PyPI官网下载：探索gym-0.4.10.tar.gz的Python库与分布式技术》 在Python的世界里，PyPI（Python Package Index）是开发者们获取和分享开源软件包的重要平台。本文将围绕“PyPI官网下载”的主题，深入探讨名为...
2020ccpc秦皇岛站部分题解
qq_57404161的博客
05-11 1192
2020ccpc秦皇岛站部分题解
O - 鸣人和佐助(BFS)
Mr.RainsdRop的胡思乱想
10-17 503
O - 鸣人和佐助 题目 佐助被大蛇丸诱骗走了，鸣人在多少时间内能追上他呢？ ![在这里插入图片描述](https://img-blog.csdnimg.cn/20201016202650989.jpg?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQ1NzUwMjk2,size_16,color_FFFFFF,t_70#pic_center)
GYM 101173 L.Lost Logic（构造）
ZSQ
03-20 526
Description 有n个变量，三种赋值，要求写出至多500个约束条件使得满足这些约束条件的变量赋值只有这三种 Input 第一行一整数n表示变量数，之后三行每行n个整数表示一组赋值(2<=n<=50) Output 存在合法约束使得满足这些约束的变量赋值只有这三种那么输出这些约束，否则输出-1 Sample Input Sample Output Solution 如
HDU6606 Distribution of books 2019杭电多校第三场
binarycopycode
07-30 656
考场上二分以后只想了贪心，没想DP 看见题解里的DP1秒想出来，卧槽巨水，完了水题也不会做了。 考虑一个答案，然后dp[i]表示分配前i本书最多可以分给几个人，然后如果无法以i结尾就不更新dp数组 然后把前缀和离散化，dp[i]=max(dp[j]+1),sum[i]-sum[j]<=mid ，那么就用sum[j]>=sum[i]-mid，然后二分出最小的离散化值，权值树状数组...
codeforces1436 E. Complicated Computations
binarycopycode
10-25 498
https://codeforces.com/contest/1436/problem/E 直接用权值树状数组暴力模拟就好了，想了一年想不出来 树状数组b[i]记录值为i上次出现的位置在哪，getmi求前缀最小位置在哪里 然后枚举右端点，不断地更新最小位置，一开始b[i]初值全部设成n+1 然后考虑直接右端点就是n，每个数字上次出现的位置就是他最后出现的位置 如果对于某个数字i，last[i]!=n && getmi(i-1)>last[i]，也就是从末尾到i最后出现的位置，
HDU3333
binarycopycode
08-27 444
8.17那天我再考场上用莫队A了这题。。。。其实是可以卡掉的。正解是线段树或树状数组。 按照询问的r排序，然后r没向右移动一次，如果这个数曾经出现了，就在这个数原来出现的地方减去这个数，再在当前位置加上这个数。 #include #include #include #include #define maxl 30010 #define maxq 100010 using namespace
gym102861E. Party Company
binarycopycode
01-26 411
https://codeforces.com/gym/102861/problem/E 听课也更不上，比赛过了一车人的题也不会，没了啊 这题很显然我们可以看出，由于最后一定是一个连通块属于同一个组，那么由组长o我们可以倍增到最上面的领导o1，使得最上面恰好逼近r 那么最显然的想法就是从这个根节点o1开始，他子树中所有>=l的都会+1，那么子树问题显然可以用dfs序来维护，按dfs序标号，然后分块对每块内部排序，然后左右边界的块暴力找，中间的块一定在子树中，就直接二分后用差分数组作区间加法 O(
2019徐州网络赛 B so easy
binarycopycode
09-07 396
https://nanti.jisuanke.com/t/41384 树状数组是天，1e6*logn*logn随便跑 离散化所有点和区间，要染色的点单点作为一个点，然后两个染色点之间的区间也作为一个点，然后二分找第一个位置在哪。 但是牛逼的标程用unorderd_map+并查集，f[x]存x右边第一个在哪，初始的是f[x]=x+1 unorderd_map的存取查询操作都近似于O（1）,m...
gym102920a. autonomous vehicle 2020
最新发布
10-28
自2020年以来，自动驾驶汽车（Autonomous Vehicle）在这一领域取得了显著进展。 自动驾驶汽车通过使用先进的传感技术，如雷达、激光雷达和摄像头，以及深度学习和人工智能算法，实现了车辆的自主驾驶。它们能够感知...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值