CodeForces 1921 G题解

本文介绍了如何利用二维前缀和解决一个涉及等腰直角三角形攻击范围的问题,通过旋转和分解为固定梯形和长方形来简化计算,并给出了C++代码示例,展示了如何计算长方形和梯形的前缀和,从而解决题目中的问题。
摘要由CSDN通过智能技术生成

Part 0 传送门

CodeForces传送门

洛谷传送门

Part 1 思路

如果攻击范围是正方形或长方形,那么最简单的方法是二维前缀和,但这道题攻击范围是三角形,也可以考虑使用二维前缀和。

因为向右上,右下,左上,左下攻击的攻击范围经过几次旋转都可以变成同一个三角形,所以我们不用一个一个去算,可以去旋转地图后再用同样的方法算。

那要怎么实现二维前缀和呢?

在下图中我们要求的是绿色三角形(即攻击范围),绿色三角形=大三角形-红色梯形-蓝色梯形+黄色长方形,那么问题就变成了怎么求梯形,长方形和三角形的前缀和。

因为攻击范围是等腰直角三角形,所以上图中的大三角形也是等腰直角三角形。如果我们用蓝色梯形右上角的顶点来记蓝色梯形,那么蓝色梯形就是固定的,用红色梯形右上角的顶点来记红色梯形,那么红色梯形就是固定的。

上面,我们说了这道题实际是让我们求梯形,长方形和三角形的前缀和,但实际上我们只要记梯形和长方形的前缀和就可以了。在上图中,我们将红色梯形右上角的顶点沿大三角形的斜边往下移,都可以构成满足条件的图形,当移到蓝色梯形右下角时,就形成了一个三角形。

Part 2 实现

求长方形前缀和,可以用公式:

s u m [ i ] [ j ] = s u m [ i − 1 ] [ j ] + s u m [ i ] [ j − 1 ] − s u m [ i − 1 ] [ j − 1 ] + a [ i ] [ j ] sum[i][j]=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1]+a[i][j] sum[i][j]=sum[i1][j]+sum[i][j1]sum[i1][j1]+a[i][j]

求梯形前缀和很难求,可以一层一层来求。先求一行的前几个的前缀和与一列的前几个的前缀和,再累加即可得到梯形前缀和。

Part 3 代码(C++)

#include<bits/stdc++.h>
using namespace std;
int t,n,m,k,ans;
inline void work(){
	ans=0;
	scanf("%d%d%d",&n,&m,&k);
	vector<string> str;
	str.resize(n+1);
	for(int i=1;i<=n;i++){
		string a;
		cin>>a;
		str[i]=' '+a;
	}
	for(int dir=1;dir<=4;dir++){
		vector<vector<int>> xd1(n+10,vector<int>(m+10,0)),xd2(n+10,vector<int>(m+10,0));
		vector<vector<int>> tx1(n+10,vector<int>(m+10,0)),tx2(n+10,vector<int>(m+10,0));
		vector<vector<int>> jx(n+10,vector<int>(m+10,0));	
		for(int i=1;i<=n;i++){
			for(int j=1;j<=m;j++){
				xd1[i][j]=xd1[i][j-1];
				xd2[i][j]=xd2[i-1][j];
				if(str[i][j]=='#'){
					xd1[i][j]++;
					xd2[i][j]++;
				}
			}
		}
		for(int i=1;i<=n;i++){
			for(int j=1;j<=m;j++){
				tx1[i][j]=tx1[i-1][min(m,j+1)]+xd1[i][j];
			}
		}
		for(int i=1;i<=m;i++){
			for(int j=1;j<=n;j++){
				tx2[j][i]=tx2[min(n,j+1)][i-1]+xd2[j][i];
			}
		}
		for(int i=1;i<=n;i++){
			for(int j=1;j<=m;j++){
				jx[i][j]=jx[i-1][j]+jx[i][j-1]-jx[i-1][j-1];
				if(str[i][j]=='#'){
					jx[i][j]++;
				}
			}
		}
		for(int i=1;i<=n;i++){
			for(int j=1;j<=m;j++){
				int sum=i+j+k-1,x,y;
				if(sum<=m){
					x=1;
					y=sum;
				}else{
					x=min(n,1+sum-m);
					y=m;
				}
				ans=max(ans,tx2[x][y]-tx1[i-1][min(m,j+k+1)]-tx2[min(n,i+k+1)][j-1]+jx[i-1][j-1]);
			}
		}
		swap(n,m);
		vector<string> st(n+1,string(m+1,' '));
		for(int i=1;i<=n;i++){
			for(int j=1;j<=m;j++){
				st[i][j]=str[j][n-i+1];
			}
		}
		str=st;
	}
	printf("%d\n",ans);
}
int main(){
	scanf("%d",&t);
	for(;t--;){
		work();
	}
}

时间复杂度为 O ( 24 N M ) O(24NM) O(24NM)

  • 18
    点赞
  • 22
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值