围度还原

#include 

#define INF 0x3f3f3f3f
#define eps 1e-6

typedef long long LL;
const double pi = acos(-1.0);
const long long mod = 1e9;
using namespace std;

const int num = 36;
const double r = 400;
const double angle = 10.0 / 180.0 * pi;

struct Point
{
    double Dis,Rad;
    double x,y;
}p1[105],p2[105];

double Length(Point a,Point b)
{
    double k = (a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y);
    return sqrt(k);
}

int main()
{
    ios_base::sync_with_stdio(false); cin.tie(0);
    freopen("int.txt","r",stdin);
    freopen("out.txt","w",stdout);
    double MAX1 = 0,MIN11 = INF,MIN12 = INF;
    for(int i = 0;i < num;i++)
    {
        cin >> p1[i].Dis;
        p1[i].Rad = r - p1[i].Dis;
        MAX1 = max(MAX1,p1[i].Rad);
        if(i < num / 2)
            MIN11 = min(MIN11,p1[i].Rad);
        else
            MIN12 = min(MIN12,p1[i].Rad);
        p1[i].x = p1[i].Rad * cos(angle * i);
        p1[i].y = p1[i].Rad * sin(angle * i);

    }

    double MAX2 = 0,MIN21 = INF,MIN22 = INF;
    for(int i = 0;i < num;i++)
    {
        cin >> p2[i].Dis;
        p2[i].Rad = r - p2[i].Dis;
        MAX2 = max(MAX2,p2[i].Rad);
        if(i < num / 2)
            MIN21 = min(MIN21,p2[i].Rad);
        else
            MIN22 = min(MIN22,p2[i].Rad);
        p2[i].x = p2[i].Rad * cos(angle * i);
        p2[i].y = p2[i].Rad * sin(angle * i);

    }

    for(int i = 0;i < num;i++)
        printf("%.2f %.2f %.2f %.2f\n",p1[i].x,p1[i].y,p2[i].x,p2[i].y);


    p1[num].Rad = r - p1[0].Dis;
    double sum1_R = 0;
    for(int i = 0;i <= num;i++)
        sum1_R += p1[i].Rad;
    double average_R1 = sum1_R / (num + 1);
    double ans11 = average_R1 * pi * 2.0;

    double ans12 = 0;
    for(int i = 1;i < num;i++)
        ans12 += Length(p1[i],p1[i - 1]);
    ans12 += Length(p1[num - 1],p1[0]);

    double ans13 = 0;
    for(int i = 1;i < num;i++)
    {

        double dis1 = Length(p1[i],p1[i - 1]);
        double r1 = dis1 / 2.0 / sin(angle / 2.0);
        ans13 += angle * r1;
        //if(i == 6 || i == 7 || i == 16 || i == 17 || i == 18)
            //ans3 -= angle * r;
    }
    double dis1 = Length(p1[num - 1],p1[0]);
    double r1 = dis1 / 2.0 / sin(angle / 2.0);
    ans13 += angle * r1;

    //double kkk1 = angle * r1 - Length(p1[num - 1],p1[0]);

    p2[num].Rad = r - p2[0].Dis;
    double sum2_R = 0;
    for(int i = 0;i <= num;i++)
        sum2_R += p2[i].Rad;
    double average_R2 = sum2_R / (num + 1);
    double ans21 = average_R2 * pi * 2.0;

    double ans22 = 0;
    for(int i = 1;i < num;i++)
        ans22 += Length(p2[i],p2[i - 1]);
    ans22 += Length(p2[num - 1],p2[0]);

    double ans23 = 0;
    for(int i = 1;i < num;i++)
    {

        double dis2 = Length(p2[i],p2[i - 1]);
        double r2 = dis2 / 2.0 / sin(angle / 2.0);
        ans23 += angle * r2;
    }
    double dis2 = Length(p2[num - 1],p2[0]);
    double r2 = dis2 / 2.0 / sin(angle / 2.0);
    ans23 += angle * r2;
    //double kkk2 = angle * r2 - Length(p2[num - 1],p2[0]);

    //printf("%.10f %.10f\n",kkk1,kkk2);

    printf("测量点之间的角度差：%.3f度\n",angle * 180 / pi);
    printf("测量点的个数：%d\n",num);
    printf("测量部位：无干扰腰\n");
    printf("单位：mm \n");
    printf("人体模型半径平均值%.3f\n",average_R1);
    printf("围度计算值1（平均半径，当成圆来计算）：%.3f\n",ans11);
    printf("围度计算值2（两点间距离公式算弦长）：%.3f\n",ans12);
    printf("围度计算值3（用弦长拟合圆弧）: %.3f\n",ans13);
    printf("围度实测值：800\n");
    printf("两个椭圆半周长拟合值：%.3f\n", pi / 2 * (sqrt(2 * MAX1 * MAX1 + 2 * MIN11 * MIN11) +  sqrt(2 * MAX1 * MAX1 + 2 * MIN12 * MIN12)));
    puts("");

    printf("测量部位：20.5度水平间距下测量值\n");
    printf("人体模型半径平均值%.3f\n",average_R2);
    printf("围度计算值1（平均半径，当成圆来计算）：%.3f\n",ans21);
    printf("围度计算值2（两点间距离公式算弦长）：%.3f\n",ans22);
    printf("围度计算值3（用弦长拟合圆弧）: %.3f\n",ans23);
    printf("围度实测值：800\n");
    printf("两个椭圆半周长拟合值：%.3f\n", pi / 2 * (sqrt(2 * MAX2 * MAX2 + 2 * MIN21 * MIN21) +  sqrt(2 * MAX2 * MAX2 + 2 * MIN22 * MIN22)));


    return 0;
}