LA 3695 Distant Galaxy

题目

LA 3695 Distant Galaxy

题解

没啥好说的，先离散， 我是枚举两行，固定两行之后就成了序列上的问题了。sum[i]表示两行一共的横向前缀和，f[i]表示第i列在这两行中间的有多少，这样对于确定的两列i,j有ans=f[j]+f[i]+sum[j-1]-sum[i],然后发现对于确定的i，找到一个位置j使得sum[j-1]+f[j]最大即可，所以倒着枚举左端点，可以做到O(1)最大值。(S[i][j]表示第i列，前j行的前缀和，cnt[i][j]表示第i列第j行行星个数)

注：要考虑这个矩形为一条线的情况，即所有行星x（或y）相等时答案为n。

代码

//QWsin
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=100+10;

int n,k1,k2;
int x[maxn],y[maxn],t1[maxn],t2[maxn];

int sum[maxn],f[maxn],S[maxn][maxn],cnt[maxn][maxn];
inline int work(int p1,int p2)
{
    for(int i=1;i<=k1;++i)
    {
        sum[i]=sum[i-1]+cnt[i][p1]+cnt[i][p2];
        f[i]=S[i][p2]-S[i][p1-1];
    }

    int pos=k1,ret=0;
    for(int i=k1-1;i>=1;--i)
    {
        ret=max(ret,f[pos]+sum[pos-1]+f[i]-sum[i]);
        if(f[i]+sum[i-1] > f[pos]+sum[pos-1]) pos=i;
    }
    return ret;
}

int kase=0;
inline void print(int ans){printf("Case %d: %d\n",++kase,ans);}
inline void solve()
{
    memset(cnt,0,sizeof(cnt));

    for(int i=1;i<=n;++i) 
    {
        scanf("%d%d",x+i,y+i);
        t1[i]=x[i],t2[i]=y[i];
    }
    sort(t1+1,t1+n+1);
    sort(t2+1,t2+n+1);

    k1=unique(t1+1,t1+n+1)-t1-1;
    k2=unique(t2+1,t2+n+1)-t2-1;

    if(k1==1||k2==1) {print(n);return ;}

    for(int i=1;i<=n;++i)
    {
        x[i]=lower_bound(t1+1,t1+k1+1,x[i])-t1;
        y[i]=lower_bound(t2+1,t2+k2+1,y[i])-t2;
        ++cnt[x[i]][y[i]];
    }

    for(int i=1;i<=k1;++i)
        for(int j=1;j<=k2;++j)
            S[i][j]=S[i][j-1]+cnt[i][j];

    int ans=0;
    for(int p1=1;p1<=k2;++p1)
        for(int p2=p1+1;p2<=k2;++p2)
            ans=max(ans,work(p1,p2));
    print(ans);
}

int main()
{
    while(scanf("%d",&n)==1&&n) solve();
    return 0;
}