Problem:
かぐや

いよいよ今月末に「はやぶさ２」が打ち上げられます。４年前に「はやぶさ」が帰ってきたときは、日本中が盛り上がったのを覚えている人も多いと思います。７年前には、「かぐや」が打ち上げられ、月の周りを回りながらたくさんの鮮明な画像を地球に送ってくれました。

上の図は、地球を原点とした空間座標（z軸は紙面から垂直に下から上に向いているものとします）に、月の軌道といくつかの月の位置、月の周りを回るかぐやの軌道を描いたものです。月の軌道は、x軸とy軸を通る平面上にある、原点を中心とする円とします。月の周りを回るかぐやの軌道は、x軸とz軸を通る平面に平行な平面上にある円とし、その中心は月の中心と一致します。月は、その軌道に沿って描かれた矢印の方向に回っているものとします。

右の図で、月の位置はＡ，Ｂ，Ｃの３つが描かれています。月を横切る直線がかぐやの軌道です。かぐやは月の周りを回っているので軌道は円ですが、z軸の正の方向から見ているので、図のx軸に平行な直線に見えます（月の位置が変わっても、つねにx軸に平行であることに注意してください）。かぐやは、その軌道上に描かれた矢印の方向に回っているものとします。

かぐやが地球から見て月の裏側に隠れてしまうと、地球と直接通信ができなくなります。かぐやの管制を担当するあなたは、かぐやが月の裏側に隠れる時間が、与えられた時間の中でどれだけになるかをプログラムで求めようとしています。

月の地球に対する位置と分単位での時間tが与えられたとき、その位置から始めてt分後までの間に、かぐやが月の裏側に隠れる時間を求めるプログラムを作成してください。ただし、地球とかぐやは点、月は半径 1800 km の球と考えます。月は半径 380000 km の軌道を 2500000 秒で一周し、かぐやは月の表面から高度 100 km の円上を２時間で一周するものとします。かぐやの最初の位置は、かぐやの軌道がz座標で最大の値をとる位置とします。

入力

入力は以下の形式で与えられる。

m t

m (0 ≤ m < 360) は月の位置を上の図のx軸の正の部分からy軸の正の部分に向かって反時計回りに測った角度を整数で表したものである。t (1 ≤ t ≤ 10000) は分で測った時間を表す整数である。

出力

はじめの位置から t 分経過するまでの間に、かぐやが月の裏側に隠れる時間（分）を実数で出力する。ただし、誤差がプラスマイナス 1.0 分を超えてはならない。この条件を満たせば小数点以下何桁表示してもよい。

入出力例

入力例１

90 10

出力例１

0.0

入力例２

0 120

出力例２

47.73

Solution in C (judge verdict: Accepted)

// AOJ 0311: Kaguya
// 2017.10.31

#include <stdio.h>
#include <math.h>

#define M_PI 3.1415926535897932384626433832795

#define dt 0.01

int main()
{
	double t, m, t_end;
	double dR, dTm, dTk, dWm0, dWm, dWk, ans;
	double cm, sm, ck, sk, xm, ym, zm, xk, yk, zk, xk_leng2, dot, d, s; 

	scanf("%lf %lf", &m, &t_end);
	dR   = 1800.0; 		//????????????
	dTm  = 380000.0; 		//????????¬??¢??????
	dTk  = dR + 100.0; 	//??????????????¬??¢??????
	dWm0 = m*M_PI/180.0; 
	dWm  = 2.0*M_PI/(2500000/60.0);
	dWk  = 2.0*M_PI/(120.0); 
	for (ans = 0, t = 0; t + 0.5*dt < t_end; t += dt) {
		cm = cos(dWm * t + dWm0);
		sm = sin(dWm * t + dWm0);
		ck = cos(dWk * t);
		sk = sin(dWk * t);
		xm = dTm*cm;
		ym = dTm*sm;
		zm = 0.0;
		xk = xm - dTk*sk; 
		yk = ym; 
		zk = zm + dTk*ck; 
		xk_leng2 = xk*xk + yk*yk + zk*zk; 
		dot = xk*xm + yk*ym + zk*zm; 
		d = dot*dot - xk_leng2*(dTm*dTm - dR*dR); 
		if (d > 0.0) {
			s = (dot + sqrt(d))/xk_leng2;
			if (s < 1.0) ans += dt;
		}
	}
	printf("%lf\n", ans);
}