ARC032 D - アットコーダーモンスターズ

https://arc032.contest.atcoder.jp/tasks/arc032_4#

切比雪夫距离，放在3000*3000的平面上，

一个集合就是恰好包含这个集合的矩形，价值是矩形长、宽较大的一个。

先处理价值最小：

可以二分。

也可以不二分。

f[i][j]表示，(j,i)为左上角的最小价值。（也就是正方形的边长）

显然有：f[i][j]<=f[i][j+1]+1

所以每次+1，合法的话，不断往下即可。

方案数：

还是枚举(j,i)，处理(j,i)为左上角，且集合的矩形左边界是j，上边界是i的所有集合。

卡边界？容斥一下就好了。

注意，f[i][z]=i而不是1，矩形可以超出3000*3000的方阵，出去部分特判（反正都是0）。

#include<bits/stdc++.h>
#define reg register int
#define il inline
#define fi first
#define se second
#define mk(a,b) make_pair(a,b)
#define numb (ch^'0')
#define pb push_back
#define solid const auto &
#define enter cout<<endl
#define pii pair<int,int>
using namespace std;
typedef long long ll;
template<class T>il void rd(T &x){
    char ch;x=;bool fl=false;while(!isdigit(ch=getchar()))(ch=='-')&&(fl=true);
    for(x=numb;isdigit(ch=getchar());x=x*+numb);(fl==true)&&(x=-x);}
template<class T>il void output(T x){if(x/)output(x/);putchar(x%+'');}
template<class T>il void ot(T x){if(x<) putchar('-'),x=-x;output(x);putchar(' ');}
template<class T>il void prt(T a[],int st,int nd){for(reg i=st;i<=nd;++i) ot(a[i]);putchar('\n');}
namespace Modulo{
const int mod=1e9+;
il int ad(int x,int y){return x+y>=mod?x+y-mod:x+y;}
il int sub(int x,int y){return ad(x,mod-y);}
il int mul(int x,int y){return (ll)x*y%mod;}
il void inc(int &x,int y){x=ad(x,y);}
il void inc2(int &x,int y){x=mul(x,y);}
il int qm(int x,int y=mod-){int ret=;while(y){if(y&) ret=mul(x,ret);x=mul(x,x);y>>=;}return ret;}
template<class ...Args>il int ad(const int a,const int b,const Args &...args) {return ad(ad(a,b),args...);}
template<class ...Args>il int mul(const int a,const int b,const Args &...args) {return mul(mul(a,b),args...);}
}
using namespace Modulo;
namespace Miracle{
const int N=;
const int M=1e5+;
int n,K;
int len;
int f[N][N];
int sum[N][N];
int jie[M],inv[M];
int h,z;
int C(int n,int m){
    if(n<||m<||n<m) return ;
    return mul(jie[n],inv[m],inv[n-m]);
}
int s(int x,int y){
    if(x<=||y<=) return ;
    x=min(x,z);y=min(y,h);
    return sum[x][y];
}
int calc(int x1,int y1,int x2,int y2){
    // return sum[x2][y2]-sum[x1-1][y2]-sum[x2][y1-1]+sum[x1-1][y1-1];
    return s(x2,y2)-s(x1-,y2)-s(x2,y1-)+s(x1-,y1-);
}
int main(){
    rd(n);rd(K);
    int x,y;
    jie[]=;
    for(reg i=;i<=n;++i) jie[i]=mul(jie[i-],i);
    inv[n]=qm(jie[n]);
    for(reg i=n-;i>=;--i) inv[i]=mul(inv[i+],i+);

    for(reg i=;i<=n;++i){
        rd(x);rd(y);
        ++x;++y;
        sum[x][y]+=;
        h=max(h,y);z=max(z,x);
    }
    for(reg i=;i<=z;++i){
        for(reg j=;j<=h;++j){
            sum[i][j]+=sum[i-][j]+sum[i][j-]-sum[i-][j-];
            // cout<<" i "<<i<<" j "<<j<<" : "<<sum[i][j]<<endl;
        }

    }
    len=0x3f3f3f3f;
    for(reg i=;i<=h;++i){
        for(reg j=z;j>=;--j){
            // f[i][j]=min(f[i][j+1]+1,min(i,z-j+1));
            if(j==z) f[i][j]=i;
            else f[i][j]=f[i][j+]+;
            // cout<<" sum "<<calc(j,i-f[i][j]+1,j+f[i][j]-1,i)<<endl;
            if(calc(j,i-f[i][j]+,j+f[i][j]-,i)>=K){
                // cout<<" ok "<<endl;
                while(calc(j,i-f[i][j]+,j+f[i][j]-,i)>=K){
                    len=min(len,f[i][j]);--f[i][j];
                }
                ++f[i][j];
            }
            // cout<<" i "<<i<<" j "<<j<<" : "<<f[i][j]<<endl;
        }
    }
    int ans=;
    for(reg i=;i<=h;++i){
        for(reg j=z;j>=;--j){
            if(f[i][j]==len){
                int now=;
                int S=calc(j,i-f[i][j]+,j+f[i][j]-,i);
                if(S<K) continue;
                int s1=calc(j,i-len+,j+len-,i-);
                int s2=calc(j+,i-len+,j+len-,i);
                int s3=calc(j+,i-len+,j+len-,i-);
                now=ad(C(S,K),mod-C(s1,K),mod-C(s2,K),C(s3,K));
                inc(ans,now);
            }
        }
    }

    printf("%d\n",len-);
    printf("%d\n",ans);
    return ;
}

}
signed main(){
    // freopen("data.in","r",stdin);
    // freopen("my.out","w",stdout);
    Miracle::main();
    return ;
}

/*
   Author: *Miracle*

*/