BC 2015年百度之星程序设计大赛 - 初赛(1)(矩形面积-旋转卡壳)

矩形面积

Accepts: 717

Submissions: 1619

Time Limit: 2000/1000 MS (Java/Others)

Memory Limit: 32768/32768 K (Java/Others)

Problem Description

小度熊有一个桌面,小度熊剪了很多矩形放在桌面上,小度熊想知道能把这些矩形包围起来的面积最小的矩形的面积是多少。

Input

第一行一个正整数 T,代表测试数据组数(1≤T≤20 ),接下来
T 组测试数据。

每组测试数据占若干行,第一行一个正整数 N(1≤N<≤1000) ,代表矩形的数量。接下来
N 行,每行 8 个整数x 1 ,y 1 ,x 2 ,y 2 ,x 3 ,y 3 ,x 4 ,y 4  ,代表矩形的四个点坐标,坐标绝对值不会超过10000。

Output

对于每组测试数据,输出两行:

第一行输出"Case #i:",i 代表第 i 组测试数据。 第二行包含1 个数字,代表面积最小的矩形的面积,结果保留到整数位。

Sample Input

2
2
5 10 5 8 3 10 3 8
8 8 8 6 7 8 7 6
1
0 0 2 2 2 0 0 2

Sample Output

Case #1:
17
Case #2:
4

旋转卡壳,和bzoj上那题基本相同

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<functional>
#include<cmath>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=next[p])
#define MAXN (10000+10)
#define INF (1000000000)
#define eps 1e-6
struct P
{
    double x,y;
    P(){}
    P(double _x,double _y):x(_x),y(_y){}
    friend bool operator<(P a,P b){return (fabs(a.y-b.y)<eps)?a.x<b.x:a.y<b.y;  }
    friend bool operator==(P a,P b){return fabs(a.x-b.x)<eps&&fabs(a.y-b.y)<eps;}
    friend bool operator!=(P a,P b){return !(a==b);}  

}a[MAXN],s[MAXN],ansp[5];
int size=0;
double ans=INF;
double dis2(P a,P b){return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);}
struct V
{
    double x,y;
    V(){}
    V(double _x,double _y):x(_x),y(_y){}
    V(P a,P b):x(b.x-a.x),y(b.y-a.y){}
    friend double operator*(V a,V b){return a.x*b.y-a.y*b.x;}
    friend V operator*(double a,V b){return V(a*b.x,a*b.y);}
    friend double operator/(V a,V b){return a.x*b.x+a.y*b.y;}
    friend P operator+(P a,V b){return P(a.x+b.x,a.y+b.y);}
    friend P operator-(P a,V b){return P(a.x-b.x,a.y-b.y);}
    friend V operator~(V a){return V(a.y,-a.x);}
    double dis2(){return x*x+y*y;   }
};
int cmp(P A,P B)
{
    double tmp=V(a[1],A)*V(a[1],B);
    if (tmp>0) return 1;
    else if (fabs(tmp)<eps) return (-dis2(A,a[1])-dis2(B,a[1])>0);
    return 0;
}
int n;  

int main()
{
//	freopen("F.in","r",stdin);

	int T;
	cin>>T;
	For(kcase,T)
	{
		ans=INF;
	    scanf("%d",&n);
	    n*=4;
	    for(int i=1;i<=n;i++) scanf("%lf%lf",&a[i].x,&a[i].y);
	    for (int i=2;i<=n;i++) if (a[i]<a[1]) swap(a[1],a[i]);
	    sort(a+2,a+1+n,cmp);
	    s[1]=a[1];size=1;
	    for (int i=2;i<=n;)
	        if (size<2||V(s[size-1],s[size])*V(s[size],a[i])>eps) s[++size]=a[i++];
	        else size--;
	    s[0]=s[size];  

	    int l=1,r=1,t=1;
	    for (int i=0;i<size;i++)
	    {
	        while (V(s[i],s[i+1])*V(s[i],s[t+1])-V(s[i],s[i+1])*V(s[i],s[t])>-eps) t=(t+1)%size;
	        while (V(s[i],s[i+1])/V(s[i],s[r+1])-V(s[i],s[i+1])/V(s[i],s[r])>-eps) r=(r+1)%size;
	        if (i==0) l=r;
	        while (V(s[i],s[i+1])/V(s[i],s[l+1])-V(s[i],s[i+1])/V(s[i],s[l])<eps) l=(l+1)%size;
	        double Dis2=dis2(s[i],s[i+1]),wlxdis=V(s[i],s[i+1])/V(s[i],s[l]),wrxdis=V(s[i],s[i+1])/V(s[i],s[r]),hxdis=V(s[i],s[i+1])*V(s[i],s[t]);
	        double tmp=hxdis*(wrxdis-wlxdis)/Dis2;
	        if (tmp<0) tmp=-tmp;
	        if (ans>tmp)
	        {
	            ans=tmp;
	            ansp[0]=s[i]-(wlxdis/Dis2)*V(s[i+1],s[i]);
	            ansp[1]=s[i]+(wrxdis/Dis2)*V(s[i],s[i+1]);
	            ansp[2]=ansp[1]+(hxdis/Dis2)*(~V(s[i+1],s[i]));
	            ansp[3]=ansp[0]+(hxdis/Dis2)*(~V(s[i+1],s[i]));
	        }
	    }
	    int p=0;
	    for (int i=1;i<4;i++) if (ansp[i]<ansp[p]) p=i;//p=0;
		printf("Case #%d:\n",kcase);
	    printf("%.0lf\n",ans);    

	}
    return 0;
}  
时间: 2024-10-09 07:28:48

BC 2015年百度之星程序设计大赛 - 初赛(1)(矩形面积-旋转卡壳)的相关文章

BC 2015年百度之星程序设计大赛 - 初赛(1)(序列变换-二分答案贪心)

序列变换 Accepts: 816 Submissions: 3578 Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Problem Description 给定序列A={A 1 ,A 2 ,...,A n } , 要求改变序列A中的某些元素,形成一个严格单调的序列B(严格单调的定义为:B i <B i+1 ,1≤i<N ). 我们定义从序列A到序列B变换的代价为cost(A,B

BC 2015年百度之星程序设计大赛 - 初赛(1)(KPI-主席树)

KPI Accepts: 517 Submissions: 2185 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Problem Description 你工作以后, KPI 就是你的全部了. 我开发了一个服务,取得了很大的知名度.数十亿的请求被推到一个大管道后同时服务从管头拉取请求.让我们来定义每个请求都有一个重要值.我的KPI是由当前管道内请求的重要值的中间值来计算.现在给你

hdu 5253 连接的管道(kruskal)(2015年百度之星程序设计大赛 - 初赛(2))

连接的管道 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 1323    Accepted Submission(s): 519 Problem Description 老 Jack 有一片农田.以往几年都是靠天吃饭的. 可是今年老天格外的不开眼.大旱.所以老 Jack 决定用管道将他的全部相邻的农田全部都串联起来.这样他就能够从远处

2015年百度之星程序设计大赛 - 初赛(1)

今天焊收音机,没来得及做 题目出的还是比较有搞头的,也很有意思 首先第一题,没想法 然后看了一下正确率最高的题 求包围所有点的最小矩形,显然是个凸包问题,很久没写过几何,直接找了个模板套了一下 hdu 5251 题解 然后照着校友的顺序开始做题,二分无误,就是少了个零,wa了好多发,明天第二场一定要先开个小号 hdu 5248 题解 接着又回到第一题,找到规律 hdu 5246 题解 然后就24:00了,还要补番,第二天还要继续焊收音机,于是等下次再做啦

2015年百度之星程序设计大赛 - 初赛(2)

目前只会四道--> 1002:很明显的MST 把,但是注意不要每个点都连接四条边,会爆内存. 1 #include <iostream> 2 #include <cstdio> 3 #include <cstring> 4 #include<cmath> 5 #include <algorithm> 6 7 typedef long long ll; 8 using namespace std; 9 10 #define N 1011 11

2015年百度之星程序设计大赛 - 资格赛:IP聚合

1 /* 2 pro:2015年百度之星程序设计大赛 - 资格赛:IP聚合 3 ID:forever95 4 Date:2015-5-25 5 */ 6 //事实证明C++就是一坨屎 7 #include <iostream> 8 #include <string> 9 #include <cstdio> 10 #include <cstdlib> 11 #include <istream> 12 #include <algorithm&g

HDU 6119 小小粉丝度度熊 【预处理+尺取法】(2017&quot;百度之星&quot;程序设计大赛 - 初赛(B))

小小粉丝度度熊 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 1572    Accepted Submission(s): 513 Problem Description 度度熊喜欢着喵哈哈村的大明星--星星小姐. 为什么度度熊会喜欢星星小姐呢? 首先星星小姐笑起来非常动人,其次星星小姐唱歌也非常好听. 但这都不是最重要的,最重要的是

2014年百度之星程序设计大赛 - 初赛(第一轮) hdu Grids (卡特兰数 大数除法取余 扩展gcd)

题目链接 分析:打表以后就能发现时卡特兰数, 但是有除法取余. f[i] = f[i-1]*(4*i - 2)/(i+1); 看了一下网上的题解,照着题解写了下面的代码,不过还是不明白,为什么用扩展gcd, 不是用逆元吗.. 网上还有别人的解释,没看懂,贴一下: (a / b) % m = ( a % (m*b)) / b 笔者注:鉴于ACM题目特别喜欢M=1000000007,为质数: 当gcd(b,m) = 1, 有性质: (a/b)%m = (a*b^-1)%m, 其中b^-1是b模m的逆

2014年百度之星程序设计大赛 - 初赛(第二轮)

1001 暴力 #include <cstdio> #include <algorithm> #include <cstring> using namespace std; const int maxn = 100100; int ll[maxn], rr[maxn]; struct node { int x, y, bj; }e[maxn]; int main() { int cas = 1; int T; scanf("%d", &T);