Min-Max
Time Limit: 1000MS   Memory Limit: 65536K
Total Submissions: 2192   Accepted: 502

Description

Define the following function




Given the value C of F(p1, p2 ... pn), can you find the minimum and maximum value of F(q1, q2 ... qn)?

Input

The input contains several test cases. For each test case, it contains three lines.



Line 1: two integers n (1<= n <= 50000) and C.

Line 2: n integers p1, p2 ... pn (|pi| < 1000 for 1 <= i <= n).

Line 3: n integers q1, q2 ... qn (|qi| < 1000 for 1 <= i <= n).

Output

For each test case, output the minimum and maximum value in a single line with the fraction rounded to 3 decimal places.

Sample Input

2 1
3 1
0 2

Sample Output

2.000 2.000

Source

题目一看就很熟悉很熟悉很熟悉啊!

我怎么一眼看到就想起詹森不等式呢。

。。。

事实上是重心公式。C是重心的x,而y在凸包上。

所以问题就变成求凸包啦!

求完凸包就要绕着凸包上每一条边求极值并更新。

/***********************************************************
> OS : Linux 3.13.0-24-generic (Mint-17)
> Author : yaolong
> Mail : dengyaolong@yeah.net
> Time : 2014年10月14日 星期二 07时44分53秒
**********************************************************/
#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const double INF = 1e50;
const int N = 61111;
struct Point
{
int x, y;
} ;
Point p[N], stk[N];
Point minp;
int top;
double cross ( Point &o, Point &a, Point &b )
{
return ( a.x - o.x ) * ( b.y - o.y ) - ( a.y - o.y ) * ( b.x - o.x );
}
double dist ( Point &A, Point & B )
{
return hypot ( A.x - B.x, A.y - B.y );
}
bool cmp ( Point A, Point B )
{
double k = cross ( minp, A, B );
if ( k < 0 ) return 0;
if ( k > 0 ) return 1;
return dist ( minp, A ) < dist ( minp, B );
}
void Gramham ( int n )
{
int i;
for ( i = 1; i < n; i++ )
{
if ( p[i].y < p[0].y || ( p[i].y == p[0].y && p[i].x < p[0].x ) )
{
swap ( p[i], p[0] );
}
}
minp = p[0];
p[n] = p[0];
sort ( p + 1, p + n, cmp );
stk[0] = p[0];
stk[1] = p[1];
top = 1;
for ( i = 2; i < n; i++ )
{
while ( top >= 1 && cross ( stk[top - 1], stk[top ], p[i] ) <= 0 ) --top;
stk[++top] = p[i];
}
}
double mmin, mmax;
int c;
void update ( Point a, Point b )
{
if ( a.x > b.x )
{
swap ( a, b );
}
if ( a.x <= c && b.x >= c )
{
if ( a.x == c && b.x == c )
{
mmax = max ( mmax, ( double ) max ( a.y, b.y ) );
mmin = min ( mmin, ( double ) min ( a.y, b.y ) );
}
else
{
double k = ( ( double ) c - a.x ) / ( ( double ) b.x - a.x ) * ( b.y - a.y ) + a.y;
mmax = max ( mmax, k );
mmin = min ( mmin, k );
}
}
}
int main()
{
int n, i;
while ( ~scanf ( "%d%d", &n, &c ) )
{
for ( i = 0; i < n; i++ )
{
scanf ( "%d", &p[i].x );
}
for ( i = 0; i < n; i++ )
{
scanf ( "%d", &p[i].y );
}
if ( n == 1 )
{
printf ( "%.3f %.3f\n", ( double ) p[0].y, ( double ) p[0].y );
continue;
}
Gramham ( n );
stk[++top] = stk[0];
mmin = INF, mmax = -INF;
for ( i = 1; i <= top; i++ )
{
update ( stk[i - 1], stk[i] );
}
printf ( "%.3f %.3f\n", mmin, mmax );
}
return 0;
}

poj2595(凸包)的更多相关文章

  1. [poj1113][Wall] (水平序+graham算法 求凸包)

    Description Once upon a time there was a greedy King who ordered his chief Architect to build a wall ...

  2. ZOJ 3871 Convex Hull(计算几何、凸包)

    题意:给n个点,|x[i]|,|y[i]| <= 1e9.求在所有情况下的子集下(子集点数>=3),凸包的面积和. 这题主要有几个方面,一个是凸包的面积,可以直接用线段的有向面积和求得,这 ...

  3. UVALive 2453 Wall (凸包)

    题意:给你一个多边形的城堡(多个点),使用最短周长的城墙将这个城堡围起来并保证城墙的每个点到城堡上的每个点的距离都不小于l 题解:因为两点间的直线一定比折线短,所以这样做 先使用所有点求得一个凸包,接 ...

  4. UVA 11168 Airport(凸包+直线方程)

    题意:给你n[1,10000]个点,求出一条直线,让所有的点都在都在直线的一侧并且到直线的距离总和最小,输出最小平均值(最小值除以点数) 题解:根据题意可以知道任意角度画一条直线(所有点都在一边),然 ...

  5. 关于2016.12.12——T1的反思:凸包的意义与应用

    2016.12.12 T1 给n个圆,保证圆圆相离,求将圆围起来的最小周长.n<=100 就像上图.考场上,我就想用切线的角度来做凸包.以圆心x,y排序,像点凸包一样,不过用两圆之间的下切线角度 ...

  6. poj1228--稳定凸包

    题目大意:给你一个凸包上的某些点(可能在凸包内),询问是否能确定这个凸包. 思路:先求出题目给出的点的凸包,看看在凸包的每条边内(不包括端点)有没有点,若有,则这条边是确定的,若没有,则这条边不确定, ...

  7. POJ 2225 / ZOJ 1438 / UVA 1438 Asteroids --三维凸包,求多面体重心

    题意: 两个凸多面体,可以任意摆放,最多贴着,问他们重心的最短距离. 解法: 由于给出的是凸多面体,先构出两个三维凸包,再求其重心,求重心仿照求三角形重心的方式,然后再求两个多面体的重心到每个多面体的 ...

  8. HDU 4946 Area of Mushroom(构造凸包)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4946 题目大意:在一个平面上有n个点p1,p2,p3,p4....pn,每个点可以以v的速度在平面上移 ...

  9. HDU 1392 凸包模板题,求凸包周长

    1.HDU 1392 Surround the Trees 2.题意:就是求凸包周长 3.总结:第一次做计算几何,没办法,还是看了大牛的博客 #include<iostream> #inc ...

随机推荐

  1. [android开发篇]使用系统权限流程

    1 声明权限https://developer.android.com/training/permissions/declaring.html 每款 Android 应用都在访问受限的沙盒中运行.如果 ...

  2. [android 应用框架api篇] bluetooth

    bluetooth接口 android.bluetooth.jar 官网网址: 接口类: https://developer.android.com/reference/android/bluetoo ...

  3. 【Luogu】P1417烹调方案(排序01背包)

    题目链接 对食材进行排序,重载运算符代码如下: struct food{ long long a,b,c; bool operator <(const food &a)const{ re ...

  4. hihoCoder #1471 拥堵的城市

    这道题目是hihoCoder Challenge 27的C题,我考虑了5天:(. 计数问题.由于树的结构的特殊性(树具有递归结构),不难想到思路是树形DP.由于这是[计数问题]而非[优化问题],我们思 ...

  5. BZOJ 2820 YY的GCD ——莫比乌斯反演

    我们可以枚举每一个质数,那么答案就是 $\sum_{p}\sum_{d<=n}\mu(d)*\lfloor n / pd \rfloor *\lfloor m / pd \rfloor$ 直接做 ...

  6. SPOJ NSUBSTR Substrings ——后缀自动机

    建后缀自动机 然后统计次数,只需要算出right集合的大小即可, 然后更新f[l[i]]和rit[i]取个max 然后根据rit集合短的一定包含长的的性质,从后往前更新一遍即可 #include &l ...

  7. 算法复习——费用流模板(poj2135)

    题目: Farm Tour Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 16898   Accepted: 6543 De ...

  8. 陌上花开(bzoj 3262)

    Description 有n朵花,每朵花有三个属性:花形(s).颜色(c).气味(m),又三个整数表示.现要对每朵花评级,一朵花的级别是它拥有的美丽能超过的花的数量.定义一朵花A比另一朵花B要美丽,当 ...

  9. Mysql字符集与校对规则

    字符集是一套字符和编码的集合,校对规则是用于比较字符集的一套规则. 所以字符集有两部分组成字符集合和对应的编码集合.比如说,现在有这几个字符:A B a b, 假设它们对应的编码分别是00, 01, ...

  10. 【Git】Git 本地的撤销修改和删除操作

    一:撤销操作 比如我现在在readme.txt文件里面增加一行 内容为555555555555,我们先通过命令查看如下: 在我未提交之前,我发现添加5555555555555内容有误,所以我得马上恢复 ...