向量的运算

因为必修4没有学好导致计算几何挂掉了，所以来补一补TAT
1.(自由)向量的定义：
起始点任意，方向和长度固定的有向线段。
2.向量的代数表示：
向量a记作$\vec a$
3.向量的几何表示：
可以用一个有向线段来表示向量。
我们现在先要证明平面向量基本定理：
定理内容：
对于平面向量$\vec a = m\vec i + n\vec j$，其中$\vec i,\vec j$为任意两个不共线的向量。
这个是个公理啊TAT
4.向量的坐标表示：
在二维向量中，分别取与x轴y轴平行的两个单位(长度)向量$\vec i,\vec j$。
则由平面向量基本定理可以知道，任意向量都能被$\vec i,\vec j$表示出来。
设$\vec a = m\vec i + n\vec j$
那么我们就可以把向量a写成坐标形式：$\vec a(m,n)$
很好，我们已经有了向量的坐标表示了。

---

向量的运算：
数乘：$\vec a = (m,n),k\vec a = (km,kn);$
加法：$\vec a = (m,n),\vec b = (p,q) \vec c = \vec a + \vec b = (m + p,n + q)$
减法同理。
数量积(点积)的代数定义：$\vec a = (m,n),\vec b = (p,q),\vec a\cdot \vec b = mp + nq;$
数量积(点积)的几何定义：$\vec a \cdot \vec b = |a||b|cos\theta$
其中$|a|,|b|$为$\vec a,\vec b$的模长。
我们如果能证明这两个是等价的，应该就能计算一些很有用的东西了。
参见百度百科：数量积
向量(叉)积的代数定义：$a \times b = (x_a * y_b - y_a * x_b)$
可以写成行列式形式：
$|x_1 y_1|$
$|x_2 y_2|$
这样就好记多了。
实际上，n维方阵的行列式的值代表的是n个向量组成的图形的n维积。
向量（叉）积的几何定义：好像没有几何定义TAT
向量积的(大小)计算式：$\vec a\times\vec b = |a||b|sin\theta$
即a和b的平行四边形的面积。

---

我好像忘记我为什么要写向量的运算了TAT
直线求交点：
定比分点公式：
在直线AB上，对于任意一点P(x,y)，有：

$$\lambda = \vec {AP} / \vec{PB}$$
$$x = \frac {x_a + \lambda x_b}{1 + \lambda},y = \frac {y_a + \lambda y_b}{1 + \lambda}$$
证明可以直接向量证明。

---

求直线AB和CD的交点O，我们可以求OC ： OD，然后根据定比分点来计算O的坐标。
现在问题在于怎么计算比值。
$OC ： OD = S_{\Delta OAB} : S_{\Delta OAC}$
根据叉积求面积即可。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#define Rep(i,n) for(int i = 1;i <= n;i ++)
using namespace std;
const int N = 100005;
struct Point 
{
    double x,y;
    Point(double x = 0.0,double y = 0.0)
    {
        this -> x = x;
        this -> y = y;
    }
}p[N];
Point operator+(const Point &x,const Point &y){return (Point){x.x + y.x,x.y + y.y};}
Point operator-(const Point &x,const Point &y){return (Point){x.x - y.x,x.y - y.y};}
double operator/(const Point &x,const Point &y){return x.x / y.x;}
double dot(const Point &x,const Point &y){return x.x * y.x + x.y * y.y;}
double crs(const Point &x,const Point &y){return x.x * y.y - x.y * y.x;}
Point Get(const Point &x,const Point &y,const Point &p,const Point &q) //计算xy与pq的交点 
{
    double S1 = fabs(crs(x - p,y - x)),S2 = fabs(crs(x - q,y - x));
    return (Point){(p.x * S2 + S1 * q.x) / (S1 + S2),(p.y * S2 + S1 * q.y) / (S1 + S2)};
}

---

向量的运算到此就截止了，如果有其他时间我应该会补上三角函数的公式。