HDU NO.3548 Enumerate the Triangles

来源:互联网 发布:广联达预算软件价格 编辑:程序博客网 时间:2024/06/06 02:30

题目大意:

              给定1000个点,求围成三角形的最大周长。

解题思路:

             先把边保存下来,减小后面不必要的复杂度,之后有三个循环1000^3,会超时,有个重要的剪枝,就是在第二个循环里面如果这两个点的长度的两倍大于等于最小周长,那么这个解不是最优,直接减掉,注意,这个剪枝非常重要。是解决这道问题的关键所在。

解题感想:

             合理的剪枝很重要,如果if(2 * edge[i][j] >= Min)这样减就减到了500+ms,如果if(edge[i][j] >= Min)这样减就飙到了900+ms。

 

 

#include<iostream>#include<cmath>using namespace std;const int MAX = 1005;typedef struct point{    int x, y;} P;P po[MAX];double edge[MAX][MAX];inline double dealTri(int a, int b, int c){    return edge[a][b] + edge[a][c] + edge[b][c];}inline bool judge(int a, int b, int c){    if((edge[a][c] + edge[a][b] > edge[b][c]) && (edge[a][c] + edge[b][c] > edge[a][b]) && (edge[a][b] + edge[b][c] > edge[a][c]))        return true;    else        return false;}int main(void){    int cas, cas_c = 1;    scanf("%d", &cas);    while(cas--)    {        memset(edge, 0, sizeof(edge));        int n;        scanf("%d", &n);        scanf("%d %d", &po[0].x, &po[0].y);        for(int i = 1; i < n; i++)        {            scanf("%d %d", &po[i].x, &po[i].y);        }        for(int i = 0; i < n; i++)        {            for(int j = i + 1; j < n; j++)            {                edge[i][j] = edge[j][i] = sqrt((double)((po[i].x - po[j].x) * (po[i].x - po[j].x) + (po[i].y - po[j].y) * (po[i].y - po[j].y)));            }        }        double Min = 10000000000;        bool flag = false;        for(int i = 0; i < n; i++)        {            for(int j = i + 1; j < n; j++)            {                if(edge[i][j] >= Min)                    continue;                for(int k = j + 1; k < n; k++)                {                    if(judge(i, j, k))//是三角形                    {                        flag = true;                        double temp = dealTri(i, j, k);                        if(temp < Min)                            Min = temp;                    }                }            }        }        printf("Case %d: ", cas_c++);        if(flag)            printf("%.3lf\n", Min);        else            printf("No Solution\n");    }    return 0;}


 

0 0