hdu 1451 Area in Triangle(计算几何 三角形)

Given a triangle field and a rope of a certain length (Figure-1), you are required to use the rope to enclose a region within the field and make the region as large as possible.

Input

The input has several sets of test data. Each set is one line containing four numbers separated by a space. The first three indicate the lengths of the edges of the triangle field, and the fourth is the length of the rope. Each of the four numbers have exactly four digits after the decimal point. The line containing four zeros ends the input and should not be processed. You can assume each of the edges are not longer than 100.0000 and the length of the rope is not longer than the perimeter of the field.
Output

Output one line for each case in the following format: 

Case i: X 

Where i is the case number, and X is the largest area which is rounded to two digits after the decimal point. 

 Sample Input

12.0000 23.0000 17.0000 40.0000

84.0000 35.0000 91.0000 210.0000

100.0000 100.0000 100.0000 181.3800

0 0 0 0

 Sample Output

Case 1: 89.35

Case 2: 1470.00

Case 3: 2618.00

题意:给你一个三角形的三个边,再给你一条长度为len的绳子.让你用绳子在三角形内圈出来一个区域让这个区域面积最大

思路:

我可以分类思考这个问题

1.如果绳子比三角形的周长还要长的话,最大的面积显然是这个三角形的面积

2.如果绳子很短呢?显然围成一个圆的时候面积是最大的.三角形里面最大的圆就是它的内切圆了

也就是说如果绳子的长度比三角形内切圆的周长还要小的时候,让它围成一个圆

3.当绳子介于1 2 的长度之间的呢?

当自由线从内切圆那种情况继续膨胀到能与三角形的边贴近但长度小于三角形周长时，将这个已经围成的面积划分为三个部分：
能构成一个更小的内切圆的三段弧，以三段弧的中心连结起来的一个更小且与原三角形相似的三角形，与原三角形贴近的三条边
所围成的三个矩形面积.

 #include <bits/stdc++.h>

 using namespace std;
 #define eps 1e-9
 #define pi acos(-1.0)
 #define zero(x) (((x)>0?(x):-(x))<eps)
 double a,b,c,len;
 int main()
 {
     //freopen("de.txt","r",stdin);
     int casee = ;
     while (~scanf("%lf%lf%lf%lf",&a,&b,&c,&len)){
         //len为自由线的长度；p0为原三角形的周长；p1为原三角形的半周长；
         //R为原三角形的内切圆半径；r为相似三角形的内切圆半径。
         if (zero(a+b+c+len))
             break;
         printf("Case %d: ",++casee);
         double p0 = a+b+c;
         double p1 = p0/;
         double S = sqrt(p1*(p1-a)*(p1-b)*(p1-c));
         double R = *S/p0;//三角形内切圆公式S=p0*R/2; R为内切圆半径
         if (len>=p0)
         {
             printf("%.2f\n",S);
             continue;
         }
         else if (*pi*R-len>eps){
             R = len/pi/;
             S = pi*R*R;
             printf("%.2f\n",S);
             continue;
         }
         else{
             double r = (p0-len)/(p0/R-*pi);
             //利用的就是三角形相似的原理；公式；p0/R*(R-r)=len-2*pi*r;左边是通过内切圆半径与周长的关系求
             //得小三角形的周长；右边是通过自由线的长度减掉三段弧得到相似三角形的周长；
             double radio = (R-r)/R;//相似比
             a*=radio;
             b*=radio;
             c*=radio;
             double p=(a+b+c)/;
             S = pi*r*r+sqrt(p*(p-a)*(p-b)*(p-c))+r**p;
             printf("%.2f\n",S);
             continue;

         }
     }
     return ;
 }