c算法几何

① 计算几何的全部算法

1. 矢量减法

设二维矢量 P = （x1,y1） ，Q = (x2,y2)
则矢量减法定义为： P - Q = ( x1 - x2 , y1 - y2 )
显然有性质 P - Q = - ( Q - P )
如不加说明，下面所有的点都看作矢量，两点的减法就是矢量相减；

2.矢量叉积

设矢量P = （x1,y1） ，Q = (x2,y2)
则矢量叉积定义为： P × Q = x1*y2 - x2*y1 得到的是一个标量
显然有性质 P × Q = - ( Q × P ) P × ( - Q ) = - ( P × Q )
如不加说明，下面所有的点都看作矢量，点的乘法看作矢量叉积；

叉乘的重要性质：

> 若 P × Q > 0 , 则P 在Q的顺时针方向
> 若 P × Q < 0 , 则P 在Q的逆时针方向
> 若 P × Q = 0 , 则P 与Q共线，但可能同向也可能反向

3.判断点在线段上

设点为Q，线段为P1P2 ，判断点Q在该线段上的依据是：

( Q - P1 ) × ( P2 - P1 ) = 0 且 Q 在以 P1，P2为对角顶点的矩形内

4.判断两线段是否相交

我们分两步确定两条线段是否相交：

(1)． 快速排斥试验

设以线段 P1P2 为对角线的矩形为R， 设以线段 Q1Q2 为对角线的矩形为T，如果
R和T不相交，显然两线段不会相交；

(2)． 跨立试验

如果两线段相交，则两线段必然相互跨立对方，如图1所示。在图1中，P1P2跨立
Q1Q2 ，则矢量 ( P1 - Q1 ) 和( P2 - Q1 )位于矢量( Q2 - Q1 ) 的两侧，即
( P1 - Q1 ) × ( Q2 - Q1 ) * ( P2 - Q1 ) × ( Q2 - Q1 ) < 0
上式可改写成
( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) > 0
当( P1 - Q1 ) × ( Q2 - Q1 ) = 0 时，说明( P1 - Q1 ) 和 ( Q2 - Q1 )共线，
但是因为已经通过快速排斥试验，所以 P1 一定在线段 Q1Q2上；同理，
( Q2 - Q1 ) ×( P2 - Q1 ) = 0 说明 P2 一定在线段 Q1Q2上。

所以判断P1P2跨立Q1Q2的依据是：

( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) ≥ 0

同理判断Q1Q2跨立P1P2的依据是：

( Q1 - P1 ) × ( P2 - P1 ) * ( P2 - P1 ) × ( Q2 - P1 ) ≥ 0

至此已经完全解决判断线段是否相交的问题。

5.判断线段和直线是否相交

如果线段 P1P2和直线Q1Q2相交，则P1P2跨立Q1Q2，即：

( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) ≥ 0

6.判断矩形是否包含点

只要判断该点的横坐标和纵坐标是否夹在矩形的左右边和上下边之间。

6.判断线段、折线、多边形是否在矩形中

因为矩形是个凸集，所以只要判断所有端点是否都在矩形中就可以了。

7.判断矩形是否在矩形中

只要比较左右边界和上下边界就可以了。

8.判断圆是否在矩形中

圆在矩形中的充要条件是：圆心在矩形中且圆的半径小于等于圆心到矩形四边的距
离的最小值。

9.判断点是否在多边形中

以点P为端点，向左方作射线L，由于多边形是有界的，所以射线L的左端一定在多
边形外，考虑沿着L从无穷远处开始自左向右移动，遇到和多边形的第一个交点的
时候，进入到了多边形的内部，遇到第二个交点的时候，离开了多边形，……所
以很容易看出当L和多边形的交点数目C是奇数的时候，P在多边形内，是偶数的话
P在多边形外。

但是有些特殊情况要加以考虑。如果L和多边形的顶点相交，有些情况下交点只能
计算一个，有些情况下交点不应被计算（自己画个图就明白了）；如果L和多边形
的一条边重合，这条边应该被忽略不计。为了统一起见，我们在计算射线L和多边
形的交点的时候，1。对于多边形的水平边不作考虑；2。对于多边形的顶点和L相
交的情况，如果该顶点是其所属的边上纵坐标较大的顶点，则计数，否则忽略；
3。对于P在多边形边上的情形，直接可判断P属于多边行。由此得出算法的伪代码
如下：

1. count ← 0;
2. 以P为端点，作从右向左的射线L;
3. for 多边形的每条边s
4. do if P在边s上
5. then return true;
6. if s不是水平的
7. then if s的一个端点在L上且该端点是s两端点中纵坐标较大的端点
9. then count ← count+1
10. else if s和L相交
11. then count ← count+1;
12. if count mod 2 = 1
13. then return true
14. else return false;

其中做射线L的方法是：设P'的纵坐标和P相同，横坐标为正无穷大（很大的一个正
数），则P和P'就确定了射线L。这个算法的复杂度为O(n)。

10.判断线段是否在多边形内

线段在多边形内的一个必要条件是线段的两个端点都在多边形内；

如果线段和多边形的某条边内交（两线段内交是指两线段相交且交点不在两线段的
端点），因为多边形的边的左右两侧分属多边形内外不同部分，所以线段一定会有
一部分在多边形外。于是我们得到线段在多边形内的第二个必要条件：线段和多边
形的所有边都不内交；

线段和多边形交于线段的两端点并不会影响线段是否在多边形内；但是如果多边形
的某个顶点和线段相交，还必须判断两相邻交点之间的线段是否包含与多边形内部。
因此我们可以先求出所有和线段相交的多边形的顶点，然后按照X-Y坐标排序，这样
相邻的两个点就是在线段上相邻的两交点，如果任意相邻两点的中点也在多边形内，
则该线段一定在多边形内。证明如下：

命题1：

如果线段和多边形的两相邻交点P1 ，P2的中点P' 也在多边形内，则P1, P2之间的
所有点都在多边形内。

证明：

假设P1,P2之间含有不在多边形内的点，不妨设该点为Q，在P1, P'之间，因为多边
形是闭合曲线，所以其内外部之间有界，而P1属于多边行内部，Q属于多边性外部，
P'属于多边性内部，P1-Q-P'完全连续，所以P1Q和QP'一定跨越多边形的边界，因此
在P1,P'之间至少还有两个该线段和多边形的交点，这和P1P2是相邻两交点矛盾，故
命题成立。证毕

由命题1直接可得出推论：

推论2：

设多边形和线段PQ的交点依次为P1,P2,……Pn，其中Pi和Pi+1是相邻两交点，线段
PQ在多边形内的充要条件是：P，Q在多边形内且对于i =1, 2,……, n-1，Pi ,Pi+1
的中点也在多边形内。

在实际编程中，没有必要计算所有的交点，首先应判断线段和多边形的边是否内交
，倘若线段和多边形的某条边内交则线段一定在多边形外；如果线段和多边形的每
一条边都不内交，则线段和多边形的交点一定是线段的端点或者多边形的顶点，只
要判断点是否在线段上就可以了。

至此我们得出算法如下：

1. if 线端PQ的端点不都在多边形内
2. then return false;
3. 点集pointSet初始化为空;
4. for 多边形的每条边s
5. do if 线段的某个端点在s上
6. then 将该端点加入pointSet;
7. else if s的某个端点在线段PQ上
8. then 将该端点加入pointSet;
9. else if s和线段PQ相交 // 这时候可以肯定是内交
10. then return false;
11. 将pointSet中的点按照X-Y坐标排序，X坐标小的排在前面，
对于X坐标相同的点，Y坐标小的排在前面；
12. for pointSet中每两个相邻点 pointSet[i] , pointSet[ i+1]
13. do if pointSet[i] , pointSet[ i+1] 的中点不在多边形中
14. then return false;
15. return true;

这个算法的复杂度也是O(n)。其中的排序因为交点数目肯定远小于多边形的顶点数
目n，所以最多是常数级的复杂度，几乎可以忽略不计。

11.判断折线在多边形内

只要判断折线的每条线段是否都在多边形内即可。设折线有m条线段，多边形有n个
顶点，则复杂度为O(m*n)。

12.判断多边形是否在多边形内

只要判断多边形的每条边是否都在多边形内即可。判断一个有m个顶点的多边形是
否在一个有n个顶点的多边形内复杂度为O(m*n)。

13.判断矩形是否在多边形内

将矩形转化为多边形，然后再判断是否在多边形内。

14.判断圆是否在多边形内

只要计算圆心到多边形的每条边的最短距离，如果该距离大于等于圆半径则该圆在
多边形内。计算圆心到多边形每条边最短距离的算法在后文阐述。

15.判断点是否在圆内

计算圆心到该点的距离，如果小于等于半径则该点在圆内。

16.判断线段、折线、矩形、多边形是否在圆内

因为圆是凸集，所以只要判断是否每个顶点都在圆内即可。

17.判断圆是否在圆内

设两圆为O1,O2，半径分别为r1, r2，要判断O2是否在O1内。先比较r1，r2的大小
，如果r1<r2则O2不可能在O1内；否则如果两圆心的距离大于r1 - r2 ，则O2不在
O1内；否则O2在O1内。

18.计算点到线段的最近点

如果该线段平行于X轴（Y轴），则过点point作该线段所在直线的垂线，垂足很容
易求得，然后计算出垂足，如果垂足在线段上则返回垂足，否则返回离垂足近的端
点；

如果该线段不平行于X轴也不平行于Y轴，则斜率存在且不为0。设线段的两端点为
pt1和pt2，斜率为：
k = ( pt2.y - pt1. y ) / (pt2.x - pt1.x );
该直线方程为：
y = k* ( x - pt1.x) + pt1.y
其垂线的斜率为 - 1 / k，
垂线方程为：
y = (-1/k) * (x - point.x) + point.y
联立两直线方程解得：
x = ( k^2 * pt1.x + k * (point.y - pt1.y ) + point.x ) / ( k^2 + 1)
y = k * ( x - pt1.x) + pt1.y;

然后再判断垂足是否在线段上，如果在线段上则返回垂足；如果不在则计算两端点
到垂足的距离，选择距离垂足较近的端点返回。

19.计算点到折线、矩形、多边形的最近点

只要分别计算点到每条线段的最近点，记录最近距离，取其中最近距离最小的点即
可。

20.计算点到圆的最近距离

如果该点在圆心，则返回UNDEFINED
连接点P和圆心O，如果PO平行于X轴，则根据P在O的左边还是右边计算出最近点的
横坐标为centerPoint.x - radius 或 centerPoint.x + radius， 如图4 (a)所示；
如果PO平行于Y轴，则根据P在O的上边还是下边计算出最近点的纵坐标为
centerPoint.y + radius 或 centerPoint.y - radius， 如图4 (b)所示。

如果PO不平行于X轴和Y轴，则PO的斜率存在且不为0，如图4(c)所示。这时直线PO
斜率为
k = （ P.y - O.y ）/ ( P.x - O.x )
直线PO的方程为：
y = k * ( x - P.x) + P.y
设圆方程为:
(x - O.x ) ^2 + ( y - O.y ) ^2 = r ^2，
联立两方程组可以解出直线PO和圆的交点，取其中离P点较近的交点即可。

21.计算两条共线的线段的交点

对于两条共线的线段，它们之间的位置关系有图5所示的几种情况。
图5(a)中两条线段没有交点；图5 (b) 和 (d) 中两条线段有无穷焦点；图5 (c)
中两条线段有一个交点。设line1是两条线段中较长的一条，line2是较短的一条，
如果line1包含了line2的两个端点，则是图5(d)的情况，两线段有无穷交点；如
果line1只包含line2的一个端点，那么如果line1的某个端点等于被line1包含的
line2的那个端点，则是图5(c)的情况，这时两线段只有一个交点，否则就是
图5(c)的情况，两线段也是有无穷的交点；如果line1不包含line2的任何端点，
则是图5(a)的情况，这时两线段没有交点。

22.计算线段或直线与线段的交点

设一条线段为L0 = P1P2，另一条线段或直线为L1 = Q1Q2 ，要计算的就是L0和L1
的交点。

1.首先判断L0和L1是否相交（方法已在前文讨论过），如果不相交则没有交点，
否则说明L0和L1一定有交点，下面就将L0和L1都看作直线来考虑。

2.如果P1和P2横坐标相同，即L0平行于Y轴
a)若L1也平行于Y轴，
i.若P1的纵坐标和Q1的纵坐标相同，说明L0和L1共线，假如L1是直线的话他们有
无穷的交点，假如L1是线段的话可用"计算两条共线线段的交点"的算法求他们
的交点（该方法在前文已讨论过）；
ii.否则说明L0和L1平行，他们没有交点；
b)若L1不平行于Y轴，则交点横坐标为P1的横坐标，代入到L1的直线方程中可以计
算出交点纵坐标；
3.如果P1和P2横坐标不同，但是Q1和Q2横坐标相同，即L1平行于Y轴，则交点横
坐标为Q1的横坐标，代入到L0的直线方程中可以计算出交点纵坐标；
4.如果P1和P2纵坐标相同，即L0平行于X轴
a)若L1也平行于X轴，
i.若P1的横坐标和Q1的横坐标相同，说明L0和L1共线，假如L1是直线的话他们
有无穷的交点，假如L1是线段的话可用"计算两条共线线段的交点"的算法求
他们的交点（该方法在前文已讨论过）；
ii.否则说明L0和L1平行，他们没有交点；

b)若L1不平行于X轴，则交点纵坐标为P1的纵坐标，代入到L1的直线方程中可以计
算出交点横坐标；
5.如果P1和P2纵坐标不同，但是Q1和Q2纵坐标相同，即L1平行于X轴，则交点纵坐标
为Q1的纵坐标，代入到L0的直线方程中可以计算出交点横坐标；
6.剩下的情况就是L1和L0的斜率均存在且不为0的情况
a)计算出L0的斜率K0，L1的斜率K1 ；
b)如果K1 = K2
i.如果Q1在L0上，则说明L0和L1共线，假如L1是直线的话有无穷交点，假如L1
是线段的话可用"计算两条共线线段的交点"的算法求他们的交点（该方法在
前文已讨论过）；
ii.如果Q1不在L0上，则说明L0和L1平行，他们没有交点。
c)联立两直线的方程组可以解出交点来

说明：这个算法并不复杂，但是要分情况讨论清楚，尤其是当两条线段共线的情况
需要单独考虑，所以在前文将求两条共线线段的算法单独写出来。另外，一开始就
先利用矢量叉乘判断线段与线段（或直线）是否相交，如果结果是相交，那么在后
面就可以将线段全部看作直线来考虑。

23.求线段或直线与折线、矩形、多边形的交点

分别求与每条边的交点即可。

24.求线段或直线与圆的交点

设圆心为O，圆半径为r，直线（或线段）L上的两点为P1,P2。
1.如果L是线段且P1，P2都包含在圆O内，则没有交点；否则进行下一步
2.如果L平行于Y轴，
a)计算圆心到L的距离dis
b)如果dis > r 则L和圆没有交点；
c)利用勾股定理，可以求出两交点坐标，如图6(a)所示；但要注意考虑L和圆的相
切情况
3.如果L平行于X轴，做法与L平行于Y轴的情况类似；
4.如果L既不平行X轴也不平行Y轴，可以求出L的斜率K，然后列出L的点斜式方程
，和圆方程联立即可求解出L和圆的两个交点；
5.如果L是线段，对于2，3，4中求出的交点还要分别判断是否属于该线段的范围内。

② 关于C语言算法

从易到难有：
排序：例如冒泡、快排、分治等。http://ke..com/view/58783.htm
搜索：如深搜、广搜等。http://ke..com/view/3688332.htm
动态规划，
贪心，
组合数学：http://ke..com/view/44868.htm
二叉树，
图论，
数论，
计算几何。
《算法导论》博大精深，不适合初学者或者适合肯下苦功夫的初学者。

③ C语言算斜率小问题（纯数学几何运算问题）

其实这个很简单的问题
斜率为小数的可能最大吧，呵呵
但是小数的存储是有误差的，所以比较起来是不准确的，这个能理解吧
还有一个问题就是分母可能为零，也就是竖线，那么斜率就是无穷大了，这样在计算上是一个分母为零的异常吧
所以不能用小数来表示斜率，并且判断不能出现小数
所以就要同分，

你考虑的变号问题，同分时需要变号（如果是负数）在代数中是不存在的
这个你可以询问一下老师或者同学，并查一下代数课本，初中的知识吧，呵呵

④ c语言算法有哪些

这里整理c语言常用算法，主要有：
交换算法
查找最小值算法
冒泡排序
选择排序
插入排序
shell排序 (希尔排序)
归并排序
快速排序
二分查找算法
查找重复算法

⑤ c语言算法

离散数学离散数学作为计算机学科的基础是竞赛中涉及最多的数学分支，重中之重又在于图论和组合数学，尤其是图论。图论之所以运用最多是因为它的变化最多，而且可以轻易地结合基本数据结构和许多算法的基本思想，较多用到的知识包括连通性判断、DFS和BFS，关节点和关键路径、欧拉回路、最小生成树、最短路径、二部图匹配和网络流等等。虽然这部分的比重很大，但是往往也是竞赛中的难题所在，如果有初学者对于这部分的某些具体内容暂时感到力不从心，也不必着急，可以慢慢积累。组合数学竞赛中设计的组合计数问题大都需要用组合数学来解决，组合数学中的知识相比于图论要简单一些，很多知识对于小学上过奥校的同学来说已经十分熟悉，但是也有一些部分需要先对代数结构中的群论有初步了解才能进行学习。组合数学在竞赛中很少以难题的形式出现，但是如果积累不够，任何一道这方面的题目却都有可能成为难题。数论以素数判断和同余为模型构造出来的题目往往需要较多的数论知识来解决，这部分在竞赛中的比重并不大，但只要来上一道，也足以使知识不足的人冥思苦想上一阵时间。素数判断和同余最常见的是在以密码学为背景的题目中出现，在运用密码学常识确定大概的过程之后，核心算法往往要涉及数论的内容。计算几何计算几何相比于其它部分来说是比较独立的，就是说它和其它的知识点很少有过多的结合，较常用到的部分包括—线段相交的判断、多边形面积的计算、内点外点的判断、凸包等等。计算几何的题目难度不会很大，但也永远不会成为最弱的题。线性代数对线性代数的应用都是围绕矩阵展开的，一些表面上是模拟的题目往往可以借助于矩阵来找到更好的算法。 ~

⑥ 设计计算几何算法的编程语言是什么，也是c、c++之类吗

你问的这个很片面啊，如果是C语言中的 sqrt()函数，就是开方，当然是用C语言编的，C++也有很多同样的啊，像VB JAVA都有算法啊，windows就是用多种语言实现不同的功能的。。每种语言都有他自己的算法，也有自己的区别 ，就像C是面向过程，而C++是面向对象，C++有继承和多态等性质什么的，

⑦ C语言中 什么是算法 算法的表示有哪几种方式

算法（Algorithm）是一系列解决问题的清晰指令。
算法也可以理解为有基本运算及规定的运算顺序所构成的完整的解题步骤。或者看成按照要求设计好的有限的确切的计算序列，并且这样的步骤和序列可以解决一类问题。
一个算法应该具有以下五个重要的特征：
有穷性，确切性，输入，输出，可行性。
算法可以使用自然语言、伪代码、流程图，或者程序语言（比如C，C++）等多种不同的方法来描述。

⑧ C语言求圆周长和面积

一、数学公式：

圆周长=2*π*半径

面积=π*半径²

二、算法分析：

周长和面积都依赖半径，所以要先输入半径值，然后套用公式，计算周长和面积。 最终输出结果即可。

三、参考代码：

代码如下

#include"stdio.h"

#definePi3.14

voidmain()

{

floatr,c,area;

printf("请输入圆的半径：");

scanf("%f",&r);

c=2*Pi*r;

area=Pi*r*r;

printf("该圆的周长是%.2f，面积是%.2f ",c,area);

}

⑨ 详解C语言算法

完整的程序如下：
typedef int status;
#define OK 1
#define ERROR 0
#define TRUE 1
#define FALSE 0
#define MAXSIZE 10
#include<stdarg.h>
#include<iostream.h>
#include<stdlib.h>
typedef int ElemType;
#define MAXDIM 8
typedef struct{
ElemType * base;
int dim;
int *bounds;
int *constants;
}Array;

status InitArray(Array &A,int dim,...)
{va_list ap;
if(dim<1||dim>MAXDIM) return ERROR;
A.dim=dim;
int i;
int total=1;
A.bounds=(int *)malloc(dim*sizeof(int));
if(!A.bounds) return ERROR;
va_start(ap,dim);
for(i=0;i<dim;i++)
{A.bounds[i]=va_arg(ap,int);
if(A.bounds[i]<0) return ERROR;
total*=A.bounds[i];
}
va_end(ap);
A.base=(ElemType*)malloc(total*sizeof(ElemType));
if(!A.base) return ERROR;
A.constants=(int *)malloc(dim*sizeof(int));
if(!A.constants) return ERROR;
A.constants[dim-1]=1;
for(i=dim-2;i>=0;i--)
A.constants[i]=A.constants[i+1]*A.bounds[i+1];
return OK;
}

status DestoryArray(Array &A)
{if(!A.base) return ERROR;
free(A.base);
if(!A.bounds) return ERROR;
free(A.bounds);
if(!A.constants) return ERROR;
free(A.constants);
A.dim=0;
return OK;
}

status Locate(Array A,int & off,...)
{va_list ap;
va_start(ap,off);
off=0;
int ind;
int i;
for(i=0;i<A.dim;i++)
{ind=va_arg(ap,int);
if(ind>A.bounds[i]||ind<0) return ERROR;
off+=ind*A.constants[i];
}
va_end(ap);
return OK;
}

status Value(Array A,ElemType & e,...)
{int off;
va_list ap;
va_start(ap,e);
Locate(A,off,ap);
e=A.base[off];
va_end(ap);
return OK;
}

status Assign(Array A,ElemType e,...)
{int off;
va_list ap;
va_start(ap,e);
Locate(A,off,ap);
A.base[off]=e;
va_end(ap);
return OK;
}
int main()
{Array A;
InitArray(A,2,2,5);
int i;
for(i=0;i<10;i++)
{ cout<<"请输入数据"<<endl;
cin>>A.base[i];
}

ElemType e;
int off;
Locate(A,off,2,3);
cout<<"偏移是"<<off<<endl;
Value(A,e,2,4);
cout<<"第二行第四列是"<<e<<endl;
cin>>e;
Assign(A,e,2,2);
Value(A,e,2,2);
cout<<"第二行第二列是"<<e<<endl;
DestoryArray(A);
return 0;
}

执行InitArray(A,2,2,5)后
A.bounds[0]=2, A.bounds[1]=5
A.constants[1]=1, A.constants[0]=5
并且为A.base分配了10个ElemType型的元素空间

PS:关键是搞懂省略号参数列表的用法