#include #include #define eb emplace_back #define all(a) (a).begin(),(a).end() #define RD(T,...) T __VA_ARGS__;li(__VA_ARGS__) #define LL(...) RD(ll,__VA_ARGS__) #define VVL(n,m,...) vvl __VA_ARGS__;setsize({n,m},__VA_ARGS__);li(__VA_ARGS__) #define JO(a,b) a##b #define jo(a,b) JO(a,b) #define uq_symbol(a) jo(a,__LINE__) #define FE(v,e,...) for(au&&__VA_OPT__([)e __VA_OPT__(,__VA_ARGS__]):v) #define FO(n) for(ll uq_symbol(i)=n;uq_symbol(i)-->0;) #define FOR(i,...) for(au[i,i##O,i##E]=rng(0,__VA_ARGS__);i=i##O;i-=i##E) #define I template #define J class #define au auto #define rr return using std::cin,std::cout,std::begin,std::end,std::rbegin; using std::swap,std::move,std::abs,std::prev,std::next; using std::tuple,std::array,std::bitset,std::minmax,std::get; using vo=void;using bo=bool; vo solve(); Iconstexpr bo same=std::is_same_v; using is=std::istream;using os=std::ostream; using dd=long double;using ll=long long;using ull=unsigned long long; using lll=__int128_t;using ulll=__uint128_t;using str=std::string; using namespace atcoder;using ml=modint;au&operator<<(os&o,const ml&x){rr o<T zz(T x){rr x<0?-x:x;} Ill len(const T&a){rr a.size();} struct edg{ ll t,w; edg(){} edg(ll t,ll w=1):t(t),w(w){} friend os&operator<<(os&o,const edg&e){rr o<struct cp{ A a={};B b={}; cp(){} cp(A a,B b):a(a),b(b){} cp(std::pairp):a(p.first),b(p.second){} bo operator==(const cp&c)const{rr a==c.a&&b==c.b;} bo operator<(const cp&c)const{rr a!=c.a?a(const cp&c)const{rr a!=c.a?a>c.a:b>c.b;} friend is&operator>>(is&i,cp&c){rr i>>c.a>>c.b;} friend os&operator<<(os&o,const cp&c){rr o<; Istruct tr{ A a={};B b={};C c={}; tr(){} tr(A a,B b,C c):a(a),b(b),c(c){} bo operator==(const tr&t)const{rr a==t.a&&b==t.b&&c==t.c;} au operator<(const tr&t)const{rr a!=t.a?a>(is&i,tr&t){rr i>>t.a>>t.b>>t.c;} friend os&operator<<(os&o,const tr&t){rr o<; IT&srt(T&a){sort(all(a));rr a;} Idecltype(au)rb(T&a){rr*rbegin(a);} Iau pot(T&a){au r=a.top();a.pop();rr r;} Iusing map=std::map; Ios&operator<<(os&o,const map&m){fe(m,e)o<using umap=std::unordered_map; Ios&operator<<(os&o,const umap&m){fe(m,e)o<struct ve; Iconstexpr bo isv=0; Iconstexpr bo isv> =1; Iconstexpr bo isv> =1; Iau rawv(V){if constexpr(isv)rr rawv(V(1)[0]);else rr V();} Istruct ve:std::vector{ using std::vector::vector; using T=decltype(rawv(V())); Ive(const std::vector&v={}){fe(v,e)this->eb(e);} ve&operator+=(const ve&u){au&v=*this;fo(i,len(v))v[i]+=u[i];rr v;} ve&operator-=(const ve&u){au&v=*this;fo(i,len(v))v[i]-=u[i];rr v;} ve&operator^=(const ve&u){fe(u,e)this->eb(e);rr*this;} ve&operator+=(const T&x){au&v=*this;fe(v,e)e+=x;rr v;} ve&operator-=(const T&x){au&v=*this;fe(v,e)e-=x;rr v;} ve&operator*=(const T&x){au&v=*this;fe(v,e)e*=x;rr v;} Iau&operator+=(const bitset&a){fo(i,n)(*this)[i]+=a[i];rr*this;} Iau&operator-=(const bitset&a){fo(i,n)(*this)[i]-=a[i];rr*this;} ve operator+(const ve&u)const{rr ve(*this)+=u;} ve operator-(const ve&u)const{rr ve(*this)-=u;} ve operator^(const ve&u)const{rr ve(*this)^=u;} ve operator+(const T&x)const{rr ve(*this)+=x;} ve operator-(const T&x)const{rr ve(*this)-=x;} ve operator*(const T&x)const{rr ve(*this)*=x;} ve&operator++(){rr*this+=1;} ve&operator--(){rr*this-=1;} ve operator-()const{rr ve(*this)*=-1;} au lower_bound(const V&x)const{rr std::lower_bound(all(*this),x);} au upper_bound(const V&x)const{rr std::upper_bound(all(*this),x);} ll arglb(const V&x)const{rr lower_bound(x)-begin(*this);} vo fill(T x={}){fe(*this,e)if constexpr(isv)e.fill(x);else e=x;} }; Iau vec(const ll(&s)[n],T x={}){ if constexpr(n==i+1)rr ve(s[i],x); else{au X=vec(s,x);rr ve(s[i],X);} } Ivo setsize(const ll(&l)[n],A&...a){((a=vec(l,rawv(a))),...);} Iusing vve=ve>;using vl=ve;using vvl=vve; using gl=ve>;using ge=ve>; struct io{io(){cin.tie(0)->sync_with_stdio(0); cout<os&operator<<(os&o,const bitset&b){fo(i,n)o<os&operator<<(os&o,const tuple&t){ apply([&](const au&...a){ll i=0;(((o<os&operator<<(os&o,const std::vector&v){fe(v,e)o<?nl:sp);rr o;} Ivo pp(const A&...a){ll i=0;((cout<is&operator>>(is&i,std::vector&v){fe(v,e)i>>e;rr i;} Ivo li(A&...a){(cin>>...>>a);} Ill bsl(F f,ll o,ll n){while(zz(o-n)>1)(f((o+n)/2)?o:n)=(o+n)/2;rr o;} Istruct rec:F{rec(F&&f):F(std::forward(f)){}Idecltype(au)operator()(A&&...a)const{rr F::operator()(*this,std::forward(a)...);}}; templateT fill(){ if constexpr(same)return x; else if constexpr(same)return{x,x}; else if constexpr(same)return{x,x,x}; else if constexpr(same>>){T r;fe(r,e)e=x;return r;} assert(0); } templateT omax(T a,T b){return bT omin(T a,T b){return aT osum(T a,T b){return a+b;} templatestruct segmax:segtree,fill>{ using segtree,fill>::segtree; }; int main(){ll T=1;fo(T)solve();} vo solve(){ LL(N,M); VVL(N,M,a); fe(a,e)srt(e); pp(bsl([&](ll X){//崖をX以下にできるか segmaxcu(N),nx(N); fo(i,N)cu.set(i,1); fo(i,1,N){ fo(j,M)nx.set(j,0); fo(j,M){ ll l=a[i-1].arglb(a[i][j]-X); ll r=a[i-1].arglb(a[i][j]+1); nx.set(j,cu.prod(l,r)); } swap(cu,nx); } rr cu.all_prod()>=1; },oo,-2)); }