题目：
现在有一个n*m的方阵,方阵里面的数字未知,但是我们知道如下约束条件:1> 每一行的数字的和2> 每一列的数字的和3> 某些格子里的数，大小有限制。比如规定第2
行第3列的数字必须大于5(或必须小于3,或必须等于10等）求解是否存在在满足所有的约束的条件下用正数来填充该方阵的方案,若有,输出填充后的方阵,否则输出IMPOSSIBLE.这道题可以转化成容量有上下界的最大流问题,将方阵的行从1……n编号,列n+1……n+m编号,添加源点s=0和汇点t=n+m+1.1>将源点和每一个行节点相连,相连所形成的边的容量和下界置为该行所有数字的和2>将每一个列节点和汇点相连,相连所形成的边的容量和下界都置为该列所有数字的和
3>从每个行节点到每个列节点连边，容量为无穷大4> 如果u行v列的数字必须大于w,则边<u,v+n>流量的下界是w+15> 如果u行v列的数字必须小于w,则边<u,v+n>容量改为w-16> 如果u行v列的数字必须等于w,则边<u,v+n>流量的下界和容量都是w找到的可行流（也是最大流），就是问题的解本题trick:
1) W可能为负数，产生流量下界为负数的情况。应处理成0
2) 数据本身可能矛盾。比如前面说了 (2,1) =1,后面又说(2,1) = 10

 1 struct CVector///矢量表示
 2 {
 3     double x, y;
 4 };
 5 ///表示从0点到 （x,y）的矢量。对矢量只关心方向
 6 ///和长度，不关心（位置）起点终点
 7
 8 ///矢量的基本运算
 9 CVector operator + (CVector p, CVector q)///矢量和
10 {
11     return CVector(p.x + q.x, p.y + q.y);
12 }
13
14 CVector operator -(CVector p, CVector q)
15 {
16     return CVector(p.x - q.x, p.y - q.y);
17 }
18
19 CVector operator *(double k, CVector p)
20 {
21     return CVector(k * p.x, k * p.y);
22 }
23
24 ///矢量的点积
25 ///性质：p·q=|p||q|cos<p,q>
26 ///a与b的点积，就是a的模乘以b在a上投影的模。若投影与a方向相反则为负值
27 ///其值的正负表示方向
28 ///功能：求同向还是异向；求投影；求出投影后用勾股定理求点到直线距离；
29 ///若a * b = 0,则 a和b垂直
30 double operator *(CVector p, CVector q)
31 {
32     return p.x * q.x + p.y * q.y;
33 }
34
35 ///矢量模长
36 ///用矢量与自身点积求模
37 double length(CVector p)///求矢量的模
38 {
39     return sqrt(p * p);
40 }
41
42 ///矢量单位化
43 ///将矢量除以自身的长度以得到同方向的单位矢量
44 CVector unit(CVector p)
45 {
46     return 1 / length(p) * p;
47 }
48
49 ///矢量的投影长度
50 ///矢量与该方向单位矢量的点积
51 ///注意：负数表示反方向
52 double project(CVector p, CVector n)///点积
53 {
54     return dot(p, unit(n));
55 }
56 double dot(CVector p, CVector q)
57 {
58     return p.x*q.x+p.y*q.y;
59 }
60
61 ///矢量的叉积
62 ///性质：在二维情况中，|p×q|=|p||q|sin<p,q>
63 ///功能：求面积；求顺时针方向还是逆时针方向；判断是否在半平面上
64 ///用法：c = a ^ b; //a,b,c都是CVector对象
65 ///a×b 为有向面积，可正可负。
66 ///若a逆时针旋转小于180度可到b,则结果为正，否则结果为负
67 double operator ^(CVector p, CVector q)
68 {
69     return p.x * q.y – q.x * p.y;
70 }
71
72 ///两个矢量所围成的三角形的面积
73 ///两个矢量的叉积的一半
74 ///注意：得到的面积为有向面积，可能为负
75 double area(CVector p, CVector q)
76 {
77     return p^q / 2;
78 }