QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#421926 | #6763. Triangle Pendant | 36champ | WA | 3ms | 3936kb | C++20 | 3.8kb | 2024-05-26 11:22:20 | 2024-05-26 11:22:21 |
Judging History
answer
#include <bits/stdc++.h>
#define double long double
#define pb push_back
using namespace std;
double eps = -1e-9;
struct point{
double x, y, z;
point operator+ (point B) {return {x + B.x, y + B.y, z + B.z};}
point operator- (point B) {return {x - B.x, y - B.y, z - B.z};}
point operator* (double d) {return {x * d, y * d, z * d};}
point operator/ (double d) {return {x / d, y / d, z / d};}
};
ostream& operator<<(ostream& os, point p)
{
return os << "("<< p.x << "," << p.y << "," << p.z << ")";
}
double dot(point A, point B)
{
return A.x * B.x + A.y * B.y + A.z * B.z;
}
double norm(point P)
{
return sqrt(dot(P, P));
}
double sq(point A, point B)
{
return dot(A - B, A - B);
}
double dist(point A, point B)
{
return sqrt(sq(A, B));
}
point normalize(point P)
{
return P / norm(P);
}
void cycle(double &a, double &b, double &c)
{
double x = a;
a = b;
b = c;
c = x;
}
void solve()
{
double la, lb, lc, a, b, c;
cin >> la >> lb >> lc >> a >> b >> c;
for(int i=0; i<3; i++)
{
//Hang by C
point A{0, 0, 0}, B{c, 0, 0}, C, D, M;
C = {(c*c+b*b-a*a)/(2*c), 0, 0};
C.y = sqrt(b * b - C.x * C.x);
M = (A + B + C) / 3;
D = C + normalize(C - M) * lc;
if(dist(A, D) < la + eps && dist(B, D) < lb + eps)
{
//cout << "HANG C " << i << "\n";
point X = normalize(M - D);
double ha = - dot(A - D, X);
double hb = - dot(B - D, X);
double hc = - dot(C - D, X);
if(i==0) cout << ha << " " << hb << " " << hc << "\n";
if(i==1) cout << hc << " " << ha << " " << hb << "\n";
if(i==2) cout << hb << " " << hc << " " << ha << "\n";
return;
}
else
{
cycle(a, b, c);
cycle(la, lb, lc);
}
}
for(int i=0; i<3; i++)
{
//Hang by AB
point A{0, 0, 0}, B{c, 0, 0}, C, D;
C = {(c*c+b*b-a*a)/(2*c), 0, 0};
C.y = sqrt(b * b - C.x * C.x);
D = {(c*c+la*la-lb*lb)/(2*c), 0, 0};
D.y = -sqrt(la*la - D.x * D.x);
//cout << A << " " << B << " " << C << " " << D << "\n";
if(dist(C, D) < lc + eps)
{
//cout << "HANG AB " << i << "\n";
point M = (A + B + C) / 3;
point X = normalize(M - D);
double ha = - dot(A - D, X);
double hb = - dot(B - D, X);
double hc = - dot(C - D, X);
if(i==0) cout << ha << " " << hb << " " << hc << "\n";
if(i==1) cout << hc << " " << ha << " " << hb << "\n";
if(i==2) cout << hb << " " << hc << " " << ha << "\n";
return;
}
else
{
cycle(a, b, c);
cycle(la, lb, lc);
}
}
//Tetrahedron
point A{0, 0, 0}, B{c, 0, 0}, C, D;
C = {(c*c+b*b-a*a)/(2*c), 0, 0};
C.y = sqrt(b * b - C.x * C.x);
D.x = (c * c + la * la - lb * lb) / (2 * c);
D.y = (dot(C, C) + la*la - lc*lc - 2 * C.x * D.x) / (2 * C.y);
D.z = sqrt(la * la - D.x * D.x - D.y * D.y);
point M = (A + B + C) / 3;
point X = normalize(M - D);
double ha = - dot(A - D, X);
double hb = - dot(B - D, X);
double hc = - dot(C - D, X);
//cout << A << " " << B << " " << C << " " << D << " " << M << "\n";
//cout << A - D << "." << X << " = " << dot(A - D, X) << "\n";
cout << ha << " " << hb << " " << hc << "\n";
}
int main()
{
ios_base::sync_with_stdio(false);
cin.tie(nullptr);
cout.tie(nullptr);
cout << fixed << setprecision(12);
int T;
cin >> T;
while(T-->0)
{
solve();
}
}
詳細信息
Test #1:
score: 100
Accepted
time: 1ms
memory: 3936kb
input:
2 1 1 1 1 1 1 2 3 3 1 1 1
output:
-0.816496580928 -0.816496580928 -0.816496580928 -2.000000000000 -2.866025403784 -2.866025403784
result:
ok 6 numbers
Test #2:
score: -100
Wrong Answer
time: 3ms
memory: 3844kb
input:
1000 21 2 14 12 13 4 29 19 13 15 10 17 29 24 15 29 24 23 29 17 30 18 9 25 27 24 30 16 4 15 28 13 17 12 21 16 16 22 10 22 15 8 15 23 24 23 27 13 26 3 27 15 16 17 5 8 20 17 6 12 24 14 13 15 19 13 27 22 18 18 23 30 22 18 14 11 29 28 7 13 22 22 17 11 19 9 16 22 20 17 21 14 20 6 29 25 20 10 19 9 27 27 17...
output:
-2.935856727587 -2.000000000000 -13.352348999858 -22.146476218003 -16.076204705477 -12.241826659011 -28.027431433360 -18.318582636183 -8.243362186282 -27.266119651334 -13.159319584326 -26.882525397693 -26.534933366926 -22.066632132649 -29.616079568530 -26.967804477492 -12.557669395974 -14.1867221345...
result:
wrong answer 382nd numbers differ - expected: '-24.86301', found: '-24.27354', error = '0.02371'