方案一:
三角形面积=底*高/2。
方案二:
sinA=sqrt(1-cosA*cosA);
cosA根据余弦定理:CosA=(c^2+b^2-a^2)/2bc;
缺点:中间因为浮点型的转换,浮点型相除,开根号,造成误差偏大。
方案三:海伦公式
S=√[p(p-a)(p-b)(p-c)] 而公式里的p为半周长: p=(a+b+c)/2
缺点:在计算中还是存在误差,但是比方案二的方法有所改观。
方案四:叉积
在二维和三维坐标系中,已知三点就可以知道三角形面积。
利用叉积计算三角形面积:|a×b|=|a||b|sin(a,b)/2 = S△
以UVALive - 5984 为例http://livearchive.onlinejudge.org/index.php
#define N 101
struct node{
char c;
int x1,x2,x3,y1,y2,y3;
int x,y;
int l,r;
}p[100];
int dis(int x1,int y1,int x2,int y2){
return (x1-x2)*(x1-x2) + (y1-y2)*(y1-y2);
}
//利用叉积计算三角形面积:|a×b|=|a||b|sin(a,b)/2 = S
int area(int x1,int y1,int x2,int y2,int x3,int y3){
return abs((x2-x1)*(y3-y1) - (y2-y1)*(x3-x1));
}
bool chk(int i,int j,int k){
if(p[k].c == 'C'){
if(dis(i,j,p[k].x,p[k].y)<=p[k].r*p[k].r)return true;
else return false;
} else if(p[k].c == 'S'){
if(i>=p[k].x && j>=p[k].y && i-p[k].x<=p[k].l && j-p[k].y<=p[k].l)return true;
else return false;
} else {
int a = area(p[k].x1, p[k].y1, p[k].x2, p[k].y2, i,j);
int b = area(p[k].x1, p[k].y1, p[k].x3, p[k].y3, i,j);
int c = area(p[k].x2, p[k].y2, p[k].x3, p[k].y3, i,j);
int d = area(p[k].x1, p[k].y1, p[k].x2, p[k].y2, p[k].x3, p[k].y3);
if(a+b+c == d)return true;
else return false;
}
}
int main(){
int t;
scanf("%d",&t);
while(t--){
char str[3];
int n;
scanf("%d",&n);
int i,j;
for(i=0;i<n;i++){
scanf("%s",str);
p[i].c = str[0];
if(str[0] == 'C'){//圆
scanf("%d%d%d",&p[i].x,&p[i].y,&p[i].r);
} else if(str[0] == 'S'){//正方形
scanf("%d%d%d",&p[i].x,&p[i].y,&p[i].l);
} else {
scanf("%d%d%d%d%d%d",&p[i].x1,&p[i].y1,&p[i].x2,&p[i].y2,&p[i].x3,&p[i].y3);
}
}
int k;
int ans = 0;
for(i=-100;i<=100;i++){
for(j=-100;j<=100;j++){
for(k=0;k<n;k++){
if(chk(i,j,k)){
ans++;
break;
}
}
}
}
printf("%d\n",ans);
}
return 0;
}