HDU3548 Enumerate the Triangles(优化)

题 目 ：

Problem Description

Little E is doing geometry works. After drawing a lot of points on a plane, he want to enumerate all the triangles which the vertexes are three of the points to find out the one with minimum perimeter. Your task is to implement his work.

 

Input

The input contains several test cases. The first line of input contains only one integer denoting the number of test cases.
The first line of each test cases contains a single integer N, denoting the number of points. (3 <= N <= 1000)
Next N lines, each line contains two integer X and Y, denoting the coordinates of a point. (0 <= X, Y <= 1000)

 

Output

For each test cases, output the minimum perimeter, if no triangles exist, output "No Solution".

 

Sample Input

2

3

0 0

1 1

2 2

4

0 0

0 2

2 1

1 1

 

Sample Output

 

Case 1: No Solution

Case 2: 4.650

题意：

平面上有n(n<=1000)点，问组成的三角形中，周长最小是多少。

优化：

周长c=L1+L2+L3，所以推得①c > 2Li，假设Li的端点为点a、b，则又有Li>=| Xa-Xb |，故②c > 2*| Xa-Xb |。

只是用优化②即可

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
#define INF 1001*1002
struct point{
    int x, y;
}p[1005];
bool cmp(point a, point b){
    return a.x == b.x ? a.y < b.y : a.x < b.x;
}
double dis(int i, int j){
    return sqrt( 1.0*(p[i].x-p[j].x)*(p[i].x-p[j].x) + (p[i].y-p[j].y)*(p[i].y-p[j].y) );
}
int main()
{
    int i, j, t, r, n;
    double IJ, IR, JR, C, minC;
    int cas = 1;
    cin>>t;
    while(t--)
    {
        cin>>n;
        for(i = 1; i <= n; i++)
            scanf("%d%d", &p[i].x, &p[i].y);
        C = minC = INF;
        sort(p+1, p+n+1, cmp);
        for(i = 1; i <= n; i++){
            for(j = i+1; j <= n; j++){
                IJ  = dis(i, j);
                /*周长c=L1+L2+L3，所以推得c > 2Li，假设Li的端点为点a、b，
                  则又有Li>=| Xa-Xb |，故c > 2*| Xa-Xb |*/
                if(minC <= 2*(p[j].x-p[i].x))break;//优化②
                if(minC <= 2*IJ)continue;//优化①

                for(r = j+1; r <= n; r++){
                    if(minC <= 2*(p[r].x-p[j].x))break;//优化②
                    IR = dis(i, r);
                    JR = dis(j, r);
                    if(IJ + IR > JR && (fabs(IJ-IR) < JR))
                        C = IJ + IR + JR;
                    minC = min(C, minC);
                }
            }
        }
        printf("Case %d: ", cas++);
        if(minC == INF)printf("No Solution
");
        else printf("%.3lf
", minC);
    }
    return 0;
}