2021杭电多校第三场 Road Discount

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

​ 给出两个 n ∗ n n*n nn的矩阵 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(kfn,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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值