QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#209975#7110. Kuririn MIRACLEPetroTarnavskyiTL 836ms4244kbC++171.5kb2023-10-10 20:38:332023-10-10 20:38:33

Judging History

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

  • [2023-10-10 20:38:33]
  • 评测
  • 测评结果:TL
  • 用时:836ms
  • 内存:4244kb
  • [2023-10-10 20:38:33]
  • 提交

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 = 400000;
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;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 10ms
memory: 4176kb

input:

1
2.00 3 30.0

output:

8.310576414006437

result:

ok found '8.3105764', expected '8.3105799', error '0.0000004'

Test #2:

score: 0
Accepted
time: 836ms
memory: 4244kb

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.457834650909218
5.488359107464506
5.848699763593380
1.123680241327300
43.266454357736990
8.305485979043274
2.377358490566038
4.015484838745667
37.948148148148142
18.652061855670105
4.155120883729960
27.178945351900254
9.625454545454545
63.311926605504588
16.613861386138616
15.093989009710338
1.017...

result:

ok 100 numbers

Test #3:

score: -100
Time Limit Exceeded

input:

1000
7.52 6.68 80.67
5.34 6.82 45.75
3.57 6.35 22.35
9.55 3.65 31.38
9.37 5.73 40.42
1.83 8.38 82.97
4.86 6.13 9.45
3.88 5.34 84.19
6.49 4.15 85.24
6.51 7.23 42.97
3.0 5.98 57.53
9.35 2.76 69.07
6.76 9.16 91.24
2.98 9.89 49.2
9.16 3.85 66.14
1.84 3.03 37.43
8.01 5.04 50.98
4.05 9.86 4.76
5.4 5.49 60...

output:


result: