POJ 1031 Fence

题目大意：在二维坐标系中，由n个点构成的封闭图形，对于图形上某条直线微元dl，满足关系di=I₀*|cosα|*dl*h，I₀=k/r；其中k和h为一个定值，r为原点到该直线微元的距离，α为两条直线的夹角，两条直线分别为原点到直线微元所在直线的垂线，原点到直线微元的连线。最后求出di的在整个封闭图形上的积分。

思路：

（1）微分公式转化，di=I₀*|cosα|*dl*h=k/r*|cosα|*dl*h=k/r*cosα*(r*dα/sin(90-α))*h=kh*dα；转化细节画个图一目了然。所以问题就转化为求原点对该几何图形的张角。

（2）张角的求法：求出原点对首尾相接的每条线段的张角θ（-PI<=θ<=PI)，进行累加，记录每一步累加的最大值max和最小值min，max-min就是张角，由于题目要求，张角小于等于2*PI。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>

using namespace std;

struct Point
{
    double x,y;
}point[105];

int main()
{
    const double PI=acos(-1.0);
    double k,h;
    double min_angle,max_angle,sum_angle;
    int n,i;
    scanf("%lf%lf%d",&k,&h,&n);
    scanf("%lf%lf",&point[0].x,&point[0].y);
    point[n]=point[0];
    min_angle=max_angle=sum_angle=0.0;
    for (i=1;i<=n;i++)
    {
        if (i<n)
            scanf("%lf%lf",&point[i].x,&point[i].y);
        sum_angle+=fmod((atan2(point[i].y,point[i].x)-atan2(point[i-1].y,point[i-1].x)+5*PI),2*PI)-PI;
        min_angle=min(sum_angle,min_angle);
        max_angle=max(sum_angle,max_angle);
    }
    printf("%.2lf\n",min(max_angle-min_angle,2*PI)*k*h);
    return 0;
}