qsc52(三角形线性插值)

题目链接：http://qscoj.cn/problem/52/

题意：中文题诶～

思路：水题，只要知道三角形插值和判断点在三角形内就OK了

关于三角形插值：http://www.cnblogs.com/wangchengfeng/p/3453194.html

先判断插入点是否在三角形内，在的话解二元一次方程组：

P.x = (1 - u - v) * P1.x + u * P2.x + v * P3.x

P.y = (1 - u - v) * P1.y + u * P2.y + v * P3.y

再将u, v代入：P = (1 - u - v) * P1 + u * P2 + v * P3

P就是答案了

代码：

 #include <iostream>
 #include <stdio.h>
 #include <math.h>
 using namespace std;

 const double res=1e-;
 struct node{
     double x, y;
 };

 double get_area(node a, node b, node c){ //***用行列式计算三角形面积
     double area=fabs(a.x*b.y+a.y*c.x+b.x*c.y-a.x*c.y-a.y*b.x-b.y*c.x);
     return area;
 }

 bool is_ok(node p1, node p2, node p3, node p){//***判断点p是否在三角形内
     double areap=get_area(p1, p2, p);
     areap+=get_area(p1, p3, p);
     areap+=get_area(p2, p3, p);
     double area=get_area(p1, p2, p3);
     if(fabs(areap-area)<res){
         return true;
     }else{
         return false;
     }
 }

 int main(void){
     double x1, y1, w1, x2, y2, w2, x3, y3, w3, x, y;
     while(cin >> x1 >> y1 >> w1 >> x2 >> y2 >> w2 >> x3 >> y3 >> w3 >> x >> y){
         node p1, p2, p3, p;
         p1.x=x1, p1.y=y1;
         p2.x=x2, p2.y=y2;
         p3.x=x3, p3.y=y3;
         p.x=x, p.y=y;
         if(!is_ok(p1, p2, p3, p)){
             cout << - << endl;
             continue;
         }
         double a1=x2-x1, b1=x3-x1, c1=x-x1;
         double a2=y2-y1, b2=y3-y1, c2=y-y1;
         double u=(b1*c2-b2*c1)/(a2*b1-a1*b2);
         double v=(c1-a1*u)/b1;
         double ans=(-u-v)*w1+u*w2+v*w3;
         printf("%.2lf\n", ans);
     }
 }