poj ACM1031

原题翻译如：
在一些平坦的场地上有一个被围栏围住的区域。篱笆具有高度h，在平面投影中，它具有由其N个顶点的笛卡尔坐标（Xi，Yi）指定的封闭多边形线（无自相交）的形式。在坐标（0，0）的点上，一个灯站在场上。灯可以位于栅栏的外部或内部，但不在其侧面，如下图所示（图中细线所示的部分不被灯照亮）：

栅栏是完全黑色的，即它既不反射也不扩散，也不能让光通过。研究和实验表明，以下法律表达了落在该栅栏任意照明点的光强度：
I 0 = k / r

其中k是不取决于所讨论点的已知常数值，r是平面投影中该点与灯之间的距离。具有宽度d1和高度h的无穷小窄垂直板的照明
dI = I 0 * |cosα| * dl * h

其中I 0是栅栏板上的光的强度，α是在该点处的栅栏侧面的法线与灯的方向之间的平面投影中的角度。
您将编写一个程序，将找到围墙的总照度，该照明被定义为所有照明板的照明总和。
输入

输入文件的第一行包含数字k，h和N，以空格分隔。k和h是实常数。N（3 <= N <= 100）是围栏的顶点数。然后N行跟随，每行包含两个实数Xi和Yi，由空格分隔。
产量

将输入文件写入栅格的总照明四舍五入到小数点后的第二位。
样品输入

0.5 1.7 3
1.0 3.0
2.0 -1.0
-4.0 -1.0
样品输出

5.34
代码如下：

#include<stdio.h>  
#include<string.h>  
#include<algorithm>  
#include<math.h>  
//#include<queue>  
#define dist(a,b) sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y))  
#define cross(a,b,c) (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x)  
#define dot(a,b,c) (b.x-a.x)*(c.x-a.x)+(b.y-a.y)*(c.y-a.y)  
#define delt(a) fabs(a)<eps?0:a>0?1:-1  
#define pi acos(-1.0)  
//#define eps 1e-8  
//#define inf 1e20  
#define N 1005  
using namespace std;
int n, m, t, nl, ml;
struct TPoint
{
    double x, y;
}pt[N], st;
void scan()
{
    scanf("%lf%lf%d", &st.x, &st.y, &n);
    for (int i = 0; i<n; i++)
        scanf("%lf%lf", &pt[i].x, &pt[i].y);
}
double getang(TPoint a, TPoint b)
{
    double ang1 = atan2(a.y, a.x), ang2 = atan2(b.y, b.x);
    if (ang1 - ang2>pi) 
        ang2 += 2 * pi;
    if (ang2 - ang1>pi) 
        ang1 += 2 * pi;
    return ang1 - ang2;
}
void solve()
{
    double minf = 0, maxf = 0, sumf = 0;
    pt[n] = pt[0];
    for (int i = 0; i<n; i++)
    {
        sumf += getang(pt[i], pt[i + 1]);
        maxf = max(sumf, maxf);
        minf = min(sumf, minf);
        if (maxf - minf>2 * pi)
        {
            maxf = minf + 2 * pi;
            break;
        }
    }
    printf("%.2f\n", (maxf - minf)*st.x*st.y);
}
int main()
{
    scan();
    solve();
    return 0;
}