题目链接: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;
}