Ural_1215. Exactness of Projectile Hit

　　/*做了两天，终于AC了。各种参考啊，判段点是否在多边形内是参考的网上的代码；
求点到直线的距离是师兄传授的一套模板，只能说有好的模板就是好啊，自己写了好几遍
没过，用师兄的模板写了一遍过了。

思路：

1、用叉积判断点是否在多边形内。

2、求点到每个边的距离的最小值。*/

//My Code:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>

using namespace std;

const int N = 110;
const double eps = 1e-6;
const double inf = 1e+8;

struct point {
double x;
double y;
}p[N], c;

int n;

double cross_product(point a, point b){
return (c.x - a.x)*(b.y - a.y) - (b.x - a.x)*(c.y - a.y);
}

int in_circle() {
int i; p[n] = p[0]; p[n+1] = p[1];
for(i = 0; i < n; i++)
if(cross_product(p[i], p[i+1]) * cross_product(p[i+1], p[i+2]) < eps)
return false;
return true;
}

double distan(point a, point b) {
return sqrt((a.x - b.x)*(a.x - b.x) + (a.y - b.y)*(a.y - b.y));
}

//好模板啊。。。
double dis(point pa, point pb) {
double A, B, C, l, s;

    A = distan(pb, c);
if(A < eps)    return 0;

    B = distan(pa, c);
if(B < eps)    return 0;

    C = distan(pa, pb);
if(C < eps)    return A;

if(B*B >= A*A + C*C)    return A;
if(A*A >= B*B + C*C)    return B;

    l = (A+B+C)/2;
    s = sqrt(l*(l-A)*(l-B)*(l-C));
return 2*s/C;
}

int main() {
//freopen("data.in", "r", stdin);

int i;
    scanf("%lf%lf%d", &c.x, &c.y, &n);
for(i = 0; i < n; i++)
        scanf("%lf%lf", &p[i].x, &p[i].y);
if(in_circle()){
        printf("0.000\n"); return 0;
    }

double ans = inf; p[n] = p[0];
for(i = 0; i < n; i++)
if(ans - dis(p[i], p[i+1]) > eps)
            ans = dis(p[i], p[i+1]);
    printf("%.3lf\n", 2*ans);
return 0;
}

 

 

转载于:https://www.cnblogs.com/vongang/archive/2011/10/31/2230792.html