ACM训练指南
ACM入门训练指南
ACM入门训练指南目标读者：想要在ACM/ICPC里进行发展，并通过SDUTOJ进行训练的初学者。
使用语言：只要会一门程序设计语言，就可以进行ACM训练了。
通过训练，可以更好地掌握语言使用能力、程序和算法设计能力。
一般通用语言如C、C++和JAVA都可以，它们有各自的优势和缺点：1．C语言设计算法效率比较高，但输入输出的格式控制比较麻烦，而ACM 对程序进行评测时对输入输出的格式要求比较高，使用C务必要熟练掌握输入输出方法。
2．C++封装了输入输出流，方便输入输出操作，减少出错的可能性；C++提供了非常强大的标准模版库（STL），使得很多在C上实现起来比较麻烦的代码，在C++上却非常方便。
3．JAVA在大型工程和安全方面有比较独特的优势，但在ACM里面却不是一种优秀的语言，因为JAVA的执行效率要比C、C++慢很多，而ACM的题目都对程序运行时间有限制，如果题目限时比较紧的话，就不适合用JAVA，然而JAVA却提供了很方便的高精度运算（大整数运算）。
建议刚学完C就用纯C来训练，在训练过程中可以学习C++，有时间再把STL 好好学学。
输入输出：初次接触ACM训练时经常会遇到的问题，就是输入和输出问题。
如果对语言的输入输出问题不是很熟悉的话，一定要先重点研究一下，特别在输入和输出时不能有冗余信息，因为学习语言时可能习惯了使用提示信息来提高程序的交互性，但ACM不需要任何交互性。
不严格按照题目要求进行输入输出的程序是无法通过系统测试的。
在线评测系统：在线评测系统，英文叫Online Judge（简称OJ），是开放的程序自动评判系统。
只要能上网，注册并登录系统后，就可以选择题目，编写程序，提交程序代码，然后由系统自动进行编译和执行，并通过系统预设测试数据来检验程序代码的正确性。
通过使用OJ训练，可以提高编程和算法设计能力，随着训练的深入，可以参加在评测系统上举行的ACM-ICPC程序设计竞赛。
很多学校都有自己的在线评测系统，里面提供了很多题目给平时学习训练用。
ACM训练计划初步设想
2,对于活动的一些建议。
a,对于英语题目，强烈要求不要直接用翻译器，自己慢慢来。平时也可收集一些好的单词。
b,在刷题过程中，以训练自己的思维为主，少看代码，以思路为主，代码自己实现为好。
c,在刷题过程中，有些模板一定要收集要，而且要理解含义，灵活运用，而且要打印出来。
1，先刷题，后总结。
如，后写代码。
c,是好题，将解题报告写于自己blog上，注重质量，把大的问题划分在小问题。
d,刷题总结过程中，以表达自己的思路为主。越有清晰，越有条理性越好。
e,刷题分类先不固定。可能一次刷题会涉及到多种类型。
f,刷题时间为上班时间。
acm培训计划
acm培训计划导言ACM (Association for Computing Machinery) 是一个国际性的计算机学会，旨在为计算机专业人士提供交流学习和培训的平台。
ACM 培训计划旨在帮助学生提升他们的算法和编程能力，从而更好地参与 ACM 竞赛。
本培训计划将围绕算法与数据结构、编程语言、数学及竞赛技巧展开，以帮助学生提升专业知识、提高团队合作能力和竞赛技能。
一、培训目标1. 提升学生算法和数据结构基础知识，使其能够灵活运用于解决实际问题。
2. 培养学生对编程语言的深刻理解和应用能力。
3. 加强学生的数学基础，提高解决问题的抽象能力。
4. 提高学生的 ACM 竞赛技巧，培养解决问题的思考和团队合作能力。
二、培训内容1. 算法与数据结构1.1. 基本算法：递归、分治、贪心、动态规划1.2. 基本数据结构：栈、队列、链表、树、图1.3. 高级算法：最短路径、最小生成树、网络流、字符串算法1.4. 算法分析与设计：时间复杂度、空间复杂度和算法优化2. 编程语言2.1. C/C++/Java/Python 等主流编程语言的基本语法和特性2.2. 编程范例分析和练习2.3. 算法实现与调试技巧3. 数学基础3.1. 离散数学基础知识3.2. 数论、组合数学和图论基础3.3. 动态规划数学建模4. ACM 竞赛技巧4.1. ACM 竞赛规则和常见题型分析4.2. 模拟训练和解题技巧分享4.3. 队伍协作与策略分享三、培训方式1. 理论授课1.1. 定期组织专家授课，系统讲解培训内容，由资深ACM 竞赛选手分享解题技巧和经验。
1.2. 组织学习交流会，鼓励学生积极提问和讨论。
2. 实践训练2.1. 组织编程实践训练，引导学生独立完成算法实现和调试。
2.2. 选派导师进行一对一指导，帮助学生解决练习中遇到的难点问题。
3. 竞赛准备3.1. 组织模拟 ACM 竞赛，帮助学生提前适应竞赛环境和节奏。
3.2. 参与区域赛和国际赛前的模拟训练，为学生提供更加真实的竞赛体验。
ACM培训计划
ACM培训计划书制作人：xxx2008/10/21一、什么是ACMACM/ICPC ( ACM International Collegiate Programming Contest)国际大学生程序设计竞，ACM/ICPC 是由国际计算机界历史悠久、颇具权威性的组织ACM (Association for Computing Machinery ，美国计算机协会) 主办的，世界上公认的规模最大、水平最高的国际大学生程序设计竞赛，其目的旨在使大学生运用计算机来充分展示自己分析问题和解决问题的能力。
该项竞赛从1970 年举办至今已历31 届，一直受到国际各知名大学的重视，并受到全世界各著名计算机公司的高度关注，在过去十几年中，APPLE 、AT&T 、MICROSOFT 和IBM 等世界著名信息企业分别担任了竞赛的赞助商。
可以说，ACM 国际大学生程序设计竞赛已成为世界各国大学生最具影响力的国际级计算机类的赛事，是广大爱好计算机编程的大学生展示才华的舞台，是著名大学计算机教育成果的直接体现，是信息企业与世界顶尖计算机人才对话的最好机会。
该项竞赛分区域预赛和国际决赛两个阶段进行，各预赛区第一名自动获得参加世界决赛的资格，世界决赛安排在每年的 3 ～ 4 月举行，而区域预赛安排在上一年的9 ～12 月在各大洲举行。
ACM/ICPC 的区域预赛是规模很大、范围很广的赛事。
仅在2003 年参加区域预赛的队伍就有来自75 个国家（地区），1411 所大学的3150 支代表队，他们分别在127 个赛场中进行比赛，以争夺全球总决赛的73 个名额，其激烈程度可想而知。
2004 年第29 届ACM/ICPC 亚洲赛区预赛共设了北京、上海、台北、高雄、汉城、德黑兰、爱媛（日本）、达卡（孟加拉国）、马尼拉、坎普尔（印度）等10 个赛站，来自亚洲各国知名高校的各个代表队进行了激烈的角逐。
中国内地从1996 年开始参加ACM/ICPC 亚洲区预赛，至今已历九届。
算法竞赛入门经典 训练指南 第三章
第三章：数据结构在算法竞赛中，数据结构是一个非常重要的知识点，它在解决问题中起着至关重要的作用。
数据结构指的是一组数据的组织方式和管理方式，能够高效地进行数据的存储和检索。
在算法竞赛中，熟练掌握各种数据结构，并结合实际问题灵活运用，对于解决复杂的问题至关重要。
1. 数组数组是最基本的数据结构之一。
它在内存中以连续的方式存储数据，可以通过下标随机访问数组中的元素。
在算法竞赛中，数组的使用非常广泛，例如在搜索、排序和动态规划等算法中都会用到数组。
数组也有其局限性，比如插入和删除操作比较耗时。
在使用数组时需要根据具体问题进行评估。
2. 链表链表是另一种常见的数据结构，它由节点组成，每个节点包含数据和指向下一个节点的指针。
与数组不同，链表的节点在内存中不一定是连续的，这样可以更加灵活地进行插入和删除操作。
在算法竞赛中，链表常常用于实现队列、栈等数据结构，同时也可以作为其他算法的辅助数据结构进行使用。
3. 栈和队列栈和队列是两种基本的数据结构，它们分别遵循“后进先出”和“先进先出”的原则。
栈可以通过压栈和弹栈操作实现，而队列则可以通过入队和出队操作来实现。
在算法竞赛中，栈和队列常常用于解决各种实际问题，比如表达式求值、迷宫搜索等。
4. 哈希表哈希表是一种以键值对形式存储数据的数据结构，它通过哈希函数将关键字映射到表中的一个位置，从而实现快速的检索。
在算法竞赛中，哈希表常常用于统计和去重等问题，通过合理设计哈希函数和解决冲突，可以大大提高算法的效率。
5. 树和图树和图是两种复杂的数据结构，它们在算法竞赛中经常用于解决各种复杂的问题。
树是一种非线性数据结构，常见的有二叉树、平衡树等，而图则是一种用边来表示数据之间关系的数据结构，包括有向图、无向图等。
在算法竞赛中，熟练掌握树和图的遍历、最短路径、最小生成树等算法，对于解决各种实际问题至关重要。
总结回顾数据结构是算法竞赛中的重要知识点，熟练掌握各种数据结构，并结合实际问题进行灵活运用，对于解决复杂的问题至关重要。
ACM简介及使用指南
期刊、杂志和会报53种； 近300个会议，4000多卷会议录； 超过40万多篇全文，及“在线计算机文献指南”数据库 中230多万条文摘题录信息； 37种SIG时事通讯 ACM附属机构出版物 ACM口述历史访谈录
使用指南快速检索区高级检索按照出版物不同类别分类浏览computingreviews此处链接至reviewscomacm电子书系列快速检索acmdl和guide文摘库快速检索框输入要搜索的关键字如ascoustics点击searchguide文摘数据库检索结果guide文摘库检索的结果数量文章出处出版社文摘结果显示可按不同要求排序acm及附属机构出版物检索结果将快速检索的结果限定在acm和附属机构的出版物范围内点击红色字体部分acm及附属机构出版物检索结果文摘名称全文pdf格式acm及附属机构出版物检索结果41高级检索在高级检索中用户可通过单独或组合字段检索限定检索文献类型获取所需数据词组或姓名wordsorphrasesnames关键字作者所在机构keywordsaffiliations出版物publication会议conferenceissnisbndoi检索acm计算机分类体系ccscomputingclassificationsystem41高级检索可在文章的标题文摘评论中检索可进行作者编辑指导者评论者检索在作者给出的文章关键词中检索作者所属的机构检索按照出版物类别检索按照会议相关信息检索issnisbndoi检索acm计算机分类体系ccs41高级检索确定出版物或出版社限定文献出版时间限定文献类型限定检索到的文献所包含的内容通过出版社或出版物检索确定会议举办方确定会议召开地点确定会议召开时间通过会议信息检索文章通过出版物的isbnissn号或doi检索检索通过ccs计算科学分类表的类目检索42检索案例检索2005年至今基于ipv6的网络安全机制研究的文章421检索结果页通过关键词筛选检索结果通过作者评论者及他们所属的机构筛选检索结果通过出版时间出版物名称出版社文献类型筛选检索结果通过会议举办方会议相关活动会议录名称筛选检索结果421检索结果页与检索式相关的期刊杂志sig出版物会议选择检索结果的排序方式选择检索结果的显示方式421检索结果页点击文章标题链接进入引文信息页点击作者链接进入该作者的profile页面下载pdf全文该文章被下载的总次数6周内和1年内被下载的次数422文章信息文章相关信息
清华大学ACM集训队培训资料(内部使用)
return 0;
}
int stonetolb(int sts)
{
return 14 * sts;
}
下面是运行情况：
Enter the weight in sone:14
14 stone = 196 pounds.
程序通过cin语句给stone提供一个值，然后在main函数中，把这个值传递给stonetolb()函数，这个植被赋给sts之后，stonetolb()用return将14*sts返回给main()。
type functionname (argumentlist)
{
statements
}
注意，定义simon()的代码在main()函数的后面，C++中不允许将函数定义在另一个函数内。每个函数的定义都是独立的，所有的函数的创建都是平等的。
simon()函数的函数头定义如下：
void simon(int n)
清华大学ACM集训队培训资料（内部使用）
一、C++基础
基本知识
所有的C++程序都是有函数组成的，函数又叫做子程序，且每个C++程序必须包含一个main函数，编译器（能够把源代码转换成目标代码的程序）把翻译后的目标代码和一些启动代码组合起来，生成可执行文件，main函数就是可执行文件的入口，所以，每个C++程序有且只有一个main函数。
如果我们修改程序1.1中的main()函数的名称，将其改为Main()，那么，IDE就会给出错误信息，比如“[Linker error] undefined reference to `WinMain@16'”，因为编译器没有找到main函数。
ACM培训大纲
实用标准文案ACM培训大纲基础内容：数据结构——》搜索——》图论DP数论博弈中级内容数据结构网络流第一章搜索1.二分搜索三分搜索2.栈3.队列4.深搜5,广搜6.第二章数据结构1.优先队列并查集2.二叉搜索树3.线段树（单点更新）4.5.精彩文档.实用标准文案第三章图论1.图的表示1.1二维数组1.2邻接表1.3前向星2.图的遍历2.1双连通分量2. 2拓扑排序3.最短路3.1迪杰斯特拉3. 2弗洛伊德4. 3 SPFA5.匹配匈牙利算法6.生成树7.网络流简介第四章动态规划1.状态转移方程2.引入3. 1 0-1背包4.2硬币问题5. 3矩阵链乘6.区间DP7.按位DP8.树形DP9.状压DP第五章数论1.欧几里得扩展欧几里得2.因数分解3. 费马小定理4.欧拉定理5.6.1筛法6. 2素数判定6. 2,1 0(Jn)方法精彩文档.实用标准文案6. 2. 2 Mi I ler-rabin 测试第六章博弈1.Nim 和2.SG函数第七章中级数据结构1.树状数组RMO 2.KMP3.AC自动机4.线段树（区间更新）5.第八章图论进阶1.网络流问题精彩文档.实用标准文案综述在很多人眼里，东北大学秦皇岛分校不算是985高校。
所以我们要用自己的能力证明我们有985 的实力。
ACM是计算机界认可度最高的一个比赛，可以说只要区域赛有过奖牌，国内任何IT公司没有理由不要。
同时，在高校之中，对一个大学计算机专业的评价，大部分人也会首先看ACM 的水平。
将ACM打出学校，在国内打出一定成绩，对扩大我校影响力很有帮助。
考虑到本校暂时没有进行专题训练的出题能力，专题训练的题目主要从UESTC 2014年集训队专题训练中获取，再加上从别的0J上找一些题目。
训练的平台设置在华中科技大学的vertual judge上面。
本人将在毕业之前承担培训任务。
在2015学年开始之前，培训计划为每两周一次，中间空闲的时间由大二或者大一熟悉C++的同学给不熟悉C++的同学进行基础的讲解。
acm训练资料(浙大模板)
1、几何1.1 注意 (4)1.2 几何公式 (4)1.3 多边形 (6)1.4 多边形切割 (9)1.5 浮点函数 (10)1.6 面积 (15)1.7 球面 (16)1.8 三角形 (17)1.9 三维几何 (19)1.10 凸包 (26)1.11 网格 (28)1.12 圆 (29)1.13 整数函数 (30)2、组合2.1 组合公式 (33)2.2 排列组合生成 (33)2.3 生成gray码 (35)2.4 置换(polya) (35)2.5 字典序全排列 (36)2.6 字典序组合 (36)3、结构3.1 并查集 (37)3.2 堆 (38)3.3 线段树 (39)3.4 子段和 (44)3.5 子阵和 (44)4、数论4.1 阶乘最后非0位 (45)4.2 模线性方程组 (46)4.3 素数 (47)4.4 欧拉函数 (48)5、数值计算5.1 定积分计算(Romberg) (49)5.2 多项式求根(牛顿法) (51)5.3 周期性方程(追赶法) (52)6、图论—NP搜索6.1 最大团 (53)6.2 最大团(n<64)(faster) (54)7、图论—连通性7.1 无向图关键点(dfs邻接阵) (56)7.2 无向图关键边(dfs邻接阵) (57)7.3 无向图的块(bfs邻接阵) (58)7.4 无向图连通分支(dfs/bfs邻接阵) (59)7.5 有向图强连通分支(dfs/bfs邻接阵) (60)7.6 有向图最小点基(邻接阵) (61)8、图论—匹配8.1 二分图最大匹配(hungary邻接表) (62)8.2 二分图最大匹配(hungary邻接阵) (63)8.3 二分图最大匹配(hungary正向表) (63)8.4二分图最佳匹配(kuhn_munkras邻接阵) (64)8.5 一般图匹配(邻接表) (65)8.6 一般图匹配(邻接阵) (66)8.7 一般图匹配(正向表) (66)9、图论—网络流9.1 最大流(邻接阵) (67)9.2 上下界最大流(邻接阵) (68)9.3 上下界最小流(邻接阵) (69)9.4 最大流无流量(邻接阵) (70)9.5 最小费用最大流(邻接阵) (70)10、图论—应用10.1 欧拉回路(邻接阵) (71)10.2 树的前序表转化 (72)10.3 树的优化算法 (73)10.4 拓扑排序(邻接阵) (74)10.5 最佳边割集 (75)10.6 最佳点割集 (76)10.7 最小边割集 (77)10.8 最小点割集 (78)10.9 最小路径覆盖 (80)11、图论—支撑树11.1 最小生成树(kruskal邻接表) (80)11.2 最小生成树(kruskal正向表) (82)11.3 最小生成树(prim+binary_heap邻接表) (83)11.4 最小生成树(prim+binary_heap正向表) (84)11.5 最小生成树(prim+mapped_heap邻接表) (85)11.6 最小生成树(prim+mapped_heap正向表) (87)11.7 最小生成树(prim邻接阵) (88)11.8 最小树形图(邻接阵) (88)12、图论—最短路径12.1 最短路径(单源bellman_ford邻接阵) (90)12.2 最短路径(单源dijkstra+bfs邻接表) (90)12.3 最短路径(单源dijkstra+bfs正向表) (91)12.4 最短路径(单源dijkstra+binary_heap邻接表) (92)12.5 最短路径(单源dijkstra+binary_heap正向表) (93)12.6 最短路径(单源dijkstra+mapped_heap邻接表) (94)12.7 最短路径(单源dijkstra+mapped_heap正向表) (95)12.8 最短路径(单源dijkstra邻接阵) (96)12.9 最短路径(多源floyd_warshall邻接阵) (97)13、应用13.1 Joseph问题 (97)13.2 N皇后构造解 (98)13.3 布尔母函数 (99)13.4 第k元素 (99)13.5 幻方构造 (100)13.6 模式匹配(kmp) (101)13.7 逆序对数 (102)13.8 字符串最小表示 (102)13.9 最长公共单调子序列 (103)13.10 最长子序列 (104)13.11 最大子串匹配 (105)13.12 最大子段和 (106)13.13 最大子阵和 (106)14、其它14.1 大数(只能处理正数) (107)14.2 分数 (113)14.3 矩阵 (115)14.4 线性方程组 (117)14.5 线性相关 (119)14.6 日期 (120)1、几何1.1注意1. 注意舍入方式(0.5的舍入方向);防止输出-0.2. 几何题注意多测试不对称数据.3. 整数几何注意xmult和dmult是否会出界;符点几何注意eps的使用.4. 避免使用斜率;注意除数是否会为0.5. 公式一定要化简后再代入.6. 判断同一个2*PI域内两角度差应该是abs(a1-a2)<beta||abs(a1-a2)>pi+pi-beta;相等应该是abs(a1-a2)<eps||abs(a1-a2)>pi+pi-eps;7. 需要的话尽量使用atan2,注意:atan2(0,0)=0,atan2(1,0)=pi/2,atan2(-1,0)=-pi/2,atan2(0,1)=0,atan2(0,-1)=pi.8. cross product = |u|*|v|*sin(a)dot product = |u|*|v|*cos(a)9. (P1-P0)x(P2-P0)结果的意义:正: <P0,P1>在<P0,P2>顺时针(0,pi)内负: <P0,P1>在<P0,P2>逆时针(0,pi)内0 : <P0,P1>,<P0,P2>共线,夹角为0或pi10. 误差限缺省使用1e-8!1.2几何公式三角形:1. 半周长P=(a+b+c)/22. 面积S=aHa/2=absin(C)/2=sqrt(P(P-a)(P-b)(P-c))3. 中线Ma=sqrt(2(b^2+c^2)-a^2)/2=sqrt(b^2+c^2+2bccos(A))/24. 角平分线Ta=sqrt(bc((b+c)^2-a^2))/(b+c)=2bccos(A/2)/(b+c)5. 高线Ha=bsin(C)=csin(B)=sqrt(b^2-((a^2+b^2-c^2)/(2a))^2)6. 内切圆半径r=S/P=asin(B/2)sin(C/2)/sin((B+C)/2)=4Rsin(A/2)sin(B/2)sin(C/2)=sqrt((P-a)(P-b)(P-c)/P)=Ptan(A/2)tan(B/2)tan(C/2)7. 外接圆半径R=abc/(4S)=a/(2sin(A))=b/(2sin(B))=c/(2sin(C))四边形:D1,D2为对角线,M对角线中点连线,A为对角线夹角1. a^2+b^2+c^2+d^2=D1^2+D2^2+4M^22. S=D1D2sin(A)/2(以下对圆的内接四边形)3. ac+bd=D1D24. S=sqrt((P-a)(P-b)(P-c)(P-d)),P为半周长正n边形:R为外接圆半径,r为内切圆半径1. 中心角A=2PI/n2. 内角C=(n-2)PI/n3. 边长a=2sqrt(R^2-r^2)=2Rsin(A/2)=2rtan(A/2)4. 面积S=nar/2=nr^2tan(A/2)=nR^2sin(A)/2=na^2/(4tan(A/2))圆:1. 弧长l=rA2. 弦长a=2sqrt(2hr-h^2)=2rsin(A/2)3. 弓形高h=r-sqrt(r^2-a^2/4)=r(1-cos(A/2))=atan(A/4)/24. 扇形面积S1=rl/2=r^2A/25. 弓形面积S2=(rl-a(r-h))/2=r^2(A-sin(A))/2棱柱:1. 体积V=Ah,A为底面积,h为高2. 侧面积S=lp,l为棱长,p为直截面周长3. 全面积T=S+2A棱锥:1. 体积V=Ah/3,A为底面积,h为高(以下对正棱锥)2. 侧面积S=lp/2,l为斜高,p为底面周长3. 全面积T=S+A棱台:1. 体积V=(A1+A2+sqrt(A1A2))h/3,A1.A2为上下底面积,h为高(以下为正棱台)2. 侧面积S=(p1+p2)l/2,p1.p2为上下底面周长,l为斜高3. 全面积T=S+A1+A2圆柱:1. 侧面积S=2PIrh2. 全面积T=2PIr(h+r)3. 体积V=PIr^2h圆锥:1. 母线l=sqrt(h^2+r^2)2. 侧面积S=PIrl3. 全面积T=PIr(l+r)4. 体积V=PIr^2h/3圆台:1. 母线l=sqrt(h^2+(r1-r2)^2)2. 侧面积S=PI(r1+r2)l3. 全面积T=PIr1(l+r1)+PIr2(l+r2)4. 体积V=PI(r1^2+r2^2+r1r2)h/3球:1. 全面积T=4PIr^22. 体积V=4PIr^3/3球台:1. 侧面积S=2PIrh2. 全面积T=PI(2rh+r1^2+r2^2)3. 体积V=PIh(3(r1^2+r2^2)+h^2)/6球扇形:1. 全面积T=PIr(2h+r0),h为球冠高,r0为球冠底面半径2. 体积V=2PIr^2h/31.3多边形#include <stdlib.h>#include <math.h>#define MAXN 1000#define offset 10000#define eps 1e-8#define zero(x) (((x)>0?(x):-(x))<eps)#define _sign(x) ((x)>eps?1:((x)<-eps?2:0))struct point{double x,y;};struct line{point a,b;};double xmult(point p1,point p2,point p0){return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);}//判定凸多边形,顶点按顺时针或逆时针给出,允许相邻边共线int is_convex(int n,point* p){int i,s[3]={1,1,1};for (i=0;i<n&&s[1]|s[2];i++)s[_sign(xmult(p[(i+1)%n],p[(i+2)%n],p[i]))]=0;return s[1]|s[2];}//判定凸多边形,顶点按顺时针或逆时针给出,不允许相邻边共线int is_convex_v2(int n,point* p){int i,s[3]={1,1,1};for (i=0;i<n&&s[0]&&s[1]|s[2];i++)s[_sign(xmult(p[(i+1)%n],p[(i+2)%n],p[i]))]=0;return s[0]&&s[1]|s[2];}//判点在凸多边形内或多边形边上,顶点按顺时针或逆时针给出int inside_convex(point q,int n,point* p){int i,s[3]={1,1,1};for (i=0;i<n&&s[1]|s[2];i++)s[_sign(xmult(p[(i+1)%n],q,p[i]))]=0;return s[1]|s[2];}//判点在凸多边形内,顶点按顺时针或逆时针给出,在多边形边上返回0int inside_convex_v2(point q,int n,point* p){int i,s[3]={1,1,1};for (i=0;i<n&&s[0]&&s[1]|s[2];i++)s[_sign(xmult(p[(i+1)%n],q,p[i]))]=0;return s[0]&&s[1]|s[2];}//判点在任意多边形内,顶点按顺时针或逆时针给出//on_edge表示点在多边形边上时的返回值,offset为多边形坐标上限int inside_polygon(point q,int n,point* p,int on_edge=1){point q2;int i=0,count;while (i<n)for (count=i=0,q2.x=rand()+offset,q2.y=rand()+offset;i<n;i++)if(zero(xmult(q,p[i],p[(i+1)%n]))&&(p[i].x-q.x)*(p[(i+1)%n].x-q.x)<eps&&(p[i].y-q.y)*(p[(i+1)% n].y-q.y)<eps)return on_edge;else if (zero(xmult(q,q2,p[i])))break;else if (xmult(q,p[i],q2)*xmult(q,p[(i+1)%n],q2)<-eps&&xmult(p[i],q,p[(i+1)%n])*xmult(p[i],q2,p[(i+1) %n])<-eps)count++;return count&1;}inline int opposite_side(point p1,point p2,point l1,point l2){return xmult(l1,p1,l2)*xmult(l1,p2,l2)<-eps;}inline int dot_online_in(point p,point l1,point l2){return zero(xmult(p,l1,l2))&&(l1.x-p.x)*(l2.x-p.x)<eps&&(l1.y-p.y)*(l2.y-p.y)<eps;}//判线段在任意多边形内,顶点按顺时针或逆时针给出,与边界相交返回1int inside_polygon(point l1,point l2,int n,point* p){point t[MAXN],tt;int i,j,k=0;if (!inside_polygon(l1,n,p)||!inside_polygon(l2,n,p))return 0;for (i=0;i<n;i++)if (opposite_side(l1,l2,p[i],p[(i+1)%n])&&opposite_side(p[i],p[(i+1)%n],l1,l2)) return 0;else if (dot_online_in(l1,p[i],p[(i+1)%n]))t[k++]=l1;else if (dot_online_in(l2,p[i],p[(i+1)%n]))t[k++]=l2;else if (dot_online_in(p[i],l1,l2))t[k++]=p[i];for (i=0;i<k;i++)for (j=i+1;j<k;j++){tt.x=(t[i].x+t[j].x)/2;tt.y=(t[i].y+t[j].y)/2;if (!inside_polygon(tt,n,p))return 0;}return 1;}point intersection(line u,line v){point ret=u.a;double t=((u.a.x-v.a.x)*(v.a.y-v.b.y)-(u.a.y-v.a.y)*(v.a.x-v.b.x)) /((u.a.x-u.b.x)*(v.a.y-v.b.y)-(u.a.y-u.b.y)*(v.a.x-v.b.x));ret.x+=(u.b.x-u.a.x)*t;ret.y+=(u.b.y-u.a.y)*t;return ret;}point barycenter(point a,point b,point c){line u,v;u.a.x=(a.x+b.x)/2;u.a.y=(a.y+b.y)/2;u.b=c;v.a.x=(a.x+c.x)/2;v.a.y=(a.y+c.y)/2;v.b=b;return intersection(u,v);}//多边形重心point barycenter(int n,point* p){point ret,t;double t1=0,t2;int i;ret.x=ret.y=0;for (i=1;i<n-1;i++)if (fabs(t2=xmult(p[0],p[i],p[i+1]))>eps){t=barycenter(p[0],p[i],p[i+1]);ret.x+=t.x*t2;ret.y+=t.y*t2;t1+=t2;}if (fabs(t1)>eps)ret.x/=t1,ret.y/=t1;return ret;}1.4多边形切割//多边形切割//可用于半平面交#define MAXN 100#define eps 1e-8#define zero(x) (((x)>0?(x):-(x))<eps)struct point{double x,y;};double xmult(point p1,point p2,point p0){return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);}int same_side(point p1,point p2,point l1,point l2){return xmult(l1,p1,l2)*xmult(l1,p2,l2)>eps;}point intersection(point u1,point u2,point v1,point v2){point ret=u1;double t=((u1.x-v1.x)*(v1.y-v2.y)-(u1.y-v1.y)*(v1.x-v2.x))/((u1.x-u2.x)*(v1.y-v2.y)-(u1.y-u2.y)*(v1.x-v2.x));ret.x+=(u2.x-u1.x)*t;ret.y+=(u2.y-u1.y)*t;return ret;}//将多边形沿l1,l2确定的直线切割在side侧切割,保证l1,l2,side不共线void polygon_cut(int& n,point* p,point l1,point l2,point side){point pp[100];int m=0,i;for (i=0;i<n;i++){if (same_side(p[i],side,l1,l2))pp[m++]=p[i];if(!same_side(p[i],p[(i+1)%n],l1,l2)&&!(zero(xmult(p[i],l1,l2))&&zero(xmult(p[(i+1)%n],l1,l2)))) pp[m++]=intersection(p[i],p[(i+1)%n],l1,l2);}for (n=i=0;i<m;i++)if (!i||!zero(pp[i].x-pp[i-1].x)||!zero(pp[i].y-pp[i-1].y))p[n++]=pp[i];if (zero(p[n-1].x-p[0].x)&&zero(p[n-1].y-p[0].y))n--;if (n<3)n=0;}1.5浮点函数//浮点几何函数库#include <math.h>#define eps 1e-8#define zero(x) (((x)>0?(x):-(x))<eps)struct point{double x,y;};struct line{point a,b;};//计算cross product (P1-P0)x(P2-P0)double xmult(point p1,point p2,point p0){return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);}double xmult(double x1,double y1,double x2,double y2,double x0,double y0){ return (x1-x0)*(y2-y0)-(x2-x0)*(y1-y0);}//计算dot product (P1-P0).(P2-P0)double dmult(point p1,point p2,point p0){return (p1.x-p0.x)*(p2.x-p0.x)+(p1.y-p0.y)*(p2.y-p0.y);}double dmult(double x1,double y1,double x2,double y2,double x0,double y0){ return (x1-x0)*(x2-x0)+(y1-y0)*(y2-y0);}//两点距离double distance(point p1,point p2){return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));}double distance(double x1,double y1,double x2,double y2){return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));}//判三点共线int dots_inline(point p1,point p2,point p3){return zero(xmult(p1,p2,p3));}int dots_inline(double x1,double y1,double x2,double y2,double x3,double y3){ return zero(xmult(x1,y1,x2,y2,x3,y3));}//判点是否在线段上,包括端点int dot_online_in(point p,line l){return zero(xmult(p,l.a,l.b))&&(l.a.x-p.x)*(l.b.x-p.x)<eps&&(l.a.y-p.y)*(l.b.y-p.y)<eps; }int dot_online_in(point p,point l1,point l2){return zero(xmult(p,l1,l2))&&(l1.x-p.x)*(l2.x-p.x)<eps&&(l1.y-p.y)*(l2.y-p.y)<eps;}int dot_online_in(double x,double y,double x1,double y1,double x2,double y2){return zero(xmult(x,y,x1,y1,x2,y2))&&(x1-x)*(x2-x)<eps&&(y1-y)*(y2-y)<eps;}//判点是否在线段上,不包括端点int dot_online_ex(point p,line l){returndot_online_in(p,l)&&(!zero(p.x-l.a.x)||!zero(p.y-l.a.y))&&(!zero(p.x-l.b.x)||!zero(p.y-l.b.y)); }int dot_online_ex(point p,point l1,point l2){returndot_online_in(p,l1,l2)&&(!zero(p.x-l1.x)||!zero(p.y-l1.y))&&(!zero(p.x-l2.x)||!zero(p.y-l2.y)); }int dot_online_ex(double x,double y,double x1,double y1,double x2,double y2){ returndot_online_in(x,y,x1,y1,x2,y2)&&(!zero(x-x1)||!zero(y-y1))&&(!zero(x-x2)||!zero(y-y2));}//判两点在线段同侧,点在线段上返回0int same_side(point p1,point p2,line l){return xmult(l.a,p1,l.b)*xmult(l.a,p2,l.b)>eps;}int same_side(point p1,point p2,point l1,point l2){return xmult(l1,p1,l2)*xmult(l1,p2,l2)>eps;}//判两点在线段异侧,点在线段上返回0int opposite_side(point p1,point p2,line l){return xmult(l.a,p1,l.b)*xmult(l.a,p2,l.b)<-eps;}int opposite_side(point p1,point p2,point l1,point l2){return xmult(l1,p1,l2)*xmult(l1,p2,l2)<-eps;}//判两直线平行int parallel(line u,line v){return zero((u.a.x-u.b.x)*(v.a.y-v.b.y)-(v.a.x-v.b.x)*(u.a.y-u.b.y));}int parallel(point u1,point u2,point v1,point v2){return zero((u1.x-u2.x)*(v1.y-v2.y)-(v1.x-v2.x)*(u1.y-u2.y));}//判两直线垂直int perpendicular(line u,line v){return zero((u.a.x-u.b.x)*(v.a.x-v.b.x)+(u.a.y-u.b.y)*(v.a.y-v.b.y));int perpendicular(point u1,point u2,point v1,point v2){return zero((u1.x-u2.x)*(v1.x-v2.x)+(u1.y-u2.y)*(v1.y-v2.y));}//判两线段相交,包括端点和部分重合int intersect_in(line u,line v){if (!dots_inline(u.a,u.b,v.a)||!dots_inline(u.a,u.b,v.b))return !same_side(u.a,u.b,v)&&!same_side(v.a,v.b,u);return dot_online_in(u.a,v)||dot_online_in(u.b,v)||dot_online_in(v.a,u)||dot_online_in(v.b,u); }int intersect_in(point u1,point u2,point v1,point v2){if (!dots_inline(u1,u2,v1)||!dots_inline(u1,u2,v2))return !same_side(u1,u2,v1,v2)&&!same_side(v1,v2,u1,u2);returndot_online_in(u1,v1,v2)||dot_online_in(u2,v1,v2)||dot_online_in(v1,u1,u2)||dot_online_in(v2,u1,u 2);}//判两线段相交,不包括端点和部分重合int intersect_ex(line u,line v){return opposite_side(u.a,u.b,v)&&opposite_side(v.a,v.b,u);}int intersect_ex(point u1,point u2,point v1,point v2){return opposite_side(u1,u2,v1,v2)&&opposite_side(v1,v2,u1,u2);}//计算两直线交点,注意事先判断直线是否平行!//线段交点请另外判线段相交(同时还是要判断是否平行!)point intersection(line u,line v){point ret=u.a;double t=((u.a.x-v.a.x)*(v.a.y-v.b.y)-(u.a.y-v.a.y)*(v.a.x-v.b.x))/((u.a.x-u.b.x)*(v.a.y-v.b.y)-(u.a.y-u.b.y)*(v.a.x-v.b.x));ret.x+=(u.b.x-u.a.x)*t;ret.y+=(u.b.y-u.a.y)*t;return ret;}point intersection(point u1,point u2,point v1,point v2){point ret=u1;double t=((u1.x-v1.x)*(v1.y-v2.y)-(u1.y-v1.y)*(v1.x-v2.x))/((u1.x-u2.x)*(v1.y-v2.y)-(u1.y-u2.y)*(v1.x-v2.x));ret.x+=(u2.x-u1.x)*t;ret.y+=(u2.y-u1.y)*t;return ret;//点到直线上的最近点point ptoline(point p,line l){point t=p;t.x+=l.a.y-l.b.y,t.y+=l.b.x-l.a.x;return intersection(p,t,l.a,l.b);}point ptoline(point p,point l1,point l2){point t=p;t.x+=l1.y-l2.y,t.y+=l2.x-l1.x;return intersection(p,t,l1,l2);}//点到直线距离double disptoline(point p,line l){return fabs(xmult(p,l.a,l.b))/distance(l.a,l.b);}double disptoline(point p,point l1,point l2){return fabs(xmult(p,l1,l2))/distance(l1,l2);}double disptoline(double x,double y,double x1,double y1,double x2,double y2){ return fabs(xmult(x,y,x1,y1,x2,y2))/distance(x1,y1,x2,y2);}//点到线段上的最近点point ptoseg(point p,line l){point t=p;t.x+=l.a.y-l.b.y,t.y+=l.b.x-l.a.x;if (xmult(l.a,t,p)*xmult(l.b,t,p)>eps)return distance(p,l.a)<distance(p,l.b)?l.a:l.b;return intersection(p,t,l.a,l.b);}point ptoseg(point p,point l1,point l2){point t=p;t.x+=l1.y-l2.y,t.y+=l2.x-l1.x;if (xmult(l1,t,p)*xmult(l2,t,p)>eps)return distance(p,l1)<distance(p,l2)?l1:l2;return intersection(p,t,l1,l2);}//点到线段距离double disptoseg(point p,line l){point t=p;t.x+=l.a.y-l.b.y,t.y+=l.b.x-l.a.x;if (xmult(l.a,t,p)*xmult(l.b,t,p)>eps)return distance(p,l.a)<distance(p,l.b)?distance(p,l.a):distance(p,l.b);return fabs(xmult(p,l.a,l.b))/distance(l.a,l.b);}double disptoseg(point p,point l1,point l2){point t=p;t.x+=l1.y-l2.y,t.y+=l2.x-l1.x;if (xmult(l1,t,p)*xmult(l2,t,p)>eps)return distance(p,l1)<distance(p,l2)?distance(p,l1):distance(p,l2);return fabs(xmult(p,l1,l2))/distance(l1,l2);}//矢量V以P为顶点逆时针旋转angle并放大scale倍point rotate(point v,point p,double angle,double scale){point ret=p;v.x-=p.x,v.y-=p.y;p.x=scale*cos(angle);p.y=scale*sin(angle);ret.x+=v.x*p.x-v.y*p.y;ret.y+=v.x*p.y+v.y*p.x;return ret;}1.6面积#include <math.h>struct point{double x,y;};//计算cross product (P1-P0)x(P2-P0)double xmult(point p1,point p2,point p0){return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);}double xmult(double x1,double y1,double x2,double y2,double x0,double y0){ return (x1-x0)*(y2-y0)-(x2-x0)*(y1-y0);}//计算三角形面积,输入三顶点double area_triangle(point p1,point p2,point p3){return fabs(xmult(p1,p2,p3))/2;}double area_triangle(double x1,double y1,double x2,double y2,double x3,double y3){ return fabs(xmult(x1,y1,x2,y2,x3,y3))/2;}//计算三角形面积,输入三边长double area_triangle(double a,double b,double c){double s=(a+b+c)/2;return sqrt(s*(s-a)*(s-b)*(s-c));}//计算多边形面积,顶点按顺时针或逆时针给出double area_polygon(int n,point* p){double s1=0,s2=0;int i;for (i=0;i<n;i++)s1+=p[(i+1)%n].y*p[i].x,s2+=p[(i+1)%n].y*p[(i+2)%n].x;return fabs(s1-s2)/2;}1.7球面#include <math.h>const double pi=acos(-1);//计算圆心角lat表示纬度,-90<=w<=90,lng表示经度//返回两点所在大圆劣弧对应圆心角,0<=angle<=pidouble angle(double lng1,double lat1,double lng2,double lat2){ double dlng=fabs(lng1-lng2)*pi/180;while (dlng>=pi+pi)dlng-=pi+pi;if (dlng>pi)dlng=pi+pi-dlng;lat1*=pi/180,lat2*=pi/180;return acos(cos(lat1)*cos(lat2)*cos(dlng)+sin(lat1)*sin(lat2));}//计算距离,r为球半径double line_dist(double r,double lng1,double lat1,double lng2,double lat2){ double dlng=fabs(lng1-lng2)*pi/180;while (dlng>=pi+pi)dlng-=pi+pi;if (dlng>pi)dlng=pi+pi-dlng;lat1*=pi/180,lat2*=pi/180;return r*sqrt(2-2*(cos(lat1)*cos(lat2)*cos(dlng)+sin(lat1)*sin(lat2))); }//计算球面距离,r为球半径inline double sphere_dist(double r,double lng1,double lat1,double lng2,double lat2){ return r*angle(lng1,lat1,lng2,lat2);}1.8三角形#include <math.h>struct point{double x,y;};struct line{point a,b;};double distance(point p1,point p2){return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));}point intersection(line u,line v){point ret=u.a;double t=((u.a.x-v.a.x)*(v.a.y-v.b.y)-(u.a.y-v.a.y)*(v.a.x-v.b.x))/((u.a.x-u.b.x)*(v.a.y-v.b.y)-(u.a.y-u.b.y)*(v.a.x-v.b.x));ret.x+=(u.b.x-u.a.x)*t;ret.y+=(u.b.y-u.a.y)*t;return ret;}//外心point circumcenter(point a,point b,point c){line u,v;u.a.x=(a.x+b.x)/2;u.a.y=(a.y+b.y)/2;u.b.x=u.a.x-a.y+b.y;u.b.y=u.a.y+a.x-b.x;v.a.x=(a.x+c.x)/2;v.a.y=(a.y+c.y)/2;v.b.x=v.a.x-a.y+c.y;v.b.y=v.a.y+a.x-c.x;return intersection(u,v);}//内心point incenter(point a,point b,point c){line u,v;double m,n;u.a=a;m=atan2(b.y-a.y,b.x-a.x);n=atan2(c.y-a.y,c.x-a.x);u.b.x=u.a.x+cos((m+n)/2);u.b.y=u.a.y+sin((m+n)/2);v.a=b;m=atan2(a.y-b.y,a.x-b.x);n=atan2(c.y-b.y,c.x-b.x);v.b.x=v.a.x+cos((m+n)/2);v.b.y=v.a.y+sin((m+n)/2);return intersection(u,v);}//垂心point perpencenter(point a,point b,point c){line u,v;u.a=c;u.b.x=u.a.x-a.y+b.y;u.b.y=u.a.y+a.x-b.x;v.a=b;v.b.x=v.a.x-a.y+c.y;v.b.y=v.a.y+a.x-c.x;return intersection(u,v);}//重心//到三角形三顶点距离的平方和最小的点//三角形内到三边距离之积最大的点point barycenter(point a,point b,point c){line u,v;u.a.x=(a.x+b.x)/2;u.a.y=(a.y+b.y)/2;u.b=c;v.a.x=(a.x+c.x)/2;v.a.y=(a.y+c.y)/2;v.b=b;return intersection(u,v);}//费马点//到三角形三顶点距离之和最小的点point fermentpoint(point a,point b,point c){point u,v;double step=fabs(a.x)+fabs(a.y)+fabs(b.x)+fabs(b.y)+fabs(c.x)+fabs(c.y);int i,j,k;u.x=(a.x+b.x+c.x)/3;u.y=(a.y+b.y+c.y)/3;while (step>1e-10)for (k=0;k<10;step/=2,k++)for (i=-1;i<=1;i++)for (j=-1;j<=1;j++){v.x=u.x+step*i;v.y=u.y+step*j;if(distance(u,a)+distance(u,b)+distance(u,c)>distance(v,a)+distance(v,b)+distance(v,c))u=v;}return u;}1.9三维几何//三维几何函数库#include <math.h>#define eps 1e-8#define zero(x) (((x)>0?(x):-(x))<eps)struct point3{double x,y,z;};struct line3{point3 a,b;};struct plane3{point3 a,b,c;};//计算cross product U x Vpoint3 xmult(point3 u,point3 v){point3 ret;ret.x=u.y*v.z-v.y*u.z;ret.y=u.z*v.x-u.x*v.z;ret.z=u.x*v.y-u.y*v.x;return ret;}//计算dot product U . Vdouble dmult(point3 u,point3 v){return u.x*v.x+u.y*v.y+u.z*v.z;}//矢量差U - Vpoint3 subt(point3 u,point3 v){point3 ret;ret.x=u.x-v.x;ret.y=u.y-v.y;ret.z=u.z-v.z;return ret;}//取平面法向量point3 pvec(plane3 s){return xmult(subt(s.a,s.b),subt(s.b,s.c));}point3 pvec(point3 s1,point3 s2,point3 s3){return xmult(subt(s1,s2),subt(s2,s3));}//两点距离,单参数取向量大小double distance(point3 p1,point3 p2){return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y)+(p1.z-p2.z)*(p1.z-p2.z)); }//向量大小double vlen(point3 p){return sqrt(p.x*p.x+p.y*p.y+p.z*p.z);}//判三点共线int dots_inline(point3 p1,point3 p2,point3 p3){return vlen(xmult(subt(p1,p2),subt(p2,p3)))<eps;}//判四点共面int dots_onplane(point3 a,point3 b,point3 c,point3 d){return zero(dmult(pvec(a,b,c),subt(d,a)));}//判点是否在线段上,包括端点和共线int dot_online_in(point3 p,line3 l){return zero(vlen(xmult(subt(p,l.a),subt(p,l.b))))&&(l.a.x-p.x)*(l.b.x-p.x)<eps&& (l.a.y-p.y)*(l.b.y-p.y)<eps&&(l.a.z-p.z)*(l.b.z-p.z)<eps;}int dot_online_in(point3 p,point3 l1,point3 l2){return zero(vlen(xmult(subt(p,l1),subt(p,l2))))&&(l1.x-p.x)*(l2.x-p.x)<eps&& (l1.y-p.y)*(l2.y-p.y)<eps&&(l1.z-p.z)*(l2.z-p.z)<eps;}//判点是否在线段上,不包括端点int dot_online_ex(point3 p,line3 l){return dot_online_in(p,l)&&(!zero(p.x-l.a.x)||!zero(p.y-l.a.y)||!zero(p.z-l.a.z))&&(!zero(p.x-l.b.x)||!zero(p.y-l.b.y)||!zero(p.z-l.b.z));}int dot_online_ex(point3 p,point3 l1,point3 l2){return dot_online_in(p,l1,l2)&&(!zero(p.x-l1.x)||!zero(p.y-l1.y)||!zero(p.z-l1.z))&& (!zero(p.x-l2.x)||!zero(p.y-l2.y)||!zero(p.z-l2.z));}//判点是否在空间三角形上,包括边界,三点共线无意义int dot_inplane_in(point3 p,plane3 s){return zero(vlen(xmult(subt(s.a,s.b),subt(s.a,s.c)))-vlen(xmult(subt(p,s.a),subt(p,s.b)))- vlen(xmult(subt(p,s.b),subt(p,s.c)))-vlen(xmult(subt(p,s.c),subt(p,s.a))));}int dot_inplane_in(point3 p,point3 s1,point3 s2,point3 s3){return zero(vlen(xmult(subt(s1,s2),subt(s1,s3)))-vlen(xmult(subt(p,s1),subt(p,s2)))- vlen(xmult(subt(p,s2),subt(p,s3)))-vlen(xmult(subt(p,s3),subt(p,s1))));}//判点是否在空间三角形上,不包括边界,三点共线无意义int dot_inplane_ex(point3 p,plane3 s){return dot_inplane_in(p,s)&&vlen(xmult(subt(p,s.a),subt(p,s.b)))>eps&&vlen(xmult(subt(p,s.b),subt(p,s.c)))>eps&&vlen(xmult(subt(p,s.c),subt(p,s.a)))>eps; }int dot_inplane_ex(point3 p,point3 s1,point3 s2,point3 s3){return dot_inplane_in(p,s1,s2,s3)&&vlen(xmult(subt(p,s1),subt(p,s2)))>eps&& vlen(xmult(subt(p,s2),subt(p,s3)))>eps&&vlen(xmult(subt(p,s3),subt(p,s1)))>eps; }//判两点在线段同侧,点在线段上返回0,不共面无意义int same_side(point3 p1,point3 p2,line3 l){return dmult(xmult(subt(l.a,l.b),subt(p1,l.b)),xmult(subt(l.a,l.b),subt(p2,l.b)))>eps;}int same_side(point3 p1,point3 p2,point3 l1,point3 l2){return dmult(xmult(subt(l1,l2),subt(p1,l2)),xmult(subt(l1,l2),subt(p2,l2)))>eps;}//判两点在线段异侧,点在线段上返回0,不共面无意义int opposite_side(point3 p1,point3 p2,line3 l){return dmult(xmult(subt(l.a,l.b),subt(p1,l.b)),xmult(subt(l.a,l.b),subt(p2,l.b)))<-eps;}int opposite_side(point3 p1,point3 p2,point3 l1,point3 l2){return dmult(xmult(subt(l1,l2),subt(p1,l2)),xmult(subt(l1,l2),subt(p2,l2)))<-eps;}//判两点在平面同侧,点在平面上返回0int same_side(point3 p1,point3 p2,plane3 s){return dmult(pvec(s),subt(p1,s.a))*dmult(pvec(s),subt(p2,s.a))>eps;}int same_side(point3 p1,point3 p2,point3 s1,point3 s2,point3 s3){return dmult(pvec(s1,s2,s3),subt(p1,s1))*dmult(pvec(s1,s2,s3),subt(p2,s1))>eps; }//判两点在平面异侧,点在平面上返回0int opposite_side(point3 p1,point3 p2,plane3 s){return dmult(pvec(s),subt(p1,s.a))*dmult(pvec(s),subt(p2,s.a))<-eps;}int opposite_side(point3 p1,point3 p2,point3 s1,point3 s2,point3 s3){return dmult(pvec(s1,s2,s3),subt(p1,s1))*dmult(pvec(s1,s2,s3),subt(p2,s1))<-eps; }//判两直线平行int parallel(line3 u,line3 v){return vlen(xmult(subt(u.a,u.b),subt(v.a,v.b)))<eps;}int parallel(point3 u1,point3 u2,point3 v1,point3 v2){return vlen(xmult(subt(u1,u2),subt(v1,v2)))<eps;}//判两平面平行int parallel(plane3 u,plane3 v){return vlen(xmult(pvec(u),pvec(v)))<eps;}int parallel(point3 u1,point3 u2,point3 u3,point3 v1,point3 v2,point3 v3){ return vlen(xmult(pvec(u1,u2,u3),pvec(v1,v2,v3)))<eps;}//判直线与平面平行int parallel(line3 l,plane3 s){return zero(dmult(subt(l.a,l.b),pvec(s)));}int parallel(point3 l1,point3 l2,point3 s1,point3 s2,point3 s3){return zero(dmult(subt(l1,l2),pvec(s1,s2,s3)));}//判两直线垂直int perpendicular(line3 u,line3 v){return zero(dmult(subt(u.a,u.b),subt(v.a,v.b)));}int perpendicular(point3 u1,point3 u2,point3 v1,point3 v2){return zero(dmult(subt(u1,u2),subt(v1,v2)));}//判两平面垂直int perpendicular(plane3 u,plane3 v){return zero(dmult(pvec(u),pvec(v)));}int perpendicular(point3 u1,point3 u2,point3 u3,point3 v1,point3 v2,point3 v3){ return zero(dmult(pvec(u1,u2,u3),pvec(v1,v2,v3)));}//判直线与平面平行int perpendicular(line3 l,plane3 s){return vlen(xmult(subt(l.a,l.b),pvec(s)))<eps;}int perpendicular(point3 l1,point3 l2,point3 s1,point3 s2,point3 s3){return vlen(xmult(subt(l1,l2),pvec(s1,s2,s3)))<eps;}//判两线段相交,包括端点和部分重合int intersect_in(line3 u,line3 v){if (!dots_onplane(u.a,u.b,v.a,v.b))return 0;if (!dots_inline(u.a,u.b,v.a)||!dots_inline(u.a,u.b,v.b))return !same_side(u.a,u.b,v)&&!same_side(v.a,v.b,u);return dot_online_in(u.a,v)||dot_online_in(u.b,v)||dot_online_in(v.a,u)||dot_online_in(v.b,u); }int intersect_in(point3 u1,point3 u2,point3 v1,point3 v2){if (!dots_onplane(u1,u2,v1,v2))return 0;if (!dots_inline(u1,u2,v1)||!dots_inline(u1,u2,v2))return !same_side(u1,u2,v1,v2)&&!same_side(v1,v2,u1,u2);returndot_online_in(u1,v1,v2)||dot_online_in(u2,v1,v2)||dot_online_in(v1,u1,u2)||dot_online_in(v2,u1,u 2);}//判两线段相交,不包括端点和部分重合int intersect_ex(line3 u,line3 v){return dots_onplane(u.a,u.b,v.a,v.b)&&opposite_side(u.a,u.b,v)&&opposite_side(v.a,v.b,u); }int intersect_ex(point3 u1,point3 u2,point3 v1,point3 v2){returndots_onplane(u1,u2,v1,v2)&&opposite_side(u1,u2,v1,v2)&&opposite_side(v1,v2,u1,u2);}//判线段与空间三角形相交,包括交于边界和(部分)包含int intersect_in(line3 l,plane3 s){return !same_side(l.a,l.b,s)&&!same_side(s.a,s.b,l.a,l.b,s.c)&& !same_side(s.b,s.c,l.a,l.b,s.a)&&!same_side(s.c,s.a,l.a,l.b,s.b);}int intersect_in(point3 l1,point3 l2,point3 s1,point3 s2,point3 s3){ return !same_side(l1,l2,s1,s2,s3)&&!same_side(s1,s2,l1,l2,s3)&& !same_side(s2,s3,l1,l2,s1)&&!same_side(s3,s1,l1,l2,s2);}//判线段与空间三角形相交,不包括交于边界和(部分)包含int intersect_ex(line3 l,plane3 s){return opposite_side(l.a,l.b,s)&&opposite_side(s.a,s.b,l.a,l.b,s.c)&& opposite_side(s.b,s.c,l.a,l.b,s.a)&&opposite_side(s.c,s.a,l.a,l.b,s.b); }int intersect_ex(point3 l1,point3 l2,point3 s1,point3 s2,point3 s3){ return opposite_side(l1,l2,s1,s2,s3)&&opposite_side(s1,s2,l1,l2,s3)&& opposite_side(s2,s3,l1,l2,s1)&&opposite_side(s3,s1,l1,l2,s2);}//计算两直线交点,注意事先判断直线是否共面和平行!//线段交点请另外判线段相交(同时还是要判断是否平行!)point3 intersection(line3 u,line3 v){point3 ret=u.a;double t=((u.a.x-v.a.x)*(v.a.y-v.b.y)-(u.a.y-v.a.y)*(v.a.x-v.b.x))/((u.a.x-u.b.x)*(v.a.y-v.b.y)-(u.a.y-u.b.y)*(v.a.x-v.b.x));ret.x+=(u.b.x-u.a.x)*t;ret.y+=(u.b.y-u.a.y)*t;ret.z+=(u.b.z-u.a.z)*t;return ret;}point3 intersection(point3 u1,point3 u2,point3 v1,point3 v2){point3 ret=u1;double t=((u1.x-v1.x)*(v1.y-v2.y)-(u1.y-v1.y)*(v1.x-v2.x))/((u1.x-u2.x)*(v1.y-v2.y)-(u1.y-u2.y)*(v1.x-v2.x));ret.x+=(u2.x-u1.x)*t;ret.y+=(u2.y-u1.y)*t;ret.z+=(u2.z-u1.z)*t;return ret;}//计算直线与平面交点,注意事先判断是否平行,并保证三点不共线!。
ACM新手入门指南
ACM/ICPC新手入门指南前言：这篇指南不对ACM/ICPC国际大学生程序设计竞赛进行介绍，计算机学子如果不了解的可以在百度上进行搜索查询，这里介绍的只是一个计算机学生想要在ACM/ICPC里进行发展的初学者。
内容比较简单通俗，完全是给新接触的人看的，已经接触过的请飘过，该干嘛的干嘛去。
语言关:要进行程序设计，也就必然要熟悉编程语言，只要掌握了一门语言，就可以进行ACM训练了。
一般通用语言如C、C++、JAVA都可以，这三种语言都有自己的优势和缺点，C在效率方面比较好；但C++封装了输入输出流，方便了我们的操作也减少出错的可能性，而且C++提供了非常强大的标准模版库（STL），使得很多在C上实现起来比较麻烦的代码，在C++上却非常方便；JAVA在大型工程和安全方面都有比较独特的优势，但在ACM里面却不是一种优秀的语言，因为JAVA的执行效率要比C、C++慢很多，如果题目限时比较紧的话，就不适合用JAVA，当然JAVA为我们提供了很方便的高精度运算（大整数运算），所以个人认为，刚学完C的可以用纯C来写训练，在训练过程中可以学学C++，有时间的把STL也好好学学，这样可以减少很多不必要的劳动。
初次接触ACM训练的同学经常会遇到问题，就是输入和输出问题，所以如果对语言的输入输出问题不是很熟悉的话，要抽几天时间重点看看,特别有些初学者在输出时总会输出冗余信息，可能认为有交互性吧，但这是ACM不允许的，它不需要任何交互性。
不严格按照题目要求进行输入输出的程序是无法通过系统测试的。
熟悉在线评测系统在线评测系统，英文叫Online Judge，（简称OJ）里面提供了很多题目给我们平时训练之用。
这里以浙江大学的在线评测系统为例，网址是 先在上面进行注册，注册完后就可以进行题目的训练了，点击主页上的“Problems”，就可以看到里面的题库，可以选任何一个题来做，里面的题目不是由易到难进行排列，而初学者要选择比较简单的题目来做。
ACM培训大纲
实用标准文案ACM培训大纲基础内容：数据结构——》搜索——》图论DP数论博弈中级内容数据结构网络流第一章搜索1.二分搜索三分搜索2.栈 3.队列 4.深搜 5.广搜 6.第二章数据结构1.优先队列并查集 2.二叉搜索树3.线段树（单点更新） 4.Trie5.精彩文档．实用标准文案第三章图论1.图的表示1.1二维数组1.2邻接表1.3前向星2.图的遍历2.1双连通分量2.2拓扑排序3.最短路3.1迪杰斯特拉3.2弗洛伊德3.3SPFA4.匹配匈牙利算法5.生成树6.网络流简介第四章动态规划1.状态转移方程2.引入2.10-1背包2.2硬币问题2.3矩阵链乘3.区间DP4.按位DP5.树形DP6.状压DP第五章数论1.欧几里得扩展欧几里得 2.因数分解3.费马小定理 4.欧拉定理 5.素数6.6.1筛法6.2素数判定6.2.1O(√n)方法精彩文档．实用标准文案6.2.2Miller-rabin测试第六章博弈1.Nim和2.SG函数第七章中级数据结构1.树状数组RMQ 2.KMP3.AC自动机4.线段树（区间更新）5.第八章图论进阶1.网络流问题精彩文档．实用标准文案综述在很多人眼里，东北大学秦皇岛分校不算是985高校。
所以我们要用自己的能力证明我们有985的实力。
ACM是计算机界认可度最高的一个比赛，可以说只要区域赛有过奖牌，国内任何IT公司没有理由不要。
同时，在高校之中，对一个大学计算机专业的评价，大部分人也会首先看ACM 的水平。
将ACM打出学校，在国内打出一定成绩，对扩大我校影响力很有帮助。
考虑到本校暂时没有进行专题训练的出题能力，专题训练的题目主要从UESTC 2014年集训队专题训练中获取，再加上从别的OJ上找一些题目。
训练的平台设置在华中科技大学的vertual judge上面。
本人将在毕业之前承担培训任务。
在2015学年开始之前，培训计划为每两周一次，中间空闲的时间由大二或者大一熟悉C++的同学给不熟悉C++的同学进行基础的讲解。
ACM训练计划建议（转）
ACM训练计划建议（转）前⾔：⽼师要我们整理⼀份训练计划给下⼀届的学弟学妹们，整理出来了，费了不少笔墨，就也将它放到博客园上供⼤家参考。
菜鸟之作，⼤⽜勿喷，如有不当或补充之处，欢迎指出。
本建议书分为三个阶段，⼤⼀、⼤⼆、⼤三。
⼤四暂没整理，⼀⽅⾯是⼤四要⾯临考验和找⼯作的问题，坚持继续acm的很少，另⼀⽅⾯，本⼈还没⼤四……下⾯以个⼈经验分析⼀下这三个阶段建议学习的内容和具体的训练计划。
正⽂：⼤⼀（第⼀阶段）： ⼤⼀是时间最充裕的⼀段时间，也是可塑性最⾼的⼀个阶段。
⼤⼀你有很多⾃由时间可以⾃⼰分配，建议这段时间先打好c/c++基础，或者是任何⼀门语⾔的基础，尽量做到“半精通”。
因为像c++这种语⾔，⼏年内达到精通都是不可能的。
⽽我这⾥所说的“半精通”，实际上是将课本完全搞透的基础上，再在课外拓展⼀些内容，加深对语⾔本⾝的理解。
为什么我这⾥强调语⾔的重要性呢？⼀⽅⾯是为了以后搞acm的需要，语⾔通畅了，可以保证你实现⼤部分算法没有障碍，⽽⽐赛时，时间是很重要的。
另⼀⽅⾯，也是你⾝为学计算机的学⽣的⼀个必须要学习的能⼒，语⾔就是你⼿中的剑，以后能披荆斩棘⾛多远，剑有多锋利占很⼤因素。
另⼀⽅⾯，建议打好数学的基础。
学长⾝为过来⼈，深受数学烂的苦。
acm越到后期的时候，其实⽤到的数学知识就越多。
像有些题⽬，⾚裸裸的就是求积分，还有⼀些题⽬，将求期望嵌⼊到了DP的题⽬⾥…… 其实这些还好说，都是显式的题⽬，还有⼀种隐式的东西，对acm帮助最⼤，那就是数学思维。
有数学思维和没有数学思维的区别，就是⼀道题有N种思路和N*N种思路的区别。
这越到后期越明显，思路很重要，⽐赛时⾯对⼀道题团队⾥⾸先要有多种思路可供分析，然后⼤家讨论哪种思路是最可⾏的，确定之后就是最擅长这个思路的⼈实现算法。
如果分析⼀道题⽬的时候，产⽣的思路很少，那么⽆疑会降低这道题的AC命中率。
啰嗦了这么多，其实我就想给⼤⼀的新⽣们强调两点，语⾔和数学。
ACM训练指南
13级暑期训练指南1.以上《入门经典》均指《算法竞赛入门经典》（第2版），《训练指南》均指《算法竞赛入门经典训练指南》。
2.学习内容为各章节，检查内容为各章最后一节所附uva上的习题。
3.检查方式为三项：3.1每周训练完成后，时间不会早于本周六23:59，统计本周训练章节所附uva习题数（不统计例题，仅统计练习）。
a.《训练指南》上的章节有统计两次的情况，第二次为增量统计。
b.统计和排名方式：把各周训练题数汇总，按汇总题数排名。
3.2每周训练最后阶段为练习赛。
a.时间和场次安排：最少为2场1个半小时（3～4题）的短赛和1场3小时的长赛（5～6题），最多为4场短赛和3场长赛，占用周六一天或周五、周六两天。
b.比赛内容：一般为本周训练章节所附习题中没人做过的。
如果所剩题目不足比赛数量和难度要求，有两种选择题目的方法：一是已结束章节中没人做过的习题，二是从其他OJ和来源选择合适难度的题目。
c.比赛地点：VirtualJudge。
每人在VJ上注册帐号，并且把nickname修改为自己的姓名，方便统计。
d.统计和排名方式：把所有比赛场次比赛结果的做题数和总时间汇总，按汇总的做题数排名，同题数按汇总时间数排名。
3.3各章节做题完成度汇总。
统计范围：a.《入门经典》上的每章为一个单位：包含例题和练习。
b.《训练指南》上各章的每小节的例题和各个难度均为一个单位：即例题、初级、中级、高级。
计算方法：按照各个单位完成题数百分比统计，并汇总，百分比按如下4种情况乘以系数计算。
c.完成度小于19%，系数为0；d.完成度小于41%，但大于等于19%，系数为0.6；e.完成度小于83%，但大于等于41%，系数为1.0；f.完成度大于等于83%，系数为1.2。
统计时间：集训结束之时，组队之前，3.4因训练第一周，甚至第二周我都在北京出差，可能耽误第一二周的检查项（3.1和3.2）。
因此如果我不能回来，那么各人根据训练计划的内容安排自动提前即可。
ACM在线训练方法09
ACM 在线训练方法（请仔细阅读）
一、ACM 训练目标
1、寒假前，有没有小组成员 AC 题目 20+。 2、开学时，有没有小组成员 AC 题目 50+，甚至出现 100+。
二、练习步骤
1、注册注意事项 pku Online Judge 网址：/JudgeOnline/ 在 中的 Online Judge 系统注册，并在学校名中注明“cfz08”。注册名今 后尽量不要改变。一方面便于管理，另一方面由于每个用户的历史记录非常完整，今后可 以作为一个非正式的自我证明材料。 2、做哪些题？ 在线题库中题目数量众多，且题目的难易程度、题目类型都不按顺序来。例如，有特 别简单的，也有超级难的(提交次数与 AC 次数的比值可以一定程度上反映难度)。 (1)先搞定把 1000 号题目(也就是第一个题目)，此题用于确定你掌握了该系统的使用方 法以及对 I/O 的要求。 (2)POJ 题山题海，比赛前能踏实解决 200 个问题的已属牛人。所以，请按照本文后“附 1”的题目分类表，根据你的学习计划做题。 (3)绝不要轻易 copy 网络现成的代码后，就算自己过了此题——毫无意义。切记！ 3、看一道题->独立思考、尝试解决
poj1768,poj1184,poj1872,poj1324,poj204 6,poj1482
我的ACM入门指南
ACM/ICPC入门指南ACM简介：ACM:即国际大学生程序设计竞赛。
ACM国际大学生程序设计竞赛由国际计算机学界著名的ACM学会主办，是世界上规模最大、水平最高的国际大学生程序竞赛，每年举办一次。
可以说，ACM国际大学生程序设计竞赛已成为世界各国大学生最具影响力的国际级计算机类的赛事。
入门篇每个ACMer平时都是靠做题来训练的，现在就以杭州电子科技大学的OJ为例简单介绍一下如何使用OJ系统，在浏览器上输入杭电的网址首先，我们需要注册一个ID，点击Register new ID填写好自己的个人信息之后，点击submit，就可以完成注册了，然后在主页右上角的登陆框输入ID，密码便可以登陆，这个时候点击主页上的 Problem Archive 打开OJ的题库点击题号为 1000 的题目 A + B ProblemProblem Description 一栏写的就是题目的描述，很多初学者一开始对全英文的题目会很不习惯，这是很正常的，多查查字典，很快就可以习惯了，而且也可以提高自己的英语阅读水平、、、、、、、、Input 一栏写的是题目的输入要求，这里提示：每行将会有两个整数 A 和 B，输入进程一直延续到文件读完而结束（也就是要循环输入）Output 一栏写的是题目的输出要求，这里提示：对于每组数据，在一行中输入A+B的值接下来的 Sample Input 就是示例输入数据，而 Sample Output 就是相应的输出结果在读完题目之后，我们就可以在自己电脑上的编译器上编写代码，等编写好代码之后，点击页面中的Submit，接着把代码复制上去，再选择相应的语言完成之后点击Submit，就可以看到OJ系统对我们的程序的判题结果了这里返回了一个 Accepted，也就是说通过了OJ的判断，题目答对了！1000A +B Problem Time Limit: 2000/1000 MS (Java/Others)Memory Limit: 65536/32768 K (Java/Others) Total Submission(s): 158161 Accepted Submission(s): 50186 Problem Description Calculate A + B.Input Each line will contain two integers A and B. Process to end of file. Output For each case, output A + B in one line.Sample Input 1 1Sample Output 2Author HDOJStatistic | Submit | Discuss | Note#include<stdio.h>int main() {int a,b;while(scanf("%d %d",&a,&b)!=EOF)printf("%d\n",a+b);}而OJ可能会返回的结果总共有：Accepted 通过 Wrong Answer 错误 Compilation Error 编译错误Runtime Error 运行错误 Time Limit Exceeded 时间超出限制Presentation Error 输出格式错误 Memory Limit Exceeded 内存超出限制Output Limit Exceeded 输出了多余的东西OJ的使用基本上就是这样，十分的简单，想要成为一个强大的ACMer的话，赶紧多做题吧 ^_^，在这里同时推荐一些OJ给大家： 杭州电子科技大学的OJ，里面有很多简单的入门题，非常适合初学者/JudgeOnline 北京大学的OJ，是目前国内最具权威的OJ 官方的OJ，里面收录了历届亚洲地区赛和世界总决赛的题目/usacogate2美国信息学奥林匹克竞赛提高篇ACM比赛是涉及到很多算法的，大概包括图论、动态规划、计算几何、数论、搜索等等，如何提高自己相关的知识是十分关键的，当然前提是你首先必须得精通一门编程语言，一般ACM正规的比赛上可以选用的语言有C、C++和Java，但由于ACM比赛很看重时间空间等效率，所以大部分情况下Java是不占优势的，而80%以上的参赛选手使用的语言都是C++，因此在这里向大家推荐一本C++的经典书籍《C++ Primer》（第4版），里面详细地介绍了C++的各种语法和STL等知识，什么时候整本书看完了，理解了，你的C++才算达到基本的水平掌握好一门语言之后，就要开始系统地学习算法知识了，这里推荐《算法设计与分析基础》这本算法书设计的算法知识十分的全面，而且讲的知识比较基础，很容易就可以看懂，适合初学者，By the way，如果你是计算机学院的学生的话，学校会在大三（大四）教授算法课的时候发这本书下来。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
ACM练习建议一位高手对我的建议：一般要做到50行以内的程序不用调试、100行以内的二分钟内调试成功.acm 主要是考算法的，主要时间是花在思考算法上，不是花在写程序与debug上。
下面给个计划你练练：第一阶段：练经典常用算法，下面的每个算法给我打上十到二十遍，同时自己精简代码，因为太常用，所以要练到写时不用想，10-15分钟内打完，甚至关掉显示器都可以把程序打出来.1.最短路(Floyd、Dijstra,BellmanFord)2.最小生成树(先写个prim,kruscal要用并查集，不好写)3.大数（高精度）加减乘除4.二分查找. (代码可在五行以内)5.叉乘、判线段相交、然后写个凸包.6.BFS、DFS,同时熟练hash表(要熟，要灵活,代码要简)7.数学上的有：辗转相除（两行内），线段交点、多角形面积公式.8. 调用系统的qsort, 技巧很多，慢慢掌握.9. 任意进制间的转换第二阶段：练习复杂一点，但也较常用的算法。
如：1. 二分图匹配（匈牙利），最小路径覆盖2. 网络流，最小费用流。
3. 线段树.4. 并查集。
5. 熟悉动态规划的各个典型：LCS、最长递增子串、三角剖分、记忆化dp6.博弈类算法。
博弈树，二进制法等。
7.最大团，最大独立集。
8.判断点在多边形内。
9. 差分约束系统.10. 双向广度搜索、A*算法，最小耗散优先.第三阶段：前两个阶段是打基础，第三阶段是锻炼在比赛中可以快速建立模型、想新算法。
这就要平时多做做综合的题型了。
1. 把oibh上的论文看看（大概几百篇的，我只看了一点点，呵呵）。
2. 平时扫扫zoj上的难题啦，别老做那些不用想的题.(中大acm的版主经常说我挑简单的来做:-P )3. 多参加网上的比赛，感受一下比赛的气氛，评估自己的实力.4. 一道题不要过了就算，问一下人，有更好的算法也打一下。
5. 做过的题要记好:-)50题第一类搜索(至少4题)1011 1033 1129 2049 2056 2488 2492 (稍难，也可并查集)第二类最短路(至少3题)1062 1125 1797 2253 2679 Bellman-Ford (难)第三类动态规划(至少6题，2479 and 2593必做)2479 and 2593 1015 1042 (也可贪心) 1141 1050 1080 1221 1260 2411 (稍难) 1276 第四类贪心(至少2题)1065 2054 (难) 1521 2709第五类并查集(至少2题)1861 1182 (难) 1308 2524第六类最小生成树(至少2题, 而且Prim 和Kruskal 至少各用一次)1251 1258 1789 2485第七类二分图(至少3题)1325 1469 2195 (KM 算法或最小费用最大流) (难) 2446 1422 and 2594第八类最大流(至少2题)1087 1459 1149 2516 (最小费用最大流) (难)第九类快速查找(B-Search, Hash and so on) (至少3题)2503 2513 (+Euler回路的判定) 1035 1200 2002第十类数论(至少2题)1061 1142 2262 2407 1811(难) 2447 (难)第十一类线段树(无最少题数要求)2352 (可用简单方法) 2528第十二类计算几何(至少2题，1113凸包算法必做)1113 1292 2148 (难) 2653 1584第十三类高精度(至少3题，1001必做)1001 1047 1131 1503 1504 1060 and 1996 (多项式)SCU1002, 1003, 1004 (/soj)第十四类模拟(至少5题)1029 and 1013 1083 and 2028 2234 and 1067 1012 1026 1068 1120 2271 2632第十五类数学(至少4题)2249 1023 2506 1079 1019 and 1095 1905 and 1064 (二分)POJ分类1、排序1423, 1694, 1723, 1727, 1763, 1788, 1828, 1838, 1840, 2201, 2376, 2377, 2380, 1318, 1877, 1928, 1971, 1974, 1990, 2001, 2002, 2092, 2379,1002（需要字符处理，排序用快排即可）1007（稳定的排序）2159（题意较难懂）2231 2371（简单排序）2388（顺序统计算法）2418（二叉排序树）2、搜索、回溯、遍历1022 1111 1118 1129 1190 1562 1564 1573 1655 2184 2225 2243 2312 2362 2378 23861010,1011,1018,1020,1054,1062,1256,1321,1363,1501，1650,1659,1664,1753,2078,2083,2303,2310,2329简单：1128, 1166, 1176, 1231, 1256, 1270, 1321, 1543, 1606, 1664, 1731, 1742, 1745, 1847,1915, 1950, 2038, 2157, 2182, 2183, 2381, 2386, 2426,不易：1024, 1054, 1117, 1167, 1708, 1746, 1775, 1878, 1903, 1966, 2046, 2197, 2349,推荐：1011, 1190, 1191, 1416, 1579, 1632, 1639, 1659, 1680, 1683, 1691, 1709, 1714, 1753, 1771, 1826, 1855, 1856, 1890, 1924, 1935, 1948, 1979, 1980, 2170, 2288, 2331, 2339,2340,1979（和迷宫类似）1980（对剪枝要求较高）3、历法1008 2080 （这种题要小心）4、枚举1012，1046，1387，1411，2245，2326，2363，2381，1054（剪枝要求较高），1650 （小数的精度问题）5、数据结构的典型算法容易：1182, 1656, 2021, 2023, 2051, 2153, 2227, 2236, 2247, 2352, 2395,不易：1145, 1177, 1195, 1227, 1661, 1834,推荐：1330, 1338, 1451, 1470, 1634, 1689, 1693, 1703, 1724, 1988, 2004, 2010, 2119, 2274, 1125(弗洛伊德算法) ，2421（图的最小生成树）6、动态规划1037 A decorative fence、1050 To the Max、1088 滑雪、1125 Stockbroker Grapevine、1141 Brackets Sequence、1159 Palindrome、1160 Post Office、1163 The Triangle、1458 Common Subsequence、1579 Function Run Fun、1887 Testing the CATCHER、1953 World Cup Noise、2386 Lake Counting7、贪心1042, 1065, 1230, 1323, 1477, 1716, 1784,1328 1755（或用单纯形方法），2054，1017，1328，1862，1922 ，2054，2209，2313，2325，2370。
8、模拟容易：1006, 1008, 1013, 1016, 1017, 1169, 1298, 1326, 1350, 1363, 1676, 1786, 1791, 1835, 1970, 2317, 2325, 2390,不易：1012, 1082, 1099, 1114, 1642, 1677, 1684, 1886,1281 1928 2083 2141 20159、递归166410、字符串处理1488, 1598, 1686, 1706, 1747, 1748, 1750, 1760, 1782, 1790, 1866, 1888, 1896, 1951, 2003,2121, 2141, 2145, 2159, 2337, 2359, 2372, 2406, 2408, 1016 1051 1126 1318 1572 1917 19362039 2083 2136 2271 2317 2330，2121 240311、数论1006,1014,1023,1061,1152,1183,1730,226212、几何有关的题目凸包：1113, 1228, 1794, 2007, 2187,1113 wall，2187 beauty contest容易：1319, 1654, 1673, 1675, 1836, 2074, 2137, 2318,不易：1685, 1687, 1696, 1873, 1901, 2172, 2333,13、任意精度运算、数字游戏、高精度计算1001 1023 1047 1060 1079 1131 1140 1142 1207 1220 1284 1289 1306 1316 1338 1405 1454 15031504 1519 1565 1650 1969 2000 2006 2081 2247 2262 2305 2316 23891001, 1220, 1405, 1503,1001（高精度乘法）2413(高精度加法，还有二分查找)14、概率统计1037,105015、小费用最大流、最大流2195 going home，2400 supervisor, supervisee，1087 a plug for UNIX，1149 PIGS，1273 drainageditches，1274 the perfect stall，1325 machine schedule，1459 power network，2239 selectingcourses16、压缩存储的DP1038 bugs integrated inc，1185 炮兵阵地，2430 lazy cow17、最长公共子串（LCS）1080 human gene functions，1159 palindrome，1458 common subsequence，2192 zipper 18、图论及组合数学2421 Constructing Roads、2369 Permutations、2234 Matches Game、2243 Knight Moves、2249 Binomial Showdown、2255 Tree Recovery、2084 Game of Connections、1906 Three powers、1833 排列、1850 Code、1562 Oil Deposits、1496 Word Index、1306 Combinations、1125 Stockbroker Grapevine、1129 Channel Allocation、1146 ID Codes、1095 Trees Made to Order、找规律2247 Humble Numbers、2309 BST、2346 Lucky tickets、2370 Democracy in danger、2365 Rope、2101 Honey and Milk Land 2028 When Can We Meet?、2084 Game of Connections、1915 Knight Moves、1922 Ride to School、1941 The Sierpinski Fractal、1953 World Cup Noise、1958 Strange Towers of Hanoi、1969 Count on Canton、1806 Manhattan 2025、1809 Regetni、1844 Sum、1870 Bee Breeding、1702 Eva\'s Balance、1728 A flea on a chessboard、1604 Just the Facts、1642 Stacking Cubes、1656 Counting Black、1657 Distance on Chessboard、1662 CoIns、1663 Number Steps、1313 Booklet Printing、1316 Self Numbers、1320 Street Numbers、1323 Game Prediction、1338 Ugly Numbers、1244 Slots of Fun、1250 Tanning Salon、1102 LC-Display、1147 Binary codes、1013 Counterfeit Dollar、19、博弈类1067 取石子游戏、1740 A New Stone Game、2234 Matches Game、1082 Calendar Game 、2348 Euclid\'s Game、2413 How many Fibs?、2419 Forest20、简单、模拟题1001 Exponentiation 、1002 487-3279、1003 Hangover 、1701 Dissatisfying Lift、2301 Beat the Spread!、2304 Combination Lock、2328 Guessing Game、2403 Hay Points 、2406 Power Strings、2339 Rock, Scissors, Paper、2350 Above Average、2218 Does This Make Me Look Fat?、2260 Error Correction、2262 Goldbach\'s Conjecture、2272 Bullseye、2136 Vertical Histogram、2174 Decoding Task、2183 Bovine Math Geniuses、2000 Gold Coins、2014 Flow Layout、2051 Argus、2081 Calendar、1918 Ranking List、1922 Ride to School、1970 The Game、1972 Dice Stacking、1974 The Happy Worm、1978 Hanafuda Shuffle、1979 Red and Black、1617 Crypto Columns、1666 Candy Sharing Game、1674 Sorting by Swapping、1503 Integer Inquiry、1504 Adding Reversed Numbers、1528 Perfection、1546 Basically Speaking、1547 Clay Bully、1573 Robot Motion、1575 Easier Done Than Said?、1581 A Contesting Decision、1590 Palindromes、1454 Factorial Frequencies、1363 Rails、1218 THE DRUNK JAILER、1281 MANAGER、1132 Border、1028 Web Navigation、21、初等数学1003 Hangover、1045 Bode Plot、1254 Hansel and Grethel、1269 Intersecting Lines、1401 Factorial、1410 Intersection、2363 Blocks 、2365 Rope、2242 The Circumference of the Circle、2291 Rotten Ropes、2295 A DP Problem、2126 Factoring a Polynomial、2191 Mersenne Composite Numbers、2196 Specialized Four-Digit Numbers、1914 Cramer\'s Rule、1835 宇航员、1799 Yeehaa!、1607 Deck、1244 Slots of Fun、1269 Intersecting Lines、1299 Polar Explorer、1183 反正切函数的应用、22、匹配1274, 1422, 1469, 1719, 2060, 2239,-------------------------------------------------------------------------------------------经典1011（搜索好题）1012（学会打表）10131019（它体现了很多此类问题的特点）1050（绝对经典的dp）1088（dp好题）1157（花店，经典的dp）1163（怎么经典的dp那么多呀？？？）1328（贪心）1458（最长公共子序列）1647（很好的真题，考临场分析准确和下手迅速）1654（学会多边形面积的三角形求法）1655（一类无根树的dp问题）1804（逆序对）2084（经典组合数学问题）2187（用凸包求最远点对，求出凸包后应该有O(N)的求法，可我就是调不出来）2195（二分图的最佳匹配）2242（计算几何经典）2295（等式处理）2353（dp，但要记录最佳路径）2354（立体解析几何）2362（搜索好题）2410（读懂题是关键）2411（经典dp）趣味1067（很难的数学，但仔细研究，是一片广阔的领域）1147（有O(n)的算法，需要思考）1240（直到一棵树的先序和后序遍历，那么有几种中序遍历呢？dp）1426（是数论吗？错，是图论！）1648（别用计算几何，用整点这个特点绕过精度的障碍吧）1833（找规律）1844（貌似dp或是搜索，其实是道有趣的数学题）1922（贪心，哈哈）22312305（不需要高精度噢）2328（要仔细噢）2356（数论知识）2359（约瑟夫问题变种）2392（有趣的问题）很繁的题100110081087（构图很烦，还有二分图的最大匹配）1128（USACO）124513291550（考的是读题和理解能力）1649（dp）2200（字符串处理+枚举）2358（枚举和避免重复都很烦）2361（仔细仔细再仔细）难题1014（数学证明比较难，但有那种想法更重要）1037（比较难的dp）1405（高精度算法也分有等级之分，不断改进吧）2002（不知道有没有比O(n^2*logn)更有的算法？）2054（极难，很强的思考能力）2085（组合数学）2414（dp，但要剪枝）2415（搜索）2423（计算几何+统计）多解题1002（可以用排序，也可以用统计的方法）1338（搜索和dp都可以）1664（搜索和dp都练一练吧）2082（这可是我讲的题噢）2352（桶排和二叉树都行）Note:1011: 很经典的剪支1014: 难在数学上1017: 严格的数学证明貌似不容易1021: 有点繁,考察对图形进行各种旋转的处理1083: 巧妙的思考角度1150: 分奇偶讨论,lg(n)算法1218: 三行就够了,虽然简单,但也有优劣之别1505: 二分加贪心1654: 做法也许很多吧,本人用有向面积做的1674: 计算圈的个数(算是graph 吧)1700: 数学证明不容易1742: O(m*n)的算法1863: 要耐心地慢慢写…^_^1988: 并查集2051: 堆2078: 不难，但剪支可以做到很好2082::O(n),你想到了吗？2084: 卡特兰数2182: 线段树2195: 最小费用最大流2234: 经典博弈算法2236: 并查集2299: 二分思想2395: Kruskal 最小生成树的拓展2406: KMP2411: 用二进制串来表示状态。