POJ1927 Area in Triangle

 

Time Limit: 1000MS		Memory Limit: 30000K
Total Submissions: 1458		Accepted: 759

Description

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

Source

Beijing 2004 Preliminary@POJ

数学问题 计算几何

（物理题？）

画图太麻烦了，随便贴个别处的题解吧233 http://blog.csdn.net/xuechelingxiao/article/details/40707691

没加case傻傻WA了两次

/*by SilverN*/
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<vector>
using namespace std;
const int mxn=100010;
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
double a,b,c,d;
inline double Sin(double a,double b,double c){
    double Cos=(a*a+b*b-c*c)/(2*a*b);
    return sqrt(1-Cos*Cos);
}
double calc(double a,double b,double c,double d){
    double Pi=acos(-1.0);
    double S=0.5*a*b*Sin(a,b,c);
    if(d-(a+b+c)>=0){return S;}
    double h=S*2/(a+b+c);
    if(h*2*Pi>=d){return d*d/4/Pi;}
    double Len=a+b+c;
    double T=(Len-d)/(Len-2*Pi*h);
    return S-S*T*T+(h*T)*(h*T)*Pi;
}
int main(){
    int i,j,cas=0;
    while(scanf("%lf%lf%lf%lf",&a,&b,&c,&d)!=EOF){
        if(a==0 && b==0 && c==0 && d==0)break;
        printf("Case %d: %.2f
",++cas,calc(a,b,c,d));
    }
    return 0;
}