题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4404
算是模板题吧,但是这个模板不好搞,也就是没有积累这个模板通常没有
这个题目分两种情况,在上升过程中爆炸,在下降过程中爆炸,求出圆心,然后直接套模板就好了!
#include <iostream>
#include <string.h>
#include <algorithm>
#include <stdio.h>
#include <cmath>
using namespace std;
#define eps 1e-8
#define PI acos(-1.0)
struct Point
{
double x,y;
}po[200];
double x_mult(Point sp, Point ep, Point op)
{
return (sp.x-op.x)*(ep.y-op.y)-(sp.y-op.y)*(ep.x-op.x);
}
double cross(Point a,Point b,Point c)
{
return (a.x-c.x)*(b.x-c.x)+(a.y-c.y)*(b.y-c.y);
}
double dist(Point a,Point b)
{
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double cal_area(Point a,Point b,Point c,double r)
{
double A,B,C,x,y,tS;
A=dist(b,c);
B=dist(a,c);
C=dist(b,a);
if(A<r&&B<r)
return x_mult(a,b,c)/2;
else if(A<r&&B>=r)
{
x=(cross(a,c,b)+sqrt(r*r*C*C-x_mult(a,c,b)*x_mult(a,c,b)))/C;
tS=x_mult(a,b,c)/2;
return asin(tS*(1-x/C)*2/r/B*(1-eps))*r*r/2+tS*x/C;
}
else if(A>=r&&B<r)
{
y=(cross(b,c,a)+sqrt(r*r*C*C-x_mult(b,c,a)*x_mult(b,c,a)))/C;
tS=x_mult(a,b,c)/2;
return asin(tS*(1-y/C)*2/r/A*(1-eps))*r*r/2+tS*y/C;
}
else if(fabs(x_mult(a,b,c))>=r*C||cross(b,c,a)<=0||cross(a,c,b)<=0)
{
if(cross(a,b,c)<0)
if(x_mult(a,b,c)<0)
return (-acos(-1.0)-asin(x_mult(a,b,c)/A/B*(1-eps)))*r*r/2;
else return (acos(-1.0)-asin(x_mult(a,b,c)/A/B*(1-eps)))*r*r/2;
else return asin(x_mult(a,b,c)/A/B*(1-eps))*r*r/2;
}
else
{
x=(cross(a,c,b)+sqrt(r*r*C*C-x_mult(a,c,b)*x_mult(a,c,b)))/C;
y=(cross(b,c,a)+sqrt(r*r*C*C-x_mult(b,c,a)*x_mult(b,c,a)))/C;
tS=x_mult(a,b,c)/2;
return (asin(tS*(1-x/C)*2/r/B*(1-eps))+asin(tS*(1-y/C)*2/r/A*(1-eps)))*r*r/2+tS*((y+x)/C-1);
}
}
int main()
{
double x0,y0,v0,ang,t,g,R,temp,now_t,ans;
Point center;
int n,i,j;
while(scanf("%lf%lf%lf%lf%lf%lf%lf",&x0,&y0,&v0,&ang,&t,&g,&R))
{
if(x0==0 && y0==0 && v0==0 && ang==0 && t==0 && g==0 && R==0)
return 0;
temp=sin(ang/180.0*PI);
center.x=x0,center.y=y0;
center.x+=v0*t*cos(ang/180.0*PI);
if(g*t <= v0)//在上升过程中爆炸
{
center.y+=(v0*temp+(v0*temp-g*t))/2*t;
}
else
{
now_t=(v0*temp)/g;
center.y+=((v0*temp)/2) *now_t;
now_t=t-now_t;
center.y-=(g*now_t*now_t)/2;
}
scanf("%d",&n);
for(i=0;i<n;i++)
scanf("%lf%lf",&po[i].x,&po[i].y);
ans=0;
for(i=0;i<n;i++)
ans+=cal_area(po[i],po[(i+1)%n],center,R);
printf("%.2lf\n",fabs(ans));
}
return 0;
}