[HDU]1006 Tick and Tick

这是HDU，Online Judge上的一道题，http://acm.hdu.edu.cn/showproblem.php?pid=1006闲来没事，研究的做做，开始思路是模拟每一秒，然后计算满足条件占所有时间的比例，最后发现计算结果不够准确，然后采用了找交集的思路，先用相对角速度，求出两两之间多长时间后会出现一次夹角循环，然后找出在12小时内他们相互满足条件的交集。我这里给出的参考也是参考了晒代码网上的参考答案。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
using namespace std;
#define CIRCLE 43200.000001
//三个指针的相对角速度
const double s_h = 719. / 120, s_m = 59. / 10, m_h = 11. /120;
//将角速度换为周期，即两指针多长时间后出现一次夹角循环，43200秒即12小时，12小时三个指针重合一次
const double tsh = 43200. / 719, tsm = 3600. / 59, tmh = 43200. / 11;
inline double max(double a, double b, double c)
{
double max = a;
if(b > max)
max = b;
if(c > max)
max = c;
return max;
}
inline double min(double a, double b, double c)
{
double min = a;
if(b < min)
min = b;
if(c < min)
min = c;
return min;
}
int main(int argc, const char *argv[]){
double d;
while(scanf("%lf", &d), d != -1){
double bsm, bsh, bmh, esm, esh, emh, begin, end, total = 0;
//两两之间第一次满足条件的时间
bsm = d / s_m;
bsh = d / s_h;
bmh = d / m_h;
//两两之间第一次满足条件结束的时间
esm = (360 - d) / s_m;
esh = (360 - d) / s_h;
emh = (360 - d) / m_h;
//通过循环找出，三者之间共同的交集,
// |____|    |____|    |____|     |____|
//   |____|    |____|    |____|    |____|
//     |____|    |____|     |____|     |____|
//因为CIRCLE为一个周期，所以满足范围必定不会跨越CIRCLE，所以et3 <= CIRCLE合理
for(double bt3 = bsh, et3 = esh; et3 <= CIRCLE; bt3 += tsh, et3 += tsh){
for(double bt2 = bmh, et2 = emh; et2 < CIRCLE; bt2 += tmh, et2 += tmh){
if(et2 < bt3)	//        et2  |______|	//et2还需要增长
continue;	//	|______|  	bt3
if(bt2 > et3)	//  |______|    bt2
break;		//        et3	 |______|	//bt2再增长也不会有交集了，break
for(double bt1 = bsm, et1 = esm; et2 <= CIRCLE; bt1 += tsm, et1 += tsm){
if(et1 < bt3 || et1 < bt2)
continue;
if(bt1 > et3 || bt1 > et2)
break;
begin = max(bt1, bt2, bt3);
end = min(et1, et2, et3);
total += (end - begin);
}
}
}
printf("%.3lf\n", total / 432);
}
return 0;
}