QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#209976#7110. Kuririn MIRACLEPetroTarnavskyiWA 418ms4208kbC++171.5kb2023-10-10 20:38:522023-10-10 20:38:53

Judging History

你现在查看的是最新测评结果

  • [2023-10-10 20:38:53]
  • 评测
  • 测评结果:WA
  • 用时:418ms
  • 内存:4208kb
  • [2023-10-10 20:38:52]
  • 提交

answer

#include <bits/stdc++.h>

using namespace std;

#define FOR(i, a, b) for(int i = (a); i < (b); i++)
#define RFOR(i, a, b) for(int i = (a) - 1; i >= (b); i--)
#define SZ(a) int(a.size())
#define ALL(a) a.begin(), a.end()
#define PB push_back
#define MP make_pair
#define F first	
#define S second

typedef long long LL;
typedef vector<int> VI;
typedef pair<int, int> PII;
typedef double db;

const int M = 200000;
const db PI = acos(-1.0);
const db EPS = 1e-7;

db v, r, d;

db f(db phi)
{
	db si = sin(phi);
	return -v / (2 * r) * (sqrt(3 + si * si) - si);
}

void solve()
{
	cin >> v >> r >> d;
	db ans = d / v;
	db phi = PI;
	const db dt = ans / M;
	FOR(i, 1, M + 1)
	{
		db t = i * dt;
		phi += dt * f(phi);
		db si = sin(phi), co = cos(phi);
		db a = d - 2 * r - v * t, b = a - 2 * r * co;
		if (abs(b) < EPS && phi < EPS)
		{
			ans = t;
			break;
		}
		if (b < -EPS)
			break;
		if (phi < PI / 2)
		{
			db l = sqrt(4 * r * r + a * a - 4 * r * a * cos(phi));
			db siBeta = 2 * r * si / l;
			assert(0 <= siBeta && siBeta < 1 + EPS);
			db beta = asin(min(1.0, siBeta));
			if (phi + beta <= PI / 2 && 2 * cos(phi + beta) >= cos(phi))
			{
				ans = min(ans, t + l / (2 * v));
				break;
			}
		}
		else if (phi < EPS)
		{
			ans = t + a / (2 * v);
			break;
		}
	}
	cout << ans << "\n";
	//cerr << ans - 8.310579933902352 << "\n";
}

int main()
{
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout << fixed << setprecision(15);
	int t;
	cin >> t;
	while (t--)
		solve();
	return 0;
}

詳細信息

Test #1:

score: 100
Accepted
time: 5ms
memory: 4208kb

input:

1
2.00 3 30.0

output:

8.310572894044750

result:

ok found '8.3105729', expected '8.3105799', error '0.0000008'

Test #2:

score: -100
Wrong Answer
time: 418ms
memory: 4024kb

input:

100
8.66 6.05 71.59
8.44 8.67 82.55
4.23 9.68 24.74
6.63 9.95 7.45
1.2 7.27 97.55
2.08 2.94 31.61
8.48 5.56 20.16
7.35 5.55 52.48
1.35 8.15 51.23
3.88 9.02 72.37
7.03 1.97 57.03
1.78 4.44 93.44
5.5 6.11 52.94
1.09 9.64 69.01
1.01 3.54 16.78
1.63 4.55 44.06
5.82 5.51 5.92
5.62 2.93 76.01
4.75 4.43 73...

output:

4.457832178995005
5.488357405882677
5.848699763593380
1.123680241327300
43.266428662029696
8.305481819823621
2.377358490566038
4.015483707940933
37.948148148148142
18.652061855670105
4.155118123835539
27.178927696135439
9.625454545454545
63.311926605504588
16.613861386138616
15.093983558138465
1.017...

result:

wrong answer 1st numbers differ - expected: '4.4578371', found: '4.4578322', error = '0.0000011'