题目：Problem - 6981 (dingbacode.com)

        题目意思是给出两个n*n的方阵，分别存每个格子的钻石个数和增加的单价，从（1,1）走到（n,n)最后得到了钻石总价是多少，每次只能向右或向下走。

        这题容易想成dp，但是钻石个数和单价对答案贡献的性质不同，很难在同一个转移方程中顾及到，所以这题其实就是一个模拟。

        记录模拟每一个格子的所有情况，即可能的钻石个数和对应的单价，然后从中删去显然不可能发展为最优解的（即一种情况的个数和单价都比另一种小），最后再从（n,n）一格中的所有情况中找最优解。

        怎样删去多余情况？一种方法是直接从（i-1,j)和（i,j-1)两个格子中发展出所有情况，然后按数量从大到小排序，然后从前往后取，若当前的单价小于已知的最大值就不取，大于就取并更新最大值。但这样会超时。所以我们可以提前保证每一个格子的情况都是有序的（按数量从大到小），然后每次更新（i,j)点情况时直接对（i-1,j)和（i,j-1)的情况进行归并排序，同时判断取舍情况，最后得到的情况也是有序的，可以被用作下一次的更新。

代码：

#include<iostream>
#include<string.h>
using namespace std;
typedef long long int ll;
int a[101][101], b[101][101];
struct point { ll x, y; };
struct node {
    int cnt;
    point d[4000];
}map[101];
point dd[4000];
int main()
{
    int tt, n;
    int i, j, k,cnt,t,t1,t2;
    ll x, y, mx,ans;
    cin >> tt;
    while (tt--)
    {
        memset(map, 0, sizeof(map));
        cin >> n;
        for (i = 1;i <= n;i++)
            for (j = 1;j <= n;j++)scanf("%d", &a[i][j]);
        for (i = 1;i <= n;i++)
            for (j = 1;j <= n;j++)scanf("%d", &b[i][j]);
        map[0].cnt = 1;
        for(i=1;i<=n;i++)
            for (j = 1;j <= n;j++)
            {
                t1 = t2 = 1;t = 0;mx = 0;
                while (t1 <= map[j].cnt && t2 <= map[j - 1].cnt)
                {
                    if (map[j].d[t1].x > map[j - 1].d[t2].x)
                    {
                        if (map[j].d[t1].y + b[i][j] > mx)
                        {
                            mx = map[j].d[t1].y + b[i][j];
                            dd[++t].x = map[j].d[t1].x + a[i][j];
                            dd[t].y = mx;
                        }
                        t1++;
                    }
                    else
                    {
                        if (map[j - 1].d[t2].y + b[i][j] > mx)
                        {
                            mx = map[j-1].d[t2].y + b[i][j];
                            dd[++t].x = map[j-1].d[t2].x + a[i][j];
                            dd[t].y = mx;
                        }
                        t2++;
                    }
                }
                while (t1 <= map[j].cnt)
                    dd[++t].x = map[j].d[t1].x + a[i][j],
                    dd[t].y = map[j].d[t1].y + b[i][j], t1++;
                while (t2 <= map[j - 1].cnt)
                    dd[++t].x = map[j - 1].d[t2].x + a[i][j],
                    dd[t].y = map[j - 1].d[t2].y + b[i][j], t2++;
                map[j].cnt = t;
                for (k = 1;k <= t;k++)
                    map[j].d[k].x = dd[k].x, map[j].d[k].y = dd[k].y;
            }
        ans = 0;
        for (i = 1;i <= map[n].cnt;i++)
            ans = max(ans, map[n].d[i].x * map[n].d[i].y);
        printf("%lld\n", ans);
    }
}