HDU 1006 Tick and Tick

题目链接:Tick and Tick

解题思路:这道题要不是看了别人的博客还是做不出来的，刚开始我想得是按秒来暴力，但是小数点后面几位一直不对，后来知道了要按照分钟来找，把每分钟happy的秒数求出来，这个秒数不一定是整数，应为求出的区间不一定前后都是整数。

设时间为h:m:s，角度为rh:rm:rs，那么rs = 6 * s;    rm = 6 * m + 0.1 * s;    rh = 30 * h + 0.5 * m + s / 120;

两两相减可以得到   d  <= |ax + b| <= 360 - d  的等式  按照a为负数或者非负数算出两个区间，将这两个区间去并集， 在对三个区间取交集。

其实这里求出6个区间后 把他们每三个一组取交集，再把每一个结果整合在一起就是这一分钟的答案。每分钟求出来后在去求百分比。

#include<stdio.h>

struct A{
double l, r;
};
typedef struct A node;

double d;

node calc(double a, double b){//不等式求解
node tem;
if(a > 0){
tem.l = (d - b) / a;
tem.r = (360 - d - b) / a;
}
else{
tem.r = (d - b) / a;
tem.l = (360 - d - b) / a;
}
if(tem.l < 0)	tem.l = 0;
if(tem.r > 60)	tem.r = 60;
if(tem.l >= tem.r)	tem.r = tem.l = 0;
return tem;
}

node com(node a, node b){//交集运算
node tem;
tem.l = a.l > b.l ? a.l : b.l;
tem.r = a.r < b.r ? a.r : b.r;
if(tem.l >= tem.r)	tem.r = tem.l = 0;
return tem;
}

double slove(int h, int m){
double a, b, tot = 0;
int i, j, k;
node ans[3][2], tem;
// 1
a = 0.1 - 6;
b = 6 * m;
ans[0][0] = calc(a, b);
ans[0][1] = calc(-a, -b);
//2
a = 1.0 / 120 - 6;
b = 30 * h + 0.5 * m;
ans[1][0] = calc(a, b);
ans[1][1] = calc(-a, -b);
//3
a = 1.0 / 120 - 0.1;
b = 30 * h + (0.5 - 6) * m;
ans[2][0] = calc(a, b);
ans[2][1] = calc(-a, -b);
//这里做每三个一组的交集运算
for(i = 0; i < 2; i++){
for(j = 0; j < 2; j++){
for(k = 0; k < 2; k++){
tem = com(ans[0][i], com(ans[1][j], ans[2][k]));
tot += tem.r - tem.l;
}
}
}
return tot;
}

int main(){
int i, j, k;
double sum;
while(scanf("%lf", &d)){
if(d == -1)		return 0;
sum = 0;
for(i = 0; i < 12; i++){
for(j = 0; j < 60; j++){
sum += slove(i, j);
}
}
printf("%.3lf\n", sum * 100.0 / 43200);
}
return 0;
}