题目:https://acm.dingbacode.com/showproblem.php?pid=6981
题意:
给出两个 n ∗ n n*n n∗n的矩阵 a a a和 b b b,求从 ( 1 , 1 ) (1,1) (1,1) 到 ( n , n ) (n,n) (n,n) 经过的路径上,矩阵 a i , j a_{i,j} ai,j的和与 b i , j b_{i,j} bi,j的乘积最大值。数据范围 1 < = n < = 100 1<=n<=100 1<=n<=100, 1 < = a i , j , b i , j < = 1 e 6 1<=a_{i,j},b_{i,j}<=1e6 1<=ai,j,bi,j<=1e6。T组样例, 1 < = T < = 10 1<=T<=10 1<=T<=10。输入数据随机。
思路:
具体思路来自标程。考虑 f [ i ] [ j ] [ k ] f[i][j][k] f[i][j][k]表示到 i , j i,j i,j这个位置,在矩阵 a a a中获得的值为 k k k时,在矩阵 b b b中最高能获得多少。 a n s = m a x ( k ∗ f n , n , k ) ans=max(k*f_{n,n,k}) ans=max(k∗fn,n,k)
对 f i , j f_{i,j} fi,j来说, f i , j , x f_{i,j,x} fi,j,x和 f i , j , y f_{i,j,y} fi,j,y,当 x < y x<y x<y且 f i , j , x < f i , j , y f_{i,j,x}<f_{i,j,y} fi,j,x<fi,j,y,那么 f i , j , x f_{i,j,x} fi,j,x就是一个无效状态,不可能发展成最优解。
在随机数据下,
n
=
100
n=100
n=100,单个
(
i
,
j
)
(i,j)
(i,j)的有效解峰值约为几千(来自随机的玄学,据说在随机下,数列的LCS是
l
o
g
log
log级的),因此我们就可以 开始瞎搞 维护每一个
(
i
,
j
)
(i,j)
(i,j)的有效值。
代码:
#include<bits/stdc++.h>
#include<stdio.h>
#define pb push_back
#define INF 0x3f3f3f3f
#define MAXN 100005
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
typedef pair<double,double>pd;
vector<P>f[105][105];
int a[105][105],b[105][105];
map<int,int>mp;
int m;
P sk[10000005];
void add(const P &t)
{
while(m&&sk[m].second<=t.second)m--;//将无效解pop出去
if(m==0||sk[m].first<t.first)sk[++m]=t;
}
void update(const vector<P>&x,const vector<P>&y,vector<P>&now)
{
int sz1=x.size(),sz2=y.size(),i=0,j=0;
m=0;
while(i<sz1&&j<sz2)add(x[i].first<y[j].first?x[i++]:y[j++]);//按照a的值升序维护
while(i<sz1)add(x[i++]);
while(j<sz2)add(y[j++]);
for(int i=1;i<=m;i++)now.pb(sk[i]);
}
int main() {
int t;scanf("%d",&t);
while(t--)
{
int n;scanf("%d",&n);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
{
scanf("%d",&a[i][j]);
f[i][j].clear();
}
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)scanf("%d",&b[i][j]);
f[1][1].pb(P(a[1][1],b[1][1]));
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
{
//printf("i = %d , j = %d\n",i,j);
if(i==1&&j==1)continue;
if(i==1)f[i][j]=f[i][j-1];
else if(j==1)f[i][j]=f[i-1][j];
else update(f[i-1][j],f[i][j-1],f[i][j]);//根据(i-1,j)和(i,j-1)维护(i,j)的有效解
for(P &p:f[i][j])p.first+=a[i][j],p.second+=b[i][j];
}
ll ans=0;
for(P p:f[n][n])ans=max(ans,1ll*p.first*p.second);
printf("%lld\n",ans);
}
return 0;
}