QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#209903#7110. Kuririn MIRACLEPetroTarnavskyi#TL 18ms3712kbC++171.3kb2023-10-10 19:14:532023-10-10 19:14:53

Judging History

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

  • [2023-10-10 19:14:53]
  • 评测
  • 测评结果:TL
  • 用时:18ms
  • 内存:3712kb
  • [2023-10-10 19:14:53]
  • 提交

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 long double db;

const int M = 400000;
const db PI = acos(-1.0);

void solve()
{
	db v, r, d;
	cin >> v >> r >> d;
	db ans = d / v;
	const db dt = ans / M;
	db phiP = PI - v / (2 * r) * sqrt(3) * dt, phiPP = PI;
	FOR(i, 2, M + 1)
	{
		db t = i * dt;
		db siP = sin(phiP);
		db der = -v / (2 * r) * (sqrt(3 + siP * siP) - siP);
		db phi = phiPP + 2 * dt * der;
		if (phi < PI / 2)
		{
			db a = d - 2 * r - v * t;
			if (a < 0)
				break;
			db si = sin(phi), co = cos(phi);
			db l = sqrt(4 * r * r + a * a - 4 * r * a * co);
			db siBeta = 2 * r * si / l;
			db beta = asin(siBeta);
			if (2 * cos(phi + beta) >= co)
			{
				ans = min(ans, t + l / (2 * v));
				break;	
			}
		}
		phiPP = phiP;
		phiP = phi;
	}
	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: 18ms
memory: 3712kb

input:

1
2.00 3 30.0

output:

8.310579934056235

result:

ok found '8.3105799', expected '8.3105799', error '0.0000000'

Test #2:

score: -100
Time Limit Exceeded

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:


result: