[JOI2012春季合宿]Constellation (凸包)

2024-02-15 15:38

本文主要是介绍[JOI2012春季合宿]Constellation (凸包),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题意

1449777-20191022193129979-160536879.png

1449777-20191022193140518-820091949.png

题解

神仙结论题。
结论: 一个点集合法当且仅当其凸包上的两种颜色点分别连续。
证明:
必要性显然。
充分性: 考虑对于一个不同色三角形\(ABC\),不妨设点\(A\)为白点,点\(B,C\)为黑点。若形内无白点,则随便连,显然成立。若形内有白点,则任取一白点\(S\), 对三角形\(SBC,BAS,CAS\)内部的点分别连边(递归构造),最后连接\(SA\).
再考虑一个凸包,设在逆时针方向上最后一个白点是\(U\), 最后一个黑点是\(V\), 则连接\(UV\), 把凸包剖成两个部分,一部分只有\(V\)是黑,一部分只有\(U\)是白。然后进行如图所示三角剖分,每个三角形内连上即可。
1449777-20191022194452336-582015186.png

知道了结论,这题就很简单了。
勿忘非空。

代码

#include<bits/stdc++.h>
#define llong long long
using namespace std;const int N = 1e5;
const int P = 1e9+7;
struct Point
{llong x,y; int c,id;Point() {}Point(int _x,int _y) {x = _x,y = _y;}
} a[N+3],ch[N+3];
bool f[N+3];
Point operator -(Point x,Point y) {return Point(x.x-y.x,x.y-y.y);}
llong Cross(Point x,Point y) {return x.x*y.y-x.y*y.x;}
int n,sz;bool cmp_ang(Point x,Point y) {return Cross(x-a[1],y-a[1])>0;}int main()
{scanf("%d",&n);for(int i=1; i<=n; i++) scanf("%lld%lld%d",&a[i].x,&a[i].y,&a[i].c);for(int i=2; i<=n; i++) {if(a[i].x<a[1].x||(a[i].x==a[1].x&&a[i].y<a[1].y)) {swap(a[i],a[1]);}}sort(a+2,a+n+1,cmp_ang);for(int i=1; i<=n; i++) a[i].id = i;sz = 1; ch[1] = a[1];for(int i=2; i<=n; i++){while(sz>=2 && Cross(ch[sz]-ch[sz-1],a[i]-ch[sz-1])<0) {sz--;}sz++; ch[sz] = a[i];}
//  printf("ch: "); for(int i=1; i<=sz; i++) printf("(%lld,%lld) ",ch[i].x,ch[i].y); puts("");int mx1 = 0,mn1 = n,mx2 = 0,mn2 = n; llong ans;for(int i=1; i<=sz; i++){if(ch[i].c==1) {mx1 = max(mx1,i),mn1 = min(mn1,i);}else if(ch[i].c==2) {mx2 = max(mx2,i),mn2 = min(mn2,i);}}
//  printf("sz=%d mn1=%d mx1=%d mn2=%d mx2=%d\n",sz,mn1,mx1,mn2,mx2);if(mx1<mn1 && mx2<mn2){ans = (1ll*sz*sz-sz+2ll)%P;}else if(mx1<mn1) //only 2{ans = 1ll;for(int i=mn2+1,j=0; i<=mn2+sz; i++){int ii = i>sz?i-sz:i;if(ch[ii].c==0) {j++;}else {ans += 1ll*j*(j+1ll)/2ll; j = 0;}}ans %= P;}else if(mx2<mn2) //only 1{ans = 1ll;for(int i=mn1+1,j=0; i<=mn1+sz; i++){int ii = i>sz?i-sz:i;if(ch[ii].c==0) {j++;}else {ans += 1ll*j*(j+1ll)/2ll; j = 0;}}ans %= P;}else{if(mn1>mn2){for(int i=mn1; i<=mx1; i++){if(ch[i].c==2) {puts("0"); return 0;}}int l = -1,r = -1;for(int i=mx1; ch[i].c!=2; i=(i==sz?1:i+1)){l++;}for(int i=mn1; ch[i].c!=2; i=(i==1?sz:i-1)){r++;}ans = 1ll*(l+1ll)*(r+1ll)%P;}else{for(int i=mn2; i<=mx2; i++){if(ch[i].c==1) {puts("0"); return 0;}}int l = -1,r = -1;for(int i=mx2; ch[i].c!=1; i=(i==sz?1:i+1)){l++;}for(int i=mn2; ch[i].c!=1; i=(i==1?sz:i-1)){r++;}ans = 1ll*(l+1ll)*(r+1ll)%P;}}for(int i=1; i<=sz; i++){f[ch[i].id] = true;}int cnt1 = 0,cnt2 = 0;for(int i=1; i<=n; i++){if(f[i]==false && a[i].c==0) {ans=(ans<<1)%P;}if(a[i].c==1) cnt1++;if(a[i].c==2) cnt2++;}if(cnt1==0) {ans--;} if(cnt2==0) {ans--;}ans = (ans+P)%P;printf("%lld\n",ans);return 0;
}

这篇关于[JOI2012春季合宿]Constellation (凸包)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/711801

相关文章

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

poj 2187 凸包or旋转qia壳法

题意: 给n(50000)个点,求这些点与点之间距离最大的距离。 解析: 先求凸包然后暴力。 或者旋转卡壳大法。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <s

Codeforces Beta Round #47 C凸包 (最终写法)

题意慢慢看。 typedef long long LL ;int cmp(double x){if(fabs(x) < 1e-8) return 0 ;return x > 0 ? 1 : -1 ;}struct point{double x , y ;point(){}point(double _x , double _y):x(_x) , y(_y){}point op

HDU 1392 HDU 1348 凸包

求凸包的周长,  注意n=1 , 2时特殊情况 int cmp(double x){if(fabs(x) < 1e-8) return 0 ;if(x > 0) return 1 ;return -1 ;}struct point{double x , y ;point(){}point(double _x , double _y):x(_x) , y(_y){}frien

【UVA】10652-Board Wrapping(凸包问题)

又增加了2个模板。 #include <iostream>#include <cstdlib>#include <cstdio>#include <string>#include <cstring>#include <cmath>#include <vector>#include <queue>#include <stack>#include <algorithm>usi

百度之星初赛1006(计算几何:能包含凸包的最小矩形面积)

矩形面积    Accepts: 717    Submissions: 1619  Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Problem Description 小度熊有一个桌面,小度熊剪了很多矩形放在桌面上,小度熊想知道能把这些

OpenCV结构分析与形状描述符(8)点集凸包计算函数convexHull()的使用

操作系统:ubuntu22.04 OpenCV版本:OpenCV4.9 IDE:Visual Studio Code 编程语言:C++11 算法描述 查找一个点集的凸包。 函数 cv::convexHull 使用斯克拉斯基算法(Sklansky’s algorithm)来查找一个二维点集的凸包,在当前实现中该算法的时间复杂度为 O(N logN)。 函数 cv::convexHull 是

OpenCV结构分析与形状描述符(9)检测轮廓相对于其凸包的凹陷缺陷函数convexityDefects()的使用

操作系统:ubuntu22.04 OpenCV版本:OpenCV4.9 IDE:Visual Studio Code 编程语言:C++11 算法描述 查找一个轮廓的凸性缺陷。 下图显示了一个手部轮廓的凸性缺陷: convexityDefects 是 OpenCV 库中的一个函数,用于检测轮廓相对于其凸包的凹陷缺陷。这个函数可以帮助识别轮廓中的凹进去的部分,通常被用来分析手部或其他物体的形状

分治算法与凸包问题

1. 什么是凸包问题? 凸包问题是计算几何中的经典问题。给定二维平面上的点集,凸包是一个最小的凸多边形,它包含了点集中所有的点。你可以把凸包想象成一根松紧带将所有点紧紧包裹住的样子,凸包的边缘仅沿着最外面的点延伸。 2. 分治法简介 分治算法是解决复杂问题的强大策略,它的思想是将问题分解为多个子问题,分别解决这些子问题后再合并得到最终解。凸包问题可以通过分治算法高效地解决,时间复杂度可以达到

多边形快速凸包算法(Melkman‘s Algorithm)

前言 平面点集的凸包算法一文介绍了如何计算平面点集或者任意多边形的凸包。对于随机的平面点集,Graham scan和Andraw's 单调链算法已经是最快的算法了。但是对于没有自相交的封闭的简单多边形,存在线性复杂度的算法。下面介绍这一优雅高效的算法。 一般的2D凸包算法,首先将点进行排序(时间复杂度),然后利用栈操作在O(n)的时间复杂度内计算凸包。初始的排序决定了最终的时间复杂度。但是本文