题意

题解

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

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

代码

#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 (凸包)相关推荐

  1. BZOJ 4388 [JOI2012春季合宿]Invitation (线段树、二叉堆、最小生成树)

    题目链接 https://www.lydsy.com/JudgeOnline/problem.php?id=4388 题解 模拟Prim算法? 原题所述的过程就是Prim算法求最大生成树的过程.于是我 ...

  2. BZOJ 4221 [JOI2012春季合宿]Kangaroo (DP)

    题目链接 https://www.lydsy.com/JudgeOnline/problem.php?id=4221 题解 orz WYC 爆切神仙DP 首先将所有袋鼠按大小排序.考虑从前往后DP, ...

  3. [JOI2012春季合宿]Rotate (链表)

    题意 题解 又是一道神仙题-- 显然的做法是大力splay,时间复杂度\(O((N+Q)N\log N)\), 可以卡掉. 正解: 使用十字链表维护矩阵,在周围增加第\(0\)行/列和第\((n+1) ...

  4. JOI-2016/17 春季合宿 切题记

    17年的合宿好难啊...感觉是我做过的最难的一套题(没有之一)了...但是可能也是价值最高的? Day1: T1 Cultivation:给你一个H*W的网格,有N<=300棵仙人掌.每一年可以 ...

  5. 【JOI2017春季合宿】Port Facility

    http://uoj.ac/problem/356 题解 思路和\(NOIP\)双栈排序差不多. 对于两个元素,若\(l_1<l_2<r_1<r_2\)那么它们不能在一个栈里,我们连 ...

  6. LOJ #2731 [JOI2016春季合宿]Solitaire (DP、组合计数)

    题目链接 https://loj.ac/problem/2731 题解 首先一个很自然的思路是,设\(dp[i][j]\)表示选了前\(i\)列,第\(2\)行第\(i\)列的格子是第\(j\)个被填 ...

  7. LOJ #2733 [JOI2016春季合宿]Sandwiches (DP)

    题目链接 https://loj.ac/problem/2733 题解 神仙题-- 首先可以观察到一个结论: 目标块的两块小三明治一定分别是最后和倒数第二个被吃的. 由此我们可以考虑这两块谁先被吃.这 ...

  8. LOJ #2734 Luogu P3615 [JOI2016春季合宿]Toilets (结论、贪心)

    题目链接 (loj) https://loj.ac/problem/2734 (luogu) https://www.luogu.org/problem/P3615 题解 嗯,考场上肝了\(3h\)然 ...

  9. JOI2017 春季合宿:Railway Trip

    自己的AC做法似乎离正解偏了十万八千里而且复杂了不少--不管怎样还是记录下来吧. 题意: 题目链接: JOISC2017 F - AtCoder JOISC2017 F - LOJ \(N\)个车站排 ...

最新文章

  1. Android问题汇总
  2. MongoDB Wiredtiger存储引擎实现原理——Copy on write的方式管理修改操作,Btree cache...
  3. easyUI的combobox选中无法显示
  4. dat文件打开乱码_5.2 实战1:解决在Linux下打开Windows汉字文本的乱码问题
  5. ng-repeat循环出来的部分调用同一个函数并且实现每个模块之间不能相互干扰
  6. oracle与sqlserver比较
  7. MySQL 数据库系统
  8. php抽象类初始化方法,php – 抽象类方法声明
  9. CF932E Team Work
  10. C++中static关键字用法详解
  11. JAAS(Java 认证和授权服务)
  12. Java主流开发框架
  13. Strategy(策略)
  14. java汉字获取拼音的方法
  15. Js中slite,splice,slice,join方法的用法区别
  16. 淹没之城(Submerged).PC单机游戏 免费下载.虚幻4 打造
  17. 比较快得Maven镜像
  18. Qt编写地图综合应用60-覆盖物坐标和搜索
  19. 超融合云数据中心,泽塔云为北京城市副中心建设添砖加瓦
  20. 数据挖掘 第四篇:OLS回归分析

热门文章

  1. 帝国CMS后台前台上传大小和类型限制的解决方法!
  2. 推荐!计算机视觉最适合入门的 8 本教程,算法与实战兼备
  3. oracle恢复某个表的数据库,如何从rman的全备份中恢复一张表?
  4. Ubuntu访问Windows共享文件夹
  5. wireshrk中的名词说明
  6. JAVA——Map基础知识
  7. 任务和特权级保护(一)——《x86汇编语言:从实模式到保护模式》读书笔记27
  8. C++ 11 多线程--线程管理
  9. es Root mapping definition has unsupported parameters解决方法
  10. DM8168学习--USB的over-current 问题总结