[xsy1100]东舰停战不可避
没有三点共线
这题的思想来源于JOI2011-2012春季训练合宿Day2T2,原题是个大毒瘤题(p.s.场上有人A,真的可怕),这题作为原题要用到的的一个结论而存在
点有两种颜色,先考虑对所有点做凸包,假如凸包上的颜色段大于$2$段,那么无解
因为连线一定在凸包内部,所以上图中红色和蓝色不可能不相交而连接起来
为了做这题,首先我们描述一个过程:$solve(A,B,C)$,其中$B,C$颜色相同且已经直接或间接地被连起来,$A$是另一种颜色,这个过程可以把$\triangle ABC$内的所有点按对应颜色连起来且连线不相交
solve($A$,$B$,$C$){
if($\exists P$ in $\triangle ABC$ and $P.color$==$A.color$){
link($A$,$P$);
solve($B$,$A$,$P$);
solve($C$,$A$,$P$);
solve($P$,$B$,$C$);
}else{
link all $P$ in $\triangle ABC$ with $B$ or $C$
}
}
代码应该挺容易懂的,如果找到和$A$同色的点就可以递归做,否则把所有三角形内的点用某种方式与$B,C$连接起来,我用的方法是把所有点极角排序,然后把每个点和左下角的点连起来
要处理两种情况
①凸包上所有点同色
把凸包上连一圈,然后再内部随便找一个不同颜色的点$P$,以它为中心绕着调用$solve(P,hull_i),hull_{i+1}$
②凸包上有两段颜色
把凸包上相同颜色的连起来,然后如图所示
就是找到每段颜色的第一个点和另一边的每一条边组成的三角形调用$solve$就可以了
因为没有三点共线,所以这样做一定可以得出答案
东舰共荣还行==
#include<stdio.h> #include<algorithm> using namespace std; #define ll long long #define inf 2147483647 struct point{ int x,y,bl,id; double ang; void get(){ if(x==0) ang=(y==0)?-inf:inf; else ang=y/(double)(x); } point(int a=0,int b=0){x=a;y=b;} }p[3010],h[3010],t[3010]; int N,n,m,stk[3010]; bool cmp1(point a,point b){return a.ang<b.ang;} bool cmp2(point a,point b){return a.id<b.id;} int nex(int x){ if(x<N)return x+1; return 1; } point operator-(point a,point b){ return point(a.x-b.x,a.y-b.y); } ll cross(point a,point b){ return a.x*(ll)b.y-a.y*(ll)b.x; } int conv(){ sort(p+1,p+n+m+1,cmp1); int i,top; top=2; stk[1]=1; stk[2]=2; for(i=3;i<=n+m;i++){ while(top>1&&cross(p[stk[top]]-p[stk[top-1]],p[i]-p[stk[top]])<0)top--; top++; stk[top]=i; } for(i=1;i<=top;i++)h[i]=p[stk[i]]; sort(p+1,p+n+m+1,cmp2); return top; } struct edge{ int x,y; edge(int a=0,int b=0){x=a;y=b;} }e[3010]; int tot; void add(int x,int y){ tot++; e[tot]=edge(x,y); } bool intri(point a,point b,point c,point d){ if(cross(b-a,c-b)<0)swap(b,c); return cross(d-b,a-d)>0&&cross(d-c,b-d)>0&&cross(d-a,c-d)>0; } void gao(int a,int b,int c){ int i,e,mx,my; for(i=1;i<=n+m;i++){ if(p[i].bl==p[a].bl&&intri(p[a],p[b],p[c],p[i])){ add(a,i); gao(b,a,i); gao(c,a,i); gao(i,b,c); return; } } e=1; t[1]=p[b]; for(i=1;i<=n+m;i++){ if(intri(p[a],p[b],p[c],p[i])){ e++; t[e]=p[i]; } } mx=inf; for(i=1;i<=e;i++){ if(t[i].x<mx||(t[i].x==mx&&t[i].y<my)){ mx=t[i].x; my=t[i].y; } } for(i=1;i<=e;i++){ t[i].x-=mx; t[i].y-=my; t[i].get(); } sort(t+1,t+e+1,cmp1); for(i=2;i<=e;i++)add(t[1].id,t[i].id); } int main(){ int i,mx,my,ct; scanf("%d%d",&n,&m); mx=inf; for(i=1;i<=n+m;i++){ scanf("%d%d",&p[i].x,&p[i].y); if(p[i].x<mx||(p[i].x==mx&&p[i].y<my)){ mx=p[i].x; my=p[i].y; } p[i].id=i; if(i>n)p[i].bl=1; } for(i=1;i<=n+m;i++){ p[i].x-=mx; p[i].y-=my; p[i].get(); } N=conv(); ct=0; for(i=1;i<=N;i++){ if(h[i].bl!=h[nex(i)].bl)ct++; } if(ct>2){ puts("touhou-kancolle war can not be avoided!"); return 0; } if(ct==0){ for(i=1;i<N;i++)add(h[i].id,h[i+1].id); for(i=1;i<=n+m;i++){ if(p[i].bl!=h[1].bl){ mx=i; break; } } for(i=1;i<=N;i++)gao(mx,h[i].id,h[nex(i)].id); }else{ for(i=1;h[i].bl==h[nex(i)].bl;i=nex(i)); for(i=nex(i);h[i].bl==h[nex(i)].bl;i=nex(i))add(h[i].id,h[nex(i)].id); for(i=nex(i);h[i].bl==h[nex(i)].bl;i=nex(i))add(h[i].id,h[nex(i)].id); mx=i; for(i=nex(i);h[i].bl==h[nex(i)].bl;i=nex(i))gao(h[mx].id,h[i].id,h[nex(i)].id); mx=i; for(i=nex(i);h[i].bl==h[nex(i)].bl;i=nex(i))gao(h[mx].id,h[i].id,h[nex(i)].id); } for(i=1;i<=tot;i++){ if(p[e[i].x].bl==0)printf("%d %d\n",e[i].x,e[i].y); } for(i=1;i<=tot;i++){ if(p[e[i].x].bl==1)printf("%d %d\n",e[i].x-n,e[i].y-n); } }
[xsy1100]东舰停战不可避的更多相关文章
- 珍爱生命,远离JS=>JS避坑记
JavaScript避坑记 转载请注明源地址: http://www.cnblogs.com/funnyzpc/p/8407952.html 上图=> 有意思的漫画,不知大家看懂了没,这里我想说 ...
- 黑马毕向东Java基础知识总结
Java基础知识总结(超级经典) 转自:百度文库 黑马毕向东JAVA基础总结笔记 侵删! 写代码: 1,明确需求.我要做什么? 2,分析思路.我要怎么做?1,2,3. 3,确定步骤.每一个思路部 ...
- 中大东校小米路由器mini实现inode上网,ipv6 wifi【中大】【东校】【inode】【ipv6】
还有不到4个月就要毕业了,前几天半夜没事捣鼓小米路由没想到竟然实现了wifi的ipv6. 正好又安利了同学一台小米路由mini,从刷机到inode到ipv6全搞了一遍. 这里将教程写出来,服务学弟妹. ...
- 东哥读书小记 之 《MacTalk人生元编程》
一直以来的自我感觉:自己是个记性偏弱的人.反正从小读书就喜欢做笔记(可自己的字写得巨丑无比,尼玛不科学呀),抄书这事儿真的就常发生俺的身上. 因为那时经常要背诵课文之类,反正为了怕自己忘记, ...
- json是个啥东东
xml 不用说 只要是搞web开发的 没听说谁不知道的 一种类似数据传输格式定义的语言 但是他却不是一个真正的轻量级的东西 其他的不说 只要传输一点很少的数据 经过他那左括号右括号 还有什么属性 一包 ...
- 长见识了,知道了collected和Graphite 这两个东东
今天下午的讨论会议中,听到了两个名词collected和Graphite这是神马东东,以前在bingo的时候也没听说过,开完会下去查了下.原来他两是监控系统的啊.以前也从来没做过系统监控方面的项目,这 ...
- Win7普通版-X86-SP1-ios-旗舰版安装
------------------------------- 自己的电脑是X64Win7系统.Acer电脑----->>现在改装成 Win7普通版-X86-SP1-ios-旗舰版安装 ...
- BPEL是个什么东东
研究团队有个做智能服务组合的,其中用到叫BPEL的东西,因为全称是Business Process Execution Language,译成中文就是商业执行过程语言,这个东东的是整合SOA的一个执行 ...
- 黑马程序员----java基础笔记中(毕向东)
<p>------<a href="http://www.itheima.com" target="blank">Java培训.Andr ...
随机推荐
- 如何使用Photoshop制作真实的尺子
前言: 日常生活中经常性的偶尔需要测量一些东西的尺寸,但刚好手头上缺乏尺子等必要的测量工具,这时候其实我们可以利用Photoshop,临时制作一把基于现实物理单位(如:厘米)的虚拟尺子. 难点: 像素 ...
- [链接] Linux下常见的~/.bashrc、/etc/profile、/etc/ld.so.config小科普以及caffe编译遇到的相关问题解决
由于博主设置禁止转载,这里贴一个链接,http://blog.csdn.net/u014266895/article/details/61928602,内容很有用,linux下很多软件问题都是各种路径 ...
- 2016广东工业大学校赛 E题 GDUT-oj1173
Problem E: 积木积水 Description 现有一堆边长为1的已经放置好的积木,小明(对的,你没看错,的确是陪伴我们成长的那个小明)想知道当下雨天来时会有多少积水.小明又是如此地喜欢二次元 ...
- lwIP内存管理机制
lwip的内存管理机制,我们以enet_lwip这个例程为例. 在使用lwip的时候,我们可以使用两种形式的内存,一种是heap(mem.c文件-mem_malloc()),一种是pool(memp. ...
- mysql 存在update不存在insert
如果在INSERT语句末尾指定了ON DUPLICATE KEY UPDATE,并且插入行后会导致在一个UNIQUE索引或PRIMARY KEY中出现重复值,则在出现重复值的行执行UPDATE:如果不 ...
- tomcat发布web项目的三种方式
tomcat发布web项目的三种方式 方式一: 配置tomcat 安装目录下的conf/server.xml <Host name="loaclhost">标签里面添加 ...
- jquery 的相关 width 和 height 方法辨析
width() 设置或返回元素的宽度(不包括内边距.边框或外边距). height() 设置或返回元素的高度(不包括内边距.边框或外边距). innerWidth() 返回元素的宽度(包括内边距). ...
- C++ Review
#include "iostream" #include "iomanip" #include "cstdio" using namespa ...
- TCP(二)
TCP半连接和全连接问题 TCP握手过程详解 如上图所示,关键部分:syns queue(半连接队列)和accept queue(全连接队列) 正常情况下的处理过程如下: 1)当server端收到cl ...
- 【洛谷 P2485】 [SDOI2011]计算器 (BSGS)
题目链接 第一问:快速幂 第二问:扩欧解线性同余方程 第三问:\(BSGS\) 三个模板 #include <cstdio> #include <cmath> #include ...