Chtholly Nota Seniorious

题目背景

大样例下发链接: https://pan.baidu.com/s/1nuVpRS1 密码: sfxg

こんなにも、たくさんの幸せをあの人に分けてもらった

だから、きっと

今の、私は

谁が何と言おうと

世界一、幸せな女の子だ

题目描述

——“假如……我是说假如喔。

万一我再过五天就会死，你能不能对我温柔一点？”

巨大的六号兽五天后将袭击浮游大陆。

无数次计算得到的残酷数据表明，只有圣剑瑟尼欧尼斯的适格精灵——珂朵莉·诺塔·瑟尼欧尼斯(Chtholly Nota Seniorious)开启妖精乡之门，才可以以生命为代价守住浮游岛。

“至少，我也希望自己不用消失，也想让别人记住。我也想留下羁绊啊。”

留给妖精少女珂朵莉的时间似乎已经不多了。

年轻的二等技官，妖精仓库的管理员，世界上最后一个人类——威廉·克梅修，数百年前曾经是一名准勇者，掌握着成为一名勇者所需要的所有知识。

大战在即，调整圣剑的状态成为了一项重要的任务。

瑟尼欧里斯（セニオリス）
圣剑的其中之一，在现存的遗迹兵装中，拥有最强大的力量。
拥有非常特殊的资质，只有极少一部分的人才能使用。
由四十一个护符组成。能将所有事物包含不死者都回归「死亡」。

威廉需要调整圣剑的状态，因此他将瑟尼欧尼斯拆分护符，组成了一个nn行mm列的矩阵。

每一个护符都有自己的魔力值。现在为了测试圣剑，你需要将这些护符分成 A,B两部分。

要求如下：

1. 圣剑的所有护符，恰好都属于两部分中的一部分。

2. 每个部分内部的方块之间，可以通过上下左右相互到达，而且每个内部的方块之间互相到达，最多允许拐一次弯。

例如

ＡＡＡＡＡ　　ＡＡＡＡＡ　　ＡＡＡＡＡ
ＡＡＢＡＡ　　ＢａＡＡＡ　　ＡＡＡＢＢ
ＡＢＢＢＡ　　ＢＢＡＡＡ　　ＡＡＡＢＢ
ＡＡＢＡＡ　　ＢａＡＡＡ　　ＡＢＢＢＢ
ＡＡＡＡＡ　　ＡＡＡＡＡ　　ＢＢＢＢＢ

　　(1)　　　　　(2)　　　　　(3)　　　　　　

其中(1)(2)是不允许的分法，(3)是允许的分法。在(2)中，a属于A区域，这两个a元素之间互相到达，没有办法最多只拐一次弯。

现在要问，所有合法的分法中，A区域的极差与B区域的极差 中间较大的一个的 最小值 是多少？

好心而可爱的在一旁默默观察奈芙莲悄悄地告诉你，极差就是区域内最大值减去最小值。

夜晚的风吹拂着，68号岛上的景色竟与地上的森林无异。转念又想，黄金妖精本身就是与森林之中出现，成长，消亡的神秘存在啊。

时间不早了，早上训练中落败的珂朵莉即将回来了。您要尽快和威廉一起调整好圣剑，千万不能迟哟。

输入输出格式

输入格式：

第一行两个自然数n,mn,m

接下来nn行，每行mm个自然数A_{i,j}Ai,j​表示权值

输出格式：

一个整数表示答案。

输入输出样例

输入样例#1：

4 4
1 12 6 11
11 4 2 14
10 1 9 20
4 17 13 10

输出样例#1：

11

说明

样例解释

1  12 6        11
11 4  2        14
10 1  9        20
4        17 13 10

分法不唯一，如图是一种合法的分法。左边部分极差12-1=11，右边一块极差20-10=10，所以答案取这两个中较大者11。没有别的分法，可以使答案更小。

数据范围与约定

测试点	n	m
#1-2	\le 10≤10	\le 10≤10
#3-4	1	\le 2000≤2000
#5-7	\le 200≤200	\le 200≤200
#8-10	\le 2000≤2000	\le 2000≤2000

对于所有的权值1\le A_{i,j} \le 10^91≤Ai,j​≤109

我们需要按行扫这个矩阵两遍，第一遍从上往下扫，第二遍从下往上扫。

基于贪心的原则，对于每一行，我们从右往左扫。当遇到这一行第一个（最靠右的一个）满足的格子的时候停止，然后下一行从同一个位置开始往左扫。这样贪心可以保证满足条件的同时给下一行更多的选择。

把左边的部分记为A部分，右边的部分记为B部分，则扫的时候顺便维护下每一部分扫到上一行为止的最大最小值，在扫这行的过程中就判断下最大值减最小值会不会爆……

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 int a[][],premax[][],premin[][];
 int sufmax[][],sufmin[][],n,m,ans;
 bool check1(int k)
 {int i,j;
     int lim=m,maxa=,mina=2e9,maxb=,minb=2e9;
     for (i=;i<=n;i++)
     {
         int flag=;
         for (j=lim;j>=&&flag==;j--)
         {
             if (max(max(maxa,premax[i][j])-min(mina,premin[i][j]),max(maxb,sufmax[i][j+])-min(minb,sufmin[i][j+]))<=k)
             {
                 lim=j;
                 maxa=max(premax[i][j],maxa);
                 if (j+<=m)
                 maxb=max(sufmax[i][j+],maxb);
                 mina=min(premin[i][j],mina);
                 if (j+<=m)
                 minb=min(sufmin[i][j+],minb);
                 flag=;
             }
         }
         if (flag==) return ;
     }
     return ;
 }
 bool check2(int k)
 {int i,j;
     int lim=m,maxa=,mina=2e9,maxb=,minb=2e9;
     for (i=n;i>=;i--)
     {
         int flag=;
         for (j=lim;j>=&&flag==;j--)
         {
             if (max(max(maxa,premax[i][j])-min(mina,premin[i][j]),max(maxb,sufmax[i][j+])-min(minb,sufmin[i][j+]))<=k)
             {
                 lim=j;
                 maxa=max(premax[i][j],maxa);
                 if (j+<=m)
                 maxb=max(sufmax[i][j+],maxb);
                 mina=min(premin[i][j],mina);
                 if (j+<=m)
                 minb=min(sufmin[i][j+],minb);
                 flag=;
             }
         }
         if (flag==) return ;
     }
     return ;
 }
 int main()
 {int i,j;
   cin>>n>>m;
   for (i=;i<=n;i++)
   for (j=;j<=m;j++)
   scanf("%d",&a[i][j]);
   for (i=;i<=n;i++)
   {
       premax[i][]=premin[i][]=a[i][];
       for (j=;j<=m;j++)
       {
           premax[i][j]=max(premax[i][j-],a[i][j]);
           premin[i][j]=min(premin[i][j-],a[i][j]);
       }
       sufmax[i][m]=sufmin[i][m]=a[i][m];
       for (j=m-;j>=;j--)
       {
           sufmax[i][j]=max(sufmax[i][j+],a[i][j]);
           sufmin[i][j]=min(sufmin[i][j+],a[i][j]);
       }
       sufmin[i][m+]=premin[i][]=2e9;
   }
       int l=;int r=1e9;
       while (l<=r)
       {
           int mid=(l+r)/;
           if (check1(mid)||check2(mid)) ans=mid,r=mid-;
           else l=mid+;
       }
       cout<<ans;
 }