HDU 4404 Worms(任意多边形与圆相交部分的面积)

题目链接：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;
}