给出平面上的n个点,找一个矩形,使得边界上包含尽量多的点。

【输入格式】

输入的第一行为数据组数T。每组数据的第一行为整数n(1≤n≤100);以下n行每行两个整数,即各个点的坐标(坐标均为绝对值不超过109的整数)。输入结束标志为n=0。

【输出格式】

对于每组数据,输出边界点个数的最大值。

先贴题解的精妙预处理:

不难发现,除非所有输入点都在同一行或者同一列上(此时答案为n),最优矩形的4条边都至少有一个点(一个角上的点同时算在两条边上)。这样,我们可以枚举4条边界所穿过的点,然后统计点数。这样做的时间复杂度为O(n5)(统计点数还需要O(n)时间),无法承受。

和《子序列》一题类似,可以考虑部分枚举,即只枚举矩形的上下边界,用其他方法确定左右边界,过程如图1-31所示。

对于竖线i,我们用left[i]表示竖线左边位于上下边界上的点数(注意,不统计位于该竖线上的点),on[i]和on2[i]表示竖线上位于上下边界之间的点数(区别在于on[i]不统计位于上下边界上的点数,而on2[i]要统计)。这样,给定左右边界ij时,矩形边界上的点数为left[j]-left[i]+on[i]+on2[j]。当右边界j确定时,on[i]-left[i]应最大。

枚举完上下边界后,我们先花O(n)时间按照从左到右的顺序扫描一遍所有点,计算left、on[i]和on2[i]数组,然后枚举右边界j,同时维护on[i]-left[i](i<j)的最大值。这一步本质上等价于例题《开放式学分制》。

再讲解一下这题的离散化处理

给每个点一个新的坐标 x坐标 为其原x坐标的排序后位置

y坐标 为其原y坐标的排序后位置

对于对距离不敏感的点集问题 这样十分方便处理问题

贴代码:

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <ctime>
#include <algorithm>
#include <iostream>
#include <sstream>
#include <string>
#define oo 0x13131313
const int maxn=100+5;
using namespace std;
struct point
{
int x,y;
};
point P[maxn];
int N;
int YN,XN;
int ans;
int FF[maxn][maxn];
int CE[maxn][maxn];
int UP[maxn][maxn];
void input()
{
YN=0;XN=0;
memset(FF,0,sizeof(FF));
memset(CE,0,sizeof(CE));
memset(UP,0,sizeof(UP));
for(int i=1;i<=N;i++)
scanf("%d%d",&P[i].x,&P[i].y);
}
int cmp1(const void *i,const void *j)
{
point *ii=(point *)i,*jj=(point *)j;
return ii->x-jj->x;
}
void getXN()
{
int i;
qsort(P+1,N,sizeof(P[1]),cmp1);
for(i=1;i<=N;i++)
{
if(P[i].x!=P[i-1].x)
{
XN++;
P[i-1].x=XN-1;
}
else P[i-1].x=XN;
}
P[N].x=XN;
}
int cmp2(const void *i,const void *j)
{
point *ii=(point *)i,*jj=(point *)j;
return ii->y-jj->y;
}
void getYN_FF()
{
int i;
qsort(P+1,N,sizeof(P[1]),cmp2);
for(i=1;i<=N;i++)
{
if(P[i].y!=P[i-1].y)
{
YN++;
P[i-1].y=YN-1;
}
else P[i-1].y=YN;
FF[P[i-1].x][P[i-1].y]++;
}
P[N].y=YN;
FF[P[N].x][P[N].y]++;
}
void get_UPCE()
{
for(int i=1;i<=XN;i++)
for(int j=1;j<=YN;j++)
{
CE[i][j]=CE[i][j-1]+FF[i][j];
UP[i][j]=UP[i-1][j]+FF[i][j];
}
}
int On[maxn],Left[maxn],On2[maxn];
void get_On_Left_On2(int a,int b)
{
memset(On,0,sizeof(On));
memset(Left,0,sizeof(Left));
memset(On2,0,sizeof(On2));
for(int i=1;i<=YN;i++)
{
On[i]=UP[b-1][i]-UP[a][i];
On2[i]=UP[b][i]-UP[a-1][i];
Left[i]=CE[a][i-1]+CE[b][i-1];
}
}
void solve()
{
ans=0;
int tempans=0,tempi=0;
for(int i=1;i<=XN;i++)
for(int i2=i+1;i2<=XN;i2++)
{
get_On_Left_On2(i,i2);
tempans=0;
for(int j=1;j<=YN;j++)
{
if(ans<Left[j]+On2[j]+tempans)
ans=Left[j]+On2[j]+tempans;
if(tempans<On[j]-Left[j])
tempans=On[j]-Left[j];
}
}
}
int main()
{
freopen("a.in","r",stdin);
freopen("a.out","w",stdout);
int Case=0;
while(cin>>N&&N)
{
Case++;
input();
if(N<=3) {
printf("Case %d: %d\n",Case,N);
continue;
}
getXN();
getYN_FF();
get_UPCE();
solve();
printf("Case %d: %d\n",Case,ans);
}
}

【巧妙预处理系列+离散化处理】【uva1382】Distant Galaxy的更多相关文章

  1. 【巧妙预处理系列】【UVA1330】City game

    最大子矩阵(City Game, SEERC 2004, LA 3029) 给定一个m×n的矩阵,其中一些格子是空地(F),其他是障碍(R).找出一个全部由F组成的面积最大的子矩阵,输出其面积乘以3后 ...

  2. hdu Distant Galaxy(遥远的银河)

    Distant Galaxy Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)To ...

  3. LA3695 Distant Galaxy

    Distant Galaxy https://vjudge.net/problem/UVALive-3695 You are observing a distant galaxy using a te ...

  4. 【poj3141】 Distant Galaxy

    http://poj.org/problem?id=3141 (题目链接) 题意 给出平面上n个点,找出一个矩形,使边界上包含尽量多的点. solution 不难发现,除非所有输入点都在同一行或同一列 ...

  5. uva 1382 - Distant Galaxy

    题目连接:http://acm.hust.edu.cn/vjudge/contest/view.action?cid=91208#problem/G 题意:  给出平面上的n个点,找出一个矩形,使得边 ...

  6. LA 3695 Distant Galaxy

    给出n个点的坐标(坐标均为正数),求最多有多少点能同在一个矩形的边界上. 题解里是构造了这样的几个数组,图中表示的很明白了. 首先枚举两条水平线,然后left[i]表示竖线i左边位于水平线上的点,on ...

  7. 【巧妙算法系列】【UVA 11384】 Help is needed for Dexter 正整数序列

    Help is needed for Dexter Time Limit: 3 Second Dexter is tired of Dee Dee. So he decided to keep Dee ...

  8. 【巧妙算法系列】【Uva 11464】 - Even Parity 偶数矩阵

    偶数矩阵(Even Parity, UVa 11464) 给你一个n×n的01矩阵(每个元素非0即1),你的任务是把尽量少的0变成1,使得每个元素的上.下.左.右的元素(如果存在的话)之和均为偶数.比 ...

  9. UVa LA 3695 - Distant Galaxy 前缀和,状态拆分,动态规划 难度: 2

    题目 https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_pr ...

随机推荐

  1. 期望dp-hdu-4336-Card Collector

    题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=4336 题目大意: 有n种卡片,每包中至多有一种卡片,概率分别为p1,p2,...pn,可能有的没有卡 ...

  2. Xcode5和6上新建工程如何本地化启动页面

    建议阅读本篇文章前先具备iOS本地化的基本知识,Google中搜索“iOS本地化”,有成片的教程~~ 最近有个app需要支持英语.简体中文.繁体中文,由于启动页面上有文字,所以也不得不做下本地化处理. ...

  3. html5 音频

    目前,web页面上没有标准的方式来播放音频文件,大多数的音频文件是使用插件来播放,而众多的浏览器使用了不同的插件.而html5的到来,给我们提供了一个标准的方式来播放web中音频文件,用户不再为浏览器 ...

  4. Oracle更改字符集

    更改oracle的字符集: sqlplus / as sysdba SQL> shutdown immediate; Database closed. Database dismounted. ...

  5. VS2010程序打包操作(超详细的)

    1.  在vs2010 选择“新建项目”----“其他项目类型”----“Visual Studio Installerà“安装项目”: 命名为:Setup1 . 这是在VS2010中将有三个文件夹, ...

  6. CSS 规划

    1)在拿到psd效果图后,我会将所有的效果图看一遍,然后开始打草稿,在纸上画,在PS中标注都可以.一般都是做个大概的统计,太详细不太现实,既费时间也费精力. 2)网站的浏览器兼容性要顾及到哪些,IE6 ...

  7. Emacs入门快捷键

    打开emacs开始一个程序最基本操作: 1.打开Emacs,执行 $ emacs 2.建立一个新的程序文件. 按C-x C-f 然后在屏幕的底部出现minibuffer,光标提示你输入文件名称, 文件 ...

  8. QT TCP/IP

    QT 网络通信(TCP/IP) 服务端: 一.监听新的客户端接入(QTcpServer) 重写函数 incomingConnection(qintptr socketDescriptor) 二.服务端 ...

  9. [C++程序设计]返回指针值的函数

    定义指针函数的一般形式为 类型名 *函数名(参数表列); 例如 int *a(int x,int y);

  10. Python爬虫学习:四、headers和data的获取

    之前在学习爬虫时,偶尔会遇到一些问题是有些网站需要登录后才能爬取内容,有的网站会识别是否是由浏览器发出的请求. 一.headers的获取 就以博客园的首页为例:http://www.cnblogs.c ...