1009

思路:用dp[i][j].first 来表示到这个点收集了多少的钻石 ,dp[i][j].second 来表示到这个点钻石的价格,对于每个点按first从小到大排序,每次加入数据X的时候,如果X的first和second都比已有的某一个小则可以丢弃,整个过程用列表来维护。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
#include<bits/stdc++.h> 
#define pint std::pair<int,int>
typedef long long ll;
const ll mod=1e9+7;
const ll maxn=105;
const double pi=acos(-1);
using namespace std;
int _,n,m;
int a[maxn][maxn],b[maxn][maxn];
vector< pint > vec[maxn][maxn];
int top = 0;
pint rem[maxn*maxn*maxn];
void check(pint x)
{
//按第一个值从小到大排序的,X的first一定大,这里就是去掉不可能解的情况
while(top && rem[top].second <= x.second) top--;
if(!top || rem[top].first < x.first) rem[++top]=x;
}
void solve(vector<pint> a, vector<pint> b, vector<pint>& ans)
{
int sza=a.size();
int szb=b.size();

top=0;
int topa =0,topb=0;
//先比较first,拿小的去比较 也就是小的先进rem
while (topa < sza && topb < szb) check((a[topa].first < b[topb].first) ? a[topa++] : b[topb++]);
while (topa < sza) check(a[topa++]);
while (topb < szb) check(b[topb++]);
for (int i = 1; i <= top; i++) ans.push_back(rem[i]);
}
int main()
{
ios::sync_with_stdio(false);
//freopen("1009.in","r",stdin);
//freopen("mine.out","w",stdout);
for(cin>>_;_;_--)
{
cin>>n;
for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)cin>>a[i][j],vec[i][j].clear();
for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)cin>>b[i][j];
vec[1][1].push_back(make_pair(a[1][1], b[1][1]));
//特别处理行列
for(int j=2;j<=n;j++)
{
vec[1][j].push_back(make_pair(vec[1][j-1][0].first+a[1][j],vec[1][j-1][0].second+b[1][j]));
}
for(int i=2;i<=n;i++)
{
vec[i][1].push_back(make_pair(vec[i-1][1][0].first+a[i][1],vec[i-1][1][0].second+b[i][1]));
}
for(int i=2;i<=n;i++)
{
for(int j=2;j<=n;j++)
{
solve(vec[i - 1][j], vec[i][j - 1], vec[i][j]);
for (auto& v : vec[i][j])v.first += a[i][j], v.second += b[i][j];
}
}
ll ans = 0;
for (auto& v : vec[n][n]) ans = max(ans, 1ll * v.first * v.second);
printf("%lld\n", ans);
}
return 0;
}