2021杭电多校第三场 Road Discount

题目：https://acm.dingbacode.com/showproblem.php?pid=6981

题意：

​ 给出两个$n\ast n$的矩阵$a$和$b$，求从$(1,1)$ 到$(n,n)$ 经过的路径上，矩阵$a_{i,j}$的和与$b_{i,j}$的乘积最大值。数据范围$1<=n<=100$，$1<=a_{i,j},b_{i,j}<=1e6$。T组样例，$1<=T<=10$。输入数据随机。

思路：

​ 具体思路来自标程。考虑$f[i][j][k]$表示到$i,j$这个位置，在矩阵$a$中获得的值为$k$时，在矩阵$b$中最高能获得多少。$ans=max(k\ast f_{n,n,k})$

​ 对$f_{i,j}$来说，$f_{i,j,x}$和$f_{i,j,y}$，当$x<y$且$f_{i,j,x}<f_{i,j,y}$，那么$f_{i,j,x}$就是一个无效状态，不可能发展成最优解。

​ 在随机数据下，$n=100$，单个$(i,j)$的有效解峰值约为几千（来自随机的玄学，据说在随机下，数列的LCS是$log$级的），因此我们就可以 开始瞎搞 维护每一个$(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;
}