2020牛客暑期多校第二场 F - Fake Maxpooling(dp/单调队列)
发布日期:2021-05-07 02:28:40 浏览次数:26 分类:精选文章

本文共 3102 字,大约阅读时间需要 10 分钟。


方法一:dp

这个做法还是很神奇的。当 k > 1 k>1 k>1时,考虑一个 k ∗ k k*k kk的子矩阵,我们将最大值假定为右下角的值,那么我们统计答案时,只需要统计矩阵 [ ( k , k ) , ( n , m ) ] [(k,k),(n,m)] [(k,k),(n,m)],而 k ∗ k k*k kk的子矩阵一定可以由若干个重叠的 2 ∗ 2 2*2 22的矩阵构成,那么我们就只需要每次将左边和上面的最值传递到当前位置,这样最后的答案和最后的矩阵没有关系,只和统计开始的行列有关系

#include 
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;#define fi first#define se second#define pb push_back#define ins insert#define lowbit(x) (x&(-x))#define mkp(x,y) make_pair(x,y)#define mem(a,x) memset(a,x,sizeof a);typedef long long ll;typedef long double ld;typedef unsigned long long ull;typedef pair
P;const double eps=1e-8;const double pi=acos(-1.0);const int inf=0x3f3f3f3f;const ll INF=1e18;const int Mod=1e9+7;const int maxn=2e5+10;int gcd[5005][5005],lcm[5005][5005];inline int max(int x,int y){ return x>y?x:y;}int main(){ //freopen("in.txt","r",stdin); //freopen("out.txt","w",stdout); //ios::sync_with_stdio(0),cin.tie(0),cout.tie(0); int n,m,K; scanf("%d%d%d",&n,&m,&K); for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) if(!gcd[i][j]) for(int k=1;i*k<=n && j*k<=m;k++) gcd[i*k][j*k]=k,lcm[i*k][j*k]=i*j*k; for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) if(K!=1) lcm[i][j]=max(lcm[i][j],max(lcm[i-1][j],lcm[i][j-1])); ll ans=0; for(int i=K;i<=n;i++) for(int j=K;j<=m;j++) ans+=lcm[i][j]; printf("%lld\n",ans); return 0;}

方法二:单调队列

实际上对每一行和每一列来说,我们只需想办法求出每 k k k长度的一段的最大值。这个时候可能会想到线段树,ST表等等,但是内存显然不允许这么做,实际上这不就是滑动窗口的最值问题吗?那么使用单调队列,具体做法为:

首先操作每一行,将每个行中 1 − k 1-k 1k的最值存到第 k k k个位置,然后操作每一列,将每列中 1 − k 1-k 1k的最值存到第 k k k个位置,那么最后只需要查询每个 k ∗ k k*k kk矩阵的右下角求和

#include 
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;#define fi first#define se second#define pb push_back#define ins insert#define Vector Point#define lowbit(x) (x&(-x))#define mkp(x,y) make_pair(x,y)#define mem(a,x) memset(a,x,sizeof a);typedef long long ll;typedef long double ld;typedef unsigned long long ull;typedef pair
pii;typedef pair
pdd;const double eps=1e-8;const double pi=acos(-1.0);const int inf=0x3f3f3f3f;const double dinf=1e300;const ll INF=1e18;const int Mod=1e9+7;const int maxn=2e5+10;int ans[5005][5005],a[5005][5005],q[5005];int main(){ //freopen("in.txt","r",stdin); //freopen("out.txt","w",stdout); ios::sync_with_stdio(0),cin.tie(0),cout.tie(0); int n,m,k; cin>>n>>m>>k; for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) if(!ans[i][j]) for(int k=1;k*i<=n && k*j<=m;k++) ans[k*i][k*j]=k,a[k*i][k*j]=i*j*k; memset(ans,0,sizeof ans); int l,r; for(int i=1;i<=n;i++){ l=r=1; for(int j=1;j<=m;j++){ while(r-l>0 && q[l]<=j-k) l++; while(r-l>0 && a[i][q[r-1]]<=a[i][j]) r--; q[r++]=j; ans[i][j]=a[i][q[l]]; } } memset(q,0,sizeof q); for(int j=k;j<=m;j++){ l=r=1; for(int i=1;i<=n;i++){ while(r-l>0 && q[l]<=i-k) l++; while(r-l>0 && ans[q[r-1]][j]<=ans[i][j]) r--; q[r++]=i; ans[i][j]=ans[q[l]][j]; } } ll sum=0; for(int i=k;i<=n;i++) for(int j=k;j<=m;j++) sum+=ans[i][j]; cout<
<
上一篇:字符串哈希
下一篇:2020牛客暑期多校第二场 B - Boundary(简单计算几何)

发表评论

最新留言

哈哈,博客排版真的漂亮呢~
[***.90.31.176]2025年04月03日 23时54分25秒